Polycation-π interactions are a driving force for molecular recognition by an intrinsically disordered oncoprotein family

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.

Atomistic modeling of liquid-liquid phase equilibrium explains dependence of critical temperature on γ-crystallin sequence

Liquid-liquid phase separation of protein solutions has regained heightened attention for its biological importance and pathogenic relevance. Coarse-grained models are limited when explaining residue-level effects on phase equilibrium. Here we report phase diagrams for γ-crystallins using atomistic modeling. The calculations were made possible by combining our FMAP method for computing chemical potentials and Brownian dynamics simulations for configurational sampling of dense protein solutions, yielding the binodal and critic temperature (Tc). We obtain a higher Tc for a known high-Tc γ-crystallin, γF, than for a low-Tc paralog, γB. The difference in Tc is corroborated by a gap in second virial coefficient. Decomposition of inter-protein interactions reveals one amino-acid substitution between γB and γF, from Ser to Trp at position 130, as the major contributor to the difference in Tc. This type of analysis enables us to link phase equilibrium to amino-acid sequence and to design mutations for altering phase equilibrium.

Barrier properties of Nup98 FG phases ruled by FG motif identity and inter-FG spacer length

Nup98 FG repeat domains comprise hydrophobic FG motifs linked through uncharged spacers. FG motifs capture nuclear transport receptors (NTRs) during nuclear pore complex (NPC) passage, confer inter-repeat cohesion, and condense the domains into a selective phase with NPC-typical barrier properties. We show that shortening inter-FG spacers enhances cohesion, increases phase density, and tightens such barrier - all consistent with a sieve-like phase. Phase separation tolerates mutating the Nup98-typical GLFG motifs, provided domain-hydrophobicity remains preserved. NTR-entry, however, is sensitive to (certain) deviations from canonical FG motifs, suggesting co-evolutionary adaptation. Unexpectedly, we observed that arginines promote FG-phase-entry apparently also by hydrophobic interactions/ hydrogen-bonding and not just through cation-π interactions. Although incompatible with NTR·cargo complexes, a YG phase displays remarkable transport selectivity, particularly for engineered GFP^NTR-variants. GLFG to FSFG mutations make the FG phase hypercohesive, precluding NTR-entry. Extending spacers relaxes this hypercohesion. Thus, antagonism between cohesion and NTR·FG interactions is key to transport selectivity.

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.

Analytical Formulation and Field-Theoretic Simulation of Sequence-Specific Phase Separation of Protein-Like Heteropolymers with Short- and Long-Spatial-Range Interactions

A theory for sequence-dependent liquid-liquid phase separation (LLPS) of intrinsically disordered proteins (IDPs) in the study of biomolecular condensates is formulated by extending the random phase approximation (RPA) and field-theoretic simulation (FTS) of heteropolymers with spatially long-range Coulomb interactions to include the fundamental effects of short-range, hydrophobic-like interactions between amino acid residues. To this end, short-range effects are modeled by Yukawa interactions between multiple nonelectrostatic charges derived from an eigenvalue decomposition of pairwise residue-residue contact energies. Chain excluded volume is afforded by incompressibility constraints. A mean-field approximation leads to an effective Flory-Huggins χ parameter, which, in conjunction with RPA, accounts for the contact-interaction effects of amino acid composition and the sequence-pattern effects of long-range electrostatics in IDP LLPS, whereas FTS based on the formulation provides full sequence dependence for both short- and long-range interactions. This general approach is illustrated here by applications to variants of a natural IDP in the context of several different amino-acid interaction schemes as well as a set of different model hydrophobic-polar sequences sharing the same composition. Effectiveness of the methodology is verified by coarse-grained explicit-chain molecular dynamics simulations.

A simple thermodynamic description of phase separation of Nup98 FG domains

The permeability barrier of nuclear pore complexes (NPCs) controls nucleocytoplasmic transport. It retains inert macromolecules but allows facilitated passage of nuclear transport receptors that shuttle cargoes into or out of nuclei. The barrier can be described as a condensed phase assembled from cohesive FG repeat domains, including foremost the charge-depleted FG domain of Nup98. We found that Nup98 FG domains show an LCST-type phase separation, and we provide comprehensive and orthogonal experimental datasets for a quantitative description of this behaviour. A derived thermodynamic model correlates saturation concentration with repeat number, temperature, and ionic strength. It allows estimating the enthalpy, entropy, and ΔG (0.2 kJ/mol, 0.1 k_B·T) contributions per repeat to phase separation and inter-repeat cohesion. While changing the cohesion strength strongly impacts the strictness of barrier, these numbers provide boundary conditions for in-depth modelling not only of barrier assembly but also of NPC passage.

Minimalist Design of an Intrinsically Disordered Protein-Mimicking Scaffold for an Artificial Membraneless Organelle

Liquid-liquid phase separation (LLPS) is an emerging and universal mechanism for intracellular organization, particularly, by forming membraneless organelles (MLOs) hosting intrinsically disordered proteins (IDPs) as scaffolds. Genetic engineering is generally applied to reconstruct IDPs harboring over 100 amino acid residues. Here, we report the first design of synthetic hybrids consisting of short oligopeptides of fewer than 10 residues as "stickers" and dextran as a "spacer" to recapitulate the characteristics of IDPs, as exemplified by the multivalent FUS protein. Hybrids undergo LLPS into micron-sized liquid droplets resembling LLPS in vitro and in living cells. Moreover, the droplets formed are capable of recruiting proteins and RNAs and providing a favorable environment for a biochemical reaction with highly enriched components, thereby mimicking the function of natural MLOs. This simple yet versatile model system can help elucidate the molecular interactions implicated in MLOs and pave ways to a new type of biomimetic materials.

Atomic resolution dynamics of cohesive interactions in phase-separated Nup98 FG domains

Cohesive FG domains assemble into a condensed phase forming the selective permeability barrier of nuclear pore complexes. Nanoscopic insight into fundamental cohesive interactions has long been hampered by the sequence heterogeneity of native FG domains. We overcome this challenge by utilizing an engineered perfectly repetitive sequence and a combination of solution and magic angle spinning NMR spectroscopy. We map the dynamics of cohesive interactions in both phase-separated and soluble states at atomic resolution using TROSY for rotational correlation time (TRACT) measurements. We find that FG repeats exhibit nanosecond-range rotational correlation times and remain disordered in both states, although FRAP measurements show slow translation of phase-separated FG domains. NOESY measurements enable the direct detection of contacts involved in cohesive interactions. Finally, increasing salt concentration and temperature enhance phase separation and decrease local mobility of FG repeats. This lower critical solution temperature (LCST) behaviour indicates that cohesive interactions are driven by entropy.

The permeability barrier of nuclear pores is formed by disordered and yet self-interacting FG repeat domains, whose sequence heterogeneity is a challenge for mechanistic insights. Here the authors overcome this challenge and characterize the protein’s dynamics by applying NMR techniques to an FG phase system that has been simplified to its essentials.

Matrin3: Disorder and ALS Pathogenesis

Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disorder characterized by the degeneration of both upper and lower motor neurons in the brain and spinal cord. ALS is associated with protein misfolding and inclusion formation involving RNA-binding proteins, including TAR DNA-binding protein (TDP-43) and fused in sarcoma (FUS). The 125-kDa Matrin3 is a highly conserved nuclear DNA/RNA-binding protein that is implicated in many cellular processes, including binding and stabilizing mRNA, regulating mRNA nuclear export, modulating alternative splicing, and managing chromosomal distribution. Mutations in MATR3, the gene encoding Matrin3, have been identified as causal in familial ALS (fALS). Matrin3 lacks a prion-like domain that characterizes many other ALS-associated RNA-binding proteins, including TDP-43 and FUS, however, our bioinformatics analyses and preliminary studies document that Matrin3 contains long intrinsically disordered regions that may facilitate promiscuous interactions with many proteins and may contribute to its misfolding. In addition, these disordered regions in Matrin3 undergo numerous post-translational modifications, including phosphorylation, ubiquitination and acetylation that modulate the function and misfolding of the protein. Here we discuss the disordered nature of Matrin3 and review the factors that may promote its misfolding and aggregation, two elements that might explain its role in ALS pathogenesis.

Physics-driven coarse-grained model for biomolecular phase separation with near-quantitative accuracy

Various physics- and data-driven sequence-dependent protein coarse-grained models have been developed to study biomolecular phase separation and elucidate the dominant physicochemical driving forces. Here, we present Mpipi, a multiscale coarse-grained model that describes almost quantitatively the change in protein critical temperatures as a function of amino-acid sequence. The model is parameterised from both atomistic simulations and bioinformatics data and accounts for the dominant role of π-π and hybrid cation-π/π-π interactions and the much stronger attractive contacts established by arginines than lysines. We provide a comprehensive set of benchmarks for Mpipi and seven other residue-level coarse-grained models against experimental radii of gyration and quantitative in-vitro phase diagrams; Mpipi predictions agree well with experiment on both fronts. Moreover, it can account for protein-RNA interactions, correctly predicts the multiphase behaviour of a charge-matched poly-arginine/poly-lysine/RNA system, and recapitulates experimental LLPS trends for sequence mutations on FUS, DDX4 and LAF-1 proteins.

Small-Angle X-ray Scattering Signatures of Conformational Heterogeneity and Homogeneity of Disordered Protein Ensembles

An accurate account of disordered protein conformations is of central importance to deciphering the physicochemical basis of biological functions of intrinsically disordered proteins and the folding-unfolding energetics of globular proteins. Physically, disordered ensembles of nonhomopolymeric polypeptides are expected to be heterogeneous, i.e., they should differ from those homogeneous ensembles of homopolymers that harbor an essentially unique relationship between average values of end-to-end distance R_EE and radius of gyration R_g. It was posited recently, however, that small-angle X-ray scattering (SAXS) data on conformational dimensions of disordered proteins can be rationalized almost exclusively by homopolymer ensembles. Assessing this perspective, chain-model simulations are used to evaluate the discriminatory power of SAXS-determined molecular form factors (MFFs) with regard to homogeneous versus heterogeneous ensembles. The general approach adopted here is not bound by any assumption about ensemble encodability, in that the postulated heterogeneous ensembles we evaluated are not restricted to those entailed by simple interaction schemes. Our analysis of MFFs for certain heterogeneous ensembles with more narrowly distributed R_EE and R_g indicates that while they deviate from MFFs of homogeneous ensembles, the differences can be rather small. Remarkably, some heterogeneous ensembles with asphericity and R_EE drastically different from those of homogeneous ensembles can nonetheless exhibit practically identical MFFs, demonstrating that SAXS MFFs do not afford unique characterizations of basic properties of conformational ensembles in general. In other words, the ensemble to MFF mapping is practically many-to-one and likely nonsmooth. Heteropolymeric variations of the R_EE-R_g relationship were further showcased using an analytical perturbation theory developed here for flexible heteropolymers. Ramifications of our findings for interpretation of experimental data are discussed.

Subcompartmentalization of polyampholyte species in organelle-like condensates is promoted by charge-pattern mismatch and strong excluded-volume interaction

Polyampholyte field theory and explicit-chain molecular dynamics models of sequence-specific phase separation of a system with two intrinsically disordered protein (IDP) species indicate consistently that a substantial polymer excluded volume and a significant mismatch of the IDP sequence charge patterns can act in concert, but not in isolation, to demix the two IDP species upon condensation. This finding reveals an energetic-geometric interplay in a stochastic, "fuzzy" molecular recognition mechanism that may facilitate subcompartmentalization of membraneless organelles.

Phosphorylation-dependent regulation of messenger RNA transcription, processing and translation within biomolecular condensates

Regulation of messenger RNA (mRNA) transcription, processing and translation occurs in the context of biomolecular condensates. How the physical properties of condensates connect with their biological regulatory functions is an ongoing area of interest, particularly for RNA metabolic pathways. Phosphorylation has emerged as an important mechanism for regulating protein phase separation propensities and localization patterns into different condensates, affecting compositions and dynamics. Key factors in transcription, mRNA processing and translation exhibit such phosphorylation-dependent changes in their roles within condensates, including their catalytic activities. Phosphorylation is increasingly understood to regulate the exchange of proteins through functionally linked condensates to fulfil their mRNA metabolic functions.

Molecular Details of Protein Condensates Probed by Microsecond Long Atomistic Simulations

The formation of membraneless organelles in cells commonly occurs via liquid-liquid phase separation (LLPS) and is in many cases driven by multivalent interactions between intrinsically disordered proteins (IDPs). Investigating the nature of these interactions, and their effect on dynamics within the condensed phase, is therefore of critical importance but very challenging for either simulation or experiment. Here, we study these interactions and their dynamics by pairing a novel multiscale simulation strategy with microsecond all-atom MD simulations of a condensed, IDP-rich phase. We simulate two IDPs this way, the low complexity domain of FUS and the N-terminal disordered domain of LAF-1, and find good agreement with experimental information about average density, water content, and residue-residue contacts. We go significantly beyond what is known from experiments by showing that ion partitioning within the condensed phase is largely driven by the charge distribution of the proteins and-in the cases considered-shows little evidence of preferential interactions of the ions with the proteins. Furthermore, we can probe the microscopic diffusive dynamics within the condensed phase, showing that water and ions are in dynamic equilibrium between dense and dilute phases, and their diffusion is reduced in the dense phase. Despite their high concentration in the condensate, the protein molecules also remain mobile, explaining the observed liquid-like properties of this phase. We finally show that IDP self-association is driven by a combination of nonspecific hydrophobic interactions as well as hydrogen bonds, salt bridges, and π-π and cation-π interactions. The simulation approach presented here allows the structural and dynamical properties of biomolecular condensates to be studied in microscopic detail and is generally applicable to single- and multicomponent systems of proteins and nucleic acids involved in LLPS.

Comparative roles of charge, π, and hydrophobic interactions in sequence-dependent phase separation of intrinsically disordered proteins

Endeavoring toward a transferable, predictive coarse-grained explicit-chain model for biomolecular condensates underlain by liquid-liquid phase separation (LLPS) of proteins, we conducted multiple-chain simulations of the N-terminal intrinsically disordered region (IDR) of DEAD-box helicase Ddx4, as a test case, to assess roles of electrostatic, hydrophobic, cation-π, and aromatic interactions in amino acid sequence-dependent LLPS. We evaluated three different residue-residue interaction schemes with a shared electrostatic potential. Neither a common hydrophobicity scheme nor one augmented with arginine/lysine-aromatic cation-π interactions consistently accounted for available experimental LLPS data on the wild-type, a charge-scrambled, a phenylalanine-to-alanine (FtoA), and an arginine-to-lysine (RtoK) mutant of Ddx4 IDR. In contrast, interactions based on contact statistics among folded globular protein structures reproduce the overall experimental trend, including that the RtoK mutant has a much diminished LLPS propensity. Consistency between simulation and experiment was also found for RtoK mutants of P-granule protein LAF-1, underscoring that, to a degree, important LLPS-driving π-related interactions are embodied in classical statistical potentials. Further elucidation is necessary, however, especially of phenylalanine's role in condensate assembly because experiments on FtoA and tyrosine-to-phenylalanine mutants suggest that LLPS-driving phenylalanine interactions are significantly weaker than posited by common statistical potentials. Protein-protein electrostatic interactions are modulated by relative permittivity, which in general depends on aqueous protein concentration. Analytical theory suggests that this dependence entails enhanced interprotein interactions in the condensed phase but more favorable protein-solvent interactions in the dilute phase. The opposing trends lead to only a modest overall impact on LLPS.

Unraveling Molecular Interactions in Liquid-Liquid Phase Separation of Disordered Proteins by Atomistic Simulations

Membraneless organelles are dynamical cellular condensates formed via biomolecular liquid-liquid phase separation of proteins and RNA molecules. Multiple evidence suggests that in several cases disordered proteins are structural scaffolds that drive the condensation by forming a dynamic network of inter- and intramolecular contacts. Despite the blooming research activity in this field, the structural characterization of these entities is very limited, and we still do not understand how the phase behavior is encoded in the amino acid sequences of the scaffolding proteins. Here we exploited explicit-solvent atomistic simulations to investigate the N-terminal disordered region of DEAD-box helicase 4 (NDDX4), which is a well-established model for phase separation. Notably, we determined NDDX4 conformational ensemble at the single-molecule level, and we relied on a "divide-and-conquer" strategy, based on simulations of various protein fragments at high concentration, to probe intermolecular interactions in conditions mimicking real condensates. Our results provide a high-resolution picture of the molecular mechanisms underlying phase separation in agreement with NMR and mutagenesis data and suggest that clusters of arginine and aromatic residues may stabilize the assembly of several condensates.

Analytical Theory for Sequence-Specific Binary Fuzzy Complexes of Charged Intrinsically Disordered Proteins

Intrinsically disordered proteins (IDPs) are important for biological functions. In contrast to folded proteins, molecular recognition among certain IDPs is "fuzzy" in that their binding and/or phase separation are stochastically governed by the interacting IDPs' amino acid sequences, while their assembled conformations remain largely disordered. To help elucidate a basic aspect of this fascinating yet poorly understood phenomenon, the binding of a homo or heterodimeric pair of polyampholytic IDPs is modeled statistical mechanically using cluster expansion. We find that the binding affinities of binary fuzzy complexes in the model correlate strongly with a newly derived simple "joint sequence charge decoration" parameter readily calculable from the pair of IDPs' sequence charge patterns. Predictions by our analytical theory are in essential agreement with coarse-grained explicit-chain simulations. This computationally efficient theoretical framework is expected to be broadly applicable to rationalizing and predicting sequence-specific IDP-IDP polyelectrostatic interactions.

Identifying sequence perturbations to an intrinsically disordered protein that determine its phase-separation behavior

Phase separation of intrinsically disordered proteins (IDPs) commonly underlies the formation of membraneless organelles, which compartmentalize molecules intracellularly in the absence of a lipid membrane. Identifying the protein sequence features responsible for IDP phase separation is critical for understanding physiological roles and pathological consequences of biomolecular condensation, as well as for harnessing phase separation for applications in bioinspired materials design. To expand our knowledge of sequence determinants of IDP phase separation, we characterized variants of the intrinsically disordered RGG domain from LAF-1, a model protein involved in phase separation and a key component of P granules. Based on a predictive coarse-grained IDP model, we identified a region of the RGG domain that has high contact probability and is highly conserved between species; deletion of this region significantly disrupts phase separation in vitro and in vivo. We determined the effects of charge patterning on phase behavior through sequence shuffling. We designed sequences with significantly increased phase separation propensity by shuffling the wild-type sequence, which contains well-mixed charged residues, to increase charge segregation. This result indicates the natural sequence is under negative selection to moderate this mode of interaction. We measured the contributions of tyrosine and arginine residues to phase separation experimentally through mutagenesis studies and computationally through direct interrogation of different modes of interaction using all-atom simulations. Finally, we show that despite these sequence perturbations, the RGG-derived condensates remain liquid-like. Together, these studies advance our fundamental understanding of key biophysical principles and sequence features important to phase separation.

Implications of kappa-casein evolutionary diversity for the self-assembly and aggregation of casein micelles

Milk alpha-, beta- and kappa-casein proteins assemble into casein micelles in breast epithelial cells. The glycomacropeptide (GMP) tails of kappa-casein that extend from the surface of the micelle are key to assembly and aggregation. Aggregation is triggered by stomach pepsin cleavage of GMP from para-kappa-casein (PKC). While one casein micelle model emphasizes the importance of hydrophobic interactions, another focuses on polar residues. We performed an evolutionary analysis of kappa-casein primary sequence and predicted features that potentially impact on protein interactions. We noted more rapid change in the earlier period (166 to 60 Ma). Pepsin and plasmin cleavage sites were avoided in the GMP, which may partly explain its amino acid composition. Short tandem repeats have led to modest expansions of PKC, and to large GMP expansions, suggesting the GMP is less length constrained. Amino acid compositional constraints were assessed across species. Polarity and hydrophobicity properties were insufficient to explain differences between PKC and GMP. Among polar residues, threonine dominates the GMP, compared to serine, probably reflecting its preference for O-glycosylation over phosphorylation. Glutamine, enriched in the bovine PQ-rich region, is not positionally conserved in other species. Among hydrophobic residues, isoleucine is clearly preferred over leucine in the GMP, and patches of hydrophobicity are not markedly positionally conserved. PKC tyrosine and charged residues showed stronger conservation of position, suggesting a role for pi-interactions, seen in other structurally dynamic protein membraneless assemblies. Independent acquisitions of cysteines are consistent with a trend of increasing stabilization of multimers by covalent disulphide bonds, over evolutionary time. In conclusion, kappa-casein compositional and positional constraints appear to be influenced by modification preferences, protease evasion and protein-protein interactions.

Temperature, Hydrostatic Pressure, and Osmolyte Effects on Liquid-Liquid Phase Separation in Protein Condensates: Physical Chemistry and Biological Implications

Liquid-liquid phase separation (LLPS) of proteins and other biomolecules play a critical role in the organization of extracellular materials and membrane-less compartmentalization of intra-organismal spaces through the formation of condensates. Structural properties of such mesoscopic droplet-like states were studied by spectroscopy, microscopy, and other biophysical techniques. The temperature dependence of biomolecular LLPS has been studied extensively, indicating that phase-separated condensed states of proteins can be stabilized or destabilized by increasing temperature. In contrast, the physical and biological significance of hydrostatic pressure on LLPS is less appreciated. Summarized here are recent investigations of protein LLPS under pressures up to the kbar-regime. Strikingly, for the cases studied thus far, LLPSs of both globular proteins and intrinsically disordered proteins/regions are typically more sensitive to pressure than the folding of proteins, suggesting that organisms inhabiting the deep sea and sub-seafloor sediments, under pressures up to 1 kbar and beyond, have to mitigate this pressure-sensitivity to avoid unwanted destabilization of their functional biomolecular condensates. Interestingly, we found that trimethylamine-N-oxide (TMAO), an osmolyte upregulated in deep-sea fish, can significantly stabilize protein droplets under pressure, pointing to another adaptive advantage for increased TMAO concentrations in deep-sea organisms besides the osmolyte's stabilizing effect against protein unfolding. As life on Earth might have originated in the deep sea, pressure-dependent LLPS is pertinent to questions regarding prebiotic proto-cells. Herein, we offer a conceptual framework for rationalizing the recent experimental findings and present an outline of the basic thermodynamics of temperature-, pressure-, and osmolyte-dependent LLPS as well as a molecular-level statistical mechanics picture in terms of solvent-mediated interactions and void volumes.

Intrinsically Disordered Linkers Impart Processivity on Enzymes by Spatial Confinement of Binding Domains

(1) Background: Processivity is common among enzymes and mechanochemical motors that synthesize, degrade, modify or move along polymeric substrates, such as DNA, RNA, polysaccharides or proteins. Processive enzymes can make multiple rounds of modification without releasing the substrate/partner, making their operation extremely effective and economical. The molecular mechanism of processivity is rather well understood in cases when the enzyme structurally confines the substrate, such as the DNA replication factor PCNA, and also when ATP energy is used to confine the succession of molecular events, such as with mechanochemical motors. Processivity may also result from the kinetic bias of binding imposed by spatial confinement of two binding elements connected by an intrinsically disordered (ID) linker. (2) Method: By statistical physical modeling, we show that this arrangement results in processive systems, in which the linker ensures an optimized effective concentration around novel binding site(s), favoring rebinding over full release of the polymeric partner. (3) Results: By analyzing 12 such proteins, such as cellulase, and RNAse-H, we illustrate that in these proteins linker length and flexibility, and the kinetic parameters of binding elements, are fine-tuned for optimizing processivity. We also report a conservation of structural disorder, special amino acid composition of linkers, and the correlation of their length with step size. (4) Conclusion: These observations suggest a unique type of entropic chain function of ID proteins, that may impart functional advantages on diverse enzymes in a variety of biological contexts.

Coarse-grained residue-based models of disordered protein condensates: utility and limitations of simple charge pattern parameters

Biomolecular condensates undergirded by phase separations of proteins and nucleic acids serve crucial biological functions. To gain physical insights into their genetic basis, we study how liquid-liquid phase separation (LLPS) of intrinsically disordered proteins (IDPs) depends on their sequence charge patterns using a continuum Langevin chain model wherein each amino acid residue is represented by a single bead. Charge patterns are characterized by the "blockiness" measure κ and the "sequence charge decoration" (SCD) parameter. Consistent with random phase approximation (RPA) theory and lattice simulations, LLPS propensity as characterized by critical temperature Tcr* increases with increasingly negative SCD for a set of sequences showing a positive correlation between κ and -SCD. Relative to RPA, the simulated sequence-dependent variation in Tcr* is often-though not always-smaller, whereas the simulated critical volume fractions are higher. However, for a set of sequences exhibiting an anti-correlation between κ and -SCD, the simulated Tcr*'s are quite insensitive to either parameter. Additionally, we find that blocky sequences that allow for strong electrostatic repulsion can lead to coexistence curves with upward concavity as stipulated by RPA, but the LLPS propensity of a strictly alternating charge sequence was likely overestimated by RPA and lattice models because interchain stabilization of this sequence requires spatial alignments that are difficult to achieve in real space. These results help delineate the utility and limitations of the charge pattern parameters and of RPA, pointing to further efforts necessary for rationalizing the newly observed subtleties.

Differential interaction of PRMT1 with RGG-boxes of the FET family proteins EWS and TAF15

The FET sub-family (FUS/TLS, EWS, TAF15) of RNA-binding proteins have remarkably similar overall structure but diverse biological and pathological roles. The molecular basis for FET protein specialization is largely unknown. Gly-Arg-Rich regions (RGG-boxes) within FET proteins are targets for methylation by Protein-Arginine-Methyl-Transferase-1 (PRMT1) and substrate capture is thought to involve electrostatic attraction between positively charged polyRGG substrates and negatively charged surface channels of PRMT1. Unlike FUS and EWS, a high proportion of TAF15 RGG-boxes are embedded within neutrally charged YGGDR(S/G)G repeats, suggesting that they might not bind well to PRMT1. This notion runs contrary however to a report that YGGDR(S/G)G repeats are methylated by PRMT1. Using peptide-based polyRGG substrates and a novel 2-hybrid binding assay, we find that the Asp residue in YGGDR(S/G)G repeats confers poor binding to PRMT1. Our results therefore indicate that YGGDR(S/G)G repeats may contribute to TAF15 specialization by enabling differential interactions with PRMT1 and reduced overall levels of TAF15 methylation compared with other FET proteins. By analogy with molecular recognition of other disordered polyvalent ligands by globular protein partners, we also propose a dynamic polyelectrostatic model for substrate capture by PRMT1.

Conformational Heterogeneity and FRET Data Interpretation for Dimensions of Unfolded Proteins

A mathematico-physically valid formulation is required to infer properties of disordered protein conformations from single-molecule Förster resonance energy transfer (smFRET). Conformational dimensions inferred by conventional approaches that presume a homogeneous conformational ensemble can be unphysical. When all possible-heterogeneous as well as homogeneous-conformational distributions are taken into account without prejudgment, a single value of average transfer efficiency 〈E〉 between dyes at two chain ends is generally consistent with highly diverse, multiple values of the average radius of gyration 〈R_g〉. Here we utilize unbiased conformational statistics from a coarse-grained explicit-chain model to establish a general logical framework to quantify this fundamental ambiguity in smFRET inference. As an application, we address the long-standing controversy regarding the denaturant dependence of 〈R_g〉 of unfolded proteins, focusing on Protein L as an example. Conventional smFRET inference concluded that 〈R_g〉 of unfolded Protein L is highly sensitive to [GuHCl], but data from SAXS suggested a near-constant 〈R_g〉 irrespective of [GuHCl]. Strikingly, our analysis indicates that although the reported 〈E〉 values for Protein L at [GuHCl] = 1 and 7 M are very different at 0.75 and 0.45, respectively, the Bayesian R_g² distributions consistent with these two 〈E〉 values overlap by as much as 75%. Our findings suggest, in general, that the smFRET-SAXS discrepancy regarding unfolded protein dimensions likely arise from highly heterogeneous conformational ensembles at low or zero denaturant, and that additional experimental probes are needed to ascertain the nature of this heterogeneity.

Dimension conversion and scaling of disordered protein chains

To extract protein dimension and energetics information from single-molecule fluorescence resonance energy transfer spectroscopy (smFRET) data, it is essential to establish the relationship between the distributions of the radius of gyration (Rg) and the end-to-end (donor-to-acceptor) distance (Ree). Here, we performed a coarse-grained molecular dynamics simulation to obtain a conformational ensemble of denatured proteins and intrinsically disordered proteins. For any disordered chain with fixed length, there is an excellent linear correlation between the average values of Rg and Ree under various solvent conditions, but the relationship deviates from the prediction of a Gaussian chain. A modified conversion formula was proposed to analyze smFRET data. The formula reduces the discrepancy between the results obtained from FRET and small-angle X-ray scattering (SAXS). The scaling law in a coil-globule transition process was examined where a significant finite-size effect was revealed, i.e., the scaling exponent may exceed the theoretical critical boundary [1/3, 3/5] and the prefactor changes notably during the transition. The Sanchez chain model was also tested and it was shown that the mean-field approximation works well for expanded chains.

Eukaryotic transcription factors: paradigms of protein intrinsic disorder

Gene-specific transcription factors (TFs) are key regulatory components of signaling pathways, controlling, for example, cell growth, development, and stress responses. Their biological functions are determined by their molecular structures, as exemplified by their structured DNA-binding domains targeting specific cis-acting elements in genes, and by the significant lack of fixed tertiary structure in their extensive intrinsically disordered regions. Recent research in protein intrinsic disorder (ID) has changed our understanding of transcriptional activation domains from ‘negative noodles’ to ID regions with function-related, short sequence motifs and molecular recognition features with structural propensities. This review focuses on molecular aspects of TFs, which represent paradigms of ID-related features. Through specific examples, we review how the ID-associated flexibility of TFs enables them to participate in large interactomes, how they use only a few hydrophobic residues, short sequence motifs, prestructured motifs, and coupled folding and binding for their interactions with co-activators, and how their accessibility to post-translational modification affects their interactions. It is furthermore emphasized how classic biochemical concepts like allostery, conformational selection, induced fit, and feedback regulation are undergoing a revival with the appreciation of ID. The review also describes the most recent advances based on computational simulations of ID-based interaction mechanisms and structural analysis of ID in the context of full-length TFs and suggests future directions for research in TF ID.

Opposing Intermolecular Tuning of Ca2+ Affinity for Calmodulin by Neurogranin and CaMKII Peptides

We investigated the impact of bound calmodulin (CaM)-target compound structure on the affinity of calcium (Ca²⁺) by integrating coarse-grained models and all-atomistic simulations with nonequilibrium physics. We focused on binding between CaM and two specific targets, Ca²⁺/CaM-dependent protein kinase II (CaMKII) and neurogranin (Ng), as they both regulate CaM-dependent Ca²⁺ signaling pathways in neurons. It was shown experimentally that Ca²⁺/CaM (holoCaM) binds to the CaMKII peptide with overwhelmingly higher affinity than Ca²⁺-free CaM (apoCaM); the binding of CaMKII peptide to CaM in return increases the Ca²⁺ affinity for CaM. However, this reciprocal relation was not observed in the Ng peptide (Ng_13–49), which binds to apoCaM or holoCaM with binding affinities of the same order of magnitude. Unlike the holoCaM-CaMKII peptide, whose structure can be determined by crystallography, the structural description of the apoCaM-Ng_13–49 is unknown due to low binding affinity, therefore we computationally generated an ensemble of apoCaM-Ng_13–49 structures by matching the changes in the chemical shifts of CaM upon Ng_13–49 binding from nuclear magnetic resonance experiments. Next, we computed the changes in Ca²⁺ affinity for CaM with and without binding targets in atomistic models using Jarzynski’s equality. We discovered the molecular underpinnings of lowered affinity of Ca²⁺ for CaM in the presence of Ng_13–49 by showing that the N-terminal acidic region of Ng peptide pries open the β-sheet structure between the Ca²⁺ binding loops particularly at C-domain of CaM, enabling Ca²⁺ release. In contrast, CaMKII peptide increases Ca²⁺ affinity for the C-domain of CaM by stabilizing the two Ca²⁺ binding loops. We speculate that the distinctive structural difference in the bound complexes of apoCaM-Ng_13–49 and holoCaM-CaMKII delineates the importance of CaM’s progressive mechanism of target binding on its Ca²⁺ binding affinities.

Phase Separation and Single-Chain Compactness of Charged Disordered Proteins Are Strongly Correlated

Liquid-liquid phase separation of intrinsically disordered proteins (IDPs) is a major undergirding factor in the regulated formation of membraneless organelles in the cell. The phase behavior of an IDP is sensitive to its amino acid sequence. Here we apply a recent random-phase-approximation polymer theory to investigate how the tendency for multiple chains of a protein to phase-separate, as characterized by the critical temperature T^∗_cr, is related to the protein's single-chain average radius of gyration 〈R_g〉. For a set of sequences containing different permutations of an equal number of positively and negatively charged residues, we found a striking correlation T^∗_cr ∼ 〈R_g〉^-γ with γ as large as ∼6.0, indicating that electrostatic effects have similarly significant impact on promoting single-chain conformational compactness and phase separation. Moreover, T^∗_cr ∝ -SCD, where SCD is a recently proposed "sequence charge decoration" parameter determined solely by sequence information. Ramifications of our findings for deciphering the sequence dependence of IDP phase separation are discussed.

Computational and theoretical advances in studies of intrinsically disordered proteins

Intrinsically disordered proteins (IDPs) are increasingly realized to play diverse biological roles, ranging from molecular signaling to the formation of membraneless organelles. Their high degree of disorder makes them more challenging to study using the techniques of conventional structural biology, because any observable will be averaged over a heterogeneous ensemble of structures. Molecular simulations and theory are therefore a natural complement to experiment for studying the structure, dynamics and function of IDPs. The diverse time and length scales relevant to the roles played by IDPs require flexibility in the techniques applied. Here, I summarize some of the developments in simulation and theory in recent years, which have been driven by the desire to better capture IDP properties at different time- and length-scales. I also provide an outlook for how methods can be improved in the future and emerging problems which may be addressed by theory and simulation.

An allosteric conduit facilitates dynamic multisite substrate recognition by the SCFCdc4 ubiquitin ligase

The ubiquitin ligase SCF^Cdc4 mediates phosphorylation-dependent elimination of numerous substrates by binding one or more Cdc4 phosphodegrons (CPDs). Methyl-based NMR analysis of the Cdc4 WD40 domain demonstrates that Cyclin E, Sic1 and Ash1 degrons have variable effects on the primary Cdc4^WD40 binding pocket. Unexpectedly, a Sic1-derived multi-CPD substrate (pSic1) perturbs methyls around a previously documented allosteric binding site for the chemical inhibitor SCF-I2. NMR cross-saturation experiments confirm direct contact between pSic1 and the allosteric pocket. Phosphopeptide affinity measurements reveal negative allosteric communication between the primary CPD and allosteric pockets. Mathematical modelling indicates that the allosteric pocket may enhance ultrasensitivity by tethering pSic1 to Cdc4. These results suggest negative allosteric interaction between two distinct binding pockets on the Cdc4^WD40 domain may facilitate dynamic exchange of multiple CPD sites to confer ultrasensitive dependence on substrate phosphorylation.

Structures and Short Linear Motif of Disordered Transcription Factor Regions Provide Clues to the Interactome of the Cellular Hub Protein Radical-induced Cell Death1

Intrinsically disordered protein regions (IDRs) lack a well defined three-dimensional structure but often facilitate key protein functions. Some interactions between IDRs and folded protein domains rely on short linear motifs (SLiMs). These motifs are challenging to identify, but once found they can point to larger networks of interactions, such as with proteins that serve as hubs for essential cellular functions. The stress-associated plant protein radical-induced cell death1 (RCD1) is one such hub, interacting with many transcription factors via their flexible IDRs. To identify the SLiM bound by RCD1, we analyzed the IDRs in three protein partners, DREB2A (dehydration-responsive element-binding protein 2A), ANAC013, and ANAC046, considering parameters such as disorder, context, charges, and pI. Using a combined bioinformatics and experimental approach, we have identified the bipartite RCD1-binding SLiM as (DE)X(1,2)(YF)X(1,4)(DE)L, with essential contributions from conserved aromatic, acidic, and leucine residues. Detailed thermodynamic analysis revealed both favorable and unfavorable contributions from the IDRs surrounding the SLiM to the interactions with RCD1, and the SLiM affinities ranged from low nanomolar to 50 times higher K_d values. Specifically, although the SLiM was surrounded by IDRs, individual intrinsic α-helix propensities varied as shown by CD spectroscopy. NMR spectroscopy further demonstrated that DREB2A underwent coupled folding and binding with α-helix formation upon interaction with RCD1, whereas peptides from ANAC013 and ANAC046 formed different structures or were fuzzy in the complexes. These findings allow us to present a model of the stress-associated RCD1-transcription factor interactome and to contribute to the emerging understanding of the interactions between folded hubs and their intrinsically disordered partners.

Sequence-Specific Polyampholyte Phase Separation in Membraneless Organelles

Liquid-liquid phase separation of charge- and/or aromatic-enriched intrinsically disordered proteins (IDPs) is critical in the biological function of membraneless organelles. Much of the physics of this recent discovery remains to be elucidated. Here, we present a theory in the random phase approximation to account for electrostatic effects in polyampholyte phase separations, yielding predictions consistent with recent experiments on the IDP Ddx4. The theory is applicable to any charge pattern and thus provides a general analytical framework for studying sequence dependence of IDP phase separation.

RGG boxes within the TET/FET family of RNA-binding proteins are functionally distinct

The multi-functional TET (TAF15/EWS/TLS) or FET (FUS/EWS/TLS) protein family of higher organisms harbor a transcriptional-activation domain (EAD) and an RNA-binding domain (RBD). The transcriptional activation function is, however, only revealed in oncogenic TET-fusion proteins because in native TET proteins it is auto-repressed by RGG-boxes within the TET RBD. Auto-repression is suggested to involve direct cation-pi interactions between multiple Arg residues within RGG boxes and EAD aromatics. Via analysis of TET transcriptional activity in different organisms, we report herein that repression is not autonomous but instead requires additional trans-acting factors. This finding is not supportive of a proposed model whereby repression occurs via a simple intramolecular EAD/RGG-box interaction. We also show that RGG-boxes present within reiterated YGGDRGG repeats that are unique to TAF15, are defective for repression due to the conserved Asp residue. Thus, RGG boxes within TET proteins can be functionally distinguished. While our results show that YGGDRGG repeats are not involved in TAF15 auto-repression, their remarkable number and conservation strongly suggest that they may confer specialized properties to TAF15 and thus contribute to functional differentiation within the TET/FET protein family.

Simple biophysics underpins collective conformations of the intrinsically disordered proteins of the Nuclear Pore Complex

Nuclear Pore Complexes (NPCs) are key cellular transporter that control nucleocytoplasmic transport in eukaryotic cells, but its transport mechanism is still not understood. The centerpiece of NPC transport is the assembly of intrinsically disordered polypeptides, known as FG nucleoporins, lining its passageway. Their conformations and collective dynamics during transport are difficult to assess in vivo. In vitro investigations provide partially conflicting results, lending support to different models of transport, which invoke various conformational transitions of the FG nucleoporins induced by the cargo-carrying transport proteins. We show that the spatial organization of FG nucleoporin assemblies with the transport proteins can be understood within a first principles biophysical model with a minimal number of key physical variables, such as the average protein interaction strengths and spatial densities. These results address some of the outstanding controversies and suggest how molecularly divergent NPCs in different species can perform essentially the same function.

DOI:http://dx.doi.org/10.7554/eLife.10785.001

Animal, plant and fungal cells contain a structure called the nucleus, inside which the genetic material of the cell is stored. For the cell to work properly, certain proteins and other molecules need to be able to enter and exit the nucleus. This transport is carried out by pore-like molecular “devices” known as Nuclear Pore Complexes, whose architecture and mode of operation are unique among cellular transporters.

Nuclear Pore Complexes are charged with a daunting task of deciding which of the hundreds of molecules it conducts per second should go through and which should not. Small molecules can pass freely through Nuclear Pore Complexes. However, larger molecules can only pass through the pore efficiently if they are bound to specialized transport proteins that interact with the proteins – called FG nucleoporins – that line the pore. A unique feature of the FG nucleoporins is that, unlike typical proteins, they do not have a defined three-dimensional structure. Instead, they form a soft and pliable lining inside the Nuclear Pore Complex passageway.

Exactly how interacting with transport proteins affects the structure and spatial arrangements of the FG nucleoporins in a way that allows them to control transport is not well understood. This is in part because existing experimental techniques are unable to study the structures of the FG nucleoporins in enough detail to track how they change during transport. The complexity and the diversity of the FG nucleoporins also make them difficult to model in detail.

Vovk, Gu et al. have developed a theoretical model that is based on just three basic physical properties of the FG nucleoporins – their flexibility, their ability to interact with each other, and their binding with the transport proteins. Future work can refine the model by incorporating further molecular details about the interactions between FG nucleoporins and transport proteins.

The predictions made by this simple model agree well with experimental results in a wide range of situations – from single molecules to complex spatial assemblies. They also explain why some of the experimental results appear to contradict each other and suggest how several outstanding controversies in the field can be reconciled. Because the model invokes only fundamental physical principles of FG nucleoporin assemblies, it shows that some of their general properties do not depend on the exact conditions. In particular, this might shed light on why Nuclear Pore Complexes in different organisms perform essentially the same function, although the details of their molecular structure may differ. This also suggests how the FG nucleoporins can be manipulated to build artificial devices based on the same principles.

DOI:http://dx.doi.org/10.7554/eLife.10785.002

Fuzzy complexes: Specific binding without complete folding

Specific molecular recognition is assumed to require a well-defined set of contacts and devoid of conformational and interaction ambiguities. Growing experimental evidence demonstrates however, that structural multiplicity or dynamic disorder can be retained in protein complexes, termed as fuzziness. Fuzzy regions establish alternative contacts between specific partners usually via transient interactions. Nature often tailors the dynamic properties of these segments via post-translational modifications or alternative splicing to fine-tune affinity. Most experimentally characterized fuzzy complexes are involved in regulation of gene-expression, signal transduction and cell-cycle regulation. Fuzziness is also characteristic to viral protein complexes, cytoskeleton structure, and surprisingly in a few metabolic enzymes. A plausible role of fuzzy complexes in increasing half-life of intrinsically disordered proteins is also discussed.

Biophysics of protein evolution and evolutionary protein biophysics

The study of molecular evolution at the level of protein-coding genes often entails comparing large datasets of sequences to infer their evolutionary relationships. Despite the importance of a protein's structure and conformational dynamics to its function and thus its fitness, common phylogenetic methods embody minimal biophysical knowledge of proteins. To underscore the biophysical constraints on natural selection, we survey effects of protein mutations, highlighting the physical basis for marginal stability of natural globular proteins and how requirement for kinetic stability and avoidance of misfolding and misinteractions might have affected protein evolution. The biophysical underpinnings of these effects have been addressed by models with an explicit coarse-grained spatial representation of the polypeptide chain. Sequence-structure mappings based on such models are powerful conceptual tools that rationalize mutational robustness, evolvability, epistasis, promiscuous function performed by 'hidden' conformational states, resolution of adaptive conflicts and conformational switches in the evolution from one protein fold to another. Recently, protein biophysics has been applied to derive more accurate evolutionary accounts of sequence data. Methods have also been developed to exploit sequence-based evolutionary information to predict biophysical behaviours of proteins. The success of these approaches demonstrates a deep synergy between the fields of protein biophysics and protein evolution.

Relating sequence encoded information to form and function of intrinsically disordered proteins

Intrinsically disordered proteins (IDPs) showcase the importance of conformational plasticity and heterogeneity in protein function. We summarize recent advances that connect information encoded in IDP sequences to their conformational properties and functions. We focus on insights obtained through a combination of atomistic simulations and biophysical measurements that are synthesized into a coherent framework using polymer physics theories.

Binding cavities and druggability of intrinsically disordered proteins

To assess the potential of intrinsically disordered proteins (IDPs) as drug design targets, we have analyzed the ligand-binding cavities of two datasets of IDPs (containing 37 and 16 entries, respectively) and compared their properties with those of conventional ordered (folded) proteins. IDPs were predicted to possess more binding cavity than ordered proteins at similar length, supporting the proposed advantage of IDPs economizing genome and protein resources. The cavity number has a wide distribution within each conformation ensemble for IDPs. The geometries of the cavities of IDPs differ from the cavities of ordered proteins, for example, the cavities of IDPs have larger surface areas and volumes, and are more likely to be composed of a single segment. The druggability of the cavities was examined, and the average druggable probability is estimated to be 9% for IDPs, which is almost twice that for ordered proteins (5%). Some IDPs with druggable cavities that are associated with diseases are listed. The optimism versus obstacles for drug design for IDPs is also briefly discussed.

Structure-function based molecular relationships in Ewing's sarcoma

Ewing's Sarcoma Oncogene (ews) on chromosome 22q12 is encoding a ubiquitously expressed RNA-binding protein (EWS) with unknown function that is target of tumor-specific chromosomal translocations in Ewing's sarcoma family of tumors. A model of transcription complex was proposed in which the heterodimer Rpb4/7 binds to EAD, connecting it to Core RNA Pol II. The DNA-binding domain, provided by EFP, is bound to the promoter. Rpb4/7 binds RNA, stabilizing the transcription complex. The complex Rpb4/7 can stabilize the preinitiation complexes by converting the conformation of RNA Pol II. EWS may change its conformation, so that NTD becomes accessible. Two different mechanisms of interaction between EWS and RNA Pol II are proposed: (I) an intermolecular EWS-EWS interaction between two molecules, pushing conformation from “closed” to “open” state, or (II) an intramolecular interaction inside the molecule of EWS, pushing conformation of the molecule from “closed” to “open” state. The modified forms of EWS may interact with Pol II subunits hsRpb5 and hsRpb7. The EWS and EFPs binding partners are described schematically in a model, an attempt to link the transcription with the splicing. The proposed model helps to understand the functional molecular interactions in cancer, to find new partners and ways to treat cancer.

Protein intrinsic disorder in Arabidopsis NAC transcription factors: transcriptional activation by ANAC013 and ANAC046 and their interactions with RCD1

The regulatory domains of NAC [no apical meristem, ATAF (Arabidopsis transcription activation factor), cup-shaped cotyledon] transcription factors (TFs) are mostly disordered. The single molecular recognition feature (MoRF) in ANAC046 (Arabidopsis NAC domain containing protein 46) is a functional hot spot mediating interactions with RCD1 (radical-induced cell death 1), a stress-associated hub which exploits disorder and different mechanisms for interactions.

Theoretical perspectives on nonnative interactions and intrinsic disorder in protein folding and binding

The diverse biological functions of intrinsically disordered proteins (IDPs) have markedly raised our appreciation of protein conformational versatility, whereas the existence of energetically favorable yet functional detrimental nonnative interactions underscores the physical limitations of evolutionary optimization. Here we survey recent advances in using biophysical modeling to gain insight into experimentally observed nonnative behaviors and IDP properties. Simulations of IDP interactions to date focus mostly on coupled folding-binding, which follows essentially the same organizing principle as the local-nonlocal coupling mechanism in cooperative folding of monomeric globular proteins. By contrast, more innovative theories of electrostatic and aromatic interactions are needed for the conceptually novel but less-explored 'fuzzy' complexes in which the functionally bound IDPs remain largely disordered.

Synergy of aromatic residues and phosphoserines within the intrinsically disordered DNA-binding inhibitory elements of the Ets-1 transcription factor

The E26 transformation-specific (Ets-1) transcription factor is autoinhibited by a conformationally disordered serine-rich region (SRR) that transiently interacts with its DNA-binding ETS domain. In response to calcium signaling, autoinhibition is reinforced by calmodulin-dependent kinase II phosphorylation of serines within the SRR. Using mutagenesis and quantitative DNA-binding measurements, we demonstrate that phosphorylation-enhanced autoinhibition requires the presence of phenylalanine or tyrosine (ϕ) residues adjacent to the SRR phosphoacceptor serines. The introduction of additional phosphorylated Ser-ϕ-Asp, but not Ser-Ala-Asp, repeats within the SRR dramatically reinforces autoinhibition. NMR spectroscopic studies of phosphorylated and mutated SRR variants, both within their native context and as separate trans-acting peptides, confirmed that the aromatic residues and phosphoserines contribute to the formation of a dynamic complex with the ETS domain. Complementary NMR studies also identified the SRR-interacting surface of the ETS domain, which encompasses its positively charged DNA-recognition interface and an adjacent region of neutral polar and nonpolar residues. Collectively, these studies highlight the role of aromatic residues and their synergy with phosphoserines in an intrinsically disordered regulatory sequence that integrates cellular signaling and gene expression.

Digested disorder: Quarterly intrinsic disorder digest (July-August-September, 2013)

The current literature on intrinsically disordered proteins grows fast. To keep interested readers up to speed with this literature, we continue a "Digested Disorder" project and represent a new issue of reader's digest of the research papers and reviews on intrinsically disordered proteins. The only 2 criteria for inclusion in this digest are the publication date (a paper should be published within the covered time frame) and topic (a paper should be dedicated to any aspect of protein intrinsic disorder). The current digest issue covers papers published during the third quarter of 2013; i.e., during the period of June, July, and September of 2013. Similar to previous issues, the papers are grouped hierarchically by topics they cover, and for each of the included paper a short description is given on its major findings.

A simulated intermediate state for folding and aggregation provides insights into ΔN6 β2-microglobulin amyloidogenic behavior

A major component of ex vivo amyloid plaques of patients with dialysis-related amyloidosis (DRA) is a cleaved variant of β₂-microglobulin (ΔN6) lacking the first six N-terminal residues. Here we perform a computational study on ΔN6, which provides clues to understand the amyloidogenicity of the full-length β₂-microglobulin. Contrary to the wild-type form, ΔN6 is able to efficiently nucleate fibrillogenesis in vitro at physiological pH. This behavior is enhanced by a mild acidification of the medium such as that occurring in the synovial fluid of DRA patients. Results reported in this work, based on molecular simulations, indicate that deletion of the N-terminal hexapeptide triggers the formation of an intermediate state for folding and aggregation with an unstructured strand A and a native-like core. Strand A plays a pivotal role in aggregation by acting as a sticky hook in dimer assembly. This study further predicts that the detachment of strand A from the core is maximized at pH 6.2 resulting into higher aggregation efficiency. The structural mapping of the dimerization interface suggests that Tyr10, His13, Phe30 and His84 are hot-spot residues in ΔN6 amyloidogenesis.

Dialysis-related amyloidosis (DRA) is a conformational disease that affects individuals undergoing long-term haemodialysis. In DRA the progressive accumulation of protein human β₂-microglobulin (Hβ₂m) in the osteoarticular system, followed by its assembly into amyloid fibrils, eventually leads to tissue erosion and destruction. Disclosing the aggregation mechanism of Hβ₂m under physiologically relevant conditions represents a major challenge due to the inability of the protein to efficiently nucleate fibrillogenesis in vitro at physiological pH. On the other hand, ΔN6, a truncated variant of Hβ₂m, which is also a major component of ex vivo amyloid deposits extracted from DRA patients, is able to efficiently form amyloid fibrils de novo in physiological conditions. This amyloidogenic behavior is dramatically enhanced in a slightly more acidic pH (6.2) compatible with the mild acidification that occurs in the synovial fluid of DRA patients. In this work, an innovative three-stage methodological approach, relying on an array of molecular simulations, spanning different levels of resolution is used to investigate the initial stage of the de novo aggregation mechanism of ΔN6 in a physiologically relevant pH range. We identify an intermediate state for folding and aggregation, whose potential to dimerize is enhanced at pH 6.2. Our results provide rationalizations for previous experimental observations and new insights into the molecular basis of DRA.

Recognition of the HIV capsid by the TRIM5α restriction factor is mediated by a subset of pre-existing conformations of the TRIM5α SPRY domain

The binding of the TRIM5α restriction factor to the HIV capsid is mediated by the C-terminal SPRY domain of TRIM5α. Atomic-level details of this host-pathogen interaction, which involves mobile variable loops of the SPRY domain, remain unclear. Some of the key determinants of restriction are encompassed by the long and disordered v1 loop of the SPRY domain. We applied molecular modeling to elucidate the conformational repertoire of the v1 loop and its role in the interaction with the capsid. All-atom replica exchange molecular dynamics revealed multiple transient, interconverting states of the v1 loop consistent with the intrinsic disorder observed experimentally. The docking of the SPRY conformations representing 10 most populated states onto the high-resolution model of the assembled HIV-1 capsid revealed that a subset of v1 conformations produced plausible binding poses, in which the SPRY domain binds close to the pseudo-2-fold symmetry axis and the v1 loop spans the interhexamer gap. Such binding mode is well supported by the NMR binding data and known escape mutants. We speculate that the binding mode that involves interaction of the capsid with a subset of preexisting SPRY conformations arising from the intrinsic disorder of the v1 loop may explain the remarkable ability of TRIM5α to resist viral evasion by mutagenesis and to restrict divergent retroviruses.

Effects of desolvation barriers and sidechains on local–nonlocal coupling and chevron behaviors in coarse-grained models of protein folding

Coarse-grained protein chain models with desolvation barriers or sidechains lead to stronger local–nonlocal coupling and more linear chevron plots.