Refining All-Atom Protein Force Fields for Polar-Rich, Prion-like, Low-Complexity Intrinsically Disordered Proteins

Expanding the molecular language of protein liquid-liquid phase separation

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

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.

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.

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

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

Insights into Molecular Diversity within the FUS/EWS/TAF15 Protein Family: Unraveling Phase Separation of the N-Terminal Low-Complexity Domain from RNA-Binding Protein EWS

The FET protein family, comprising FUS, EWS, and TAF15, plays crucial roles in mRNA maturation, transcriptional regulation, and DNA damage response. Clinically, they are linked to Ewing family tumors and neurodegenerative diseases such as amyotrophic lateral sclerosis. The fusion protein EWS::FLI1, the causative mutation of Ewing sarcoma, arises from a genomic translocation that fuses a portion of the low-complexity domain (LCD) of EWS (EWSLCD) with the DNA binding domain of the ETS transcription factor FLI1. This fusion protein modifies transcriptional programs and disrupts native EWS functions, such as splicing. The exact role of the intrinsically disordered EWSLCD remains a topic of active investigation, but its ability to phase separate and form biomolecular condensates is believed to be central to EWS::FLI1's oncogenic properties. Here, we used paramagnetic relaxation enhancement NMR, microscopy, and all-atom molecular dynamics (MD) simulations to better understand the self-association and phase separation tendencies of the EWSLCD. Our NMR data and mutational analysis suggest that a higher density and proximity of tyrosine residues amplify the likelihood of condensate formation. MD simulations revealed that the tyrosine-rich termini exhibit compact conformations with unique contact networks and provided critical input on the relationship between contacts formed within a single molecule (intramolecular) and inside the condensed phase (intermolecular). These findings enhance our understanding of FET proteins' condensate-forming capabilities and underline differences between EWS, FUS, and TAF15.

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.

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

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

Structure-function relationships in protein homorepeats

Homorepeats (or polyX), protein segments containing repetitions of the same amino acid, are abundant in proteomes from all kingdoms of life and are involved in crucial biological functions as well as several neurodegenerative and developmental diseases. Mainly inserted in disordered segments of proteins, the structure/function relationships of homorepeats remain largely unexplored. In this review, we summarize present knowledge for the most abundant homorepeats, highlighting the role of the inherent structure and the conformational influence exerted by their flanking regions. Recent experimental and computational methods enable residue-specific investigations of these regions and promise novel structural and dynamic information for this elusive group of proteins. This information should increase our knowledge about the structural bases of phenomena such as liquid-liquid phase separation and trinucleotide repeat disorders.

"Boundary residues" between the folded RNA recognition motif and disordered RGG domains are critical for FUS-RNA binding

Fused in sarcoma (FUS) is an abundant RNA-binding protein, which drives phase separation of cellular condensates and plays multiple roles in RNA regulation. The RNA-binding ability of FUS protein is crucial to its cellular function. Here, our molecular simulation study on the FUS-RNA complex provides atomic resolution insights into the observations from biochemical studies and also illuminates our understanding of molecular driving forces that mediate the structure, stability, and interaction of the RNA recognition motif (RRM) and RGG domains of FUS with a stem-loop junction RNA. We observe clear cooperativity and division of labor among the ordered (RRM) and disordered domains (RGG1 and RGG2) of FUS that leads to an organized and tighter RNA binding. Irrespective of the length of RGG2, the RGG2-RNA interaction is confined to the stem-loop junction and the proximal stem regions. On the other hand, the RGG1 interactions are primarily with the longer RNA stem. We find that the C terminus of RRM, which make up the "boundary residues" that connect the folded RRM with the long disordered RGG2 stretch of the protein, plays a critical role in FUS-RNA binding. Our study provides high-resolution molecular insights into the FUS-RNA interactions and forms the basis for understanding the molecular origins of full-length FUS interaction with RNA.

Toward Consistent Physics-Based Modeling of Local Backbone Structures and Chirality Change of Proteins in Coarse-Grained Approaches

A reliable representation of local interactions is critical for the accuracy of modeling protein structure and dynamics at both the all-atom and coarse-grained levels. The development of local (mainly torsional) potentials was focused on careful parametrization of the predetermined (usually Fourier) formulas rather than on their physics-based derivation. In this Perspective we discuss the state-of-the-art methods for modeling local interactions, including the scale-consistent theory developed in our laboratory, which implies that the coarse-grained torsional potentials inseparably depend on the virtual-bond angles adjacent to a given dihedral and that multitorsional terms should be considered. We extend the treatment to split the residue-based torsional potentials into the site-based regular and improper torsional potentials. These considerations are illustrated with the revised torsional potentials and improper-torsional potentials involving the l-alanine residue and the improper-torsional potential corresponding to serine-residue enantiomerization. Applications of the new approach in coarse-grained modeling and revising all-atom force fields are discussed.

Complex Conformational Space of the RNA Polymerase II C-Terminal Domain upon Phosphorylation

Intrinsically disordered proteins (IDPs) have been closely studied during the past decade due to their importance in many biological processes. The disordered nature of this group of proteins makes it difficult to observe its full span of the conformational space using either experimental or computational studies. In this article, we explored the conformational space of the C-terminal domain (CTD) of RNA polymerase II (Pol II), which is also an intrinsically disordered low complexity domain, using enhanced sampling methods. We provided a detailed conformational analysis of model systems of CTD with different lengths; first with the last 44 residues of the human CTD sequence and finally the CTD model with 2-heptapeptide repeating units. We then investigated the effects of phosphorylation on CTD conformations by performing simulations at different phosphorylated states. We obtained broad conformational spaces in nonphosphorylated CTD models, and phosphorylation has complex effects on the conformations of the CTD. These complex effects depend on the length of the CTD, spacing between the multiple phosphorylation sites, ion coordination, and interactions with the nearby residues.

Insights into Molecular Diversity within the FET Family: Unraveling Phase Separation of the N-Terminal Low Complexity Domain from RNA-Binding Protein EWS

The FET family proteins, which includes FUS, EWS, and TAF15, are RNA chaperones instrumental in processes such as mRNA maturation, transcriptional regulation, and the DNA damage response. These proteins have clinical significance: chromosomal rearrangements in FET proteins are implicated in Ewing family tumors and related sarcomas. Furthermore, point mutations in FUS and TAF15 are associated with neurodegenerative conditions like amyotrophic lateral sclerosis and frontotemporal lobar dementia. The fusion protein EWS::FLI1, the causative mutation of Ewing sarcoma, arises from a genomic translocation that fuses the low-complexity domain (LCD) of EWS (EWSLCD) with the DNA binding domain of the ETS transcription factor FLI1. This fusion not only alters transcriptional programs but also hinders native EWS functions like splicing. However, the precise function of the intrinsically disordered EWSLCD is still a topic of active investigation. Due to its flexible nature, EWSLCD can form transient interactions with itself and other biomolecules, leading to the formation of biomolecular condensates through phase separation - a mechanism thought to be central to the oncogenicity of EWS::FLI1. In our study, we used paramagnetic relaxation enhancement NMR, analytical ultracentrifugation, light microscopy, and all-atom molecular dynamics (MD) simulations to better understand the self-association and phase separation tendencies of EWSLCD. Our aim was to elucidate the molecular events that underpin EWSLCD-mediated biomolecular condensation. Our NMR data suggest tyrosine residues primarily drive the interactions vital for EWSLCD phase separation. Moreover, a higher density and proximity of tyrosine residues amplify the likelihood of condensate formation. Atomistic MD simulations and hydrodynamic experiments revealed that the tyrosine-rich N and C-termini tend to populate compact conformations, establishing unique contact networks, that are connected by a predominantly extended, tyrosine-depleted, linker region. MD simulations provide critical input on the relationship between contacts formed within a single molecule (intramolecular) and inside the condensed phase (intermolecular), and changes in protein conformations upon condensation. These results offer deeper insights into the condensate-forming abilities of the FET proteins and highlights unique structural and functional nuances between EWS and its counterparts, FUS and TAF15.

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.

Clustering Heterogeneous Conformational Ensembles of Intrinsically Disordered Proteins with t-Distributed Stochastic Neighbor Embedding

Intrinsically disordered proteins (IDPs) populate a range of conformations that are best described by a heterogeneous ensemble. Grouping an IDP ensemble into "structurally similar" clusters for visualization, interpretation, and analysis purposes is a much-desired but formidable task, as the conformational space of IDPs is inherently high-dimensional and reduction techniques often result in ambiguous classifications. Here, we employ the t-distributed stochastic neighbor embedding (t-SNE) technique to generate homogeneous clusters of IDP conformations from the full heterogeneous ensemble. We illustrate the utility of t-SNE by clustering conformations of two disordered proteins, Aβ42, and α-synuclein, in their APO states and when bound to small molecule ligands. Our results shed light on ordered substates within disordered ensembles and provide structural and mechanistic insights into binding modes that confer specificity and affinity in IDP ligand binding. t-SNE projections preserve the local neighborhood information, provide interpretable visualizations of the conformational heterogeneity within each ensemble, and enable the quantification of cluster populations and their relative shifts upon ligand binding. Our approach provides a new framework for detailed investigations of the thermodynamics and kinetics of IDP ligand binding and will aid rational drug design for IDPs.

Highly Charged Proteins and Their Repulsive Interactions Antagonize Biomolecular Condensation

Biomolecular condensation is involved in various cellular processes; therefore, regulation of condensation is crucial to prevent deleterious protein aggregation and maintain a stable cellular environment. Recently, a class of highly charged proteins, known as heat-resistant obscure (Hero) proteins, was shown to protect other client proteins from pathological aggregation. However, the molecular mechanisms by which Hero proteins protect other proteins from aggregation remain unknown. In this study, we performed multiscale molecular dynamics (MD) simulations of Hero11, a Hero protein, and the C-terminal low-complexity domain (LCD) of the transactive response DNA-binding protein 43 (TDP-43), a client protein of Hero11, under various conditions to examine their interactions with each other. We found that Hero11 permeates into the condensate formed by the LCD of TDP-43 (TDP-43-LCD) and induces changes in conformation, intermolecular interactions, and dynamics of TDP-43-LCD. We also examined possible Hero11 structures in atomistic and coarse-grained MD simulations and found that Hero11 with a higher fraction of disordered region tends to assemble on the surface of the condensates. Based on the simulation results, we have proposed three possible mechanisms for Hero11's regulatory function: (i) In the dense phase, TDP-43-LCD reduces contact with each other and shows faster diffusion and decondensation due to the repulsive Hero11-Hero11 interactions. (ii) In the dilute phase, the saturation concentration of TDP-43-LCD is increased, and its conformation is relatively more extended and variant, induced by the attractive Hero11-TDP-43-LCD interactions. (iii) Hero11 on the surface of small TDP-43-LCD condensates can contribute to avoiding their fusion due to repulsive interactions. The proposed mechanisms provide new insights into the regulation of biomolecular condensation in cells under various conditions.

The structure of pathogenic huntingtin exon 1 defines the bases of its aggregation propensity

Huntington's disease is a neurodegenerative disorder caused by a CAG expansion in the first exon of the HTT gene, resulting in an extended polyglutamine (poly-Q) tract in huntingtin (httex1). The structural changes occurring to the poly-Q when increasing its length remain poorly understood due to its intrinsic flexibility and the strong compositional bias. The systematic application of site-specific isotopic labeling has enabled residue-specific NMR investigations of the poly-Q tract of pathogenic httex1 variants with 46 and 66 consecutive glutamines. Integrative data analysis reveals that the poly-Q tract adopts long α-helical conformations propagated and stabilized by glutamine side chain to backbone hydrogen bonds. We show that α-helical stability is a stronger signature in defining aggregation kinetics and the structure of the resulting fibrils than the number of glutamines. Our observations provide a structural perspective of the pathogenicity of expanded httex1 and pave the way to a deeper understanding of poly-Q-related diseases.

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

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

In Silico Studies to Predict the Role of Solvent in Guiding the Conformations of Intrinsically Disordered Peptides and Their Aggregated Protofilaments

The formation of amyloids due to the self-assembly of intrinsically disordered proteins or peptides is a hallmark for different neurodegenerative diseases. For example, amyloids formed by the amyloid beta (Aβ) peptides are responsible for the most devastating neuropathological disease, namely, Alzheimer's disease, while aggregation of α-synuclein peptides causes the etiology of another neuropathological disease, Parkinson's disease. Characterization of the intermediates and the final amyloid formed during the aggregation process is, therefore, crucial for microscopic understanding of the origin behind such diseases, as well as for the development of proper therapeutics to combat those. However, most of the research activities reported in this area have been directed toward examining the early stages of the aggregation process, including probing the conformational characteristics of the responsible protein/peptide in the monomeric state or in small oligomeric forms. This is because the small soluble oligomers have been found to be more deleterious than the final insoluble amyloids. This review discusses some of the recent findings obtained from our simulation studies on Aβ and α-synuclein monomers and small preformed Aβ aggregates. A molecular-level insight of the aggregation process with a special emphasis on the role of water in inducing the aggregation process has been provided.

Developing Bonded Potentials for a Coarse-Grained Model of Intrinsically Disordered Proteins

Recent advances in residue-level coarse-grained (CG) computational models have enabled molecular-level insights into biological condensates of intrinsically disordered proteins (IDPs), shedding light on the sequence determinants of their phase separation. The existing CG models that treat protein chains as flexible molecules connected via harmonic bonds cannot populate common secondary-structure elements. Here, we present a CG dihedral angle potential between four neighboring beads centered at Cα atoms to faithfully capture the transient helical structures of IDPs. In order to parameterize and validate our new model, we propose Cα-based helix assignment rules based on dihedral angles that succeed in reproducing the atomistic helicity results of a polyalanine peptide and folded proteins. We then introduce sequence-dependent dihedral angle potential parameters (εd) and use experimentally available helical propensities of naturally occurring 20 amino acids to find their optimal values. The single-chain helical propensities from the CG simulations for commonly studied prion-like IDPs are in excellent agreement with the NMR-based α-helix fraction, demonstrating that the new HPS-SS model can accurately produce structural features of IDPs. Furthermore, this model can be easily implemented for large-scale assembly simulations due to its simplicity.

Modified Protein-Water Interactions in CHARMM36m for Thermodynamics and Kinetics of Proteins in Dilute and Crowded Solutions

Proper balance between protein-protein and protein-water interactions is vital for atomistic molecular dynamics (MD) simulations of globular proteins as well as intrinsically disordered proteins (IDPs). The overestimation of protein-protein interactions tends to make IDPs more compact than those in experiments. Likewise, multiple proteins in crowded solutions are aggregated with each other too strongly. To optimize the balance, Lennard-Jones (LJ) interactions between protein and water are often increased about 10% (with a scaling parameter, λ = 1.1) from the existing force fields. Here, we explore the optimal scaling parameter of protein-water LJ interactions for CHARMM36m in conjunction with the modified TIP3P water model, by performing enhanced sampling MD simulations of several peptides in dilute solutions and conventional MD simulations of globular proteins in dilute and crowded solutions. In our simulations, 10% increase of protein-water LJ interaction for the CHARMM36m cannot maintain stability of a small helical peptide, (AAQAA)₃ in a dilute solution and only a small modification of protein-water LJ interaction up to the 3% increase (λ = 1.03) is allowed. The modified protein-water interactions are applicable to other peptides and globular proteins in dilute solutions without changing thermodynamic properties from the original CHARMM36m. However, it has a great impact on the diffusive properties of proteins in crowded solutions, avoiding the formation of too sticky protein-protein interactions.

Low Complexity Induces Structure in Protein Regions Predicted as Intrinsically Disordered

There is increasing evidence that many intrinsically disordered regions (IDRs) in proteins play key functional roles through interactions with other proteins or nucleic acids. These interactions often exhibit a context-dependent structural behavior. We hypothesize that low complexity regions (LCRs), often found within IDRs, could have a role in inducing local structure in IDRs. To test this, we predicted IDRs in the human proteome and analyzed their structures or those of homologous sequences in the Protein Data Bank (PDB). We then identified two types of simple LCRs within IDRs: regions with only one (polyX or homorepeats) or with only two types of amino acids (polyXY). We were able to assign structural information from the PDB more often to these LCRs than to the surrounding IDRs (polyX 61.8% > polyXY 50.5% > IDRs 39.7%). The most frequently observed polyX and polyXY within IDRs contained E (Glu) or G (Gly). Structural analyses of these sequences and of homologs indicate that polyEK regions induce helical conformations, while the other most frequent LCRs induce coil structures. Our work proposes bioinformatics methods to help in the study of the structural behavior of IDRs and provides a solid basis suggesting a structuring role of LCRs within them.

Electronic Polarization at the Interface between the p53 Transactivation Domain and Two Binding Partners

Intrinsically disordered proteins (IDPs) are an abundant class of highly charged proteins that participate in numerous crucial biological processes, often in regulatory roles. IDPs do not have one major free energy minimum with a dominant structure, instead existing as conformational ensembles of multiple semistable conformations. p53 is a prototypical protein with disordered regions and binds to many structurally diverse partners, making it a useful model for exploring the role of electrostatic interactions at IDP binding interfaces. In this study, we used the Drude-2019 force field to simulate the p53 transactivation domain with two protein partners to probe the role of electrostatic interactions in IDP protein-protein interactions. We found that the Drude-2019 polarizable force field reasonably reproduced experimental chemical shifts of the p53 transactivation domain (TAD) in one complex for which these data are available. We also found that the proteins in these complexes displayed dipole response at specific residues of each protein and that residues primarily involved in binding showed a large percent change in dipole moment between the unbound and complexed states. Probing the role of electrostatic interactions in IDP binding can allow us greater fundamental understanding of these interactions and may help with targeting p53 or its partners for drug design.

Conformational Dynamics of Intrinsically Disordered Proteins Regulate Biomolecular Condensate Chemistry

Motions in biomolecules are critical for biochemical reactions. In cells, many biochemical reactions are executed inside of biomolecular condensates formed by ultradynamic intrinsically disordered proteins. A deep understanding of the conformational dynamics of intrinsically disordered proteins in biomolecular condensates is therefore of utmost importance but is complicated by diverse obstacles. Here we review emerging data on the motions of intrinsically disordered proteins inside of liquidlike condensates. We discuss how liquid–liquid phase separation modulates internal motions across a wide range of time and length scales. We further highlight the importance of intermolecular interactions that not only drive liquid–liquid phase separation but appear as key determinants for changes in biomolecular motions and the aging of condensates in human diseases. The review provides a framework for future studies to reveal the conformational dynamics of intrinsically disordered proteins in the regulation of biomolecular condensate chemistry.

Molecular interactions contributing to FUS SYGQ LC-RGG phase separation and co-partitioning with RNA polymerase II heptads

The RNA-binding protein FUS (Fused in Sarcoma) mediates phase separation in biomolecular condensates and functions in transcription by clustering with RNA polymerase II. Specific contact residues and interaction modes formed by FUS and the C-terminal heptad repeats of RNA polymerase II (CTD) have been suggested but not probed directly. Here we show how RGG domains contribute to phase separation with the FUS N-terminal low-complexity domain (SYGQ LC) and RNA polymerase II CTD. Using NMR spectroscopy and molecular simulations, we demonstrate that many residue types, not solely arginine-tyrosine pairs, form condensed-phase contacts via several interaction modes including, but not only sp2-π and cation-π interactions. In phases also containing RNA polymerase II CTD, many residue types form contacts, including both cation-π and hydrogen-bonding interactions formed by the conserved human CTD lysines. Hence, our data suggest a surprisingly broad array of residue types and modes explain co-phase separation of FUS and RNA polymerase II.

Biophysical studies of phase separation integrating experimental and computational methods

Biomolecular phase separation that contributes to the formation of membraneless organelles and biomolecular condensates has recently gained tremendous attention because of the importance of these assemblies in physiology, disease, and engineering applications. Understanding and directing biomolecular phase separation requires a multiscale view of the biophysical properties of these phases. Yet, many classic tools to characterize biomolecular properties do not apply in these condensed phases. Here, we discuss insights obtained from spectroscopic methods, in particular nuclear magnetic resonance and optical spectroscopy, in understanding the molecular and atomic interactions that underlie the formation of protein-rich condensates. We also review approaches closely coupling nuclear magnetic resonance data with computational methods especially coarse-grained and all-atom molecular simulations, which provide insight into molecular features of phase separation. Finally, we point to future methodolical developments, particularly visualizing biophysical properties of condensates in cells.

Integrating single-molecule spectroscopy and simulations for the study of intrinsically disordered proteins

Over the last two decades, intrinsically disordered proteins and protein regions (IDRs) have emerged from a niche corner of biophysics to be recognized as essential drivers of cellular function. Various techniques have provided fundamental insight into the function and dysfunction of IDRs. Among these techniques, single-molecule fluorescence spectroscopy and molecular simulations have played a major role in shaping our modern understanding of the sequence-encoded conformational behavior of disordered proteins. While both techniques are frequently used in isolation, when combined they offer synergistic and complementary information that can help uncover complex molecular details. Here we offer an overview of single-molecule fluorescence spectroscopy and molecular simulations in the context of studying disordered proteins. We discuss the various means in which simulations and single-molecule spectroscopy can be integrated, and consider a number of studies in which this integration has uncovered biological and biophysical mechanisms.

Physics-based computational and theoretical approaches to intrinsically disordered proteins

Intrinsically disordered proteins (IDPs) are an important class of proteins that do not fold to a well-defined three-dimensional shape but rather adopt an ensemble of inter-converting conformations. This feature makes their experimental characterization challenging and invites a theoretical and computational approach to complement experimental studies. In this review, we highlight the recent progress in developing new computational and theoretical approaches to study the structure and dynamics of monomeric and order higher assemblies of IDPs, with a particular emphasis on their phase separation into protein-rich condensates.

N-terminal acetylation modestly enhances phase separation and reduces aggregation of the low-complexity domain of RNA-binding protein fused in sarcoma

The RNA-binding protein fused in sarcoma (FUS) assembles via liquid-liquid phase separation (LLPS) into functional RNA granules and aggregates in amyotrophic lateral sclerosis associated neuronal inclusions. Several studies have demonstrated that posttranslational modification (PTM) can significantly alter FUS phase separation and aggregation, particularly charge-altering phosphorylation of the nearly uncharged N-terminal low complexity domain of FUS (FUS LC). However, the occurrence and impact of N-terminal acetylation on FUS phase separation remains unexplored, even though N-terminal acetylation is the most common PTM in mammals and changes the charge at the N-terminus. First, we find that FUS is predominantly acetylated in two human cell types and stress conditions. Next, we show that recombinant FUS LC can be acetylated when co-expressed with the NatA complex in Escherichia coli. Using NMR spectroscopy, we find that N-terminal acetylated FUS LC (FUS LC Nt-Ac) does not notably alter monomeric FUS LC structure or motions. Despite no difference in structure, Nt-Ac-FUS LC phase separates more avidly than unmodified FUS LC. More importantly, N-terminal acetylation of FUS LC reduces aggregation. Our findings highlight the importance of N-terminal acetylation of proteins that undergo physiological LLPS and pathological aggregation.