A theoretical method to compute sequence dependent configurational properties in charged polymers and proteins

Zwitterionic Peptides: From Mechanism, Design Strategies to Applications

Zwitterionic peptides, as a type of peptide composed of charged residues, are electrically neutral, which combine the advantages of zwitterionic materials and biological peptides, exhibiting hydrophilicity and programmable properties. As attractive candidates for resisting nonspecific adsorption of biomacromolecules and microorganisms, zwitterionic peptides have been applied in materials science, biomedicine, and biochemistry over the past decade. In this review, the development of zwitterionic peptides has been systematically outlined and analyzed, including their mechanisms, structure-function relationships, and design strategies. Furthermore, this review emphasizes and discusses their recent applications for developing functional coatings, biosensors, drug delivery systems, and engineering proteins. Finally, future research perspectives and challenges of zwitterionic peptides are also prospected and discussed. This review is intended to provide clarity and insight into the design and applications of zwitterionic peptides.

Structural transitions of a semi-flexible polyampholyte

Polyampholytes (PAs) are charged polymers composed of positively and negatively charged monomers along their backbone. The sequence of the charged monomers and the bending of the chain significantly influence the conformation and dynamical behavior of the PA. Using coarse-grained molecular dynamics simulations, we comprehensively study the structural and dynamical properties of flexible and semi-flexible PAs. The simulation results demonstrate a flexible PA chain, displaying a transition from a coil to a globule in the parameter space of the charge sequence. In addition, the behavior of the mean-square displacement (MSD), denoted as ⟨(Δr(t))2⟩, reveals distinct dynamics, specifically for the alternating and charge-segregated sequences. The MSD follows a power-law behavior, where ⟨(Δr(t))2⟩ ∼ tβ, with β ≈ 3/5 and β ≈ 1/2 for the alternating sequence and the charge-segregated sequence in the absence of hydrodynamic interactions, respectively. However, when hydrodynamic interactions are incorporated, the exponent β shifts to ∼3/5 for the charge-segregated sequence and 2/3 for the well-mixed alternating sequence. For a semi-flexible PA chain, varying the bending rigidity and electrostatic interaction strength (Γe) leads to distinct, fascinating conformational states, including globule, bundle, and torus-like conformations. We show that PAs acquire circular and hairpin-like conformations in the intermediate bending regime. The transition between various conformations is identified in terms of the shape factor estimated from the ratios of eigenvalues of the gyration tensor.

Beyond monopole electrostatics in regulating conformations of intrinsically disordered proteins

Conformations and dynamics of an intrinsically disordered protein (IDP) depend on its composition of charged and uncharged amino acids, and their specific placement in the protein sequence. In general, the charge (positive or negative) on an amino acid residue in the protein is not a fixed quantity. Each of the ionizable groups can exist in an equilibrated distribution of fully ionized state (monopole) and an ion-pair (dipole) state formed between the ionizing group and its counterion from the background electrolyte solution. The dipole formation (counterion condensation) depends on the protein conformation, which in turn depends on the distribution of charges and dipoles on the molecule. Consequently, effective charges of ionizable groups in the IDP backbone may differ from their chemical charges in isolation-a phenomenon termed charge-regulation. Accounting for the inevitable dipolar interactions, that have so far been ignored, and using a self-consistent procedure, we present a theory of charge-regulation as a function of sequence, temperature, and ionic strength. The theory quantitatively agrees with both charge reduction and salt-dependent conformation data of Prothymosin-alpha and makes several testable predictions. We predict charged groups are less ionized in sequences where opposite charges are well mixed compared to sequences where they are strongly segregated. Emergence of dipolar interactions from charge-regulation allows spontaneous coexistence of two phases having different conformations and charge states, sensitively depending on the charge patterning. These findings highlight sequence dependent charge-regulation and its potential exploitation by biological regulators such as phosphorylation and mutations in controlling protein conformation and function.

Design of intrinsically disordered protein variants with diverse structural properties

Intrinsically disordered proteins (IDPs) perform a broad range of functions in biology, suggesting that the ability to design IDPs could help expand the repertoire of proteins with novel functions. Computational design of IDPs with specific conformational properties has, however, been difficult because of their substantial dynamics and structural complexity. We describe a general algorithm for designing IDPs with specific structural properties. We demonstrate the power of the algorithm by generating variants of naturally occurring IDPs that differ in compaction, long-range contacts, and propensity to phase separate. We experimentally tested and validated our designs and analyzed the sequence features that determine conformations. We show how our results are captured by a machine learning model, enabling us to speed up the algorithm. Our work expands the toolbox for computational protein design and will facilitate the design of proteins whose functions exploit the many properties afforded by protein disorder.

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

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

Transient Non-local Interactions Dominate the Dynamics of Measles Virus NTAIL

The RNA genome of measles virus is encapsidated by the nucleoprotein within a helical nucleocapsid that serves as template for both transcription and replication. The intrinsically disordered domain of the nucleoprotein (NTAIL), partly protruding outward from the nucleocapsid, is essential for binding the polymerase complex responsible for viral transcription and replication. As for many IDPs, binding of NTAIL occurs through a short molecular recognition element (MoRE) that folds upon binding, with the majority of NTAIL remaining disordered. Though NTAIL regions far from the MoRE influence the binding affinity, interactions between them and the MoRE have not been investigated in depth. Using an integrated approach, relying on photo-induced electron transfer (PET) experiments between tryptophan and cysteine pairs placed at different positions in the protein under varying salt and pH conditions, combined with simulations and analytical models, we identified transient interactions between two disordered regions distant in sequence, which dominate NTAIL dynamics, and regulate the conformational preferences of both the MoRE and the entire NTAIL domain. Co-evolutionary analysis corroborates our findings, and suggests an important functional role for the same intramolecular interactions. We propose mechanisms by which these non-local interactions may regulate binding to the phosphoprotein, polymerase recruitment, and ultimately viral transcription and replication. Our findings may be extended to other IDPs, where non-local intra-protein interactions affect the conformational preferences of intermolecular binding sites.

Asymmetric oligomerization state and sequence patterning can tune multiphase condensate miscibility

Endogenous biomolecular condensates, composed of a multitude of proteins and RNAs, can organize into multiphasic structures with compositionally distinct phases. This multiphasic organization is generally understood to be critical for facilitating their proper biological function. However, the biophysical principles driving multiphase formation are not completely understood. Here we use in vivo condensate reconstitution experiments and coarse-grained molecular simulations to investigate how oligomerization and sequence interactions modulate multiphase organization in biomolecular condensates. We demonstrate that increasing the oligomerization state of an intrinsically disordered protein results in enhanced immiscibility and multiphase formation. Interestingly, we find that oligomerization tunes the miscibility of intrinsically disordered proteins in an asymmetric manner, with the effect being more pronounced when the intrinsically disordered protein, exhibiting stronger homotypic interactions, is oligomerized. Our findings suggest that oligomerization is a flexible biophysical mechanism that cells can exploit to tune the internal organization of biomolecular condensates and their associated biological functions.

Sequence dependence of critical properties for two-letter chains

Histogram-reweighting grand canonical Monte Carlo simulations are used to obtain the critical properties of lattice chains composed of solvophilic and solvophobic monomers. The model is a modification of one proposed by Larson et al. [J. Chem. Phys. 83, 2411 (1985)], lowering the "contrast" between beads of different types to prevent aggregation into finite-size micelles that would mask true phase separation between bulk high- and low-density phases. Oligomeric chains of lengths between 5 and 24 beads are studied. Mixed-field finite-size scaling methods are used to obtain the critical properties with typical relative accuracies of better than 10-4 for the critical temperature and 10-3 for the critical volume fraction. Diblock chains are found to have lower critical temperatures and volume fractions relative to the corresponding homopolymers. The addition of solvophilic blocks of increasing length to a fixed-length solvophobic segment results in a decrease of both the critical temperature and the critical volume fraction, with an eventual slow asymptotic approach to the long-chain limiting behavior. Moving a single solvophobic or solvophilic bead along a chain leads to a minimum or maximum in the critical temperature, with no change in the critical volume fraction. Chains of identical length and composition have a significant spread in their critical properties, depending on their precise sequence. The present study has implications for understanding biomolecular phase separation and for developing design rules for synthetic polymers with specific phase separation properties. It also provides data potentially useful for the further development of theoretical models for polymer and surfactant phase behavior.

Charge block-driven liquid-liquid phase separation - mechanism and biological roles

Liquid-liquid phase separation (LLPS) has increasingly been found to play pivotal roles in a number of intracellular events and reactions, and has introduced a new paradigm in cell biology to explain protein-protein and enzyme-ligand interactions beyond conventional molecular and biochemical theories. LLPS is driven by the cumulative effects of weak and promiscuous interactions, including electrostatic, hydrophobic and cation-π interactions, among polypeptides containing intrinsically disordered regions (IDRs) and describes the macroscopic behaviours of IDR-containing proteins in an intracellular milieu. Recent studies have revealed that interactions between 'charge blocks' - clusters of like charges along the polypeptide chain - strongly induce LLPS and play fundamental roles in its spatiotemporal regulation. Introducing a new parameter, termed 'charge blockiness', into physicochemical models of disordered polypeptides has yielded a better understanding of how the intrinsic amino acid sequence of a polypeptide determines the spatiotemporal occurrence of LLPS within a cell. Charge blockiness might also explain why some post-translational modifications segregate within IDRs and how they regulate LLPS. In this Review, we summarise recent progress towards understanding the mechanism and biological roles of charge block-driven LLPS and discuss how this new characteristic parameter of polypeptides offers new possibilities in the fields of structural biology and cell biology.

Effects of Charge Sequence Pattern and Lysine-to-Arginine Substitution on the Structural Stability of Bioinspired Polyampholytes

A comprehensive study focusing on the combined influence of the charge sequence pattern and the type of positively charged amino acids on the formation of secondary structures in sequence-specific polyampholytes is presented. The sequences of interest consisting exclusively of ionizable amino acids (lysine, K; arginine, R; and glutamic acid, E) are (EKEK)5, (EKKE)5, (ERER)5, (ERRE)5, and (EKER)5. The stability of the secondary structure was examined at three pH values in the presence of urea and NaCl. The results presented here underscore the combined prominent effects of the charge sequence pattern and the type of positively charged monomers on secondary structure formation. Additionally, (ERRE)5 readily aggregated across a wide range of pH. In contrast, sequences with the same charge pattern, (EKKE)5, as well as the sequences with the equivalent amino acid content, (ERER)5, exhibited no aggregate formation under equivalent pH and concentration conditions.

HYPK: A marginally disordered protein sensitive to charge decoration

Intrinsically disordered proteins (IDPs) that lie close to the empirical boundary separating IDPs and folded proteins in Uversky's charge-hydropathy plot may behave as "marginal IDPs" and sensitively switch conformation upon changes in environment (temperature, crowding, and charge screening), sequence, or both. In our search for such a marginal IDP, we selected Huntingtin-interacting protein K (HYPK) near that boundary as a candidate; PKIα, also near that boundary, has lower secondary structure propensity; and Crk1, just across the boundary on the folded side, has higher secondary structure propensity. We used a qualitative Förster resonance energy transfer-based assay together with circular dichroism to simultaneously probe global and local conformation. HYPK shows several unique features indicating marginality: a cooperative transition in end-to-end distance with temperature, like Crk1 and folded proteins, but unlike PKIα; enhanced secondary structure upon crowding, in contrast to Crk1 and PKIα; and a cross-over from salt-induced expansion to compaction at high temperature, likely due to a structure-to-disorder transition not seen in Crk1 and PKIα. We then tested HYPK's sensitivity to charge patterning by designing charge-flipped variants including two specific sequences with identical amino acid composition that markedly differ in their predicted size and response to salt. The experimentally observed trends, also including mutants of PKIα, verify the predictions from sequence charge decoration metrics. Marginal proteins like HYPK show features of both folded and disordered proteins that make them sensitive to physicochemical perturbations and structural control by charge patterning.

Disordered clock protein interactions and charge blocks turn an hourglass into a persistent circadian oscillator

Organismal physiology is widely regulated by the molecular circadian clock, a feedback loop composed of protein complexes whose members are enriched in intrinsically disordered regions. These regions can mediate protein-protein interactions via SLiMs, but the contribution of these disordered regions to clock protein interactions had not been elucidated. To determine the functionality of these disordered regions, we applied a synthetic peptide microarray approach to the disordered clock protein FRQ in Neurospora crassa. We identified residues required for FRQ's interaction with its partner protein FRH, the mutation of which demonstrated FRH is necessary for persistent clock oscillations but not repression of transcriptional activity. Additionally, the microarray demonstrated an enrichment of FRH binding to FRQ peptides with a net positive charge. We found that positively charged residues occurred in significant "blocks" within the amino acid sequence of FRQ and that ablation of one of these blocks affected both core clock timing and physiological clock output. Finally, we found positive charge clusters were a commonly shared molecular feature in repressive circadian clock proteins. Overall, our study suggests a mechanistic purpose for positive charge blocks and yielded insights into repressive arm protein roles in clock function.

Multiscale Computational Framework for the Liquid-Liquid Phase Separation of Intrinsically Disordered Proteins

The reversible assembly of intrinsically disordered proteins (IDPs) to form membraneless organelles (MLOs) is a fundamental process involved in the spatiotemporal regulation in living cells. MLOs formed via liquid-liquid phase separation (LLPS) serve as molecule-enhancing hubs to regulate cell functions. Owing to the complexity and dynamic nature of the protein assembly via a network of weak inter- and intra-molecular interactions, it is challenging to describe and predict the LLPS behavior. We have developed a multiscale computational model for IDPs, using the fused in sarcoma (FUS) protein and its variants as illustrative examples. To simplify the description of protein, FUS is represented as a linear chain of stickers interspaced by spacers, as inspired by the associative polymer theory. Low-complexity aromatic-rich kinked segments (LARKS) available in FUS were identified using LARKSdb and represented as "stickers". The pairwise potential energies of each pair of stickers and their β-sheet-forming propensity were estimated via molecular docking and all atomistic molecular dynamics (AA-MD) simulations. Subsequently, FUS chains were randomly positioned in a cubic lattice as coarse-grained (CG) beads, with the bead assignment based on the Kuhn length estimation of stickers and spacers. Stochastic FUS movements were modeled by Monte Carlo (MC) simulations. In addition to the Metropolis algorithm, discretized pair potential distributions between stickers were considered in the move acceptance criteria. The chosen pair potential represents one of the possible binding energy states, with its probability determined by the frequency of the binding energy distribution histogram. The fluctuations of averaged radial distribution functions (RDFs) in successive MC trial move intervals of equilibrated lattice MC simulations were used to indicate the dynamic nature of assembly/disassembly of the protein chains. This multiscale computational framework provides an economical and efficient way of predicting and describing the LLPS behavior of IDPs.

Sequence Sensitivity in Membrane Remodeling by Polyampholyte Condensates

Intrinsically disordered peptides (IDPs) have been found to undergo liquid-liquid phase separation (LLPS) and produce complex coacervates that play numerous regulatory roles in the cell. Recent experimental studies have discovered that LLPS at or near the membrane surface helps in the biomolecular organization during signaling events and can significantly alter the membrane morphology. However, the molecular mechanism and microscopic details of such processes still remain unclear. Here we study the effect of polyampholyte and polyelectrolyte condensation on two different anionic membranes, as they represent a majority of naturally occurring IDPs. The polyampholytes are fifty-residue polymers, made of glutamate(E) and lysine(K) with different charge patterns. The polyelectrolytes are separate chains of E25 and K25. We first calibrate the MARTINI v3.0 force field and then perform long-time-scale coarse-grained molecular dynamics simulations. We find that condensates formed by all the polyampholytes get adsorbed on the membrane. However, the strong polyampholytes (i.e., blocky sequences) can remodel the membranes more prominently than the weaker ones (i.e., scrambled sequences). Condensates formed by the blocky sequences induce a significant negative curvature (∼0.1 nm-1) and local demixing of lipids, whereas those by the scrambled sequences tend to wet the membrane to a greater extent without generating significant curvature or demixing. We perform several microscopic analyses to characterize the nature of the interaction between membranes and these condensates. Our analyses of interaction energetics reveal that membrane remodeling and/or wetting are favored by enhanced interactions between polyampholytes with lipids and the counterions.

Direct prediction of intrinsically disordered protein conformational properties from sequence

Intrinsically disordered regions (IDRs) are ubiquitous across all domains of life and play a range of functional roles. While folded domains are generally well described by a stable three-dimensional structure, IDRs exist in a collection of interconverting states known as an ensemble. This structural heterogeneity means that IDRs are largely absent from the Protein Data Bank, contributing to a lack of computational approaches to predict ensemble conformational properties from sequence. Here we combine rational sequence design, large-scale molecular simulations and deep learning to develop ALBATROSS, a deep-learning model for predicting ensemble dimensions of IDRs, including the radius of gyration, end-to-end distance, polymer-scaling exponent and ensemble asphericity, directly from sequences at a proteome-wide scale. ALBATROSS is lightweight, easy to use and accessible as both a locally installable software package and a point-and-click-style interface via Google Colab notebooks. We first demonstrate the applicability of our predictors by examining the generalizability of sequence-ensemble relationships in IDRs. Then, we leverage the high-throughput nature of ALBATROSS to characterize the sequence-specific biophysical behavior of IDRs within and between proteomes.

Sequence-dependent material properties of biomolecular condensates and their relation to dilute phase conformations

Material properties of phase-separated biomolecular condensates, enriched with disordered proteins, dictate many cellular functions. Contrary to the progress made in understanding the sequence-dependent phase separation of proteins, little is known about the sequence determinants of condensate material properties. Using the hydropathy scale and Martini models, we computationally decipher these relationships for charge-rich disordered protein condensates. Our computations yield dynamical, rheological, and interfacial properties of condensates that are quantitatively comparable with experimentally characterized condensates. Interestingly, we find that the material properties of model and natural proteins respond similarly to charge segregation, despite different sequence compositions. Molecular interactions within the condensates closely resemble those within the single-chain ensembles. Consequently, the material properties strongly correlate with molecular contact dynamics and single-chain structural properties. We demonstrate the potential to harness the sequence characteristics of disordered proteins for predicting and engineering the material properties of functional condensates, with insights from the dilute phase properties.

The molecular basis for cellular function of intrinsically disordered protein regions

Intrinsically disordered protein regions exist in a collection of dynamic interconverting conformations that lack a stable 3D structure. These regions are structurally heterogeneous, ubiquitous and found across all kingdoms of life. Despite the absence of a defined 3D structure, disordered regions are essential for cellular processes ranging from transcriptional control and cell signalling to subcellular organization. Through their conformational malleability and adaptability, disordered regions extend the repertoire of macromolecular interactions and are readily tunable by their structural and chemical context, making them ideal responders to regulatory cues. Recent work has led to major advances in understanding the link between protein sequence and conformational behaviour in disordered regions, yet the link between sequence and molecular function is less well defined. Here we consider the biochemical and biophysical foundations that underlie how and why disordered regions can engage in productive cellular functions, provide examples of emerging concepts and discuss how protein disorder contributes to intracellular information processing and regulation of cellular function.

A sequence-based model for identifying proteins undergoing liquid-liquid phase separation/forming fibril aggregates via machine learning

Liquid-liquid phase separation (LLPS) and the solid aggregate (also referred to as amyloid aggregates) formation of proteins, have gained significant attention in recent years due to their associations with various physiological and pathological processes in living organisms. The systematic investigation of the differences and connections between proteins undergoing LLPS and those forming amyloid fibrils at the sequence level has not yet been explored. In this research, we aim to address this gap by comparing the two types of proteins across 36 features using collected data available currently. The statistical comparison results indicate that, 24 of the selected 36 features exhibit significant difference between the two protein groups. A LLPS-Fibrils binary classification model built on these 24 features using random forest reveals that the fraction of intrinsically disordered residues (FIDR ) is identified as the most crucial feature. While, in the further three-class LLPS-Fibrils-Background classification model built on the same screened features, the composition of cysteine and that of leucine show more significant contributions than others. Through feature ablation analysis, we finally constructed a model FLFB (Feature-based LLPS-Fibrils-Background protein predictor) using six refined features, with an average area under the receiver operating characteristics of 0.83. This work indicates using sequence features and a machine learning model, proteins undergoing LLPS or forming amyloid fibrils can be identified.

Conformational ensembles of the human intrinsically disordered proteome

Intrinsically disordered proteins and regions (collectively, IDRs) are pervasive across proteomes in all kingdoms of life, help to shape biological functions and are involved in numerous diseases. IDRs populate a diverse set of transiently formed structures and defy conventional sequence-structure-function relationships1. Developments in protein science have made it possible to predict the three-dimensional structures of folded proteins at the proteome scale2. By contrast, there is a lack of knowledge about the conformational properties of IDRs, partly because the sequences of disordered proteins are poorly conserved and also because only a few of these proteins have been characterized experimentally. The inability to predict structural properties of IDRs across the proteome has limited our understanding of the functional roles of IDRs and how evolution shapes them. As a supplement to previous structural studies of individual IDRs3, we developed an efficient molecular model to generate conformational ensembles of IDRs and thereby to predict their conformational properties from sequences4,5. Here we use this model to simulate nearly all of the IDRs in the human proteome. Examining conformational ensembles of 28,058 IDRs, we show how chain compaction is correlated with cellular function and localization. We provide insights into how sequence features relate to chain compaction and, using a machine-learning model trained on our simulation data, show the conservation of conformational properties across orthologues. Our results recapitulate observations from previous studies of individual protein systems and exemplify how to link-at the proteome scale-conformational ensembles with cellular function and localization, amino acid sequence, evolutionary conservation and disease variants. Our freely available database of conformational properties will encourage further experimental investigation and enable the generation of hypotheses about the biological roles and evolution of IDRs.

Sequence-Dependent Material Properties of Biomolecular Codensates and their Relation to Dilute Phase Conformations

Material properties of phase-separated biomolecular assemblies, enriched with disordered proteins, dictate their ability to participate in many cellular functions. Despite the significant effort dedicated to understanding how the sequence of the disordered protein drives its phase separation to form condensates, little is known about the sequence determinants of condensate material properties. Here, we computationally decipher these relationships for charged disordered proteins using model sequences comprised of glutamic acid and lysine residues as well as naturally occurring sequences of LAF1's RGG domain and DDX4's N-terminal domain. We do so by delineating how the arrangement of oppositely charged residues within these sequences influences the dynamical, rheological, and interfacial properties of the condensed phase through equilibrium and non-equilibrium molecular simulations using the hydropathy scale and Martini models. Our computations yield material properties that are quantitatively comparable with experimentally characterized condensate systems. Interestingly, we find that the material properties of both the model and natural proteins respond similarly to the segregation of charges, despite their very different sequence compositions. Condensates of the highly charge-segregated sequences exhibit slower dynamics than the uniformly charge-patterned sequences, because of their comparatively long-lived molecular contacts between oppositely charged residues. Surprisingly, the molecular interactions within the condensate are highly similar to those within a single-chain for all sequences. Consequently, the condensate material properties of charged disordered proteins are strongly correlated with their dense phase contact dynamics and their single-chain structural properties. Our findings demonstrate the potential to harness the sequence characteristics of disordered proteins for predicting and engineering the material properties of functional condensates, with insights from the dilute phase properties.

Sequence-Dependent Conformational Transitions of Disordered Proteins During Condensation

Intrinsically disordered proteins (IDPs) can form biomolecular condensates through phase separation. It is recognized that the conformation of IDPs in the dense and dilute phases as well as at the interfaces of condensates can critically impact the resulting properties associated with their functionality. However, a comprehensive understanding of the conformational transitions of IDPs during condensation remains elusive. In this study, we employ a coarse-grained polyampholyte model, comprising an equal number of oppositely charged residues-glutamic acid and lysine-whereby conformations and phase behavior can be readily tuned by altering the protein sequence. By manipulating the sequence patterns from perfectly alternating to block-like, we obtain chains with ideal-like conformations to semi-compact structures in the dilute phase, while in the dense phase, the chain conformation is approximately that of an ideal chain, irrespective of the protein sequence. By performing simulations at different concentrations, we find that the chains assemble from the dilute phase through small oligomeric clusters to the dense phase, accompanied by a gradual swelling of the individual chains. We further demonstrate that these findings are applicable to several naturally occurring proteins involved in the formation of biological condensates. Concurrently, we delve deeper into the chain conformations within the condensate, revealing that chains at the interface show a strong sequence dependence, but remain more collapsed than those in the bulk-like dense phase. This study addresses critical gaps in our knowledge of IDP conformations within condensates as a function of protein sequence.

Electrostatics and hydrophobicity in the dynamics of intrinsically disordered proteins

Internal friction is a major contribution to the dynamics of intrinsically disordered proteins (IDPs). Yet, the molecular origin of internal friction has so far been elusive. Here, we investigate whether attractive electrostatic interactions in IDPs modulate internal friction differently than the hydrophobic effect. To this end, we used nanosecond fluorescence correlation spectroscopy (nsFCS) and single-molecule Förster resonance energy transfer (FRET) to quantify the conformation and dynamics of the disordered DNA-binding domains Myc, Max and Mad at different salt concentrations. We find that internal friction effects are stronger when the chain is compacted by electrostatic attractions compared to the hydrophobic effect. Although the effect is moderate, the results show that the heteropolymeric nature of IDPs is reflected in their dynamics.

Theory of the center-of-mass diffusion and viscosity of microstructured and variable sequence copolymer liquids

Biomolecular condensates formed through the phase separation of proteins and nucleic acids are widely observed, offering a fundamental means of organizing intracellular materials in a membrane-less fashion. Traditionally, these condensates have been regarded as homogeneous isotropic liquids. However, in analogy with some synthetic copolymer systems, our recent theoretical research has demonstrated that model biomolecular condensates can exhibit a microemulsion-like internal structure, contingent upon the specific sequence, inter-chain site-site interactions, and concentrated phase polymer density. Motivated by these considerations, here we present a microscopic dynamical theory for the self-diffusion constant and viscosity of a simpler class of model systems - concentrated unentangled A/B regular multiblock copolymer solutions. Our approach integrates static equilibrium local and microdomain scale structural information obtained from PRISM integral equation theory and the time evolution of the autocorrelation function of monomer scale forces at the center-of-mass level to determine the polymer diffusion constant and viscosity in a weak caging regime far from a glass or gel transition. We focus on regular multi-block systems both for simplicity and for its relevance to synthetic macromolecular science. The impact of sequence and inter-chain attraction strength on the slowing down of copolymer mass transport and flow due to local clustering enhanced collisional friction and emergent microdomain scale ordering are established. Analytic analysis and metrics employed in the study of biomolecular condensates are employed to identify key order parameters that quantity how attractive forces, packing structure, multiblock sequence, and copolymer density determine dynamical slowing down above and below the crossover to a fluctuating polymeric microemulsion state.

Conformational Properties of Polymers at Droplet Interfaces as Model Systems for Disordered Proteins

Polymer models serve as useful tools for studying the formation and physical properties of biomolecular condensates. In recent years, the interface dividing the dense and dilute phases of condensates has been discovered to be closely related to their functionality, but the conformational preferences of the constituent proteins remain unclear. To elucidate this, we perform molecular simulations of a droplet formed by phase separation of homopolymers as a surrogate model for the prion-like low-complexity domains. By systematically analyzing the polymer conformations at different locations in the droplet, we find that the chains become compact at the droplet interface compared with the droplet interior. Further, segmental analysis revealed that the end sections of the chains are enriched at the interface to maximize conformational entropy and are more expanded than the middle sections of the chains. We find that the majority of chain segments lie tangential to the droplet surface, and only the chain ends tend to align perpendicular to the interface. These trends also hold for the natural proteins FUS LC and LAF-1 RGG, which exhibit more compact chain conformations at the interface compared to the droplet interior. Our findings provide important insights into the interfacial properties of biomolecular condensates and highlight the value of using simple polymer physics models to understand the underlying mechanisms.

Design of intrinsically disordered protein variants with diverse structural properties

Intrinsically disordered proteins (IDPs) perform a wide range of functions in biology, suggesting that the ability to design IDPs could help expand the repertoire of proteins with novel functions. Designing IDPs with specific structural or functional properties has, however, been difficult, in part because determining accurate conformational ensembles of IDPs generally requires a combination of computational modelling and experiments. Motivated by recent advancements in efficient physics-based models for simulations of IDPs, we have developed a general algorithm for designing IDPs with specific structural properties. We demonstrate the power of the algorithm by generating variants of naturally occurring IDPs with different levels of compaction and that vary more than 100 fold in their propensity to undergo phase separation, even while keeping a fixed amino acid composition. We experimentally tested designs of variants of the low-complexity domain of hnRNPA1 and find high accuracy in our computational predictions, both in terms of single-chain compaction and propensity to undergo phase separation. We analyze the sequence features that determine changes in compaction and propensity to phase separate and find an overall good agreement with previous findings for naturally occurring sequences. Our general, physics-based method enables the design of disordered sequences with specified conformational properties. Our algorithm thus expands the toolbox for protein design to include also the most flexible proteins and will enable the design of proteins whose functions exploit the many properties afforded by protein disorder.

Driving forces of the complex formation between highly charged disordered proteins

Highly disordered complexes between oppositely charged intrinsically disordered proteins present a new paradigm of biomolecular interactions. Here, we investigate the driving forces of such interactions for the example of the highly positively charged linker histone H1 and its highly negatively charged chaperone, prothymosin α (ProTα). Temperature-dependent single-molecule Förster resonance energy transfer (FRET) experiments and isothermal titration calorimetry reveal ProTα-H1 binding to be enthalpically unfavorable, and salt-dependent affinity measurements suggest counterion release entropy to be an important thermodynamic driving force. Using single-molecule FRET, we also identify ternary complexes between ProTα and H1 in addition to the heterodimer at equilibrium and show how they contribute to the thermodynamics observed in ensemble experiments. Finally, we explain the observed thermodynamics quantitatively with a mean-field polyelectrolyte theory that treats counterion release explicitly. ProTα-H1 complex formation resembles the interactions between synthetic polyelectrolytes, and the underlying principles are likely to be of broad relevance for interactions between charged biomolecules in general.

Pervasive, conserved secondary structure in highly charged protein regions

Understanding how protein sequences confer function remains a defining challenge in molecular biology. Two approaches have yielded enormous insight yet are often pursued separately: structure-based, where sequence-encoded structures mediate function, and disorder-based, where sequences dictate physicochemical and dynamical properties which determine function in the absence of stable structure. Here we study highly charged protein regions (>40% charged residues), which are routinely presumed to be disordered. Using recent advances in structure prediction and experimental structures, we show that roughly 40% of these regions form well-structured helices. Features often used to predict disorder-high charge density, low hydrophobicity, low sequence complexity, and evolutionarily varying length-are also compatible with solvated, variable-length helices. We show that a simple composition classifier predicts the existence of structure far better than well-established heuristics based on charge and hydropathy. We show that helical structure is more prevalent than previously appreciated in highly charged regions of diverse proteomes and characterize the conservation of highly charged regions. Our results underscore the importance of integrating, rather than choosing between, structure- and disorder-based approaches.

Cross-Talk of Cation-π Interactions with Electrostatic and Aromatic Interactions: A Salt-Dependent Trade-off in Biomolecular Condensates

Biomolecular condensates are essential for cellular functionality, yet the complex interplay among the diverse molecular interactions that mediate their formation remains poorly understood. Here, using coarse-grained molecular dynamics simulations, we address the contribution of cation-π interactions to the stability of condensates formed via liquid-liquid phase separation. We found greater stabilization of up to 80% via cation-π interactions in condensates formed from peptides with higher aromatic residue content or less charge clustering. The contribution of cation-π interactions to droplet stability increases with increasing ionic strength, suggesting a trade-off between cation-π and electrostatic interactions. Cation-π interactions, therefore, can compensate for reduced electrostatic interactions, such as occurs at higher salt concentrations and in sequences with less charged residue content or clustering. Designing condensates with desired biophysical characteristics therefore requires quantification not only of the individual interactions but also cross-talks involving charge-charge, π-π, and cation-π interactions.

Defining the condensate landscape of fusion oncoproteins

Fusion oncoproteins (FOs) arise from chromosomal translocations in ~17% of cancers and are often oncogenic drivers. Although some FOs can promote oncogenesis by undergoing liquid-liquid phase separation (LLPS) to form aberrant biomolecular condensates, the generality of this phenomenon is unknown. We explored this question by testing 166 FOs in HeLa cells and found that 58% formed condensates. The condensate-forming FOs displayed physicochemical features distinct from those of condensate-negative FOs and segregated into distinct feature-based groups that aligned with their sub-cellular localization and biological function. Using Machine Learning, we developed a predictor of FO condensation behavior, and discovered that 67% of ~3000 additional FOs likely form condensates, with 35% of those predicted to function by altering gene expression. 47% of the predicted condensate-negative FOs were associated with cell signaling functions, suggesting a functional dichotomy between condensate-positive and -negative FOs. Our Datasets and reagents are rich resources to interrogate FO condensation in the future.

The effect of monomer polarizability on the stability and salt partitioning in model coacervates

Coacervation of charged polymer chains has been a topic of major interest in both polymer and biological sciences, as it is a subset of a phenomenon called liquid-liquid phase separation (LLPS). In this process the polymer-rich phase separates from the polymer-lean supernatant while still maintaining its liquid-like properties. LLPS has been shown to play a crucial role in cellular homeostasis by driving the formation of membraneless organelles. It also has the potential to be harnessed to aid in novel therapeutical applications. Recent studies have demonstrated that there is no one simple mechanism which drives LLPS, which is instead a result of the combined effect of electrostatic, dipolar, hydrophobic, and other weak interactions. Using coarse-grained polymer simulations we investigate the relatively unexplored effects of monomer polarizability and spatially varying dielectric constant on LLPS propensity, and these factors affect the properties of the resulting condensates. In order to produce spatial variations in the dielectric constant, all our simulations include explicit solvent and counterions. We demonstrate that polarizability has only a minor effect on the bulk behaviour of the condensates but plays a major role when ion partitioning and microstructure are considered. We observe that the major contribution comes from the nature of the neutral blocks as endowing them with an induced dipole changes their character from hydrophobic to hydrophilic. We hypothesize that the results of this work can aid in guiding future studies concerned with LLPS by providing a general framework and by highlighting important factors which influence LLPS.

ParSe 2.0: A web tool to identify drivers of protein phase separation at the proteome level

We have developed an algorithm, ParSe, which accurately identifies from the primary sequence those protein regions likely to exhibit physiological phase separation behavior. Originally, ParSe was designed to test the hypothesis that, for flexible proteins, phase separation potential is correlated to hydrodynamic size. While our results were consistent with that idea, we also found that many different descriptors could successfully differentiate between three classes of protein regions: folded, intrinsically disordered, and phase-separating intrinsically disordered. Consequently, numerous combinations of amino acid property scales can be used to make robust predictions of protein phase separation. Built from that finding, ParSe 2.0 uses an optimal set of property scales to predict domain-level organization and compute a sequence-based prediction of phase separation potential. The algorithm is fast enough to scan the whole of the human proteome in minutes on a single computer and is equally or more accurate than other published predictors in identifying proteins and regions within proteins that drive phase separation. Here, we describe a web application for ParSe 2.0 that may be accessed through a browser by visiting https://stevewhitten.github.io/Parse_v2_FASTA to quickly identify phase-separating proteins within large sequence sets, or by visiting https://stevewhitten.github.io/Parse_v2_web to evaluate individual protein sequences.

Delineating Structural Propensities of the 4E-BP2 Protein via Integrative Modeling and Clustering

The intrinsically disordered 4E-BP2 protein regulates mRNA cap-dependent translation through interaction with the predominantly folded eukaryotic initiation factor 4E (eIF4E). Phosphorylation of 4E-BP2 dramatically reduces the level of eIF4E binding, in part by stabilizing a binding-incompatible folded domain. Here, we used a Rosetta-based sampling algorithm optimized for IDRs to generate initial ensembles for two phospho forms of 4E-BP2, non- and 5-fold phosphorylated (NP and 5P, respectively), with the 5P folded domain flanked by N- and C-terminal IDRs (N-IDR and C-IDR, respectively). We then applied an integrative Bayesian approach to obtain NP and 5P conformational ensembles that agree with experimental data from nuclear magnetic resonance, small-angle X-ray scattering, and single-molecule Förster resonance energy transfer (smFRET). For the NP state, inter-residue distance scaling and 2D maps revealed the role of charge segregation and pi interactions in driving contacts between distal regions of the chain (∼70 residues apart). The 5P ensemble shows prominent contacts of the N-IDR region with the two phosphosites in the folded domain, pT37 and pT46, and, to a lesser extent, delocalized interactions with the C-IDR region. Agglomerative hierarchical clustering led to partitioning of each of the two ensembles into four clusters with different global dimensions and contact maps. This helped delineate an NP cluster that, based on our smFRET data, is compatible with the eIF4E-bound state. 5P clusters were differentiated by interactions of C-IDR with the folded domain and of the N-IDR with the two phosphosites in the folded domain. Our study provides both a better visualization of fundamental structural poses of 4E-BP2 and a set of falsifiable insights on intrachain interactions that bias folding and binding of this protein.

Conserved molecular recognition by an intrinsically disordered region in the absence of sequence conservation

Intrinsically disordered regions (IDRs) are critical for cellular function, yet often appear to lack sequence conservation when assessed by multiple sequence alignments. This raises the question of if and how function can be encoded and preserved in these regions despite massive sequence variation. To address this question, we have applied coarse-grained molecular dynamics simulations to investigate non-specific RNA binding of coronavirus nucleocapsid proteins. Coronavirus nucleocapsid proteins consist of multiple interspersed disordered and folded domains that bind RNA. We focussed here on the first two domains of coronavirus nucleocapsid proteins, the disordered N-terminal domain (NTD) followed by the folded RNA binding domain (RBD). While the NTD is highly variable across evolution, the RBD is structurally conserved. This combination makes the NTD-RBD a convenient model system to explore the interplay between an IDR adjacent to a folded domain, and how changes in IDR sequence can influence molecular recognition of a partner. Our results reveal a surprising degree of sequence-specificity encoded by both the composition and the precise order of the amino acids in the NTD. The presence of an NTD can - depending on the sequence - either suppress or enhance RNA binding. Despite this sensitivity, large-scale variation in NTD sequences is possible while certain sequence features are retained. Consequently, a conformationally-conserved fuzzy RNA:protein complex is found across nucleocapsid protein orthologs, despite large-scale changes in both NTD sequence and RBD surface chemistry. Taken together, these insights shed light on the ability of disordered regions to preserve functional characteristics despite their sequence variability.

Conformational Properties of Polymers at Droplet Interfaces as Model Systems for Disordered Proteins

Polymer models serve as useful tools for studying the formation and physical properties of biomolecular condensates. In recent years, the interface dividing the dense and dilute phases of condensates has been discovered to be closely related to their functionality, but the conformational preferences of the constituent proteins remain unclear. To elucidate this, we perform molecular simulations of a droplet formed by liquid‒liquid phase separation of homopolymers, as a surrogate model for the prion-like low-complexity domains. By systematically analyzing the polymer conformations at different locations in the droplet, we find that the chains become compact at the droplet interface compared to the droplet interior. Further, segmental analysis revealed that the end sections of the chains are enriched at the interface to maximize conformational entropy, and are more expanded than the middle sections of the chains. We find that the majority of chain segments lie tangential to the droplet surface and only the chain ends tend to align perpendicular to the interface. These trends also hold for the natural proteins FUC LC and LAF-1 RGG, which exhibit more compact chain conformations at the interface compared with the droplet interior. Our findings provide important insights into the interfacial properties of biomolecular condensates and highlight the value of using simple polymer physics models to understand the underlying mechanisms.

Theoretical and Data-Driven Approaches for Biomolecular Condensates

Biomolecular condensation processes are increasingly recognized as a fundamental mechanism that living cells use to organize biomolecules in time and space. These processes can lead to the formation of membraneless organelles that enable cells to perform distinct biochemical processes in controlled local environments, thereby supplying them with an additional degree of spatial control relative to that achieved by membrane-bound organelles. This fundamental importance of biomolecular condensation has motivated a quest to discover and understand the molecular mechanisms and determinants that drive and control this process. Within this molecular viewpoint, computational methods can provide a unique angle to studying biomolecular condensation processes by contributing the resolution and scale that are challenging to reach with experimental techniques alone. In this Review, we focus on three types of dry-lab approaches: theoretical methods, physics-driven simulations and data-driven machine learning methods. We review recent progress in using these tools for probing biomolecular condensation across all three fields and outline the key advantages and limitations of each of the approaches. We further discuss some of the key outstanding challenges that we foresee the community addressing next in order to develop a more complete picture of the molecular driving forces behind biomolecular condensation processes and their biological roles in health and disease.

Intramolecular structural heterogeneity altered by long-range contacts in an intrinsically disordered protein

Short-range interactions and long-range contacts drive the 3D folding of structured proteins. The proteins' structure has a direct impact on their biological function. However, nearly 40% of the eukaryotes proteome is composed of intrinsically disordered proteins (IDPs) and protein regions that fluctuate between ensembles of numerous conformations. Therefore, to understand their biological function, it is critical to depict how the structural ensemble statistics correlate to the IDPs' amino acid sequence. Here, using small-angle X-ray scattering and time-resolved Förster resonance energy transfer (trFRET), we study the intramolecular structural heterogeneity of the neurofilament low intrinsically disordered tail domain (NFLt). Using theoretical results of polymer physics, we find that the Flory scaling exponent of NFLt subsegments correlates linearly with their net charge, ranging from statistics of ideal to self-avoiding chains. Surprisingly, measuring the same segments in the context of the whole NFLt protein, we find that regardless of the peptide sequence, the segments' structural statistics are more expanded than when measured independently. Our findings show that while polymer physics can, to some level, relate the IDP's sequence to its ensemble conformations, long-range contacts between distant amino acids play a crucial role in determining intramolecular structures. This emphasizes the necessity of advanced polymer theories to fully describe IDPs ensembles with the hope that it will allow us to model their biological function.

Long-lived metastable knots in polyampholyte chains

Knots in proteins and DNA are known to have significant effect on their equilibrium and dynamic properties as well as on their function. While knot dynamics and thermodynamics in electrically neutral and uniformly charged polymer chains are relatively well understood, proteins are generally polyampholytes, with varied charge distributions along their backbones. Here we use simulations of knotted polymer chains to show that variation in the charge distribution on a polyampholyte chain with zero net charge leads to significant variation in the resulting knot dynamics, with some charge distributions resulting in long-lived metastable knots that escape the (open-ended) chain on a timescale that is much longer than that for knots in electrically neutral chains. The knot dynamics in such systems can be described, quantitatively, using a simple one-dimensional model where the knot undergoes biased Brownian motion along a "reaction coordinate", equal to the knot size, in the presence of a potential of mean force. In this picture, long-lived knots result from charge sequences that create large electrostatic barriers to knot escape. This model allows us to predict knot lifetimes even when those times are not directly accessible by simulations.

The molecular grammar of protein disorder guiding genome-binding locations

Intrinsically disordered regions (IDRs) direct transcription factors (TFs) towards selected genomic occurrences of their binding motif, as exemplified by budding yeast's Msn2. However, the sequence basis of IDR-directed TF binding selectivity remains unknown. To reveal this sequence grammar, we analyze the genomic localizations of >100 designed IDR mutants, each carrying up to 122 mutations within this 567-AA region. Our data points at multivalent interactions, carried by hydrophobic-mostly aliphatic-residues dispersed within a disordered environment and independent of linear sequence motifs, as the key determinants of Msn2 genomic localization. The implications of our results for the mechanistic basis of IDR-based TF binding preferences are discussed.

Theory and Simulation of Multiphase Coexistence in Biomolecular Mixtures

Biomolecular condensates constitute a newly recognized form of spatial organization in living cells. Although many condensates are believed to form as a result of phase separation, the physicochemical properties that determine the phase behavior of heterogeneous biomolecular mixtures are only beginning to be explored. Theory and simulation provide invaluable tools for probing the relationship between molecular determinants, such as protein and RNA sequences, and the emergence of phase-separated condensates in such complex environments. This review covers recent advances in the prediction and computational design of biomolecular mixtures that phase-separate into many coexisting phases. First, we review efforts to understand the phase behavior of mixtures with hundreds or thousands of species using theoretical models and statistical approaches. We then describe progress in developing analytical theories and coarse-grained simulation models to predict multiphase condensates with the molecular detail required to make contact with biophysical experiments. We conclude by summarizing the challenges ahead for modeling the inhomogeneous spatial organization of biomolecular mixtures in living cells.

The Analytical Flory Random Coil Is a Simple-to-Use Reference Model for Unfolded and Disordered Proteins

Denatured, unfolded, and intrinsically disordered proteins (collectively referred to here as unfolded proteins) can be described using analytical polymer models. These models capture various polymeric properties and can be fit to simulation results or experimental data. However, the model parameters commonly require users' decisions, making them useful for data interpretation but less clearly applicable as stand-alone reference models. Here we use all-atom simulations of polypeptides in conjunction with polymer scaling theory to parameterize an analytical model of unfolded polypeptides that behave as ideal chains (ν = 0.50). The model, which we call the analytical Flory random coil (AFRC), requires only the amino acid sequence as input and provides direct access to probability distributions of global and local conformational order parameters. The model defines a specific reference state to which experimental and computational results can be compared and normalized. As a proof-of-concept, we use the AFRC to identify sequence-specific intramolecular interactions in simulations of disordered proteins. We also use the AFRC to contextualize a curated set of 145 different radii of gyration obtained from previously published small-angle X-ray scattering experiments of disordered proteins. The AFRC is implemented as a stand-alone software package and is also available via a Google Colab notebook. In summary, the AFRC provides a simple-to-use reference polymer model that can guide intuition and aid in interpreting experimental or simulation results.

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

Assemblies of multivalent RNA-binding protein fused in sarcoma (FUS) can exist in the functional liquid-like state as well as less dynamic and potentially toxic amyloid- and hydrogel-like states. How could then cells form liquid-like condensates while avoiding their transformation to amyloids? Here, we show how posttranslational phosphorylation can provide a "handle" that prevents liquid-solid transition of intracellular condensates containing FUS. Using residue-specific coarse-grained simulations, for 85 different mammalian FUS sequences, we show how the number of phosphorylation sites and their spatial arrangement affect intracluster dynamics preventing conversion to amyloids. All atom simulations further confirm that phosphorylation can effectively reduce the β-sheet propensity in amyloid-prone fragments of FUS. A detailed evolutionary analysis shows that mammalian FUS PLDs are enriched in amyloid-prone stretches compared to control neutrally evolved sequences, suggesting that mammalian FUS proteins evolved to self-assemble. However, in stark contrast to proteins that do not phase-separate for their function, mammalian sequences have phosphosites in close proximity to these amyloid-prone regions. These results suggest that evolution uses amyloid-prone sequences in prion-like domains to enhance phase separation of condensate proteins while enriching phosphorylation sites in close proximity to safeguard against liquid-solid transitions.

Critical Assessment of Self-Consistency Checks in the All-Atom Molecular Dynamics Simulation of Intrinsically Disordered Proteins

All atom simulations can be used to quantify conformational properties of Intrinsically Disordered Proteins (IDP). However, simulations must satisfy convergence checks to ensure observables computed from simulation are reliable and reproducible. While absolute convergence is purely a theoretical concept requiring infinitely long simulation, a more practical, yet rigorous, approach is to impose Self Consistency Checks (SCCs) to gain confidence in the simulated data. Currently there is no study of SCCs in IDPs, unlike their folded counterparts. In this paper, we introduce different criteria for self-consistency checks for IDPs. Next, we impose these SCCs to critically assess the performance of different simulation protocols using the N terminal domain of HIV Integrase and the linker region of SARS-CoV-2 Nucleoprotein as two model IDPs. All simulation protocols begin with all-atom implicit solvent Monte Carlo (MC) simulation and subsequent clustering of MC generated conformations to create the representative structures of the IDPs. These representative structures serve as the initial structure for subsequent molecular dynamics (MD) runs with explicit solvent. We conclude that generating multiple short (∼3 μs) MD simulation trajectories─all starting from the most representative MC generated conformation─and merging them is the protocol of choice due to (i) its ability to satisfy multiple SCCs, (ii) consistently reproducing experimental data, and (iii) the efficiency of running independent trajectories in parallel by harnessing multiple cores available in modern GPU clusters. Running one long trajectory (greater than 20 μs) can also satisfy the first two criteria but is less desirable due to prohibitive computation time. These findings help resolve the challenge of identifying a usable starting configuration, provide an objective measure of SCC, and establish rigorous criteria to determine the minimum length (for one long simulation) or number of trajectories needed in all-atom simulation of IDPs.

Engineering synthetic biomolecular condensates

The concept of phase-separation-mediated formation of biomolecular condensates provides a new framework to understand cellular organization and cooperativity-dependent cellular functions. With growing understanding of how biological systems drive phase separation and how cellular functions are encoded by biomolecular condensates, opportunities have emerged for cellular control through engineering of synthetic biomolecular condensates. In this Review, we discuss how to construct synthetic biomolecular condensates and how they can regulate cellular functions. We first describe the fundamental principles by which biomolecular components can drive phase separation. Next, we discuss the relationship between the properties of condensates and their cellular functions, which informs the design of components to create programmable synthetic condensates. Finally, we describe recent applications of synthetic biomolecular condensates for cellular control and discuss some of the design considerations and prospective applications.

Interpreting Transient Interactions of Intrinsically Disordered Proteins

The flexible nature of intrinsically disordered proteins (IDPs) gives rise to a conformational ensemble with a diverse set of conformations. The simplest way to describe this ensemble is through a homopolymer model without any specific interactions. However, there has been growing evidence that the conformational properties of IDPs and their relevant functions can be affected by transient interactions between specific and even nonlocal pairs of amino acids. Interpreting these interactions from experimental methods, each of which is most sensitive to a different distance regime referred to as probing length, remains a challenging and unsolved problem. Here, we first show that transient interactions can be realized between short fragments of charged amino acids by generating conformational ensembles using model disordered peptides and coarse-grained simulations. Using these ensembles, we investigate how sensitive different types of experimental measurements are to the presence of transient interactions. We find methods with shorter probing lengths to be more appropriate for detecting these transient interactions, but one experimental method is not sufficient due to the existence of other weak interactions typically seen in IDPs. Finally, we develop an adjusted polymer model with an additional short-distance peak which can robustly reproduce the distance distribution function from two experimental measurements with complementary short and long probing lengths. This new model can suggest whether a homopolymer model is insufficient for describing a specific IDP and meets the challenge of quantitatively identifying specific, transient interactions from a background of nonspecific, weak interactions.

Biophysical and Integrative Characterization of Protein Intrinsic Disorder as a Prime Target for Drug Discovery

Protein intrinsic disorder is increasingly recognized for its biological and disease-driven functions. However, it represents significant challenges for biophysical studies due to its high conformational flexibility. In addressing these challenges, we highlight the complementary and distinct capabilities of a range of experimental and computational methods and further describe integrative strategies available for combining these techniques. Integrative biophysics methods provide valuable insights into the sequence–structure–function relationship of disordered proteins, setting the stage for protein intrinsic disorder to become a promising target for drug discovery. Finally, we briefly summarize recent advances in the development of new small molecule inhibitors targeting the disordered N-terminal domains of three vital transcription factors.

The analytical Flory random coil is a simple-to-use reference model for unfolded and disordered proteins

Denatured, unfolded, and intrinsically disordered proteins (collectively referred to here as unfolded proteins) can be described using analytical polymer models. These models capture various polymeric properties and can be fit to simulation results or experimental data. However, the model parameters commonly require users' decisions, making them useful for data interpretation but less clearly applicable as stand-alone reference models. Here we use all-atom simulations of polypeptides in conjunction with polymer scaling theory to parameterize an analytical model of unfolded polypeptides that behave as ideal chains (ν = 0.50). The model, which we call the analytical Flory Random Coil (AFRC), requires only the amino acid sequence as input and provides direct access to probability distributions of global and local conformational order parameters. The model defines a specific reference state to which experimental and computational results can be compared and normalized. As a proof-of-concept, we use the AFRC to identify sequence-specific intramolecular interactions in simulations of disordered proteins. We also use the AFRC to contextualize a curated set of 145 different radii of gyration obtained from previously published small-angle X-ray scattering experiments of disordered proteins. The AFRC is implemented as a stand-alone software package and is also available via a Google colab notebook. In summary, the AFRC provides a simple-to-use reference polymer model that can guide intuition and aid in interpreting experimental or simulation results.

Pervasive, conserved secondary structure in highly charged protein regions

Understanding how protein sequences confer function remains a defining challenge in molecular biology. Two approaches have yielded enormous insight yet are often pursued separately: structure-based, where sequence-encoded structures mediate function, and disorder-based, where sequences dictate physicochemical and dynamical properties which determine function in the absence of stable structure. Here we study highly charged protein regions (>40% charged residues), which are routinely presumed to be disordered. Using recent advances in structure prediction and experimental structures, we show that roughly 40% of these regions form well-structured helices. Features often used to predict disorder-high charge density, low hydrophobicity, low sequence complexity, and evolutionarily varying length-are also compatible with solvated, variable-length helices. We show that a simple composition classifier predicts the existence of structure far better than well-established heuristics based on charge and hydropathy. We show that helical structure is more prevalent than previously appreciated in highly charged regions of diverse proteomes and characterize the conservation of highly charged regions. Our results underscore the importance of integrating, rather than choosing between, structure- and disorder-based approaches.

Effects of Sequence Composition, Patterning and Hydrodynamics on the Conformation and Dynamics of Intrinsically Disordered Proteins

Intrinsically disordered proteins (IDPs) and intrinsically disordered regions (IDRs) perform diverse functions in cellular organization, transport and signaling. Unlike the well-defined structures of the classical natively folded proteins, IDPs and IDRs dynamically span large conformational and structural ensembles. This dynamic disorder impedes the study of the relationship between the amino acid sequences of the IDPs and their spatial structures and dynamics, with different experimental techniques often offering seemingly contradictory results. Although experimental and theoretical evidence indicates that some IDP properties can be understood based on their average biophysical properties and amino acid composition, other aspects of IDP function are dictated by the specifics of the amino acid sequence. We investigate the effects of several key variables on the dimensions and the dynamics of IDPs using coarse-grained polymer models. We focus on the sequence "patchiness" informed by the sequence and biophysical properties of different classes of IDPs-and in particular FG nucleoporins of the nuclear pore complex (NPC). We show that the sequence composition and patterning are well reflected in the global conformational variables such as the radius of gyration and hydrodynamic radius, while the end-to-end distance and dynamics are highly sequence-specific. We find that in good solvent conditions highly heterogeneous sequences of IDPs can be well mapped onto averaged minimal polymer models for the purpose of prediction of the IDPs dimensions and dynamic relaxation times. The coarse-grained simulations are in a good agreement with the results of atomistic MD. We discuss the implications of these results for the interpretation of the recent experimental measurements, and for the further applications of mesoscopic models of FG nucleoporins and IDPs more broadly.

Convergent behavior of extended stalk regions from staphylococcal surface proteins with widely divergent sequence patterns

Staphylococcus epidermidis and S. aureus are highly problematic bacteria in hospital settings. This stems, at least in part, from strong abilities to form biofilms on abiotic or biotic surfaces. Biofilms are well-organized multicellular aggregates of bacteria, which, when formed on indwelling medical devices, lead to infections that are difficult to treat. Cell wall-anchored (CWA) proteins are known to be important players in biofilm formation and infection. Many of these proteins have putative stalk-like regions or regions of low complexity near the cell wall-anchoring motif. Recent work demonstrated the strong propensity of the stalk region of the S. epidermidis accumulation-associated protein (Aap) to remain highly extended under solution conditions that typically induce compaction or other significant conformational changes. This behavior is consistent with the expected function of a stalk-like region that is covalently attached to the cell wall peptidoglycan and projects the adhesive domains of Aap away from the cell surface. In this study, we evaluate whether the ability to resist compaction is a common theme among stalk regions from various staphylococcal CWA proteins. Circular dichroism spectroscopy was used to examine secondary structure changes as a function of temperature and cosolvents along with sedimentation velocity analytical ultracentrifugation and SAXS to characterize structural characteristics in solution. All stalk regions tested are intrinsically disordered, lacking secondary structure beyond random coil and polyproline type II helix, and they all sample highly extended conformations. Remarkably, the Ser-Asp dipeptide repeat region of SdrC exhibited nearly identical behavior in solution when compared to the Aap Pro/Gly-rich region, despite highly divergent sequence patterns, indicating conservation of function by various distinct staphylococcal CWA protein stalk regions.

Numerical Techniques for Applications of Analytical Theories to Sequence-Dependent Phase Separations of Intrinsically Disordered Proteins

Biomolecular condensates, physically underpinned to a significant extent by liquid-liquid phase separation (LLPS), are now widely recognized by numerous experimental studies to be of fundamental biological, biomedical, and biophysical importance. In the face of experimental discoveries, analytical formulations emerged as a powerful yet tractable tool in recent theoretical investigations of the role of LLPS in the assembly and dissociation of these condensates. The pertinent LLPS often involves, though not exclusively, intrinsically disordered proteins engaging in multivalent interactions that are governed by their amino acid sequences. For researchers interested in applying these theoretical methods, here we provide a practical guide to a set of computational techniques devised for extracting sequence-dependent LLPS properties from analytical formulations. The numerical procedures covered include those for the determination of spinodal and binodal phase boundaries from a general free energy function with examples based on the random phase approximation in polymer theory, construction of tie lines for multiple-component LLPS, and field-theoretic simulation of multiple-chain heteropolymeric systems using complex Langevin dynamics. Since a more accurate physical picture often requires comparing analytical theory against explicit-chain model predictions, a commonly utilized methodology for coarse-grained molecular dynamics simulations of sequence-specific LLPS is also briefly outlined.