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.

Intra-condensate demixing of TDP-43 inside stress granules generates pathological aggregates.

Cytosolic aggregation of the nuclear protein TDP-43 is associated with many neurodegenerative diseases, but the triggers for TDP-43 aggregation are still debated. Here, we demonstrate that TDP-43 aggregation requires a double event. One is up-concentration in stress granules beyond a threshold, and the other is oxidative stress. These two events collectively induce intra-condensate demixing, giving rise to a dynamic TDP-43 enriched phase within stress granules, which subsequently transitions into pathological aggregates. Mechanistically, intra-condensate demixing is triggered by local unfolding of the RRM1 domain for intermolecular disulfide bond formation and by increased hydrophobic patch interactions in the C-terminal domain. By engineering TDP-43 variants resistant to intra-condensate demixing, we successfully eliminate pathological TDP-43 aggregates in cells. We conclude that up-concentration inside condensates and simultaneous exposure to environmental stress could be a general pathway for protein aggregation, with intra-condensate demixing constituting a key intermediate step.

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.

Specific residues and conformational plasticity define the substrate specificity of short-chain dehydrogenases/reductases.

Short-chain dehydrogenases/reductases (SDRs) are one of the most prevalent enzyme families distributed among the sequenced microorganisms. Despite the presence of a conserved catalytic tetrad and high structural similarity, these enzymes exhibit different substrate specificities. The insufficient knowledge regarding the amino acids underlying substrate specificity hinders the understanding of the SDRs' roles in diverse and significant biological processes. Here, we performed bioinformatic analysis, molecular modeling, and mutagenesis studies to identify the key residues that regulate the substrate specificities of two homologous microbial SDRs (i.e., DesE and KduD). Further, we investigated the impact of altering the physicochemical properties of these amino acids on enzyme activity. Interestingly, molecular dynamics simulations also suggest a critical role of enzyme conformational flexibility in substrate recognition and catalysis. Overall, our findings improve the understanding of microbial SDR substrate specificity and shed light on future rational design of more efficient and effective biocatalysts.

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.

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

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

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

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

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

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

Principles Governing the Phase Separation of Multidomain Proteins.

A variety of membraneless organelles, often termed "biological condensates", play an important role in the regulation of cellular processes such as gene transcription, translation, and protein quality control. On the basis of experimental and theoretical investigations, liquid-liquid phase separation (LLPS) has been proposed as a possible mechanism for the origin of biological condensates. LLPS requires multivalent macromolecules that template the formation of long-range, intermolecular interaction networks and results in the formation of condensates with defined composition and material properties. Multivalent interactions driving LLPS exhibit a wide range of modes from highly stereospecific to nonspecific and involve both folded and disordered regions. Multidomain proteins serve as suitable macromolecules for promoting phase separation and achieving disparate functions due to their potential for multivalent interactions and regulation. Here, we aim to highlight the influence of the domain architecture and interdomain interactions on the phase separation of multidomain protein condensates. First, the general principles underlying these interactions are illustrated on the basis of examples of multidomain proteins that are predominantly associated with nucleic acid binding and protein quality control and contain both folded and disordered regions. Next, the examples showcase how LLPS properties of folded and disordered regions can be leveraged to engineer multidomain constructs that form condensates with the desired assembly and functional properties. Finally, we highlight the need for improvements in coarse-grained computational models that can provide molecular-level insights into multidomain protein condensates in conjunction with experimental efforts.

NEDD8 Deamidation Inhibits Cullin RING Ligase Dynamics.

Cullin-RING ligases (CRLs) are a significant subset of Ubiquitin E3 ligases that regulate multiple cellular substrates involved in innate immunity, cytoskeleton modeling, and cell cycle. The glutamine deamidase Cycle inhibitory factor (Cif) from enteric bacteria inactivates CRLs to modulate these processes in the host cell. The covalent attachment of a Ubiquitin-like protein NEDD8 catalytically activates CRLs by driving conformational changes in the Cullin C-terminal domain (CTD). NEDDylation results in a shift from a compact to an open CTD conformation through non-covalent interactions between NEDD8 and the WHB subdomain of CTD, eliminating the latter's inhibitory interactions with the RING E3 ligase-Rbx1/2. It is unknown whether the non-covalent interactions are sufficient to stabilize Cullin CTD's catalytic conformation. We studied the dynamics of Cullin-CTD in the presence and absence of NEDD8 using atomistic molecular dynamics (MD) simulations. We uncovered that NEDD8 engages in non-covalent interactions with 4HB/αß subdomains in Cullin-CTD to promote open conformations. Cif deamidates glutamine 40 in NEDD8 to inhibit the conformational change in CRLs by an unknown mechanism. We investigated the effect of glutamine deamidation on NEDD8 and its interaction with the WHB subdomain post-NEDDylation using MD simulations and NMR spectroscopy. Our results suggest that deamidation creates a new intramolecular salt bridge in NEDD8 to destabilize the NEDD8/WHB complex and reduce CRL activity.

Deamidation disrupts native and transient contacts to weaken the interaction between UBC13 and RING-finger E3 ligases.

The deamidase OspI from enteric bacteria Shigella flexneri deamidates a glutamine residue in the host ubiquitin-conjugating enzyme UBC13 and converts it to glutamate (Q100E). Consequently, its polyubiquitination activity in complex with the RING-finger ubiquitin ligase TRAF6 and the downstream NF-κB inflammatory response is silenced. The precise role of deamidation in silencing the UBC13/TRAF6 complex is unknown. We report that deamidation inhibits the interaction between UBC13 and TRAF6 RING-domain (TRAF6RING) by perturbing both the native and transient interactions. Deamidation creates a new intramolecular salt-bridge in UBC13 that competes with a critical intermolecular salt-bridge at the native UBC13/TRAF6RING interface. Moreover, the salt-bridge competition prevents transient interactions necessary to form a typical UBC13/RING complex. Repulsion between E100 and the negatively charged surface of RING also prevents transient interactions in the UBC13/RING complex. Our findings highlight a mechanism wherein a post-translational modification perturbs the conformation and stability of transient complexes to inhibit protein-protein association.