Temperature-dependent structural changes in intrinsically disordered proteins: formation of alpha-helices or loss of polyproline II?

Optimizing Force Fields with Experimental Data Using Ensemble Reweighting and Potential Contrasting

Despite force field improvements over the past decades, we still encounter situations where simulation results disagree with experiments due to force field inaccuracies. Such situations provide opportunities to improve force fields. In this study, we introduce a novel framework for optimizing force fields using experimental data. The unique feature of this framework is that it aims to optimize force fields to match experiments while minimizing the perturbation made to the original force field. To achieve this, we combine ensemble reweighting techniques with the potential contrasting method. Ensemble reweighting is used to reweight an ensemble of conformations generated using an existing force field to match experimental data while minimizing the perturbation to the original ensemble. Potential contrasting is then utilized to optimize force field parameters to reproduce the reweighted ensemble. We demonstrate the framework's effectiveness by optimizing a coarse-grained force field for intrinsically disordered proteins using experimental radius of gyration data.

Architectonic principles of polyproline II helix bundle protein domains

Glycine rich polyproline II helix assemblies are an emerging class of natural domains found in several proteins with different functions and diverse origins. The distinct properties of these domains relative to those composed of α-helices and β-sheets could make glycine-rich polyproline II helix assemblies a useful building block for protein design. Whereas the high population of polyproline II conformers in disordered state ensembles could facilitate glycine-rich polyproline II helix folding, the architectonic bases of these structures are not well known. Here, we compare and analyze their structures to uncover common features. These protein domains are found to be highly tolerant of distinct flanking sequences. This speaks to the robustness of this fold and strongly suggests that glycine rich polyproline II assemblies could be grafted with other protein domains to engineer new structures and functions. These domains are also well packed with few or no cavities. Moreover, a significant trend towards antiparallel helix configuration is observed in all these domains and could provide stabilizing interactions among macrodipoles. Finally, extensive networks of Cα-H···OC hydrogen bonds are detected in these domains. Despite their diverse evolutionary origins and activities, glycine-rich polyproline II helix assemblies share architectonic features which could help design novel proteins.

PACSAB Server: A Web-Based Tool for the Study of Aggregation and the Conformational Ensemble of Disordered and Folded Proteins

We present in this article the PACSAB server, which is designed to provide information about the structural ensemble and interactions of both stable and disordered proteins to researchers in the field of molecular biology. The use of this tool does not require any computational skills as the user just needs to upload the structure of the protein to be studied; the server runs a simulation with the PACSAB model, a highly accurate coarse-grained model that is much more efficient than standard molecular dynamics for the exploration of the conformational space of multiprotein systems. The trajectories generated by the simulations based on this model reveal the propensity of the protein under study for aggregation, identify the residues playing a central role in the aggregation process, and reproduce the whole conformational space of disordered proteins. All of this information is shown and can be downloaded from the web page.

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.

Biophysical characterization and insights into the oligomeric nature of CD2-associated protein

INTRODUCTION

Glomerular podocytes are specialized epithelial cells localized to the blood-urine interface of the kidney. Podocyte slit-diaphragm (SD), a size-and-charge-selective junction, is instrumental in blood ultrafiltration and the formation of protein-free urine. The SD consists of macromolecular complexes of several proteins, such as nephrin, podocin, and CD2-associated protein (CD2AP). CD2AP is an adapter protein and is considered to be crucial for the integrity of SD. Mutations in the SD proteins cause nephrotic syndrome (NS), characterized by proteinuria. SD proteins' structural features must be elucidated to understand the mechanism of proteinuria in NS. In this study, we expressed, purified, and biophysically characterized heterologously expressed human CD2AP.

METHODS

Codon-optimized human CD2AP was expressed in E. coli Rosetta cells. The recombinant protein was induced with 1 mM IPTG and purified by Ni-NTA affinity chromatography. Analytical size-exclusion chromatography, blue native-PAGE, circular dichroism, and fluorescence spectroscopy were performed to decipher the oligomeric nature, secondary structural content, and tertiary packing of CD2AP.

RESULTS

Our analysis revealed that CD2AP adopts a predominantly disordered secondary structure despite exhibiting moderate tertiary packing, characterized by low helical and β-sheet content. CD2AP readily assembles into homo-oligomers, with octamers and tetramers constituting the primary population. Interestingly, the inherent flexibility of CD2AP's secondary structural elements appears resistant to thermal denaturation. Frameshift mutation (p.K579Efs*7) that leads to loss of the coiled-coil domain promotes aberrant oligomerization of CD2AP through SH3 domains.

CONCLUSION

We successfully expressed full-length human CD2AP in a heterologous system, wherein the secondary structure of CD2AP is predominantly disordered. CD2AP can form higher-order oligomers, and the significance of these oligomers and the impact of mutations in the context of size-selective permeability of SD needs further investigation.

Study of Correlated Motions to Detect the Conformational Transitions of the Intrinsically Disordered Sheep Prion Peptide

Intrinsically disordered proteins (IDPs) are known for their random structural changes throughout their sequence based on the environment. The mechanism underlying these structural changes is difficult to explain. All biological processes are known to follow the direction through which they act. A study of the correlated motion can help to understand the direction of the change. Herein, we introduced the multivariate statistical analysis (MSA) technique to study the correlated motion of the peptide. The correlated motion of the sheep prion peptide was studied with the change in the temperature and solvent. These techniques helped to identify the contributing residual motions that helped to form the different secondary structures of the protein and also the triggering factors that drive these sorts of residual motions. The structural details match the experimentally reported data. It was found that the direction of the change of the secondary structure for this peptide shifted from the C-terminal to the N-terminal with an increase in the temperature. It was found that the involvement of the hydrophobic residues present at the C-terminal and the middle residues (residues 12-17) is responsible for forming a β-sheet at the normal temperature. Hydration water was found to play an important role in this change. Insights gained from this study can be used to design strategies for desirable structural changes in the IDPs.

Experimental methods to study the structure and dynamics of intrinsically disordered regions in proteins

Eukaryotic proteins often feature long stretches of amino acids that lack a well-defined three-dimensional structure and are referred to as intrinsically disordered proteins (IDPs) or regions (IDRs). Although these proteins challenge conventional structure-function paradigms, they play vital roles in cellular processes. Recent progress in experimental techniques, such as NMR spectroscopy, single molecule FRET, high speed AFM and SAXS, have provided valuable insights into the biophysical basis of IDP function. This review discusses the advancements made in these techniques particularly for the study of disordered regions in proteins. In NMR spectroscopy new strategies such as 13C detection, non-uniform sampling, segmental isotope labeling, and rapid data acquisition methods address the challenges posed by spectral overcrowding and low stability of IDPs. The importance of various NMR parameters, including chemical shifts, hydrogen exchange rates, and relaxation measurements, to reveal transient secondary structures within IDRs and IDPs are presented. Given the high flexibility of IDPs, the review outlines NMR methods for assessing their dynamics at both fast (ps-ns) and slow (μs-ms) timescales. IDPs exert their functions through interactions with other molecules such as proteins, DNA, or RNA. NMR-based titration experiments yield insights into the thermodynamics and kinetics of these interactions. Detailed study of IDPs requires multiple experimental techniques, and thus, several methods are described for studying disordered proteins, highlighting their respective advantages and limitations. The potential for integrating these complementary techniques, each offering unique perspectives, is explored to achieve a comprehensive understanding of IDPs.

FRETpredict: a Python package for FRET efficiency predictions using rotamer libraries

Förster resonance energy transfer (FRET) is a widely-used and versatile technique for the structural characterization of biomolecules. Here, we introduce FRETpredict, an easy-to-use Python software to predict FRET efficiencies from ensembles of protein conformations. FRETpredict uses a rotamer library approach to describe the FRET probes covalently bound to the protein. The software efficiently and flexibly operates on large conformational ensembles such as those generated by molecular dynamics simulations to facilitate the validation or refinement of molecular models and the interpretation of experimental data. We provide access to rotamer libraries for many commonly used dyes and linkers and describe a general methodology to generate new rotamer libraries for FRET probes. We demonstrate the performance and accuracy of the software for different types of systems: a rigid peptide (polyproline 11), an intrinsically disordered protein (ACTR), and three folded proteins (HiSiaP, SBD2, and MalE). FRETpredict is open source (GPLv3) and is available at github.com/KULL-Centre/FRETpredict and as a Python PyPI package at pypi.org/project/FRETpredict .

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.

Phase-Dependent Adsorption of Myelin Basic Protein to Phosphatidylcholine Lipid Bilayers

The dense packing of opposite cytoplasmic surfaces of the lipid-enriched myelin membrane, responsible for the proper saltatory conduction of nerve impulses through axons, is ensured by the adhesive properties of myelin basic protein (MBP). Although preferentially interacting with negatively charged phosphatidylserine (PS) lipids, as an intrinsically disordered protein, it can easily adapt its shape to its immediate environment and thus adsorb to domains made of zwitterionic phosphatidylcholine (PC) lipids. As the molecular-level interaction pattern between MBP and PC lipid membranes suffers from scarce characterization, an experimental and computational study of multilamellar liposomes (MLVs) composed of 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) in the presence of bovine MBP is presented here. Calorimetric and temperature-dependent UV-Vis measurements identified DPPC pretransition temperature (Tp) and calorimetric enthalpy (ΔHcal) as the physicochemical parameters most responsive to the presence of MBP. Besides suggesting an increase in β-sheet fractions of structured MBP segments as DPPC lipids undergo from the gel (20 °C) to the fluid (50 °C) phase, FTIR spectra unraveled the significant contribution of lysine (Lys) residues in the adsorption pattern, especially when DPPC is in the fluid (50 °C) phase. In addition to highlighting the importance of Lys residues in the MBP adsorption on DPPC lipid bilayer, employing salt bridges (SBs) and hydrogen bonds (HBs), MD data suggest the crucial importance of the orientation of MBP with respect to the surface of the DPPC lipid bilayer.

Charge Distribution Patterns of IA3 Impact Conformational Expansion and Hydration Diffusivity of the Disordered Ensemble

IA3 is a 68 amino acid natural peptide/protein inhibitor of yeast aspartic proteinase A (YPRA) that is intrinsically disordered in solution with induced N-terminal helicity when in the protein complex with YPRA. Based on the intrinsically disordered protein (IDP) parameters of fractional net charge (FNC), net charge density per residue (NCPR), and charge patterning (κ), the two domains of IA3 are defined to occupy different domains within conformationally based subclasses of IDPs, thus making IA3 a bimodal domain IDP. Site-directed spin labeling (SDSL) electron paramagnetic resonance (EPR) spectroscopy and low-field Overhauser dynamic nuclear polarization (ODNP) spectroscopy results show that these two domains possess different degrees of compaction and hydration diffusivity behavior. This work suggests that SDSL EPR line shapes, analyzed in terms of their local tumbling volume (VL), provide insights into the compaction of the unstructured IDP ensemble in solution and that protein sequence and net charge distribution patterns within a conformational subclass can impact bound water hydration dynamics, thus possibly offering an alternative thermodynamic property that can encode conformational binding and behavior of IDPs and liquid-liquid phase separations.

Structural characterization of an intrinsically disordered protein complex using integrated small-angle neutron scattering and computing

Characterizing structural ensembles of intrinsically disordered proteins (IDPs) and intrinsically disordered regions (IDRs) of proteins is essential for studying structure-function relationships. Due to the different neutron scattering lengths of hydrogen and deuterium, selective labeling and contrast matching in small-angle neutron scattering (SANS) becomes an effective tool to study dynamic structures of disordered systems. However, experimental timescales typically capture measurements averaged over multiple conformations, leaving complex SANS data for disentanglement. We hereby demonstrate an integrated method to elucidate the structural ensemble of a complex formed by two IDRs. We use data from both full contrast and contrast matching with residue-specific deuterium labeling SANS experiments, microsecond all-atom molecular dynamics (MD) simulations with four molecular mechanics force fields, and an autoencoder-based deep learning (DL) algorithm. From our combined approach, we show that selective deuteration provides additional information that helps characterize structural ensembles. We find that among the four force fields, a99SB-disp and CHARMM36m show the strongest agreement with SANS and NMR experiments. In addition, our DL algorithm not only complements conventional structural analysis methods but also successfully differentiates NMR and MD structures which are indistinguishable on the free energy surface. Lastly, we present an ensemble that describes experimental SANS and NMR data better than MD ensembles generated by one single force field and reveal three clusters of distinct conformations. Our results demonstrate a new integrated approach for characterizing structural ensembles of IDPs.

Accurate Simulation of Coupling between Protein Secondary Structure and Liquid-Liquid Phase Separation

Intrinsically disordered proteins (IDPs) frequently mediate liquid-liquid phase separation (LLPS) that underlies the formation of membraneless organelles. Together with theory and experiment, efficient coarse-grained (CG) simulations have been instrumental in understanding sequence-specific phase separation of IDPs. However, the widely-used Cα-only models are severely limited in capturing the peptide nature of IDPs, including backbone-mediated interactions and effects of secondary structures, in LLPS. Here, we describe a hybrid resolution (HyRes) protein model for accurate description of the backbone and transient secondary structures in LLPS. With an atomistic backbone and coarse-grained side chains, HyRes accurately predicts the residue helical propensity and chain dimension of monomeric IDPs. Using GY-23 as a model system, we show that HyRes is efficient enough for direct simulation of spontaneous phase separation, and at the same time accurate enough to resolve the effects of single mutations. HyRes simulations also successfully predict increased beta-sheet formation in the condensate, consistent with available experimental data. We further utilize HyRes to study the phase separation of TPD-43, where several disease-related mutants in the conserved region (CR) have been shown to affect residual helicities and modulate LLPS propensity. The simulations successfully recapitulate the effect of these mutants on the helicity and LLPS propensity of TDP-43 CR. Analyses reveal that the balance between backbone and sidechain-mediated interactions, but not helicity itself, actually determines LLPS propensity. We believe that the HyRes model represents an important advance in the molecular simulation of LLPS and will help elucidate the coupling between IDP transient secondary structures and phase separation.

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.

The relevance of short peptides for an understanding of unfolded and intrinsically disordered proteins

Over the last thirty years the unfolded state of proteins has attracted considerable interest owing to the discovery of intrinsically disordered proteins which perform a plethora of functions despite resembling unfolded proteins to a significant extent. Research on both, unfolded and disordered proteins has revealed that their conformational properties can deviate locally from random coil behavior. In this context results from work on short oligopeptides suggest that individual amino acid residues sample the sterically allowed fraction of the Ramachandran plot to a different extent. Alanine has been found to exhibit a peculiarity in that it has a very high propensity for adopting polyproline II like conformations. This Perspectives article reviews work on short peptides aimed at exploring the Ramachandran distributions of amino acid residues in different contexts with experimental and computational means. Based on the thus provided overview the article discussed to what extent short peptides can serve as tools for exploring unfolded and disordered proteins and as benchmarks for the development of a molecular dynamics force field.

Interpreting Transient Interactions of Intrinsically Disordered Proteins

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

Residue-specific binding of Ni(II) ions influences the structure and aggregation of amyloid beta (Aβ) peptides

Alzheimer's disease (AD) is the most common cause of dementia worldwide. AD brains display deposits of insoluble amyloid plaques consisting mainly of aggregated amyloid-β (Aβ) peptides, and Aβ oligomers are likely a toxic species in AD pathology. AD patients display altered metal homeostasis, and AD plaques show elevated concentrations of metals such as Cu, Fe, and Zn. Yet, the metal chemistry in AD pathology remains unclear. Ni(II) ions are known to interact with Aβ peptides, but the nature and effects of such interactions are unknown. Here, we use numerous biophysical methods-mainly spectroscopy and imaging techniques-to characterize Aβ/Ni(II) interactions in vitro, for different Aβ variants: Aβ(1-40), Aβ(1-40)(H6A, H13A, H14A), Aβ(4-40), and Aβ(1-42). We show for the first time that Ni(II) ions display specific binding to the N-terminal segment of full-length Aβ monomers. Equimolar amounts of Ni(II) ions retard Aβ aggregation and direct it towards non-structured aggregates. The His6, His13, and His14 residues are implicated as binding ligands, and the Ni(II)·Aβ binding affinity is in the low µM range. The redox-active Ni(II) ions induce formation of dityrosine cross-links via redox chemistry, thereby creating covalent Aβ dimers. In aqueous buffer Ni(II) ions promote formation of beta sheet structure in Aβ monomers, while in a membrane-mimicking environment (SDS micelles) coil-coil helix interactions appear to be induced. For SDS-stabilized Aβ oligomers, Ni(II) ions direct the oligomers towards larger sizes and more diverse (heterogeneous) populations. All of these structural rearrangements may be relevant for the Aβ aggregation processes that are involved in AD brain pathology.

Improved predictions of phase behaviour of intrinsically disordered proteins by tuning the interaction range

The formation and viscoelastic properties of condensates of intrinsically disordered proteins (IDPs) is dictated by amino acid sequence and solution conditions. Because of the involvement of biomolecular condensates in cell physiology and disease, advancing our understanding of the relationship between protein sequence and phase separation (PS) may have important implications in the formulation of new therapeutic hypotheses. Here, we present CALVADOS 2, a coarse-grained model of IDPs that accurately predicts conformational properties and propensities to undergo PS for diverse sequences and solution conditions. In particular, we systematically study the effect of varying the range of the nonionic interactions and use our findings to improve the temperature scale of the model. We further optimize the residue-specific model parameters against experimental data on the conformational properties of 55 proteins, while also leveraging 70 hydrophobicity scales from the literature to avoid overfitting the training data. Extensive testing shows that the model accurately predicts chain compaction and PS propensity for sequences of diverse length and charge patterning, as well as at different temperatures and salt concentrations.

Improved predictions of phase behaviour of intrinsically disordered proteins by tuning the interaction range

The formation and viscoelastic properties of condensates of intrinsically disordered proteins (IDPs) is dictated by amino acid sequence and solution conditions. Because of the involvement of biomolecular condensates in cell physiology and disease, advancing our understanding of the relationship between protein sequence and phase separation (PS) may have important implications in the formulation of new therapeutic hypotheses. Here, we present CALVADOS 2, a coarse-grained model of IDPs that accurately predicts conformational properties and propensities to undergo PS for diverse sequences and solution conditions. In particular, we systematically study the effect of varying the range of the nonionic interactions and use our findings to improve the temperature scale of the model. We further optimize the residue-specific model parameters against experimental data on the conformational properties of 55 proteins, while also leveraging 70 hydrophobicity scales from the literature to avoid overfitting the training data. Extensive testing shows that the model accurately predicts chain compaction and PS propensity for sequences of diverse length and charge patterning, as well as at different temperatures and salt concentrations.

Functional Implications of Dynamic Structures of Intrinsically Disordered Proteins Revealed by High-Speed AFM Imaging

The unique functions of intrinsically disordered proteins (IDPs) depend on their dynamic protean structure that often eludes analysis. High-speed atomic force microscopy (HS-AFM) can conduct this difficult analysis by directly visualizing individual IDP molecules in dynamic motion at sub-molecular resolution. After brief descriptions of the microscopy technique, this review first shows that the intermittent tip-sample contact does not alter the dynamic structure of IDPs and then describes how the number of amino acids contained in a fully disordered region can be estimated from its HS-AFM images. Next, the functional relevance of a dumbbell-like structure that has often been observed on IDPs is discussed. Finally, the dynamic structural information of two measles virus IDPs acquired from their HS-AFM and NMR analyses is described together with its functional implications.

A context-dependent and disordered ubiquitin-binding motif

Ubiquitin is a small, globular protein that is conjugated to other proteins as a posttranslational event. A palette of small, folded domains recognizes and binds ubiquitin to translate and effectuate this posttranslational signal. Recent computational studies have suggested that protein regions can recognize ubiquitin via a process of folding upon binding. Using peptide binding arrays, bioinformatics, and NMR spectroscopy, we have uncovered a disordered ubiquitin-binding motif that likely remains disordered when bound and thus expands the palette of ubiquitin-binding proteins. We term this motif Disordered Ubiquitin-Binding Motif (DisUBM) and find it to be present in many proteins with known or predicted functions in degradation and transcription. We decompose the determinants of the motif showing it to rely on features of aromatic and negatively charged residues, and less so on distinct sequence positions in line with its disordered nature. We show that the affinity of the motif is low and moldable by the surrounding disordered chain, allowing for an enhanced interaction surface with ubiquitin, whereby the affinity increases ~ tenfold. Further affinity optimization using peptide arrays pushed the affinity into the low micromolar range, but compromised context dependence. Finally, we find that DisUBMs can emerge from unbiased screening of randomized peptide libraries, featuring in de novo cyclic peptides selected to bind ubiquitin chains. We suggest that naturally occurring DisUBMs can recognize ubiquitin as a posttranslational signal to act as affinity enhancers in IDPs that bind to folded and ubiquitylated binding partners.

NKCC1 and KCC2: Structural insights into phospho-regulation

Inhibitory neurotransmission plays a fundamental role in the central nervous system, with about 30–50% of synaptic connections being inhibitory. The action of both inhibitory neurotransmitter, gamma-aminobutyric-acid (GABA) and glycine, mainly relies on the intracellular Cl^– concentration in neurons. This is set by the interplay of the cation chloride cotransporters NKCC1 (Na⁺, K⁺, Cl^– cotransporter), a main Cl^– uptake transporter, and KCC2 (K⁺, Cl^– cotransporter), the principle Cl^– extruder in neurons. Accordingly, their dysfunction is associated with severe neurological, psychiatric, and neurodegenerative disorders. This has triggered great interest in understanding their regulation, with a strong focus on phosphorylation. Recent structural data by cryogenic electron microscopy provide the unique possibility to gain insight into the action of these phosphorylations. Interestingly, in KCC2, six out of ten (60%) known regulatory phospho-sites reside within a region of 134 amino acid residues (12% of the total residues) between helices α8 and α9 that lacks fixed or ordered three-dimensional structures. It thus represents a so-called intrinsically disordered region. Two further phospho-sites, Tyr⁹⁰³ and Thr⁹⁰⁶, are also located in a disordered region between the ß8 strand and the α8 helix. We make the case that especially the disordered region between helices α8 and α9 acts as a platform to integrate different signaling pathways and simultaneously constitute a flexible, highly dynamic linker that can survey a wide variety of distinct conformations. As each conformation can have distinct binding affinities and specificity properties, this enables regulation of [Cl^–]_i and thus the ionic driving force in a history-dependent way. This region might thus act as a molecular processor underlying the well described phenomenon of ionic plasticity that has been ascribed to inhibitory neurotransmission. Finally, it might explain the stunning long-range effects of mutations on phospho-sites in KCC2.

Negative Thermal Expansion and Disorder-to-Order Collapse of an Intrinsically Disordered Protein under Marginally Denaturing Conditions

Recent work with intrinsically disordered proteins (IDPs) has projected a myriad of their survival instincts based mainly on the total charge content, the abundance of polar residues, and the paucity of hydrophobic amino acids. This work uses a plant IDP AtPP16-1 (Arabidopsis thaliana phloem protein class 16-1), whose solution NMR structure was determined by us recently, to show legitimate negative thermal expansion (NTE) of the native state. The thermal expansion continues to be negative even as the tertiary structure is perturbed by ultralow levels of urea up to 0.4 M. The NTE of these subdenatured states is called apparent NTE because they are prone to undergo conformational changes with temperature. Hydrodynamic shrinkage of the NTE IDP is also observed by dynamic light scattering (DLS) and NMR-measured global rotational correlation time (τ_c). The protein with denatured tertiary structure but otherwise preserved native-state secondary structure collapses to a dynamically rigid state. The data are mainly based on thermal coefficients of chemical shift and nuclear relaxation measured by heteronuclear NMR. The hydrodynamic shrinkage and collapse under marginally varying solvent compositions that may arise from unstable tertiary structure and dynamic disorder of chain segments across the backbone could be a generic property of IDPs.

A fine balance of hydrophobic-electrostatic communication pathways in a pH-switching protein

Allostery is the phenomenon of coupling between distal binding sites in a protein. Such coupling is at the crux of protein function and regulation in a myriad of scenarios, yet determining the molecular mechanisms of coupling networks in proteins remains a major challenge. Here, we report mechanisms governing pH-dependent myristoyl switching in monomeric hisactophilin, whereby the myristoyl moves between a sequestered state, i.e., buried within the core of the protein, to an accessible state, in which the myristoyl has increased accessibility for membrane binding. Measurements of the pH and temperature dependence of amide chemical shifts reveal protein local structural stability and conformational heterogeneity that accompany switching. An analysis of these measurements using a thermodynamic cycle framework shows that myristoyl-proton coupling at the single-residue level exists in a fine balance and extends throughout the protein. Strikingly, small changes in the stereochemistry or size of core and surface hydrophobic residues by point mutations readily break, restore, or tune myristoyl switch energetics. Synthesizing the experimental results with those of molecular dynamics simulations illuminates atomistic details of coupling throughout the protein, featuring a large network of hydrophobic interactions that work in concert with key electrostatic interactions. The simulations were critical for discerning which of the many ionizable residues in hisactophilin are important for switching and identifying the contributions of nonnative interactions in switching. The strategy of using temperature-dependent NMR presented here offers a powerful, widely applicable way to elucidate the molecular mechanisms of allostery in proteins at high resolution.

A modular approach to map out the conformational landscapes of unbound intrinsically disordered proteins

SignificanceIntrinsically disordered proteins have the unique ability to morph in response to multiple partners and thereby process sophisticated inputs and outputs. It is, however, a mystery whether their response is passive, that is, entirely determined by the partner, or controlled via an internal, yet unknown, folding mechanism. Here we introduce an approach to examine this key question and demonstrate its potential by dissecting the conformational properties of the partially disordered protein NCBD and obtaining important clues about how it performs its biological function.

The Disordered EZH2 Loop: Atomic Level Characterization by 1HN- and 1Hα-Detected NMR Approaches, Interaction with the Long Noncoding HOTAIR RNA

The 96-residue-long loop of EZH2 is proposed to play a role in the interaction with long non-coding RNAs (lncRNAs) and to contribute to EZH2 recruitment to the chromatin. However, molecular details of RNA recognition have not been described so far. Cellular studies have suggested that phosphorylation of the Thr345 residue localized in this loop influences RNA binding; however, no mechanistic explanation has been offered. To address these issues, a systematic NMR study was performed. As the ¹H^N-detected NMR approach presents many challenges under physiological conditions, our earlier developed, as well as improved, ¹H^α-detected experiments were used. As a result of the successful resonance assignment, the obtained chemical shift values indicate the highly disordered nature of the EZH2 loop, with some nascent helical tendency in the Ser407–Ser412 region. Further investigations conducted on the phosphomimetic mutant EZH2^T345D showed that the mutation has only a local effect, and that the loop remains disordered. On the other hand, the mutation influences the cis/trans Pro346 equilibrium. Interactions of both the wild-type and the phosphomimetic mutant with the lncRNA HOTAIR₁₄₀ (1–140 nt) highlight that the Thr367–Ser375 region is affected. This segment does not resemble any of the previously reported RNA-binding motifs, therefore the identified binding region is unique. As no structural changes occur in the EZH2 loop upon RNA binding, we can consider the protein–RNA interaction as a “fuzzy” complex.

Expression, Purification, Characterization and Cellular Uptake of MeCP2 Variants

The transcriptional regulator Methyl-CpG-binding protein 2 (MeCP2) is an intrinsically disordered protein, mutations in which, are implicated in the onset of Rett Syndrome, a severe and debilitating neurodevelopmental disorder. Delivery of this protein fused to the cell-penetrating peptide TAT could allow for the intracellular replenishment of functional MeCP2 and hence potentially serve as a prospective Rett Syndrome therapy. This work outlines the expression, purification and characterization of various TAT-MeCP2 constructs as well as their full-length and shortened eGFP fusion variants. The latter two constructs were used for intracellular uptake studies with subsequent analysis via western blotting and live-cell imaging. All purified MeCP2 samples exhibited high degree of stability and very little aggregation propensity. Full length and minimal TAT-MeCP2-eGFP were found to efficiently transduce into human dermal and murine fibroblasts and localize to cell nuclei. These findings clearly support the utility of MeCP2-based protein replacement therapy as a potential Rett Syndrome treatment option.

Residual Structure in the Denatured State of the Fast-Folding UBA(1) Domain from the Human DNA Excision Repair Protein HHR23A

The structure of the first ubiquitin-associated domain from HHR23A, UBA(1), was determined by X-ray crystallography at a 1.60 Å resolution, and its stability, folding kinetics, and residual structure under denaturing conditions have been investigated. The concentration dependence of thermal denaturation and size-exclusion chromatography indicate that UBA(1) is monomeric. Guanidine hydrochloride (GdnHCl) denaturation experiments reveal that the unfolding free energy, ΔG_u°'(H₂O), of UBA(1) is 2.4 kcal mol^-1. Stopped-flow folding kinetics indicates sub-millisecond folding with only proline isomerization phases detectable at 25 °C. The full folding kinetics are observable at 4 °C, yielding a folding rate constant, k_f, in the absence of a denaturant of 13,000 s^-1 and a Tanford β-value of 0.80, consistent with a compact transition state. Evaluation of the secondary structure via circular dichroism shows that the residual helical structure in the denatured state is replaced by polyproline II structure as the GdnHCl concentration increases. Analysis of NMR secondary chemical shifts for backbone ¹⁵NH, ¹³CO, and ¹³Cα atoms between 4 and 7 M GdnHCl shows three islands of residual helical secondary structure that align in sequence with the three native-state helices. Extrapolation of the NMR data to 0 M GdnHCl demonstrates that helical structure would populate to 17-33% in the denatured state under folding conditions. Comparison with NMR data for a peptide corresponding to helix 1 indicates that this helix is stabilized by transient tertiary interactions in the denatured state of UBA(1). The high helical content in the denatured state, which is enhanced by transient tertiary interactions, suggests a diffusion-collision folding mechanism.

Improving Martini 3 for Disordered and Multidomain Proteins

Coarse-grained molecular dynamics simulations are a useful tool to determine conformational ensembles of proteins. Here, we show that the coarse-grained force field Martini 3 underestimates the global dimensions of intrinsically disordered proteins (IDPs) and multidomain proteins when compared with small-angle X-ray scattering (SAXS) data and that increasing the strength of protein-water interactions favors more expanded conformations. We find that increasing the strength of interactions between protein and water by ca. 10% results in improved agreement with the SAXS data for IDPs and multidomain proteins. We also show that this correction results in a more accurate description of self-association of IDPs and folded proteins and better agreement with paramagnetic relaxation enhancement data for most IDPs. While simulations with this revised force field still show deviations to experiments for some systems, our results suggest that it is overall a substantial improvement for coarse-grained simulations of soluble proteins.

Intrinsically disordered proteins: Ensembles at the limits of Anfinsen's dogma

Intrinsically disordered proteins (IDPs) are proteins that lack rigid 3D structure. Hence, they are often misconceived to present a challenge to Anfinsen's dogma. However, IDPs exist as ensembles that sample a quasi-continuum of rapidly interconverting conformations and, as such, may represent proteins at the extreme limit of the Anfinsen postulate. IDPs play important biological roles and are key components of the cellular protein interaction network (PIN). Many IDPs can interconvert between disordered and ordered states as they bind to appropriate partners. Conformational dynamics of IDPs contribute to conformational noise in the cell. Thus, the dysregulation of IDPs contributes to increased noise and "promiscuous" interactions. This leads to PIN rewiring to output an appropriate response underscoring the critical role of IDPs in cellular decision making. Nonetheless, IDPs are not easily tractable experimentally. Furthermore, in the absence of a reference conformation, discerning the energy landscape representation of the weakly funneled IDPs in terms of reaction coordinates is challenging. To understand conformational dynamics in real time and decipher how IDPs recognize multiple binding partners with high specificity, several sophisticated knowledge-based and physics-based in silico sampling techniques have been developed. Here, using specific examples, we highlight recent advances in energy landscape visualization and molecular dynamics simulations to discern conformational dynamics and discuss how the conformational preferences of IDPs modulate their function, especially in phenotypic switching. Finally, we discuss recent progress in identifying small molecules targeting IDPs underscoring the potential therapeutic value of IDPs. Understanding structure and function of IDPs can not only provide new insight on cellular decision making but may also help to refine and extend Anfinsen's structure/function paradigm.

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.

Macromolecular crowding induces a binding competent transient structure in intrinsically disordered Gab1

Intrinsically disordered proteins (IDPs) are an important class of proteins which lack tertiary structure elements. Their dynamic properties can depend on reversible post-translational modifications and the complex cellular milieu, which provides a crowded environment. Both influences the thermodynamic stability and folding of globular proteins as well as the conformational plasticity of IDPs. Here we investigate the intrinsically disordered C-terminal region (amino acids 613-694) of human Grb2-associated binding protein 1 (Gab1), which binds to the disease-relevant Src homolog region2 (SH2) domain-containing protein tyrosine phosphatase SHP2 (PTPN11). This binding is mediated by phosphorylation at Tyr 627 and Tyr 659 in Gab1. We characterize induced structure in Gab1_613-694 and binding to SHP2 by NMR, CD and ITC under non-crowding and crowding conditions, employing chemical and biological crowding agents and compare the results of the non-phosphorylated and tyrosine phosphorylated C-terminal Gab1 fragment. Our results show that under crowding conditions pre-structured motifs in two distinct regions of Gab1 are formed whereas phosphorylation has no impact on the dynamics and IDP character. These structured regions are identical to the binding regions towards SHP2. Therefore, biological crowders could induce some SHP2 binding capacity. Our results therefore indicate that high concentrations of macromolecules stabilize the preformed or excited binding state in the C-terminal Gab1 region and foster the binding to the SH2 tandem motif of SHP2, even in the absence of tyrosine phosphorylation.

Molecular Characterization of an Intrinsically Disordered Chaperone Reveals Net-Charge Regulation in Chaperone Action

Molecular chaperones are diverse biomacromolecules involved in the maintenance of cellular protein homeostasis (proteostasis). Here we demonstrate that in contrast to most chaperones with defined three-dimensional structures, the acid-inducible protein Asr in Escherichia coli is intrinsically disordered and exhibits varied aggregation-preventing or aggregation-promoting activities, acting as a "conditionally active chaperone". Bioinformatics and experimental analyses of Asr showed that it is devoid of hydrophobic patches but rich in positive charges and local polyproline II backbone structures. Asr contributes to the integrity of the bacterial outer membrane under mildly acidic conditions in vivo and possesses chaperone activities toward model clients in vitro. Notably, its chaperone activity is dependent on the net charges of clients: on the one hand, it inhibits the aggregation of clients with similar net charges; on the other hand, it stimulates the aggregation of clients with opposite net charges. Mutational analysis confirmed that positively charged residues in Asr are essential for the varied effects on protein aggregation, suggesting that electrostatic interactions are the major driving forces underlying Asr's proteostasis-related activity. These findings present a unique example of an intrinsically disordered molecular chaperone with distinctive dual functions-as an aggregase or as a chaperone-depending on the net charges of clients.

Accurate Description of Protein-Protein Recognition and Protein Aggregation with the Implicit-Solvent-Based PACSAB Protein Model

We used the PACSAB protein model, based on the implicit solvation approach, to simulate protein–protein recognition and study the effect of helical structure on the association of aggregating peptides. After optimization, the PACSAB force field was able to reproduce correctly both the correct binding interface in ubiquitin dimerization and the conformational ensemble of the disordered protein activator for hormone and retinoid receptor (ACTR). The PACSAB model allowed us to predict the native binding of ACTR with its binding partner, reproducing the refolding upon binding mechanism of the disordered protein.

Accurate model of liquid-liquid phase behavior of intrinsically disordered proteins from optimization of single-chain properties

Many intrinsically disordered proteins (IDPs) may undergo liquid-liquid phase separation (LLPS) and participate in the formation of membraneless organelles in the cell, thereby contributing to the regulation and compartmentalization of intracellular biochemical reactions. The phase behavior of IDPs is sequence dependent, and its investigation through molecular simulations requires protein models that combine computational efficiency with an accurate description of intramolecular and intermolecular interactions. We developed a general coarse-grained model of IDPs, with residue-level detail, based on an extensive set of experimental data on single-chain properties. Ensemble-averaged experimental observables are predicted from molecular simulations, and a data-driven parameter-learning procedure is used to identify the residue-specific model parameters that minimize the discrepancy between predictions and experiments. The model accurately reproduces the experimentally observed conformational propensities of a set of IDPs. Through two-body as well as large-scale molecular simulations, we show that the optimization of the intramolecular interactions results in improved predictions of protein self-association and LLPS.

Accurate model of liquid-liquid phase behavior of intrinsically disordered proteins from optimization of single-chain properties

Many intrinsically disordered proteins (IDPs) may undergo liquid-liquid phase separation (LLPS) and participate in the formation of membraneless organelles in the cell, thereby contributing to the regulation and compartmentalization of intracellular biochemical reactions. The phase behavior of IDPs is sequence dependent, and its investigation through molecular simulations requires protein models that combine computational efficiency with an accurate description of intramolecular and intermolecular interactions. We developed a general coarse-grained model of IDPs, with residue-level detail, based on an extensive set of experimental data on single-chain properties. Ensemble-averaged experimental observables are predicted from molecular simulations, and a data-driven parameter-learning procedure is used to identify the residue-specific model parameters that minimize the discrepancy between predictions and experiments. The model accurately reproduces the experimentally observed conformational propensities of a set of IDPs. Through two-body as well as large-scale molecular simulations, we show that the optimization of the intramolecular interactions results in improved predictions of protein self-association and LLPS.

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.

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.

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.

Investigating the conformational dynamics of SARS-CoV-2 NSP6 protein with emphasis on non-transmembrane 91-112 & 231-290 regions

The NSP6 protein of SARS-CoV-2 is a transmembrane protein, with some regions lying outside the membrane. Besides a brief role of NSP6 in autophagosome formation, this is not studied significantly. Also, there is no structural information available to date. Based on the prediction by TMHMM server for transmembrane prediction, it is found that the N-terminal residues (1-11), middle region residues (91–112), and C-terminal residues (231–290) lies outside the membrane. Molecular Dynamics (MD) simulations showed that NSP6 consists of helical structures. In contrast, the membrane outside lying region (91–112) showed partial helicity, which was further used as a model and obtained disordered type conformation during 1.5 μs. Additionally, a 200ns simulation study of residues 231–290 have shown significant conformational changes. As compared to helical and beta-sheet conformations in its structure model, the 200ns simulation resulted in the loss of beta-sheet structures while helical regions remained intact. Further, we have experimentally characterized the residue 91–112 by using reductionist approaches. CD spectroscopy suggests that the NSP6 (91–112) is disordered-like region in isolation, which gains helical conformation in different biological mimic environmental conditions. These studies can be helpful to study NSP6 (91–112) interactions with host proteins, where different protein conformations might play a significant role. The present study adds up more information about the NSP6 protein aspect, which could be exploited for its host protein interaction and pathogenesis.

Salt dependent phase behavior of intrinsically disordered proteins from a coarse-grained model with explicit water and ions

Multivalent proteins and nucleic acids can self-assemble into biomolecular condensates that contribute to compartmentalize the cell interior. Computer simulations offer a unique view to elucidate the mechanisms and key intermolecular interactions behind the dynamic formation and dissolution of these condensates. In this work, we present a novel approach to include explicit water and salt in sequence-dependent coarse-grained (CG) models for proteins and RNA, enabling the study of biomolecular condensate formation in a salt-dependent manner. Our framework combines a reparameterized version of the HPS protein force field with the monoatomic mW water model and the mW-ion potential for NaCl. We show how our CG model qualitatively captures the experimental radius of the gyration trend of a subset of intrinsically disordered proteins and reproduces the experimental protein concentration and water percentage of the human fused in sarcoma (FUS) low-complexity-domain droplets at physiological salt concentration. Moreover, we perform seeding simulations as a function of salt concentration for two antagonist systems: the engineered peptide PR₂₅ and poly-uridine/poly-arginine mixtures, finding good agreement with their reported in vitro phase behavior with salt concentration in both cases. Taken together, our work represents a step forward towards extending sequence-dependent CG models to include water and salt, and to consider their key role in biomolecular condensate self-assembly.

Ocular disease-associated mutations diminish the mitotic chromosome retention ability of PAX6

Transcription factors play a key role in maintaining cell identity. One mechanism of such cell memory after multiple rounds of cell division cycles is through persistent mitotic chromosome binding, although how individual transcription factors achieve mitotic chromosome retention is not completely understood. Here we show that PAX6, a lineage-determining transcription factor, coats mitotic chromosomes. Using deletion and point mutants associated with human ocular diseases in live-cell imaging analysis, we identified two regions, MCR-D1 and MCR-D2, that were responsible for mitotic chromosome retention of PAX6. We also identified three nuclear localization signals (NLSs) that contributed to mitotic chromosome retention independent of their nuclear import functions. Full mitotic chromosome retention required the presence of DNA-binding domains as well as NLSs within MCR-Ds. Furthermore, disease-associated mutations and NLS mutations changed the distribution of intrinsically disordered regions (IDRs) in PAX6. Our findings not only identify PAX6 as a novel mitotic chromosome retention factor but also demonstrate that the mechanism of mitotic chromosome retention involves sequence-specific DNA binding, NLSs, and molecular conformation determined by IDRs. These findings link mitotic chromosome retention with PAX6-related pathogenesis and imply similar mechanisms for other lineage-determining factors in the PAX family.

Insight into Calcium-Binding Motifs of Intrinsically Disordered Proteins

Motifs within proteins help us categorize their functions. Intrinsically disordered proteins (IDPs) are rich in short linear motifs, conferring them many different roles. IDPs are also frequently highly charged and, therefore, likely to interact with ions. Canonical calcium-binding motifs, such as the EF-hand, often rely on the formation of stabilizing flanking helices, which are a key characteristic of folded proteins, but are absent in IDPs. In this study, we probe the existence of a calcium-binding motif relevant to IDPs. Upon screening several carefully selected IDPs using NMR spectroscopy supplemented with affinity quantification by colorimetric assays, we found calcium-binding motifs in IDPs which could be categorized into at least two groups-an Excalibur-like motif, sequentially similar to the EF-hand loop, and a condensed-charge motif carrying repetitive negative charges. The motifs show an affinity for calcium typically in the ~100 μM range relevant to regulatory functions and, while calcium binding to the condensed-charge motif had little effect on the overall compaction of the IDP chain, calcium binding to Excalibur-like motifs resulted in changes in compaction. Thus, calcium binding to IDPs may serve various structural and functional roles that have previously been underreported.

The disordered PCI-binding human proteins CSNAP and DSS1 have diverged in structure and function

Intrinsically disordered proteins (IDPs) regularly constitute components of larger protein assemblies contributing to architectural stability. Two small, highly acidic IDPs have been linked to the so-called PCI complexes carrying PCI-domain subunits, including the proteasome lid and the COP9 signalosome. These two IDPs, DSS1 and CSNAP, have been proposed to have similar structural propensities and functions, but they display differences in their interactions and interactome sizes. Here we characterized the structural properties of human DSS1 and CSNAP at the residue level using NMR spectroscopy and probed their propensities to bind ubiquitin. We find that distinct structural features present in DSS1 are completely absent in CSNAP, and vice versa, with lack of relevant ubiquitin binding to CSNAP, suggesting the two proteins to have diverged in both structure and function. Our work additionally highlights that different local features of seemingly similar IDPs, even subtle sequence variance, may endow them with different functional traits. Such traits may underlie their potential to engage in multiple interactions thereby impacting their interactome sizes.

Conformational dynamics of 13 amino acids long NSP11 of SARS-CoV-2 under membrane mimetics and different solvent conditions

The intrinsically disordered proteins/regions (IDPs/IDPRs) are known to be responsible for multiple cellular processes and are associated with many chronic diseases. In viruses, the existence of a disordered proteome is also proven and is related to its conformational dynamics inside the host. The SARS-CoV-2 has a large proteome, in which, structure and functions of all proteins are not known yet, along with non-structural protein 11 (nsp11). In this study, we have performed extensive experimentation on nsp11. Our results based on the CD spectroscopy gives characteristic disordered spectrum for IDPs. Further, we investigated the conformational behavior of nsp11 in the presence of membrane mimetic environment, α-helix inducer, and natural osmolyte. In the presence of negatively charged and neutral liposomes, nsp11 remains disordered. However, with SDS micelle, it adopted an α-helical conformation, suggesting the helical propensity of nsp11. Finally, we again confirmed the IDP behavior of nsp11 using MD simulations. In future, this conformational dynamic study could help to clarify its functional importance in SARS-CoV-2 infection.

SARS-CoV-2 NSP1 C-terminal (residues 131-180) is an intrinsically disordered region in isolation

The NSP1- C terminal structure in complex with ribosome using cryo-EM is available now, and the N-terminal region structure in isolation is also deciphered in literature. However, as a reductionist approach, the conformation of NSP1- C terminal region (NSP1-CTR; amino acids 131-180) has not been studied in isolation. We found that NSP1-CTR conformation is disordered in an aqueous solution. Further, we examined the conformational propensity towards alpha-helical structure using trifluoroethanol, we observed induction of helical structure conformation using CD spectroscopy. Additionally, in SDS, NSP1-CTR shows a conformational change from disordered to ordered, possibly gaining alpha-helix in part. But in the presence of neutral lipid DOPC, a slight change in conformation is observed, which implies the possible role of hydrophobic interaction and electrostatic interaction on the conformational changes of NSP1. Fluorescence-based studies have shown a blue shift and fluorescence quenching in the presence of SDS, TFE, and lipid vesicles. In agreement with these results, fluorescence lifetime and fluorescence anisotropy decay suggest a change in conformational dynamics. The zeta potential studies further validated that the conformational dynamics are primarily because of hydrophobic interaction. These experimental studies were complemented through Molecular Dynamics (MD) simulations, which have shown a good correlation and testifies our experiments. We believe that the intrinsically disordered nature of the NSP1-CTR will have implications for enhanced molecular recognition feature properties of this IDR, which may add disorder to order transition and disorder-based binding promiscuity with its interacting proteins.

High resolution ensemble description of metamorphic and intrinsically disordered proteins using an efficient hybrid parallel tempering scheme

Mapping free energy landscapes of complex multi-funneled metamorphic proteins and weakly-funneled intrinsically disordered proteins (IDPs) remains challenging. While rare-event sampling molecular dynamics simulations can be useful, they often need to either impose restraints or reweigh the generated data to match experiments. Here, we present a parallel-tempering method that takes advantage of accelerated water dynamics and allows efficient and accurate conformational sampling across a wide variety of proteins. We demonstrate the improved sampling efficiency by benchmarking against standard model systems such as alanine di-peptide, TRP-cage and β-hairpin. The method successfully scales to large metamorphic proteins such as RFA-H and to highly disordered IDPs such as Histatin-5. Across the diverse proteins, the calculated ensemble averages match well with the NMR, SAXS and other biophysical experiments without the need to reweigh. By allowing accurate sampling across different landscapes, the method opens doors for sampling free energy landscape of complex uncharted proteins.

Structural and Energetic Characterization of the Denatured State from the Perspectives of Peptides, the Coil Library, and Intrinsically Disordered Proteins

The α and polyproline II (PPII) basins are the two most populated regions of the Ramachandran map when constructed from the protein coil library, a widely used denatured state model built from the segments of irregular structure found in the Protein Data Bank. This indicates the α and PPII conformations are dominant components of the ensembles of denatured structures that exist in solution for biological proteins, an observation supported in part by structural studies of short, and thus unfolded, peptides. Although intrinsic conformational propensities have been determined experimentally for the common amino acids in short peptides, and estimated from surveys of the protein coil library, the ability of these intrinsic conformational propensities to quantitatively reproduce structural behavior in intrinsically disordered proteins (IDPs), an increasingly important class of proteins in cell function, has thus far proven elusive to establish. Recently, we demonstrated that the sequence dependence of the mean hydrodynamic size of IDPs in water and the impact of heat on the coil dimensions, provide access to both the sequence dependence and thermodynamic energies that are associated with biases for the α and PPII backbone conformations. Here, we compare results from peptide-based studies of intrinsic conformational propensities and surveys of the protein coil library to those of the sequence-based analysis of heat effects on IDP hydrodynamic size, showing that a common structural and thermodynamic description of the protein denatured state is obtained.

Characterization of AMBN I and II Isoforms and Study of Their Ca2+-Binding Properties

Ameloblastin (Ambn) as an intrinsically disordered protein (IDP) stands for an important role in the formation of enamel—the hardest biomineralized tissue commonly formed in vertebrates. The human ameloblastin (AMBN) is expressed in two isoforms: full-length isoform I (AMBN ISO I) and isoform II (AMBN ISO II), which is about 15 amino acid residues shorter than AMBN ISO I. The significant feature of AMBN—its oligomerization ability—is enabled due to a specific sequence encoded by exon 5 present at the N-terminal part in both known isoforms. In this study, we characterized AMBN ISO I and AMBN ISO II by biochemical and biophysical methods to determine their common features and differences. We confirmed that both AMBN ISO I and AMBN ISO II form oligomers in in vitro conditions. Due to an important role of AMBN in biomineralization, we further addressed the calcium (Ca²⁺)-binding properties of AMBN ISO I and ISO II. The binding properties of AMBN to Ca²⁺ may explain the role of AMBN in biomineralization and more generally in Ca²⁺ homeostasis processes.