Hydrodynamic Radii of Intrinsically Disordered Proteins Determined from Experimental Polyproline II Propensities

Hydrodynamic Radii of Intrinsically Disordered Proteins: Fast Prediction by Minimum Dissipation Approximation and Experimental Validation

The diffusion coefficients of globular and fully unfolded proteins can be predicted with high accuracy solely from their mass or chain length. However, this approach fails for intrinsically disordered proteins (IDPs) containing structural domains. We propose a rapid predictive methodology for estimating the diffusion coefficients of IDPs. The methodology uses accelerated conformational sampling based on self-avoiding random walks and includes hydrodynamic interactions between coarse-grained protein subunits, modeled using the generalized Rotne-Prager-Yamakawa approximation. To estimate the hydrodynamic radius, we rely on the minimum dissipation approximation recently introduced by Cichocki et al. Using a large set of experimentally measured hydrodynamic radii of IDPs over a wide range of chain lengths and domain contributions, we demonstrate that our predictions are more accurate than the Kirkwood approximation and phenomenological approaches. Our technique may prove to be valuable in predicting the hydrodynamic properties of both fully unstructured and multidomain disordered proteins.

Diversity of hydrodynamic radii of intrinsically disordered proteins

Intrinsically disordered proteins (IDPs) form an important class of biomolecules regulating biological processes in higher organisms. The lack of a fixed spatial structure facilitates them to perform their regulatory functions and allows the efficiency of biochemical reactions to be controlled by temperature and the cellular environment. From the biophysical point of view, IDPs are biopolymers with a broad configuration state space and their actual conformation depends on non-covalent interactions of its amino acid side chain groups at given temperature and chemical conditions. Thus, the hydrodynamic radius (Rh) of an IDP of a given polymer length (N) is a sequence- and environment-dependent variable. We have reviewed the literature values of hydrodynamic radii of IDPs determined experimentally by SEC, AUC, PFG NMR, DLS, and FCS, and complement them with our FCS results obtained for a series of protein fragments involved in the regulation of human gene expression. The data collected herein show that the values of hydrodynamic radii of IDPs can span the full space between the folded globular and denatured proteins in the Rh(N) diagram.

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.

Illuminating Intrinsically Disordered Proteins with Integrative Structural Biology

Intense study of intrinsically disordered proteins (IDPs) did not begin in earnest until the late 1990s when a few groups, working independently, convinced the community that these 'weird' proteins could have important functions. Over the past two decades, it has become clear that IDPs play critical roles in a multitude of biological phenomena with prominent examples including coordination in signaling hubs, enabling gene regulation, and regulating ion channels, just to name a few. One contributing factor that delayed appreciation of IDP functional significance is the experimental difficulty in characterizing their dynamic conformations. The combined application of multiple methods, termed integrative structural biology, has emerged as an essential approach to understanding IDP phenomena. Here, we review some of the recent applications of the integrative structural biology philosophy to study IDPs.

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

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

Different Forms of Disorder in NMDA-Sensitive Glutamate Receptor Cytoplasmic Domains Are Associated with Differences in Condensate Formation

The N-methyl-D-aspartate (NMDA)-sensitive glutamate receptor (NMDAR) helps assemble downstream signaling pathways through protein interactions within the postsynaptic density (PSD), which are mediated by its intracellular C-terminal domain (CTD). The most abundant NMDAR subunits in the brain are GluN2A and GluN2B, which are associated with a developmental switch in NMDAR composition. Previously, we used single molecule fluorescence resonance energy transfer (smFRET) to show that the GluN2B CTD contained an intrinsically disordered region with slow, hop-like conformational dynamics. The CTD from GluN2B also undergoes liquid-liquid phase separation (LLPS) with synaptic proteins. Here, we extend these observations to the GluN2A CTD. Sequence analysis showed that both subunits contain a form of intrinsic disorder classified as weak polyampholytes. However, only GluN2B contained matched patterning of arginine and aromatic residues, which are linked to LLPS. To examine the conformational distribution, we used discrete molecular dynamics (DMD), which revealed that GluN2A favors extended disordered states containing secondary structures while GluN2B favors disordered globular states. In contrast to GluN2B, smFRET measurements found that GluN2A lacked slow conformational dynamics. Thus, simulation and experiments found differences in the form of disorder. To understand how this affects protein interactions, we compared the ability of these two NMDAR isoforms to undergo LLPS. We found that GluN2B readily formed condensates with PSD-95 and SynGAP, while GluN2A failed to support LLPS and instead showed a propensity for colloidal aggregation. That GluN2A fails to support this same condensate formation suggests a developmental switch in LLPS propensity.

Intrinsically disordered regions that drive phase separation form a robustly distinct protein class

Protein phase separation is thought to be a primary driving force for the formation of membrane-less organelles, which control a wide range of biological functions from stress response to ribosome biogenesis. Among phase-separating (PS) proteins, many have intrinsically disordered regions (IDRs) that are needed for phase separation to occur. Accurate identification of IDRs that drive phase separation is important for testing the underlying mechanisms of phase separation, identifying biological processes that rely on phase separation, and designing sequences that modulate phase separation. To identify IDRs that drive phase separation, we first curated datasets of folded, ID, and PS ID sequences. We then used these sequence sets to examine how broadly existing amino acid property scales can be used to distinguish between the three classes of protein regions. We found that there are robust property differences between the classes and, consequently, that numerous combinations of amino acid property scales can be used to make robust predictions of protein phase separation. This result indicates that multiple, redundant mechanisms contribute to the formation of phase-separated droplets from IDRs. The top-performing scales were used to further optimize our previously developed predictor of PS IDRs, ParSe. We then modified ParSe to account for interactions between amino acids and obtained reasonable predictive power for mutations that have been designed to test the role of amino acid interactions in driving protein phase separation. Collectively, our findings provide further insight into the classification of IDRs and the elements involved in protein phase separation.

Prediction of Variable-Length B-Cell Epitopes for Antipeptide Paratopes Using the Program HAPTIC

BACKGROUND

B-cell epitope prediction for antipeptide antibody responses enables peptide-based vaccine design and related translational applications. This entails estimating epitopeparatope binding free-energy changes from antigen sequence; but attempts to do so assuming uniform epitope length (e.g., of hexapeptide sequences, each spanning a typical paratope diameter when fully extended) have neglected empirically established variation in epitope length.

OBJECTIVE

This work aimed to develop a sequence-based physicochemical approach to variablelength B-cell epitope prediction for antipeptide paratopes recognizing flexibly disordered targets.

METHODS

Said approach was developed by analogy between epitope-paratope binding and protein folding modeled as polymer collapse, treating paratope structure implicitly. Epitope-paratope binding was thus conceptually resolved into processes of epitope compaction, collapse and contact, with epitope collapse presenting the main entropic barrier limiting epitope length among nonpolyproline sequences. The resulting algorithm was implemented as a computer program, namely the Heuristic Affinity Prediction Tool for Immune Complexes (HAPTIC), which is freely accessible via an online interface (http://badong.freeshell.org/haptic.htm). This was used in conjunction with published data on representative known peptide immunogens.

RESULTS

HAPTIC predicted immunodominant epitope sequences with lengths limited by penalties for both compaction and collapse, consistent with known paratope-bound structures of flexibly disordered epitopes. In most cases, the predicted association constant was greater than its experimentally determined counterpart but below the predicted upper bound for affinity maturation in vivo.

CONCLUSION

HAPTIC provides a physicochemically plausible means for estimating the affinity of antipeptide paratopes for sterically accessible and flexibly disordered peptidic antigen sequences by explicitly considering candidate B-cell epitopes of variable length.

DHX15-independent roles for TFIP11 in U6 snRNA modification, U4/U6.U5 tri-snRNP assembly and pre-mRNA splicing fidelity

The U6 snRNA, the core catalytic component of the spliceosome, is extensively modified post-transcriptionally, with 2'-O-methylation being most common. However, how U6 2'-O-methylation is regulated remains largely unknown. Here we report that TFIP11, the human homolog of the yeast spliceosome disassembly factor Ntr1, localizes to nucleoli and Cajal Bodies and is essential for the 2'-O-methylation of U6. Mechanistically, we demonstrate that TFIP11 knockdown reduces the association of U6 snRNA with fibrillarin and associated snoRNAs, therefore altering U6 2'-O-methylation. We show U6 snRNA hypomethylation is associated with changes in assembly of the U4/U6.U5 tri-snRNP leading to defects in spliceosome assembly and alterations in splicing fidelity. Strikingly, this function of TFIP11 is independent of the RNA helicase DHX15, its known partner in yeast. In sum, our study demonstrates an unrecognized function for TFIP11 in U6 snRNP modification and U4/U6.U5 tri-snRNP assembly, identifying TFIP11 as a critical spliceosome assembly regulator.

Integrative structural dynamics probing of the conformational heterogeneity in synaptosomal-associated protein 25

SNAP-25 (synaptosomal-associated protein of 25 kDa) is a prototypical intrinsically disordered protein (IDP) that is unstructured by itself but forms coiled-coil helices in the SNARE complex. With high conformational heterogeneity, detailed structural dynamics of unbound SNAP-25 remain elusive. Here, we report an integrative method to probe the structural dynamics of SNAP-25 by combining replica-exchange discrete molecular dynamics (rxDMD) simulations and label-based experiments at ensemble and single-molecule levels. The rxDMD simulations systematically characterize the coil-to-molten globular transition and reconstruct structural ensemble consistent with prior ensemble experiments. Label-based experiments using Förster resonance energy transfer and double electron-electron resonance further probe the conformational dynamics of SNAP-25. Agreements between simulations and experiments under both ensemble and single-molecule conditions allow us to assign specific helix-coil transitions in SNAP-25 that occur in submillisecond timescales and potentially play a vital role in forming the SNARE complex. We expect that this integrative approach may help further our understanding of IDPs.

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.

Accelerating Reaction Rates of Biomolecules by Using Shear Stress in Artificial Capillary Systems

Biomimetics is a design principle within chemistry, biology, and engineering, but chemistry biomimetic approaches have been generally limited to emulating nature's chemical toolkit while emulation of nature's physical toolkit has remained largely unexplored. To begin to explore this, we designed biophysically mimetic microfluidic reactors with characteristic length scales and shear stresses observed within capillaries. We modeled the effect of shear with molecular dynamics studies and showed that this induces specific normally buried residues to become solvent accessible. We then showed using kinetics experiments that rates of reaction of these specific residues in fact increase in a shear-dependent fashion. We applied our results in the creation of a new microfluidic approach for the multidimensional study of cysteine biomarkers. Finally, we used our approach to establish dissociation of the therapeutic antibody trastuzumab in a reducing environment. Our results have implications for the efficacy of existing therapeutic antibodies in blood plasma as well as suggesting in general that biophysically mimetic chemistry is exploited in biology and should be explored as a research area.

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.

Hidden dynamic signatures drive substrate selectivity in the disordered phosphoproteome

The discovery that more than 40% of the eukaryotic proteome is intrinsically disordered, and that these disordered segments are enriched in phosphorylation sites, suggests that conformational heterogeneity may be important to kinase selectivity. Indeed, phosphorylation prediction programs reliant on classic notions of conserved sequence information (i.e., “vertical information”) are only partially effective. We find that the conformational equilibrium of the phosphorylatable site, whose information is embedded in sequence-averaged energetic and structural properties of the protein (i.e., “horizontal information”), plays a major role in distinguishing phosphorylatable versus nonphosphorylatable sites. In fact, employing both horizontal and vertical information produces a state-of-the-art phosphorylation predictor, wherein the conformational equilibrium of the disordered chain is the dominant contributor.

Phosphorylation sites are hyperabundant in the eukaryotic disordered proteome, suggesting that conformational fluctuations play a major role in determining to what extent a kinase interacts with a particular substrate. In biophysical terms, substrate selectivity may be determined not just by the structural–chemical complementarity between the kinase and its protein substrates but also by the free energy difference between the conformational ensembles that are, or are not, recognized by the kinase. To test this hypothesis, we developed a statistical-thermodynamics-based informatics framework, which allows us to probe for the contribution of equilibrium fluctuations to phosphorylation, as evaluated by the ability to predict Ser/Thr/Tyr phosphorylation sites in the disordered proteome. Essential to this framework is a decomposition of substrate sequence information into two types: vertical information encoding conserved kinase specificity motifs and horizontal information encoding substrate conformational equilibrium that is embedded, but often not apparent, within position-specific conservation patterns. We find not only that conformational fluctuations play a major role but also that they are the dominant contribution to substrate selectivity. In fact, the main substrate classifier distinguishing selectivity is the magnitude of change in local compaction of the disordered chain upon phosphorylation of these mostly singly phosphorylated sites. In addition to providing fundamental insights into the consequences of phosphorylation across the proteome, our approach provides a statistical-thermodynamic strategy for partitioning any sequence-based search into contributions from structural–chemical complementarity and those from changes in conformational equilibrium.

Identification of an α-MoRF in the Intrinsically Disordered Region of the Escargot Transcription Factor

Molecular recognition features (MoRFs) are common in intrinsically disordered proteins (IDPs) and intrinsically disordered regions (IDRs). MoRFs are in constant order-disorder structural transitions and adopt well-defined structures once they are bound to their targets. Here, we study Escargot (Esg), a transcription factor in Drosophila melanogaster that regulates multiple cellular functions, and consists of a disordered N-terminal domain and a group of zinc fingers at its C-terminal domain. We analyzed the N-terminal domain of Esg with disorder predictors and identified a region of 45 amino acids with high probability to form ordered structures, which we named S2. Through 54 μs of molecular dynamics (MD) simulations using CHARMM36 and implicit solvent (generalized Born/surface area (GBSA)), we characterized the conformational landscape of S2 and found an α-MoRF of ∼16 amino acids stabilized by key contacts within the helix. To test the importance of these contacts in the stability of the α-MoRF, we evaluated the effect of point mutations that would impair these interactions, running 24 μs of MD for each mutation. The mutations had mild effects on the MoRF, and in some cases, led to gain of residual structure through long-range contacts of the α-MoRF and the rest of the S2 region. As this could be an effect of the force field and solvent model we used, we benchmarked our simulation protocol by carrying out 32 μs of MD for the (AAQAA)₃ peptide. The results of the benchmark indicate that the global amount of helix in shorter peptides like (AAQAA)₃ is reasonably predicted. Careful analysis of the runs of S2 and its mutants suggests that the mutation to hydrophobic residues may have nucleated long-range hydrophobic and aromatic interactions that stabilize the MoRF. Finally, we have identified a set of residues that stabilize an α-MoRF in a region still without functional annotations in Esg.

Predicting Conformational Properties of Intrinsically Disordered Proteins from Sequence

Intrinsically disordered proteins (IDPs) can adopt a range of conformations from globules to swollen coils. This large range of conformational preferences for different IDPs raises the question of how conformational preferences are encoded by sequence. Global compositional features of a sequence such as the fraction of charged residues and the net charge per residue engender certain conformational biases. However, more specific sequence features such as the patterning of oppositely charged residues, expansion driving residues, or residues that can undergo posttranslational modifications can also influence the conformational ensembles of an IDP. Here, we outline how to calculate important global compositional features and patterning metrics that can be used to classify IDPs into different conformational classes and predict relative changes in conformation for sequences with the same amino acid composition. Although increased effort has been devoted to determining conformational properties of IDPs in recent years, quantitative predictions of conformation directly from sequence remain difficult and often inaccurate. Thus, if quantitative predictions of conformational properties are desired, then sequence-specific simulations must be performed.

Impact of Heat on Coil Hydrodynamic Size Yields the Energetics of Denatured State Conformational Bias

Conformational equilibria in the protein denatured state have key roles regulating folding, stability, and function. The extent of conformational bias in the protein denatured state under folding conditions, however, has thus far proven elusive to quantify, particularly with regard to its sequence dependence and energetic character. To better understand the structural preferences of the denatured state, we analyzed both the sequence dependence to the mean hydrodynamic size of disordered proteins in water and the impact of heat on the coil dimensions, showing that the sequence dependence and thermodynamic energies associated with intrinsic biases for the α and polyproline II (PPII) backbone conformations can be obtained. Experiments that evaluate how the hydrodynamic size changes with compositional changes in the protein reveal amino acid specific preferences for PPII that are in good quantitative agreement with calorimetry-measured values from unfolded peptides and those inferred by survey of the protein coil library. At temperatures above 25 °C, the denatured state follows the predictions of a PPII-dominant ensemble. Heat effects on coil hydrodynamic size indicate the α bias is comparable to the PPII bias at cold temperatures. Though historically thought to give poor resolution to structural details, the hydrodynamic size of the unfolded state is found to be an effective reporter on the extent of the biases for the α and PPII backbone conformations.

Integrating disorder in globular multidomain proteins: Fuzzy sensors and the role of SH3 domains

Intrinsically disordered proteins represent about one third of eukaryotic proteins. An additional third correspond to proteins containing folded domains as well as large intrinsically disordered regions (IDR). While IDRs may represent functionally autonomous domains, in some instances it has become clear that they provide a new layer of regulation for the activity displayed by the folded domains. The sensitivity of the conformational ensembles defining the properties of IDR to small changes in the cellular environment and the capacity to modulate this response through post-translational modifications makes IDR ideal sensors enabling continuous, integrative responses to complex cellular inputs. Folded domains (FD), on the other hand, are ideal effectors, e.g. by catalyzing enzymatic reactions or participating in binary on/off switches. In this perspective review we discuss the possible role of intramolecular fuzzy complexes to integrate the very different dynamic scales of IDR and FD, inspired on the recent observations of such dynamic complexes in Src family kinases, and we explore the possible general role of the SH3 domains connecting IDRs and FD.

Sequence Reversal Prevents Chain Collapse and Yields Heat-Sensitive Intrinsic Disorder

Sequence patterns of charge, hydrophobicity, hydrogen bonding, and other amino acid physicochemical properties contribute to mechanisms of protein folding, but how sequence composition and patterns influence the conformational dynamics of the denatured state ensemble is not fully understood. To investigate structure-sequence relationships in the denatured state, we reversed the sequence of staphylococcal nuclease and characterized its structure, thermodynamic character, and hydrodynamic radius using circular dichroism spectroscopy, dynamic light scattering, analytical ultracentrifugation, and size-exclusion chromatography as a function of temperature. The macromolecular size of "Retro-nuclease" is highly expanded in solution with characteristics similar to biological intrinsically disordered proteins. In contradistinction to a disordered state, Retro-nuclease exhibits a broad sigmoid transition of its hydrodynamic dimensions as temperature is increased, indicating a thermodynamically controlled compaction. Counterintuitively, the magnitude of these temperature-induced hydrodynamic changes exceed that observed from thermal denaturation of folded unaltered staphylococcal nuclease. Undetectable by calorimetry and intrinsic tryptophan fluorescence, the lack of heat capacity or fluorescence changes throughout the thermal transition indicate canonical hydrophobic collapse did not drive the Retro-nuclease structural transitions. Temperature-dependent circular dichroism spectroscopy performed on Retro-nuclease and computer simulations correlate to temperature sensitivity in the intrinsic sampling of backbone conformations for polyproline II and α-helix. The experimental results indicate a role for sequence direction in mediating the collapse of the polypeptide chain, whereas the simulation trends illustrate the generality of the observed heat effects on disordered protein structure.

Side chain electrostatic interactions and pH-dependent expansion of the intrinsically disordered, highly acidic carboxyl-terminus of γ-tubulin

Intramolecular electrostatic attraction and repulsion strongly influence the conformational sampling of intrinsically disordered proteins and domains (IDPs). In order to better understand this complex relationship, we have used nuclear magnetic resonance to measure side chain pK_a values and pH-dependent translational diffusion coefficients for the unstructured and highly acidic carboxyl-terminus of γ-tubulin (γ-CT), providing insight into how the net charge of an IDP relates to overall expansion or collapse of the conformational ensemble. Many of the pK_a values in the γ-CT are shifted upward by 0.3-0.4 units and exhibit negatively cooperative ionization pH profiles, likely due to the large net negative charge that accumulates on the molecule as the pH is raised. pK_a shifts of this magnitude correspond to electrostatic interaction energies between the affected residues and the rest of the charged molecule that are each on the order of 1 kcal mol^-1 . Diffusion of the γ-CT slowed with increasing net charge, indicative of an expanding hydrodynamic radius (r_H ). The degree of expansion agreed quantitatively with what has been seen from comparisons of IDPs with different charge content, yielding the general trend that every 0.1 increase in relative charge (|Q|/res) produces a roughly 5% increase in r_H . While γ-CT pH titration data followed this trend nearly perfectly, there were substantially larger deviations for the database of different IDP sequences. This suggests that other aspects of an IDP's primary amino acid sequence beyond net charge influence the sensitivity of r_H to electrostatic interactions.

Metal-binding polymorphism in late embryogenesis abundant protein AtLEA4-5, an intrinsically disordered protein

Late embryogenesis abundant (LEA) proteins accumulate in plants during adverse conditions and their main attributed function is to confer tolerance to stress. One of the deleterious effects of the adverse environment is the accumulation of metal ions to levels that generate reactive oxygen species, compromising the survival of cells. AtLEA4-5, a member of group 4 of LEAs in Arabidopsis, is an intrinsically disordered protein. It has been shown that their N-terminal region is able to undergo transitions to partially folded states and prevent the inactivation of enzymes. We have characterized metal ion binding to AtLEA4-5 by circular dichroism, electronic absorbance spectroscopy (UV–vis), electron paramagnetic resonance, dynamic light scattering, and isothermal titration calorimetry. The data shows that AtLEA4-5 contains a single binding site for Ni(II), while Zn(II) and Cu(II) have multiple binding sites and promote oligomerization. The Cu(II) interacts preferentially with histidine residues mostly located in the C-terminal region with moderate affinity and different coordination modes. These results and the lack of a stable secondary structure formation indicate that an ensemble of conformations remains accessible to the metal for binding, suggesting the formation of a fuzzy complex. Our results support the multifunctionality of LEA proteins and suggest that the C-terminal region of AtLEA4-5 could be responsible for antioxidant activity, scavenging metal ions under stress conditions while the N-terminal could function as a chaperone.

Order, Disorder, and Temperature-Driven Compaction in a Designed Elastin Protein

Artificial minielastin constructs have been designed that replicate the structure and function of natural elastins in a simpler context, allowing the NMR observation of structure and dynamics of elastin-like proteins with complete residue-specific resolution. We find that the alanine-rich cross-linking domains of elastin have a partially helical structure, but only when capped by proline-rich hydrophobic domains. We also find that the hydrophobic domains, composed of prominent 6-residue repeats VPGVGG and APGVGV found in natural elastins, appear random coil by both NMR chemical shift analysis and circular dichroism. However, these elastin hydrophobic domains exhibit structural bias for a dynamically disordered conformation that is neither helical nor β sheet with a degree of nonrandom structural bias which is dependent on residue type and position in the sequence. Another nonrandom-coil aspect of hydrophobic domain structure lies in the fact that, in contrast to other intrinsically disordered proteins, these hydrophobic domains retain a relatively condensed conformation whether attached to cross-linking domains or not. Importantly, these domains and the proteins containing them constrict with increasing temperature by up to 30% in volume without becoming more ordered. This property is often observed in nonbiological polymers and suggests that temperature-driven constriction is a new type of protein structural change that is linked to elastin's biological functions of coacervation-driven assembly and elastic recoil.

CIDER: Resources to Analyze Sequence-Ensemble Relationships of Intrinsically Disordered Proteins

Intrinsically disordered proteins and regions (IDPs) represent a large class of proteins that are defined by conformational heterogeneity and lack of persistent tertiary/secondary structure. IDPs play important roles in a range of biological functions, and their dysregulation is central to numerous diseases, including neurodegeneration and cancer. The conformational ensembles of IDPs are encoded by their amino acid sequences. Here, we present two computational tools that are designed to enable rapid and high-throughput analyses of a wide range of physicochemical properties encoded by IDP sequences. The first, CIDER, is a user-friendly webserver that enables rapid analysis of IDP sequences. The second, localCIDER, is a high-performance software package that enables a wide range of analyses relevant to IDP sequences. In addition to introducing the two packages, we demonstrate the utility of these resources using examples where sequence analysis offers biophysical insights.

The Proline/Glycine-Rich Region of the Biofilm Adhesion Protein Aap Forms an Extended Stalk that Resists Compaction

Staphylococcus epidermidis is one of the primary bacterial species responsible for healthcare-associated infections. The most significant virulence factor for S. epidermidis is its ability to form a biofilm, which renders the bacteria highly resistant to host immune responses and antibiotic action. Intercellular adhesion within the biofilm is mediated by the accumulation-associated protein (Aap), a cell wall-anchored protein that self-assembles in a zinc-dependent manner. The C-terminal portion of Aap contains a 135-aa-long, proline/glycine-rich region (PGR) that has not yet been characterized. The region contains a set of 18 nearly identical AEPGKP repeats. Analysis of the PGR using biophysical techniques demonstrated the region is a highly extended, intrinsically disordered polypeptide with unusually high polyproline type II helix propensity. In contrast to many intrinsically disordered polypeptides, there was a minimal temperature dependence of the global conformational state of PGR in solution as measured by analytical ultracentrifugation and dynamic light scattering. Furthermore, PGR was resistant to conformational collapse or α-helix formation upon the addition of the osmolyte trimethylamine N-oxide or the cosolvent 2,2,2-trifluoroethanol. Collectively, these results suggest PGR functions as a resilient, extended stalk that projects the rest of Aap outward from the bacterial cell wall, promoting intercellular adhesion between cells in the biofilm. This work sheds light on regions of low complexity often found near the attachment point of bacterial cell wall-anchored proteins.

Intrinsic α helix propensities compact hydrodynamic radii in intrinsically disordered proteins

Proteins that lack tertiary stability under normal conditions, known as intrinsically disordered, exhibit a wide range of biological activities. Molecular descriptions for the biology of intrinsically disordered proteins (IDPs) consequently rely on disordered structural models, which in turn require experiments that assess the origins to structural features observed. For example, while hydrodynamic size is mostly insensitive to sequence composition in chemically denatured proteins, IDPs show strong sequence-specific effects in the hydrodynamic radius (R_h ) when measured under normal conditions. To investigate sequence-modulation of IDP R_h , disordered ensembles generated by a hard sphere collision model modified with a structure-based parameterization of the solution energetics were used to parse the contributions of net charge, main chain dihedral angle bias, and excluded volume on hydrodynamic size. Ensembles for polypeptides 10-35 residues in length were then used to establish power-law scaling relationships for comparison to experimental R_h from 26 IDPs. Results showed the expected outcomes of increased hydrodynamic size from increases in excluded volume and net charge, and compaction from chain-solvent interactions. Chain bias representing intrinsic preferences for α helix and polyproline II (PP_II ), however, modulated R_h with intricate dependence on the simulated propensities. PP_II propensities at levels expected in IDPs correlated with heightened R_h sensitivity to even weak α helix propensities, indicating bias for common (φ, ψ) are important determinants of hydrodynamic size. Moreover, data show that IDP R_h can be predicted from sequence with good accuracy from a small set of physicochemical properties, namely intrinsic conformational propensities and net charge. Proteins 2017; 85:296-311. © 2016 Wiley Periodicals, Inc.

NMR spectroscopic studies of a TAT-derived model peptide in imidazolium-based ILs: influence on chemical shifts and the cis/trans equilibrium state

The impact of ionic liquids on the chemical shifts and the cis/trans equilibrium state of a model peptide was systematically investigated by NMR spectroscopy.