Temperature-dependent solvation modulates the dimensions of disordered proteins

Temperature-Dependent Coarse-Grained Model for Simulations of Intrinsically Disordered Protein LCST and UCST Liquid-Liquid Phase Separations

Many intrinsically disordered proteins (IDPs) can undergo a liquid-liquid phase separation (LLPS) in water, depending on solution conditions (temperature, pH, and ionic strength). There are two types of LLPS that are controlled by temperature: those occurring above a lower critical solution temperature (LCST) and those occurring below an upper critical solution temperature (UCST). IDP coarse-grained (CG) models are particularly appropriate for investigating the physical and chemical factors that govern their LLPS and supramolecular organization. However, the development of CG models allowing simulations of both LCST and UCST behavior of temperature-sensitive IDPs is still in its infancy. In this context, we present here a novel temperature-dependent (TD) CG model for IDP simulations based on the MARTINI 3 force field. The model was developed by modifying the Lennard-Jones potentials between apolar or charged solute beads and water with a TD rescaling factor. It was parametrized to fit the TD potentials of mean force (PMF) between two apolar or two charged molecules computed using all-atom (AA) simulations. We show that the TD CG model is able to reproduce the experimentally known LLPS of both LCST and UCST low-complexity sequences and to estimate phase transition temperatures comparable to experimental measurements.

Collapse and Protein Folding: Should We Be Surprised That Biothermodynamics Works So Well?

A complete understanding of protein function and dynamics requires the characterization of the multiple thermodynamic states, including the denatured state ensemble (DSE). Whereas residual structure in the DSE (as well as in partially folded states) is pertinent in many biological contexts, here we are interested in how such structure affects protein thermodynamics. We examine issues related to chain collapse in light of new developments, focusing on potential complications arising from differences in the DSE's properties under various conditions. Despite some variability in the degree of collapse and structure in the DSE, stability measurements are remarkably consistent between two standard methods, calorimetry and chemical denaturation, as well as with hydrogen-deuterium exchange. This robustness is due in part to the DSEs obtained with different perturbations being thermodynamically equivalent and hence able to serve as a common reference state. An examination of the properties of the DSE points to it as being a highly expanded ensemble with minimal amounts of stable hydrogen bonded structure. These two features are likely to be critical in the broad and successful application of thermodynamics to protein folding. Our review concludes with a discussion of the impact of these findings on folding mechanisms and pathways.

Structural Context Modulates the Conformational Ensemble of the Intrinsically Disordered Amino Terminus of α-Synuclein

Regions of intrinsic disorder play crucial roles in biological systems, yet they often elude characterization by conventional biophysical techniques. To capture conformational distributions across different time scales, we employed a freezing approach coupled with solid-state NMR analysis. Using segmentally isotopically labeled α-synuclein (α-syn), we investigated the conformational ensembles of the six alanines, three glycines, and a single site (L8) in the disordered amino terminus under three distinct conditions: in 8 M urea, as a frozen monomer in buffer, and within the disordered regions flanking the amyloid core. The experimental spectra varied significantly among these conditions and deviated from those of a statistical coil. In 8 M urea, monomeric α-syn exhibited the most restricted conformational sampling, rarely accessing chemical shifts characteristic of α-helices or β-strands. In buffer, monomeric α-syn showed a broader conformational sampling, favoring α-helical conformations and, to a lesser extent, random coil states. Notably, amino acids in the disordered regions flanking the amyloid core demonstrated the most extensive conformational sampling, with broad peaks encompassing the entire range of possible chemical shifts and a marked increase in highly extended β-strand conformations. Collectively, this work demonstrates that intrinsically disordered regions exhibit distinct conformational ensembles, which are influenced not only by the chemical environment but also by the conformations of adjacent protein sequences. The differences in the conformational ensembles of the disordered amino terminus may explain why the monomer and the amyloid form of α-syn interact with different biomolecules inside cells.

Accurate prediction of thermoresponsive phase behavior of disordered proteins

Protein responses to environmental stress, particularly temperature fluctuations, have long been a subject of investigation, with a focus on how proteins maintain homeostasis and exhibit thermoresponsive properties. While UCST-type (upper critical solution temperature) phase behavior has been studied extensively and can now be predicted reliably using computational models, LCST-type (lower critical solution temperature) phase transitions remain less explored, with a lack of computational models capable of accurate prediction. This gap limits our ability to probe fully how proteins undergo phase transitions in response to temperature changes. Here, we introduce Mpipi-T, a residue-level coarse-grained model designed to predict LCST-type phase behavior of proteins. Parametrized using both atomistic simulations and experimental data, Mpipi-T accounts for entropically driven protein phase separation that occurs upon heating. Accordingly, Mpipi-T predicts temperature-driven protein behavior quantitatively in both single- and multi-chain systems. Beyond its predictive capabilities, we demonstrate that Mpipi-T provides a framework for uncovering the molecular mechanisms underlying heat stress responses, offering new insights into how proteins sense and adapt to thermal changes in biological systems.

Optimal scaling of protein-water interactions coupled with targeted torsional refinements yields balanced force fields suitable for simulations of single-chain folded proteins, disordered polypeptides, and protein-protein complexes

All-atom molecular dynamics (MD) simulations based on physics-based force fields, serve as an essential complement to experiments for investigating protein structure, dynamics, and interactions. Despite significant advances in force field development, achieving a consistent balance of molecular interactions that stabilize folded proteins and protein-protein complexes while simultaneously capturing the conformational dynamics of intrinsically disordered polypeptides (IDPs), remains challenging. In this work, we systematically evaluated two current state-of-the-art force fields (i) AMBER ff03ws, and (ii) AMBER ff99SBws, by comprehensively assessing their performance on both folded domains and IDPs. By selectively scaling side chain-water interactions for uncharged residues, the refined AMBER ff03w-sc force field demonstrated improved conformational stability of folded proteins while maintaining accurate representations of IDPs. However, AMBER ff03w-sc failed to correct the discrepancies in NMR-derived ps-ns timescale backbone dynamics associated with flexible loops. Interestingly, AMBER ff99SBws retained its structural stability despite the application of upscaled interactions with water for both sidechain and backbone atoms and displayed robust agreement with NMR-derived backbone dynamics. Further, a targeted refinement of glutamine backbone torsion parameters, yielded AMBER ff99SBws-STQ', which effectively resolved discrepancies associated with glutamine α-helicity predictions. Extensive validation against small angle X-ray scattering (SAXS) and NMR chemical shifts, revealed that both refined force fields accurately reproduced chain dimensions and secondary structure propensities of disordered peptides and prion-like domains. Importantly, both force fields reliably maintained the stability of protein-protein complexes over microsecond timescales. Our systematic refinement strategies provide improved accuracy and transferability for simulating diverse protein systems, from folded domains to IDPs and protein complexes.

Revealing an origin of temperature-dependent structural change in intrinsically disordered proteins

Intrinsically disordered proteins (IDPs) show structural changes stimulated by changes in external conditions. This study aims to reveal the temperature dependence of the structure and the dynamics of the intrinsically disordered region of the helicase-associated endonuclease for fork-structured DNA, one of the typical IDPs, using an integrative approach. Small-angle X-ray scattering (SAXS) and circular dichroism (CD) studies revealed that the radius of gyration and ellipticity at 222 nm remained constant up to 313-323 K, followed by a decline above this temperature range. NMR studies revealed the absence of a promotion of the α helix. As a result, SAXS, CD, and NMR data strongly suggest that these temperature-dependent structural changes were primarily due to a reduction in the content of the polyproline II (PPII) helix. Moreover, quasielastic neutron scattering studies revealed a slight change in the activation energy in a similar temperature range. Considering the concept of glass transition, it is posited that dynamical cooperativity between the PPII helix and water may play a significant role in these structural changes. The findings suggest that internal dynamics are crucial for regulating the structure of IDPs, highlighting the importance of considering dynamical cooperativity in future studies of protein behavior under varying temperature conditions.

Intrinsic stiffness and Θ-solvent regime in intrinsically disordered proteins: Implications for liquid-liquid phase separation

Liquid-liquid phase separation (LLPS) exhibited by intrinsically disordered proteins (IDPs) depends on the solvation state around the Θ-regime, which separates good from poor solvent. Experimentally, the Θ-solvent regime of the finite length (N) IDPs, as probed by small angle X-ray scattering (SAXS) and single molecular fluorescence resonance energy transfer (smFRET), is in disagreement. Using computer simulations of a coarse-grained IDP model, we address the effect of chain length on the Θ-regime of IDPs with polar side chains (polyglutamine) and hydrophobic side chains (polyleucine) subject to varying concentrations of cosolvents [ C ] , urea (denaturant) or trimethylamine N-oxide (protective osmolyte) in water. Due to their intrinsic stiffness, these IDPs are always expanded on short-length scales, independent of the solvent quality. As a result, for short IDP sequences ( ≈ 10 to 25 residues), their propensity to exhibit LLPS cannot be inferred from single-chain properties. Further, for finite-size IDPs, the cosolvent concentration to attain the Θ-regime ( [ C Θ ] ) extracted from the structure factor emulating SAXS and pair distances mimicking smFRET differs. They converge to the same cosolvent concentration only at large N, indicating that finite size corrections vary for different IDP properties. We show that the radius of gyration ( R g ) of the IDPs in the Θ-solvent regime satisfies the scaling relation R g 2 = N f ( c N ) , which can be exploited to accurately extract [ C Θ ] ( c = ( [ C ] / [ C Θ ] - 1 ) ). We demonstrate the importance of finite size aspects originating from the chain stiffness and thermal blob size in analyzing IDP properties to identify the Θ-solvent regime.

Structural and Functional Relevance of Charge Based Transient Interactions inside Intrinsically Disordered Proteins

Intrinsically disordered proteins (IDPs) perform a wide range of biological functions without adopting stable, well-defined, three-dimensional structures. Instead, IDPs exist as dynamic ensembles of flexible conformations, traditionally thought to be governed by weak, nonspecific interactions, which are well described by homopolymer theory. However, recent research highlights the presence of transient, specific interactions in several IDPs, suggesting that factors beyond overall size influence their conformational behavior. In this study, we investigate how the spatial arrangement of charged amino acids within IDP sequences shapes the prevalence of transient, specific interactions. Through a series of model peptides, we establish a quantitative empirical relationship between the fraction of transient interactions and a novel sequence metric, termed effective charged patch length, which characterizes the ability of charged patches to drive these interactions. By examining IDP ensembles with varying levels of transient interactions, we further explore their heteropolymeric structural behavior in phase-separated condensates, where we observe the formation of a condensate-spanning network structure. Additionally, we perform a proteome-wide scan for charge-based transient interactions within disordered regions of the human proteome, revealing that approximately 10% of these regions exhibit such charge-driven transient interactions, leading to heteropolymeric behaviors in their conformational ensembles. Finally, we examine how these charge-based transient interactions correlate with molecular functions, identifying specific biological roles in which these interactions are enriched.

Structural dynamics of the intrinsically disordered linker region of cardiac troponin T

The cardiac troponin complex, composed of troponins I, T, and C, plays a central role in regulating the calcium-dependent interactions between myosin and the thin filament. Mutations in troponin can cause cardiomyopathies; however, it is still a major challenge to connect how changes in sequence affect troponin's function. Recent high-resolution structures of the thin filament revealed critical insights into the structure-function relationship of troponin, but there remain large, unresolved segments of troponin, including the troponin-T linker region that is a hotspot for cardiomyopathy mutations. This linker region is predicted to be intrinsically disordered, with behaviors that are not well described by traditional structural approaches; however, this proposal has not been experimentally verified. Here, we used a combination of single-molecule Förster resonance energy transfer (FRET), molecular dynamics simulations, and functional reconstitution assays to investigate the troponin-T linker region. We show that in the context of both isolated troponin and the fully regulated troponin complex, the linker behaves as a dynamic, intrinsically disordered region. This region undergoes polyampholyte expansion in the presence of high salt and distinct conformational changes during the assembly of the troponin complex. We also examine the ΔE160 hypertrophic cardiomyopathy mutation in the linker and demonstrate that it does not affect the conformational dynamics of the linker, rather it allosterically affects interactions with other troponin complex subunits, leading to increased molecular contractility. Taken together, our data clearly demonstrate the importance of disorder within the troponin-T linker and provide new insights into the molecular mechanisms driving the pathogenesis of cardiomyopathies.

The dimerization domain of SARS CoV 2 Nucleocapsid protein is partially disordered as a monomer and forms a high affinity dynamic complex

The SARS-CoV-2 Nucleocapsid (N) is a 419 amino acids protein that drives the compaction and packaging of the viral genome. This compaction is aided not only by protein-RNA interactions, but also by protein-protein interactions that contribute to increasing the valence of the nucleocapsid protein. Here, we focused on quantifying the mechanisms that control dimer formation. Single-molecule Förster Resonance Energy Transfer enabled us to investigate the conformations of the dimerization domain in the context of the full-length protein as well as the energetics associated with dimerization. Under monomeric conditions, we observed significantly expanded configurations of the dimerization domain (compared to the folded dimer structure), which are consistent with a dynamic conformational ensemble. The addition of unlabeled protein stabilizes a folded dimer configuration with a high mean transfer efficiency, in agreement with predictions based on known structures. Dimerization is characterized by a dissociation constant of ~ 12 nM at 23 °C and is driven by strong enthalpic interactions between the two protein subunits, which originate from the coupled folding and binding. Interestingly, the dimer structure retains some of the conformational heterogeneity of the monomeric units, and the addition of denaturant reveals that the dimer domain can significantly expand before being completely destabilized. Our findings suggest that the inherent flexibility of the monomer form is required to adopt the specific fold of the dimer domain, where the two subunits interlock with one another. We proposed that the retained flexibility of the dimer form may favor the capture and interactions with RNA, and that the temperature dependence of dimerization may explain some of the previous observations regarding the phase separation propensity of the N protein.

Direct prediction of intermolecular interactions driven by disordered regions

Intrinsically disordered regions (IDRs) are critical for a wide variety of cellular functions, many of which involve interactions with partner proteins. Molecular recognition is typically considered through the lens of sequence-specific binding events. However, a growing body of work has shown that IDRs often interact with partners in a manner that does not depend on the precise order of the amino acid order, instead driven by complementary chemical interactions leading to disordered bound-state complexes. Despite this emerging paradigm, we lack tools to describe, quantify, predict, and interpret these types of structurally heterogeneous interactions from the underlying amino acid sequences. Here, we repurpose the chemical physics developed originally for molecular simulations to develop an approach for predicting intermolecular interactions between IDRs and partner proteins. Our approach enables the direct prediction of phase diagrams, the identification of chemically-specific interaction hotspots on IDRs, and a route to develop and test mechanistic hypotheses regarding IDR function in the context of molecular recognition. We use our approach to examine a range of systems and questions to highlight its versatility and applicability.

How hydrophobicity, side chains, and salt affect the dimensions of disordered proteins

Despite the generally accepted role of the hydrophobic effect as the driving force for folding, many intrinsically disordered proteins (IDPs), including those with hydrophobic content typical of foldable proteins, behave nearly as self-avoiding random walks (SARWs) under physiological conditions. Here, we tested how temperature and ionic conditions influence the dimensions of the N-terminal domain of pertactin (PNt), an IDP with an amino acid composition typical of folded proteins. While PNt contracts somewhat with temperature, it nevertheless remains expanded over 10-58°C, with a Flory exponent, ν, >0.50. Both low and high ionic strength also produce contraction in PNt, but this contraction is mitigated by reducing charge segregation. With 46% glycine and low hydrophobicity, the reduced form of snow flea anti-freeze protein (red-sfAFP) is unaffected by temperature and ionic strength and persists as a near-SARW, ν ~ 0.54, arguing that the thermal contraction of PNt is due to stronger interactions between hydrophobic side chains. Additionally, red-sfAFP is a proxy for the polypeptide backbone, which has been thought to collapse in water. Increasing the glycine segregation in red-sfAFP had minimal effect on ν. Water remained a good solvent even with 21 consecutive glycine residues (ν > 0.5), and red-sfAFP variants lacked stable backbone hydrogen bonds according to hydrogen exchange. Similarly, changing glycine segregation has little impact on ν in other glycine-rich proteins. These findings underscore the generality that many disordered states can be expanded and unstructured, and that the hydrophobic effect alone is insufficient to drive significant chain collapse for typical protein sequences.

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.

Cold unfolding of heat-responsive TRPV3

The homotetrameric thermosensitive transient receptor potential vanilloid 1-4 (TRPV1-4) channels in sensory neurons are strongly responsive to heat stimuli. However, their cold activations have not been reported in line with the nonzero heat capacity difference during heat or cold unfolding transitions. Here, along with the experimental examinations of the predicted ring size changes in different domains against the central pore during channel gating at various temperatures, the K169A mutant of reduced human TRPV3 was first found to be activated and inactivated by cold below 42°C. Further thermoring analyses revealed distinct heat and cold unfolding pathways, which resulted in different protein thermostabilities. Thus, both cold and heat unfolding transitions of thermosensitive TRPV1-4 channels may exist once a mutation destabilizes the closed state.

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.

Structural biases in disordered proteins are prevalent in the cell

Intrinsically disordered proteins and protein regions (IDPs) are prevalent in all proteomes and are essential to cellular function. Unlike folded proteins, IDPs exist in an ensemble of dissimilar conformations. Despite this structural plasticity, intramolecular interactions create sequence-specific structural biases that determine an IDP ensemble's three-dimensional shape. Such structural biases can be key to IDP function and are often measured in vitro, but whether those biases are preserved inside the cell is unclear. Here we show that structural biases in IDP ensembles found in vitro are recapitulated inside human-derived cells. We further reveal that structural biases can change in a sequence-dependent manner due to changes in the intracellular milieu, subcellular localization, and intramolecular interactions with tethered well-folded domains. We propose that the structural sensitivity of IDP ensembles can be leveraged for biological function, can be the underlying cause of IDP-driven pathology or can be used to design disorder-based biosensors and actuators.

Understanding the structural basis of the binding specificity of c-di-AMP to M. smegmatis RecA using computational biology approach

Mycobacterium tuberculosis RecA (MtRecA), a protein involved in DNA repair, homologous recombination and SOS pathway, contributes to the development of multidrug resistance. ATP binding-site in RecA has been a drug target to disable RecA dependent DNA repair. For the first time, experiments have shown the existence and binding of c-di-AMP to a novel allosteric site in the C-terminal-Domain (CTD) of Mycobacterium smegmatis RecA (MsRecA), a close homolog of MtRecA. In addition, it was observed that the c-di-AMP was not binding to Escherichia coli RecA (EcRecA). This article analyses the possible interactions of the three RecA homologs with the various c-di-AMP conformations to gain insights into the structural basis of the natural preference of c-di-AMP to MsRecA and not to EcRecA, using the structural biology tools. The comparative analysis, based on amino acid composition, homology, motifs, residue types, docking, molecular dynamics simulations and binding free energy calculations, indeed, conclusively indicates strong binding of c-di-AMP to MsRecA. Having very similar results as MsRecA, it is highly plausible for c-di-AMP to strongly bind MtRecA as well. These insights from the in-silico studies adds a new therapeutic approach against TB through design and development of novel allosteric inhibitors for the first time against MtRecA.Communicated by Ramaswamy H. Sarma.

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.

Modulating hierarchical self-assembly in thermoresponsive intrinsically disordered proteins through high-temperature incubation time

The cornerstone of structural biology is the unique relationship between protein sequence and the 3D structure at equilibrium. Although intrinsically disordered proteins (IDPs) do not fold into a specific 3D structure, breaking this paradigm, some IDPs exhibit large-scale organization, such as liquid-liquid phase separation. In such cases, the structural plasticity has the potential to form numerous self-assembled structures out of thermal equilibrium. Here, we report that high-temperature incubation time is a defining parameter for micro and nanoscale self-assembly of resilin-like IDPs. Interestingly, high-resolution scanning electron microscopy micrographs reveal that an extended incubation time leads to the formation of micron-size rods and ellipsoids that depend on the amino acid sequence. More surprisingly, a prolonged incubation time also induces amino acid composition-dependent formation of short-range nanoscale order, such as periodic lamellar nanostructures. We, therefore, suggest that regulating the period of high-temperature incubation, in the one-phase regime, can serve as a unique method of controlling the hierarchical self-assembly mechanism of structurally disordered proteins.

Intrinsically disordered regions are poised to act as sensors of cellular chemistry

Intrinsically disordered proteins and protein regions (IDRs) are abundant in eukaryotic proteomes and play a wide variety of essential roles. Instead of folding into a stable structure, IDRs exist in an ensemble of interconverting conformations whose structure is biased by sequence-dependent interactions. The absence of a stable 3D structure, combined with high solvent accessibility, means that IDR conformational biases are inherently sensitive to changes in their environment. Here, we argue that IDRs are ideally poised to act as sensors and actuators of cellular physicochemistry. We review the physical principles that underlie IDR sensitivity, the molecular mechanisms that translate this sensitivity to function, and recent studies where environmental sensing by IDRs may play a key role in their downstream function.

Conformational entropy in molecular recognition of intrinsically disordered proteins

Broad conformational ensembles make intrinsically disordered proteins or regions entropically intriguing. Although methodologically challenging and understudied, emerging studies into their changes in conformational entropy (ΔS°conf) upon complex formation have provided both quantitative and qualitative insight. Recent work based on thermodynamics from isothermal titration calorimetry and NMR spectroscopy uncovers an expanded repertoire of regulatory mechanisms, where ΔS°conf plays roles in partner selection, state behavior, functional buffering, allosteric regulation, and drug design. We highlight these mechanisms to display the large entropic reservoir of IDPs for the regulation of molecular communication. We call upon the field to make efforts to contribute to this insight as more studies are needed for forwarding mechanistic decoding of intrinsically disordered proteins and their complexes.

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.

Modulating Hierarchical Self-Assembly In Thermoresponsive Intrinsically Disordered Proteins Through High-Temperature Incubation Time

The cornerstone of structural biology is the unique relationship between protein sequence and the 3D structure at equilibrium. Although intrinsically disordered proteins (IDPs) do not fold into a specific 3D structure, breaking this paradigm, some IDPs exhibit large-scale organization, such as liquid-liquid phase separation. In such cases, the structural plasticity has the potential to form numerous self-assembled structures out of thermal equilibrium. Here, we report that high-temperature incubation time is a defining parameter for micro and nanoscale self-assembly of resilin-like IDPs. Interestingly, high-resolution scanning electron microscopy micrographs reveal that an extended incubation time leads to the formation of micron-size rods and ellipsoids that depend on the amino acid sequence. More surprisingly, a prolonged incubation time also induces amino acid composition-dependent formation of short-range nanoscale order, such as periodic lamellar nanostructures. We can correlate the lamellar structures to β-sheet formation and demonstrate similarities between the observed nanoscopic structural arrangement and spider silk. We, therefore, suggest that regulating the period of high-temperature incubation, in the one-phase regime, can serve as a unique method of controlling the hierarchical self-assembly mechanism of structurally disordered proteins.

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.

Data-Driven Models for Predicting Intrinsically Disordered Protein Polymer Physics Directly from Composition or Sequence

The molecular-level understanding of intrinsically disordered proteins is challenging due to experimental characterization difficulties. Computational understanding of IDPs also requires fundamental advances, as the leading tools for predicting protein folding (e.g., AlphaFold), typically fail to describe the structural ensembles of IDPs. The focus of this paper is to 1) develop new representations for intrinsically disordered proteins and 2) pair these representations with classical machine learning and deep learning models to predict the radius of gyration and derived scaling exponent of IDPs. Here, we build a new physically-motivated feature called the bag of amino acid interactions representation, which encodes pairwise interactions explicitly into the representation. This feature essentially counts and weights all possible non-bonded interactions in a sequence and thus is, in principle, compatible with arbitrary sequence lengths. To see how well this new feature performs, both categorical and physically-motivated featurization techniques are tested on a computational dataset containing 10,000 sequences simulated at the coarse-grained level. The results indicate that this new feature outperforms the other purely categorical and physically-motivated features and possesses solid extrapolation capabilities. For future use, this feature can potentially provide physical insights into amino acid interactions, including their temperature dependence, and be applied to other protein spaces.

Decoupling of Interactions between Model-Charged Peptides Reveals Key Factors Responsible for Liquid-Liquid Phase Separation

Liquid-liquid phase separation (LLPS) by disordered proteins has been shown to govern biological processes and cause numerous diseases. Therefore, a deeper understanding of the interactions and their variation with external factors is key to modulating the LLPS behavior of different systems and protecting proteins from pathological aggregation. In this context, we have looked at interactions between similarly charged peptides to understand the molecular features that may drive or prevent condensate formation under various conditions. We have studied dimer formation for model peptides where charged and noncharged amino acids have been placed alternatively. Using arginine and glutamic acid as the charged residues and varying the other residues with glycine, alanine, and proline to alter hydrophobicity, we have obtained the free-energy surface (FES) for the dimer formation for these systems under high salt concentration at two different temperatures using all-atom molecular dynamics simulations and the well-tempered metadynamics method. Our results indicate that a combination of effects such as hydrophobicity, arginine-arginine interactions, or water release from the solvation shell makes dimerization free energy more favorable for the positively charged peptides with lower flexibility. For the negatively charged peptides, the crucial role of water has been found in governing the FES. Systems having charged residues and phenylalanine in the peptide sequence also have been studied at high salt concentrations using unbiased simulations. In this case, only the positively charged peptides were found to aggregate through temperature-dependent hydrophobic and cation-π interactions. Overall, our study indicates that the negatively charged peptides are more likely to remain in the dilute phase under various conditions compared to the positively charged systems. The findings from our study would be helpful in designing and controlling systems to obtain LLPS behavior for therapeutic usage.

Combining Experiments and Simulations to Examine the Temperature-Dependent Behavior of a Disordered Protein

Intrinsically disordered proteins are a class of proteins that lack stable folded conformations and instead adopt a range of conformations that determine their biochemical functions. The temperature-dependent behavior of such disordered proteins is complex and can vary depending on the specific protein and environment. Here, we have used molecular dynamics simulations and previously published experimental data to investigate the temperature-dependent behavior of histatin 5, a 24-residue-long polypeptide. We examined the hypothesis that histatin 5 undergoes a loss of polyproline II (PPII) structure with increasing temperature, leading to more compact conformations. We found that the conformational ensembles generated by the simulations generally agree with small-angle X-ray scattering data for histatin 5, but show some discrepancies with the hydrodynamic radius as probed by pulsed-field gradient NMR spectroscopy, and with the secondary structure information derived from circular dichroism. We attempted to reconcile these differences by reweighting the conformational ensembles against the scattering and NMR data. By doing so, we were in part able to capture the temperature-dependent behavior of histatin 5 and to link the observed decrease in hydrodynamic radius with increasing temperature to a loss of PPII structure. We were, however, unable to achieve agreement with both the scattering and NMR data within experimental errors. We discuss different possible reasons for this including inaccuracies in the force field, differences in conditions of the NMR and scattering experiments, and issues related to the calculation of the hydrodynamic radius from conformational ensembles. Our study highlights the importance of integrating multiple types of experimental data when modeling conformational ensembles of disordered proteins and how environmental factors such as the temperature influence them.

Excluded Volume and Weak Interactions in Crowded Solutions Modulate Conformations and RNA Binding of an Intrinsically Disordered Tail

The cellular milieu is a solution crowded with a significant concentration of different components (proteins, nucleic acids, metabolites, etc.). Such a crowded environment affects protein conformations, dynamics, and interactions. Intrinsically disordered proteins and regions are particularly sensitive to these effects. Here, we investigate the impact on an intrinsically disordered tail that flanks a folded domain, the N-terminal domain, and the RNA-binding domain of the SARS-CoV-2 nucleocapsid protein. We mimic the crowded environment of the cell using polyethylene glycol (PEG) and study its impact on protein conformations using single-molecule Förster resonance energy transfer. We found that high-molecular-weight PEG induces a collapse of the disordered N-terminal tail, whereas low-molecular-weight PEG induces a chain expansion. Our data can be explained by accounting for two opposing contributions: favorable interactions between the protein and crowder molecules and screening of excluded volume interactions. We further characterized the interaction between protein and RNA in the presence of crowding agents. While for all PEG molecules tested, we observed an increase in the binding affinity, the trend is not monotonic as a function of the degree of PEG polymerization. This points to the role of nonspecific protein-PEG interactions on binding in addition to the entropic effects due to crowding. To separate the enthalpic and entropic components of the effects, we investigated the temperature dependence of the association constants in the absence and presence of crowders. Finally, we compared the effects of crowding across mutations in the disordered region and found that the threefold difference in association constants for two naturally occurring variants of the SARS-CoV-2 nucleocapsid protein is reduced to almost identical affinities in the presence of crowders. Overall, our data provide new insights into understanding and modeling the contribution of crowding effects on disordered regions, including the impact of interactions between proteins and crowders and their interplay when binding a ligand.

Revisiting macromolecular hydration with HullRadSAS

Hydration of biological macromolecules is important for their stability and function. Historically, attempts have been made to describe the degree of macromolecular hydration using a single parameter over a narrow range of values. Here, we describe a method to calculate two types of hydration: surface shell water and entrained water. A consideration of these two types of hydration helps to explain the "hydration problem" in hydrodynamics. The combination of these two types of hydration allows accurate calculation of hydrodynamic volume and related macromolecular properties such as sedimentation and diffusion coefficients, intrinsic viscosities, and the concentration-dependent non-ideality identified with sedimentation velocity experiments.

Disordered proteins mitigate the temperature dependence of site-specific binding free energies

Biophysical characterization of protein-protein interactions involving disordered proteins is challenging. A common simplification is to measure the thermodynamics and kinetics of disordered site binding using peptides containing only the minimum residues necessary. We should not assume, however, that these few residues tell the whole story. Son of sevenless, a multidomain signaling protein from Drosophila melanogaster, is critical to the mitogen-activated protein kinase pathway, passing an external signal to Ras, which leads to cellular responses. The disordered 55 kDa C-terminal domain of Son of sevenless is an autoinhibitor that blocks guanidine exchange factor activity. Activation requires another protein, Downstream of receptor kinase (Drk), which contains two Src homology 3 domains. Here, we utilized NMR spectroscopy and isothermal titration calorimetry to quantify the thermodynamics and kinetics of the N-terminal Src homology 3 domain binding to the strongest sites incorporated into the flanking disordered sequences. Comparing these results to those for isolated peptides provides information about how the larger domain affects binding. The affinities of sites on the disordered domain are like those of the peptides at low temperatures but less sensitive to temperature. Our results, combined with observations showing that intrinsically disordered proteins become more compact with increasing temperature, suggest a mechanism for this effect.

Myofilament-associated proteins with intrinsic disorder (MAPIDs) and their resolution by computational modeling

The cardiac sarcomere is a cellular structure in the heart that enables muscle cells to contract. Dozens of proteins belong to the cardiac sarcomere, which work in tandem to generate force and adapt to demands on cardiac output. Intriguingly, the majority of these proteins have significant intrinsic disorder that contributes to their functions, yet the biophysics of these intrinsically disordered regions (IDRs) have been characterized in limited detail. In this review, we first enumerate these myofilament-associated proteins with intrinsic disorder (MAPIDs) and recent biophysical studies to characterize their IDRs. We secondly summarize the biophysics governing IDR properties and the state-of-the-art in computational tools toward MAPID identification and characterization of their conformation ensembles. We conclude with an overview of future computational approaches toward broadening the understanding of intrinsic disorder in the cardiac sarcomere.

The biophysics of disordered proteins from the point of view of single-molecule fluorescence spectroscopy

Intrinsically disordered proteins (IDPs) and regions (IDRs) have emerged as key players across many biological functions and diseases. Differently from structured proteins, disordered proteins lack stable structure and are particularly sensitive to changes in the surrounding environment. Investigation of disordered ensembles requires new approaches and concepts for quantifying conformations, dynamics, and interactions. Here, we provide a short description of the fundamental biophysical properties of disordered proteins as understood through the lens of single-molecule fluorescence observations. Single-molecule Förster resonance energy transfer (FRET) and fluorescence correlation spectroscopy (FCS) provides an extensive and versatile toolbox for quantifying the characteristics of conformational distributions and the dynamics of disordered proteins across many different solution conditions, both in vitro and in living cells.

pH Induced Switch in the Conformational Ensemble of Intrinsically Disordered Protein Prothymosin-α and Its Implications for Amyloid Fibril Formation

Aggregation of intrinsically disordered proteins (IDPs) can lead to neurodegenerative diseases. Although there is experimental evidence that acidic pH promotes IDP monomer compaction leading to aggregation, the general mechanism is unclear. We studied the pH effect on the conformational ensemble of prothymosin-α (proTα), which is involved in multiple essential functions, and probed its role in aggregation using computer simulations. We show that compaction in the proTα dimension at low pH is due to the protein's collapse in the intermediate region (E41-D80) rich in glutamic acid residues, enhancing its β-sheet content. We observed by performing dimer simulations that the conformations with high β-sheet content could act as aggregation-prone (N*) states and nucleate the aggregation process. The simulations initiated using N* states form dimers within a microsecond time scale, whereas the non-N* states do not form dimers within this time scale. This study contributes to understanding the general principles of pH-induced IDP aggregation.

Sticker-and-spacer model for amyloid beta condensation and fibrillation

A major pathogenic hallmark of Alzheimer's disease is the presence of neurotoxic plaques composed of amyloid beta (Aβ) peptides in patients' brains. The pathway of plaque formation remains elusive, though some clues appear to lie in the dominant presence of Aβ_{1 − 42} in these plaques despite Aβ₁₋₄₀ making up approximately 90% of the Aβ pool. We hypothesize that this asymmetry is driven by the hydrophobicity of the two extra amino acids that are incorporated in Aβ₁₋₄₂. To investigate this hypothesis at the level of single molecules, we have developed a molecular “sticker-and-spacer lattice model” of unfolded Aβ. The model protein has a single sticker that may reversibly dimerise and elongate into semi-flexible linear chains. The growth is hampered by excluded-volume interactions that are encoded by the hydrophilic spacers but are rendered cooperative by the attractive interactions of hydrophobic spacers. For sufficiently strong hydrophobicity, the chains undergo liquid-liquid phase-separation (LLPS) into condensates that facilitate the nucleation of fibers. We find that a small fraction of Aβ₁₋₄₀ in a mixture of Aβ₁₋₄₀ and Aβ₁₋₄₂ shifts the critical concentration for LLPS to lower values. This study provides theoretical support for the hypothesis that LLPS condensates act as a precursor for aggregation and provides an explanation for the Aβ₁₋₄₂-enrichment of aggregates in terms of hydrophobic interactions.

Characterization of the Conformations of Amyloid Beta 42 in Solution That May Mediate Its Initial Hydrophobic Aggregation

Intrinsically disordered peptides, such as amyloid β42 (Aβ42), lack a well-defined structure in solution. Aβ42 can undergo abnormal aggregation and amyloidogenesis in the brain, forming fibrillar plaques, a hallmark of Alzheimer's disease. The insoluble fibrillar forms of Aβ42 exhibit well-defined, cross β-sheet structures at the molecular level and are less toxic than the soluble, intermediate disordered oligomeric forms. However, the mechanism of initial interaction of monomers and subsequent oligomerization is not well understood. The structural disorder of Aβ42 adds to the challenges of determining the structural properties of its monomers, making it difficult to understand the underlying molecular mechanism of pathogenic aggregation. Certain regions of Aβ42 are known to exhibit helical propensity in different physiological conditions. NMR spectroscopy has shown that the Aβ42 monomer at lower pH can adopt an α-helical conformation and as the pH is increased, the peptide switches to β-sheet conformation and aggregation occurs. CD spectroscopy studies of aggregation have shown the presence of an initial spike in the amount of α-helical content at the start of aggregation. Such an increase in α-helical content suggests a mechanism wherein the peptide can expose critical non-polar residues for interaction, leading to hydrophobic aggregation with other interacting peptides. We have used molecular dynamics simulations to characterize in detail the conformational landscape of monomeric Aβ42 in solution to identify molecular properties that may mediate the early stages of oligomerization. We hypothesized that conformations with α-helical structure have a higher probability of initiating aggregation because they increase the hydrophobicity of the peptide. Although random coil conformations were found to be the most dominant, as expected, α-helical conformations are thermodynamically accessible, more so than β-sheet conformations. Importantly, for the first time α-helical conformations are observed to increase the exposure of aromatic and hydrophobic residues to the aqueous solvent, favoring their hydrophobically driven interaction with other monomers to initiate aggregation. These findings constitute a first step toward characterizing the mechanism of formation of disordered, low-order oligomers of Aβ42.

Calmodulin's Interdomain Linker Is Optimized for Dynamics Signal Transmission and Calcium Binding

Linkers are ubiquitous in multidomain proteins. These linkers are integral to protein functions, and accumulating evidence suggests that the linkers' versatile roles are encoded in their sequences. However, a molecular picture of how amino acid differences in the linker influence protein function is still lacking. By using extensive Gaussian-accelerated MD coupled with dynamic network analysis, we reveal the molecular bases underlying the linker's role in Calmodulin (CaM), a highly conserved Ca²⁺-signaling hub in eukaryotes. Three CaM constructs comprising a wild-type linker, a flexible linker (four glycines at position D78-S81), and a rigid linker (four prolines at position D78-S81) were simulated. We show that the flexible linker resembles the wild type in allowing CaM to sample a large ensemble of conformations while the rigid linker confines the sampling. Our simulations recapture experimental observations that target binding enhances the Ca²⁺ affinity to CaM's EF-hand sites at the N-domain. However, only the wild-type linker can both correctly capture the Ca²⁺ binding order and maintain the α-helical structure of the domain. The other two constructs either bind Ca²⁺ in an incorrect order or exhibit unfolding of an N-domain helix. We demonstrate that the wild-type linker achieves these outcomes by transmitting interdomain dynamics efficiently. This was evidenced by stronger (anti)correlations among the linker residues, decoupling of the hydrogen bonds between A1-A15 and V35-E45, and structuring of the N-domain for Ca²⁺ binding. This decoupling was not evident for the other two constructs. Lastly, we show that the wild-type linker's optimal transmission stems from its thermodynamically favorable strain and solvation relative to the other two constructs. Our results show how the linker sequence tunes CaM function, suggesting possible mechanisms for changes in linker properties such as mutations or post-translational modifications to modulate protein/substrate binding.

Influence of ion specificity and concentration on the conformational transition of intrinsically disordered sheep prion peptide

The structural sensitivity of the IDPs with the ions has been observed experimentally; however, it is still unclear how the presence of different metal ions affects structural stability. We performed atomistic molecular dynamics simulation of sheep prion peptide (142-167) in presence of different monovalent, divalent ions at various concentrations to find out the effect of the size, charge, and ionic concentration on the structure of the peptide. It is found that Li + ions have a higher survival probability compared to Na + , K + and Mg2 + affecting the solvation structure of the protein leading to the alpha-helix structure. At high concentration, due to the increase in the ion-solvent and ion-counter interactions, the effect of the ions is screened on the surface of the protein and hence no ion specificity is observed. This study demonstrates how ions can be used to regulate the protein structure and function that can help in designing drugs.

Effects of the Temperature and Salt Concentration on the Structural Characteristics of the Protein (PDB Code 1BBL)

The effect of the temperature and salt solution on the structural characteristics of the protein 1BBL was investigated by molecular dynamics simulations. The paper presents simulation results regarding the non-bonded energy and the structural stability of the protein immersed in salt solutions with different concentrations and temperatures. Our work demonstrates that the electrostatic potential energy and van der Waals energy of the system show the opposite changes with the influence of the external environment. Since the electrostatic potential energy changes more obviously, it is dominated in the non-bonding interactions. The structural parameters, such as the root mean square deviation and the radius of gyration, increased initially and decreased afterward with the increase of the salt concentration. The protein presented the loose structure with a relative low stability when it was immersed in a monovalent solution with a salt concentration of 0.8 mol/L. The salt concentration corresponding to the maximum value of structural parameters in the monovalent salt solution was double that in the divalent salt solution. It was also concluded that the protein presented a compact and stable structure when immersed in salt solutions with a high concentration of 2.3 mol/L. The analysis of the root mean square deviation and root mean square fluctuation of the protein sample also exhibited that the structural stability and chain flexibility are strongly guided by the effect of the temperature. These conclusions help us to understand the structural characteristics of the protein immersed in the salt solutions with different concentrations and temperatures.

Competing interactions give rise to two-state behavior and switch-like transitions in charge-rich intrinsically disordered proteins

The most commonly occurring intrinsically disordered proteins (IDPs) are polyampholytes, which are defined by the duality of low net charge per residue and high fractions of charged residues. Recent experiments have uncovered nuances regarding sequence–ensemble relationships of model polyampholytic IDPs. These include differences in conformational preferences for sequences with lysine vs. arginine and the suggestion that well-mixed sequences form a range of conformations, including globules, conformations with ensemble averages that are reminiscent of ideal chains, or self-avoiding walks. Here, we explain these observations by analyzing results from atomistic simulations. We find that polyampholytic IDPs generally sample two distinct stable states, namely, globules and self-avoiding walks. Globules are favored by electrostatic attractions between oppositely charged residues, whereas self-avoiding walks are favored by favorable free energies of hydration of charged residues. We find sequence-specific temperatures of bistability at which globules and self-avoiding walks can coexist. At these temperatures, ensemble averages over coexisting states give rise to statistics that resemble ideal chains without there being an actual counterbalancing of intrachain and chain-solvent interactions. At equivalent temperatures, arginine-rich sequences tilt the preference toward globular conformations whereas lysine-rich sequences tilt the preference toward self-avoiding walks. We also identify differences between aspartate- and glutamate-containing sequences, whereby the shorter aspartate side chain engenders preferences for metastable, necklace-like conformations. Finally, although segregation of oppositely charged residues within the linear sequence maintains the overall two-state behavior, compact states are highly favored by such systems.

Heterogeneity in Protein Folding and Unfolding Reactions

Proteins have dynamic structures that undergo chain motions on time scales spanning from picoseconds to seconds. Resolving the resultant conformational heterogeneity is essential for gaining accurate insight into fundamental mechanistic aspects of the protein folding reaction. The use of high-resolution structural probes, sensitive to population distributions, has begun to enable the resolution of site-specific conformational heterogeneity at different stages of the folding reaction. Different states populated during protein folding, including the unfolded state, collapsed intermediate states, and even the native state, are found to possess significant conformational heterogeneity. Heterogeneity in protein folding and unfolding reactions originates from the reduced cooperativity of various kinds of physicochemical interactions between various structural elements of a protein, and between a protein and solvent. Heterogeneity may arise because of functional or evolutionary constraints. Conformational substates within the unfolded state and the collapsed intermediates that exchange at rates slower than the subsequent folding steps give rise to heterogeneity on the protein folding pathways. Multiple folding pathways are likely to represent distinct sequences of structure formation. Insight into the nature of the energy barriers separating different conformational states populated during (un)folding can also be obtained by resolving heterogeneity.

Seeking Solvation: Exploring the Role of Protein Hydration in Silk Gelation

The mechanism by which arthropods (e.g., spiders and many insects) can produce silk fibres from an aqueous protein (fibroin) solution has remained elusive, despite much scientific investigation. In this work, we used several techniques to explore the role of a hydration shell bound to the fibroin in native silk feedstock (NSF) from Bombyx mori silkworms. Small angle X-ray and dynamic light scattering (SAXS and DLS) revealed a coil size (radius of gyration or hydrodynamic radius) around 12 nm, providing considerable scope for hydration. Aggregation in dilute aqueous solution was observed above 65 °C, matching the gelation temperature of more concentrated solutions and suggesting that the strength of interaction with the solvent (i.e., water) was the dominant factor. Infrared (IR) spectroscopy indicated decreasing hydration as the temperature was raised, with similar changes in hydration following gelation by freezing or heating. It was found that the solubility of fibroin in water or aqueous salt solutions could be described well by a relatively simple thermodynamic model for the stability of the protein hydration shell, which suggests that the affected water is enthalpically favoured but entropically penalised, due to its reduced (vibrational or translational) dynamics. Moreover, while the majority of this investigation used fibroin from B. mori, comparisons with published work on silk proteins from other silkworms and spiders, globular proteins and peptide model systems suggest that our findings may be of much wider significance.

Multiscale relaxation dynamics and diffusion of myelin basic protein in solution studied by quasielastic neutron scattering

We report an analysis of high-resolution quasielastic neutron scattering spectra from Myelin Basic Protein (MBP) in solution, comparing the spectra at three different temperatures (283, 303, and 323 K) for a pure D₂O buffer and a mixture of D₂O buffer with 30% of deuterated trifluoroethanol (TFE). Accompanying experiments with dynamic light scattering and Circular Dichroism (CD) spectroscopy have been performed to obtain, respectively, the global diffusion constant and the secondary structure content of the molecule for both buffers as a function of temperature. Modeling the decay of the neutron intermediate scattering function by the Mittag-Leffler relaxation function, ϕ(t) = E_α(-(t/τ)^α) (0 < α < 1), we find that trifluoroethanol slows down the relaxation dynamics of the protein at 283 K and leads to a broader relaxation rate spectrum. This effect vanishes with increasing temperature, and at 323 K, its relaxation dynamics is identical in both solvents. These results are coherent with the data from dynamic light scattering, which show that the hydrodynamic radius of MBP in TFE-enriched solutions does not depend on temperature and is only slightly smaller compared to the pure D₂O buffer, except for 283 K, where it is much reduced. In accordance with these observations, the CD spectra reveal that TFE induces essentially a partial transition from β-strands to α-helices, but only a weak increase in the total secondary structure content, leaving about 50% of the protein unfolded. The results show that MBP is for all temperatures and in both buffers an intrinsically disordered protein and that TFE essentially induces a reduction in its hydrodynamic radius and its relaxation dynamics at low temperatures.

Quantifying charge state heterogeneity for proteins with multiple ionizable residues

Ionizable residues can release and take up protons and this has an influence on protein structure and function. The extent of protonation is linked to the overall pH of the solution and the local environments of ionizable residues. Binding or unbinding of a single proton generates a distinct charge microstate defined by a specific pattern of charges. Accordingly, the overall partition function is a sum over all charge microstates and Boltzmann weights of all conformations associated with each of the charge microstates. This ensemble-of-ensembles description recast as a q-canonical ensemble allows us to analyze and interpret potentiometric titrations that provide information regarding net charge as a function of pH. In the q-canonical ensemble, charge microstates are grouped into mesostates where each mesostate is a collection of microstates of the same net charge. Here, we show that leveraging the structure of the q-canonical ensemble allows us to decouple contributions of net proton binding and release from proton arrangement and conformational considerations. Through application of the q-canonical formalism to analyze potentiometric measurements of net charge in proteins with repetitive patterns of Lys and Glu residues, we determine the underlying mesostate pKa values and, more importantly, we estimate relative mesostate populations as a function of pH. This is a strength of using the q-canonical approach that cannot be replicated using purely site-specific analyses. Overall, our work shows how measurements of charge equilibria, decoupled from measurements of conformational equilibria, and analyzed using the framework of the q-canonical ensemble, provide protein-specific quantitative descriptions of pH-dependent populations of mesostates. This method is of direct relevance for measuring and understanding how different charge states contribute to conformational, binding, and phase equilibria of proteins.

Beta turn propensity and a model polymer scaling exponent identify intrinsically disordered phase-separating proteins

The complex cellular milieu can spontaneously demix, or phase separate, in a process controlled in part by intrinsically disordered (ID) proteins. A protein's propensity to phase separate is thought to be driven by a preference for protein-protein over protein-solvent interactions. The hydrodynamic size of monomeric proteins, as quantified by the polymer scaling exponent (v), is driven by a similar balance. We hypothesized that mean v, as predicted by protein sequence, would be smaller for proteins with a strong propensity to phase separate. To test this hypothesis, we analyzed protein databases containing subsets of proteins that are folded, disordered, or disordered and known to spontaneously phase separate. We find that the phase-separating disordered proteins, on average, had lower calculated values of v compared with their non-phase-separating counterparts. Moreover, these proteins had a higher sequence-predicted propensity for β-turns. Using a simple, surface area-based model, we propose a physical mechanism for this difference: transient β-turn structures reduce the desolvation penalty of forming a protein-rich phase and increase exposure of atoms involved in π/sp2 valence electron interactions. By this mechanism, β-turns could act as energetically favored nucleation points, which may explain the increased propensity for turns in ID regions (IDRs) utilized biologically for phase separation. Phase-separating IDRs, non-phase-separating IDRs, and folded regions could be distinguished by combining v and β-turn propensity. Finally, we propose a new algorithm, ParSe (partition sequence), for predicting phase-separating protein regions, and which is able to accurately identify folded, disordered, and phase-separating protein regions based on the primary sequence.

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

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

When Order Meets Disorder: Modeling and Function of the Protein Interface in Fuzzy Complexes

The degree of proteins structural organization ranges from highly structured, compact folding to intrinsic disorder, where each degree of self-organization corresponds to specific functions: well-organized structural motifs in enzymes offer a proper environment for precisely positioned functional groups to participate in catalytic reactions; at the other end of the self-organization spectrum, intrinsically disordered proteins act as binding hubs via the formation of multiple, transient and often non-specific interactions. This review focusses on cases where structurally organized proteins or domains associate with highly disordered protein chains, leading to the formation of interfaces with varying degrees of fuzziness. We present a review of the computational methods developed to provide us with information on such fuzzy interfaces, and how they integrate experimental information. The discussion focusses on two specific cases, microtubules and homologous recombination nucleoprotein filaments, where a network of intrinsically disordered tails exerts regulatory function in recruiting partner macromolecules, proteins or DNA and tuning the atomic level association. Notably, we show how computational approaches such as molecular dynamics simulations can bring new knowledge to help bridging the gap between experimental analysis, that mostly concerns ensemble properties, and the behavior of individual disordered protein chains that contribute to regulation functions.

The Protein Folding Problem: The Role of Theory

The protein folding problem was first articulated as question of how order arose from disorder in proteins: How did the various native structures of proteins arise from interatomic driving forces encoded within their amino acid sequences, and how did they fold so fast? These matters have now been largely resolved by theory and statistical mechanics combined with experiments. There are general principles. Chain randomness is overcome by solvation-based codes. And in the needle-in-a-haystack metaphor, native states are found efficiently because protein haystacks (conformational ensembles) are funnel-shaped. Order-disorder theory has now grown to encompass a large swath of protein physical science across biology.

Single-Molecule FRET of Membrane Transport Proteins

Uncovering the structure and function of biomolecules is a fundamental goal in structural biology. Membrane-embedded transport proteins are ubiquitous in all kingdoms of life. Despite structural flexibility, their mechanisms are typically studied by ensemble biochemical methods or by static high-resolution structures, which complicate a detailed understanding of their dynamics. Here, we review the recent progress of single molecule Förster Resonance Energy Transfer (smFRET) in determining mechanisms and timescales of substrate transport across membranes. These studies do not only demonstrate the versatility and suitability of state-of-the-art smFRET tools for studying membrane transport proteins but they also highlight the importance of membrane mimicking environments in preserving the function of these proteins. The current achievements advance our understanding of transport mechanisms and have the potential to facilitate future progress in drug design.

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

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

Atomic view of cosolute-induced protein denaturation probed by NMR solvent paramagnetic relaxation enhancement

The cosolvent effect arises from the interaction of cosolute molecules with a protein and alters the equilibrium between native and unfolded states. Denaturants shift the equilibrium toward the latter, while osmolytes stabilize the former. The molecular mechanism whereby cosolutes perturb protein stability is still the subject of considerable debate. Probing the molecular details of the cosolvent effect is experimentally challenging as the interactions are very weak and transient, rendering them invisible to most conventional biophysical techniques. Here, we probe cosolute-protein interactions by means of NMR solvent paramagnetic relaxation enhancement together with a formalism we recently developed to quantitatively describe, at atomic resolution, the energetics and dynamics of cosolute-protein interactions in terms of a concentration normalized equilibrium average of the interspin distance, [Formula: see text], and an effective correlation time, τ_c The system studied is the metastable drkN SH3 domain, which exists in dynamic equilibrium between native and unfolded states, thereby permitting us to probe the interactions of cosolutes with both states simultaneously under the same conditions. Two paramagnetic cosolute denaturants were investigated, one neutral and the other negatively charged, differing in the presence of a carboxyamide group versus a carboxylate. Our results demonstrate that attractive cosolute-protein backbone interactions occur largely in the unfolded state and some loop regions in the native state, electrostatic interactions reduce the [Formula: see text] values, and temperature predominantly impacts interactions with the unfolded state. Thus, destabilization of the native state in this instance arises predominantly as a consequence of interactions of the cosolutes with the unfolded state.