Rigidity,Secondary Structure,and the Universality of the Boson Peak in Proteins

Complementary neutron- and light-scattering results on nine proteins and amino acids reveal the role of rigidity and secondary structure in determining the time- and lengthscales of low-frequency collective vibrational dynamics in proteins. These dynamics manifest in a spectral feature, known as the boson peak (BP), which is common to all disordered materials. We demonstrate that BP position scales systematically with structural motifs, reflecting local rigidity: disordered proteins appear softer than α-helical proteins; which are softer than β-sheet proteins. Our analysis also reveals a universal spectral shape of the BP in proteins and amino acid mixtures; superimposable on the shape observed in typical glasses. Uniformity in the underlying physical mechanism, independent of the specific chemical composition, connects the BP vibrations to nanometer-scale heterogeneities, providing an experimental benchmark for coarse-grained simulations, structure/rigidity relationships, and engineering of proteins for novel applications.     

H and C NMR studies on the temperature-dependent water and protein dynamics in hydrated elastin,myoglobin and collagen

²H NMR spin-lattice relaxation and line-shape analyses are performed to study the temperature-dependent dynamics of water in the hydration shells of myoglobin, elastin, and collagen. The results show that the dynamical behaviors of the hydration waters are similar for these proteins when using comparable hydration levels of h = 0.25–0.43. Since water dynamics is characterized by strongly nonexponential correlation functions, we use a Cole–Cole spectral density for spin-lattice relaxation analysis, leading to correlation times, which are in nice agreement with results for the main dielectric relaxation process observed for various proteins in the literature. The temperature dependence can roughly be described by an Arrhenius law, with the possibility of a weak crossover in the vicinity of 220 K. Near ambient temperatures, the results substantially depend on the exact shape of the spectral density so that deviations from an Arrhenius behavior cannot be excluded in the high-temperature regime. However, for the studied proteins, the data give no evidence for the existence of a sharp fragile-to-strong transition reported for lysozyme at about 220 K. Line-shape analysis reveals that the mechanism for the rotational motion of hydration waters changes in the vicinity of 220 K. For myoglobin, we observe an isotropic motion at high temperatures and an anisotropic large-amplitude motion at low temperatures. Both mechanisms coexist in the vicinity of 220 K. ¹³C CP MAS spectra show that hydration results in enhanced elastin dynamics at ambient temperatures, where the enhancement varies among different amino acids. Upon cooling, the enhanced mobility decreases. Comparison of ²H and ¹³C NMR data reveals that the observed protein dynamics is slower than the water dynamics.     

DNP-Enhanced MAS NMR: A Tool to Snapshot Conformational Ensembles of α-Synuclein in Different States

Intrinsically disordered proteins dynamically sample a wide conformational space and therefore do not adopt a stable and defined three-dimensional conformation. The structural heterogeneity is related to their proper functioning in physiological processes. Knowledge of the conformational ensemble is crucial for a complete comprehension of this kind of proteins. We here present an approach that utilizes dynamic nuclear polarization-enhanced solid-state NMR spectroscopy of sparsely isotope-labeled proteins in frozen solution to take snapshots of the complete structural ensembles by exploiting the inhomogeneously broadened line-shapes. We investigated the intrinsically disordered protein α-synuclein (α-syn), which plays a key role in the etiology of Parkinson’s disease, in three different physiologically relevant states. For the free monomer in frozen solution we could see that the so-called “random coil conformation” consists of α-helical and β-sheet-like conformations, and that secondary chemical shifts of neighboring amino acids tend to be correlated, indicative of frequent formation of secondary structure elements. Based on these results, we could estimate the number of disordered regions in fibrillar α-syn as well as in α-syn bound to membranes in different protein-to-lipid ratios. Our approach thus provides quantitative information on the propensity to sample transient secondary structures in different functional states. Molecular dynamics simulations rationalize the results.     

Solution Structures,Dynamics, and Ice Growth Inhibitory Activity of Peptide Fragments Derived from an Antarctic Yeast Protein

Exotic functions of antifreeze proteins (AFP) and antifreeze glycopeptides (AFGP) have recently been attracted with much interest to develop them as commercial products. AFPs and AFGPs inhibit ice crystal growth by lowering the water freezing point without changing the water melting point. Our group isolated the Antarctic yeast Glaciozyma antarctica that expresses antifreeze protein to assist it in its survival mechanism at sub-zero temperatures. The protein is unique and novel, indicated by its low sequence homology compared to those of other AFPs. We explore the structure-function relationship of G. antarctica AFP using various approaches ranging from protein structure prediction, peptide design and antifreeze activity assays, nuclear magnetic resonance (NMR) studies and molecular dynamics simulation. The predicted secondary structure of G. antarctica AFP shows several α-helices, assumed to be responsible for its antifreeze activity. We designed several peptide fragments derived from the amino acid sequences of α-helical regions of the parent AFP and they also showed substantial antifreeze activities, below that of the original AFP. The relationship between peptide structure and activity was explored by NMR spectroscopy and molecular dynamics simulation. NMR results show that the antifreeze activity of the peptides correlates with their helicity and geometrical straightforwardness. Furthermore, molecular dynamics simulation also suggests that the activity of the designed peptides can be explained in terms of the structural rigidity/flexibility, i.e., the most active peptide demonstrates higher structural stability, lower flexibility than that of the other peptides with lower activities, and of lower rigidity. This report represents the first detailed report of downsizing a yeast AFP into its peptide fragments with measurable antifreeze activities.     

Differences in the number of intrinsically disordered regions between yeast duplicated proteins, and their relationship with functional divergence

Background

Intrinsically disordered regions are enriched in short interaction motifs that play a critical role in many protein-protein interactions. Since new short interaction motifs may easily evolve, they have the potential to rapidly change protein interactions and cellular signaling. In this work we examined the dynamics of gain and loss of intrinsically disordered regions in duplicated proteins to inspect if changes after genome duplication can create functional divergence. For this purpose we used Saccharomyces cerevisiae and the outgroup species Lachancea kluyveri.

Principal Findings

We find that genes duplicated as part of a genome duplication (ohnologs) are significantly more intrinsically disordered than singletons (p<2.2^e-16, Wilcoxon), reflecting a preference for retaining intrinsically disordered proteins in duplicate. In addition, there have been marked changes in the extent of intrinsic disorder following duplication. A large number of duplicated genes have more intrinsic disorder than their L. kluyveri ortholog (29% for duplicates versus 25% for singletons) and an even greater number have less intrinsic disorder than the L. kluyveri ortholog (37% for duplicates versus 25% for singletons). Finally, we show that the number of physical interactions is significantly greater in the more intrinsically disordered ohnolog of a pair (p = 0.003, Wilcoxon).

Conclusion

This work shows that intrinsic disorder gain and loss in a protein is a mechanism by which a genome can also diverge and innovate. The higher number of interactors for proteins that have gained intrinsic disorder compared with their duplicates may reflect the acquisition of new interaction partners or new functional roles.     

The Middle Region of an HP1-binding Protein, HP1-BP74, Associates with Linker DNA at the Entry/Exit Site of Nucleosomal DNA

In higher eukaryotic cells, DNA molecules are present as chromatin fibers, complexes of DNA with various types of proteins; chromatin fibers are highly condensed in metaphase chromosomes during mitosis. Although the formation of the metaphase chromosome structure is essential for the equal segregation of replicated chromosomal DNA into the daughter cells, the mechanism involved in the organization of metaphase chromosomes is poorly understood. To identify proteins involved in the formation and/or maintenance of metaphase chromosomes, we examined proteins that dissociated from isolated human metaphase chromosomes by 0.4 m NaCl treatment; this treatment led to significant chromosome decondensation, but the structure retained the core histones. One of the proteins identified, HP1-BP74 (heterochromatin protein 1-binding protein 74), composed of 553 amino acid residues, was further characterized. HP1-BP74 middle region (BP74Md), composed of 178 amino acid residues (Lys⁹⁷–Lys²⁷⁴), formed a chromatosome-like structure with reconstituted mononucleosomes and protected the linker DNA from micrococcal nuclease digestion by ∼25 bp. The solution structure determined by NMR revealed that the globular domain (Met¹⁵³–Thr²³⁷) located within BP74Md possesses a structure similar to that of the globular domain of linker histones, which underlies its nucleosome binding properties. Moreover, we confirmed that BP74Md and full-length HP1-BP74 directly binds to HP1 (heterochromatin protein 1) and identified the exact sites responsible for this interaction. Thus, we discovered that HP1-BP74 directly binds to HP1, and its middle region associates with linker DNA at the entry/exit site of nucleosomal DNA in vitro.     

Highly Disordered Amyloid-β Monomer Probed by Single-Molecule FRET and MD Simulation

Monomers of amyloid-β (Aβ) protein are known to be disordered, but there is considerable controversy over the existence of residual or transient conformations that can potentially promote oligomerization and fibril formation. We employed single-molecule Förster resonance energy transfer (FRET) spectroscopy with site-specific dye labeling using an unnatural amino acid and molecular dynamics simulations to investigate conformations and dynamics of Aβ isoforms with 40 (Aβ40) and 42 residues (Aβ42). The FRET efficiency distributions of both proteins measured in phosphate-buffered saline at room temperature show a single peak with very similar FRET efficiencies, indicating there is apparently only one state. 2D FRET efficiency-donor lifetime analysis reveals, however, that there is a broad distribution of rapidly interconverting conformations. Using nanosecond fluorescence correlation spectroscopy, we measured the timescale of the fluctuations between these conformations to be ～35 ns, similar to that of disordered proteins. These results suggest that both Aβ40 and Aβ42 populate an ensemble of rapidly reconfiguring unfolded states, with no long-lived conformational state distinguishable from that of the disordered ensemble. To gain molecular-level insights into these observations, we performed molecular dynamics simulations with a force field optimized to describe disordered proteins. We find, as in experiments, that both peptides populate configurations consistent with random polymer chains, with the vast majority of conformations lacking significant secondary structure, giving rise to very similar ensemble-averaged FRET efficiencies.     

Predicting mostly disordered proteins by using structure-unknown protein data

Background  

Predicting intrinsically disordered proteins is important in structural biology because they are thought to carry out various cellular functions even though they have no stable three-dimensional structure. We know the structures of far more ordered proteins than disordered proteins. The structural distribution of proteins in nature can therefore be inferred to differ from that of proteins whose structures have been determined experimentally. We know many more protein sequences than we do protein structures, and many of the known sequences can be expected to be those of disordered proteins. Thus it would be efficient to use the information of structure-unknown proteins in order to avoid training data sparseness. We propose a novel method for predicting which proteins are mostly disordered by using spectral graph transducer and training with a huge amount of structure-unknown sequences as well as structure-known sequences.     

Dynamics of the intrinsically disordered protein NUPR1 in isolation and in its fuzzy complexes with DNA and prothymosin α

Intrinsically disordered proteins (IDPs) explore diverse conformations in their free states and, a few of them, also in their molecular complexes. This functional plasticity is essential for the function of IDPs, although their dynamics in both free and bound states is poorly understood. NUPR1 is a protumoral multifunctional IDP, activated during the acute phases of pancreatitis. It interacts with DNA and other IDPs, such as prothymosin α (ProTα), with dissociation constants of ~0.5 μM, and a 1:1 stoichiometry. We studied the structure and picosecond-to-nanosecond (ps-ns) dynamics by using both NMR and SAXS in: (i) isolated NUPR1; (ii) the NUPR1/ProTα complex; and (iii) the NUPR1/double stranded (ds) GGGCGCGCCC complex. Our SAXS findings show that NUPR1 remained disordered when bound to either partner, adopting a worm-like conformation; the fuzziness of bound NUPR1 was also pinpointed by NMR. Residues with the largest values of the relaxation rates (R₁, R_1ρ, R₂ and η_xy), in the free and bound species, were mainly clustered around the 30s region of the sequence, which agree with one of the protein hot-spots already identified by site-directed mutagenesis. Not only residues in this region had larger relaxation rates, but they also moved slower than the rest of the molecule, as indicated by the reduced spectral density approach (RSDA). Upon binding, the energy landscape of NUPR1 was not funneled down to a specific, well-folded conformation, but rather its backbone flexibility was kept, with distinct motions occurring at the hot-spot region.     

Nervous system myelin in the electric ray, Torpedo marmorata: Morphological characterization of the membrane and biochemical analysis of its protein components

In Torpedo, PNS as well as CNS myelines are characterized by clearly separated double intraperiod lines. CNS myelin of Torpedo contains two glycosylated hydrophobic proteins labelled T1 (25,800 Da1) and T2 (29,700 Da1), and two basic proteins BP1 and BP2, migrating like mammalian large basic protein (BP2) and pre-small basic protein (BP1) (Barbarese et al., 1977). PNS myelin of Torpedo carries only BP1 and is characterized by a closely spaced doublet of the glycosylated hydrophobic proteins Con A+ (29,700 Da1) and Con A? (31,000 Da1); the latter does not bind Concanavalin A. These glycosylated proteins (T1, T2, Con A+, Con A?) contain mannose, N-acetylglucosamine and galactose, but lack fucose and sialic acids. They have isoleucine at their amino terminus. They bind anti-rat PNS myelin P₀ antibodies but do not react with anti-rat CNS myelin PLP antibodies. Limited proteolyses of isolated proteins suggest sequence homologies between T1 and T2, and possibly between Con A+ and Con A?. The two basic proteins BP1 and BP2 bind antibodies directed against human myelin basic protein. All Torpedo myelin proteins electrofocus in pH regions characteristic of their mammalian counterparts.     

Markov models of amino acid substitution to study proteins with intrinsically disordered regions

Background

Intrinsically disordered proteins (IDPs) or proteins with disordered regions (IDRs) do not have a well-defined tertiary structure, but perform a multitude of functions, often relying on their native disorder to achieve the binding flexibility through changing to alternative conformations. Intrinsic disorder is frequently found in all three kingdoms of life, and may occur in short stretches or span whole proteins. To date most studies contrasting the differences between ordered and disordered proteins focused on simple summary statistics. Here, we propose an evolutionary approach to study IDPs, and contrast patterns specific to ordered protein regions and the corresponding IDRs.

Results

Two empirical Markov models of amino acid substitutions were estimated, based on a large set of multiple sequence alignments with experimentally verified annotations of disordered regions from the DisProt database of IDPs. We applied new methods to detect differences in Markovian evolution and evolutionary rates between IDRs and the corresponding ordered protein regions. Further, we investigated the distribution of IDPs among functional categories, biochemical pathways and their preponderance to contain tandem repeats.

Conclusions

We find significant differences in the evolution between ordered and disordered regions of proteins. Most importantly we find that disorder promoting amino acids are more conserved in IDRs, indicating that in some cases not only amino acid composition but the specific sequence is important for function. This conjecture is also reinforced by the observation that for of our data set IDRs evolve more slowly than the ordered parts of the proteins, while we still support the common view that IDRs in general evolve more quickly. The improvement in model fit indicates a possible improvement for various types of analyses e.g. de novo disorder prediction using a phylogenetic Hidden Markov Model based on our matrices showed a performance similar to other disorder predictors.     

Intrinsically Disordered Regions May Lower the Hydration Free Energy in Proteins: A Case Study of Nudix Hydrolase in the Bacterium Deinococcus radiodurans

The proteome of the radiation- and desiccation-resistant bacterium D. radiodurans features a group of proteins that contain significant intrinsically disordered regions that are not present in non-extremophile homologues. Interestingly, this group includes a number of housekeeping and repair proteins such as DNA polymerase III, nudix hydrolase and rotamase. Here, we focus on a member of the nudix hydrolase family from D. radiodurans possessing low-complexity N- and C-terminal tails, which exhibit sequence signatures of intrinsic disorder and have unknown function. The enzyme catalyzes the hydrolysis of oxidatively damaged and mutagenic nucleotides, and it is thought to play an important role in D. radiodurans during the recovery phase after exposure to ionizing radiation or desiccation. We use molecular dynamics simulations to study the dynamics of the protein, and study its hydration free energy using the GB/SA formalism. We show that the presence of disordered tails significantly decreases the hydration free energy of the whole protein. We hypothesize that the tails increase the chances of the protein to be located in the remaining water patches in the desiccated cell, where it is protected from the desiccation effects and can function normally. We extrapolate this to other intrinsically disordered regions in proteins, and propose a novel function for them: intrinsically disordered regions increase the “surface-properties” of the folded domains they are attached to, making them on the whole more hydrophilic and potentially influencing, in this way, their localization and cellular activity.     

Short peptides as predictors for the structure of polyarginine sequences in disordered proteins

Intrinsically disordered proteins and intrinsically disordered regions are frequently enriched in charged amino acids. Intrinsically disordered regions are regularly involved in important biological processes in which one or more charged residues is the driving force behind a protein-biomolecule interaction. Several lines of experimental and computational evidence suggest that polypeptides and proteins that carry high net charges have a high preference for extended conformations with average end-to-end distances exceeding expectations for self-avoiding random coils. Here, we show that charged arginine residues even in short glycine-capped model peptides (GRRG and GRRRG) significantly affect the conformational propensities of each other when compared with the intrinsic propensities of a mostly unperturbed arginine in the tripeptide GRG. A conformational analysis based on experimentally determined J-coupling constants from heteronuclear NMR spectroscopy and amide I′ band profiles from vibrational spectroscopy reveals that nearest-neighbor interactions stabilize extended β-strand conformations at the expense of polyproline II and turn conformations. The results from molecular dynamics simulations with a CHARMM36m force field and TIP3P water reproduce our results only to a limited extent. The use of the Ramachandran distribution of the central residue of GRRRG in a calculation of end-to-end distances of polyarginines of different length yielded the expected power law behavior. The scaling coefficient of 0.66 suggests that such peptides would be more extended than predicted by a self-avoiding random walk. Our findings thus support in principle theoretical predictions.     

Does Lack of Secondary Structure Imply Intrinsic Disorder in Proteins? A Sequence Analysis

Intrinsically disordered proteins (IDPs)/protein regions (IDPRs) lack unique three-dimensional structure at the level of secondary and/or tertiary structure and are represented as an ensemble of interchanging conformations. To investigate the role of presence/absence of secondary structures in promoting intrinsic disorder in proteins, a comparative sequence analysis of IDPs, IDPRs and proteins with minimal secondary structures (less than 5%) is required. A sequence analysis reveals proteins with minimal secondary structure content have high mean net positive charge, low mean net hydrophobicity and low sequence complexity. Interestingly, analysis of the relative local electrostatic interactions reveal that an increase in the relative repulsive interactions between amino acids separated by three or four residues lead to either loss of secondary structure or intrinsic disorder. IDPRs show increase in both local negative-negative and positive-positive repulsive interactions. While IDPs show a marked increase in the local negative-negative interactions, proteins with minimal secondary structure depict an increase in the local positive-positive interactions. IDPs and IDPRs are enriched in D, E and Q residues, while proteins with minimal secondary structure are depleted of these residues. Proteins with minimal secondary structures have higher content of G and C, while IDPs and IDPRs are depleted of these residues. These results confirm that proteins with minimal secondary structure have a distinctly different propensity for charge, hydrophobicity, specific amino acids and local electrostatic interactions as compared to IDPs/IDPRs. Thus we conclude that lack of secondary structure may be a necessary but not a sufficient condition for intrinsic disorder in proteins.     

On the Role of Aggregation Prone Regions in Protein Evolution,Stability, and Enzymatic Catalysis: Insights from Diverse Analyses

The various roles that aggregation prone regions (APRs) are capable of playing in proteins are investigated here via comprehensive analyses of multiple non-redundant datasets containing randomly generated amino acid sequences, monomeric proteins, intrinsically disordered proteins (IDPs) and catalytic residues. Results from this study indicate that the aggregation propensities of monomeric protein sequences have been minimized compared to random sequences with uniform and natural amino acid compositions, as observed by a lower average aggregation propensity and fewer APRs that are shorter in length and more often punctuated by gate-keeper residues. However, evidence for evolutionary selective pressure to disrupt these sequence regions among homologous proteins is inconsistent. APRs are less conserved than average sequence identity among closely related homologues (≥80% sequence identity with a parent) but APRs are more conserved than average sequence identity among homologues that have at least 50% sequence identity with a parent. Structural analyses of APRs indicate that APRs are three times more likely to contain ordered versus disordered residues and that APRs frequently contribute more towards stabilizing proteins than equal length segments from the same protein. Catalytic residues and APRs were also found to be in structural contact significantly more often than expected by random chance. Our findings suggest that proteins have evolved by optimizing their risk of aggregation for cellular environments by both minimizing aggregation prone regions and by conserving those that are important for folding and function. In many cases, these sequence optimizations are insufficient to develop recombinant proteins into commercial products. Rational design strategies aimed at improving protein solubility for biotechnological purposes should carefully evaluate the contributions made by candidate APRs, targeted for disruption, towards protein structure and activity.     

Combination of 15N reverse labeling and afterglow spectroscopy for assigning membrane protein spectra by magic-angle-spinning solid-state NMR: application to the multidrug resistance protein EmrE

Magic-angle-spinning (MAS) solid-state NMR spectroscopy has emerged as a viable method to characterize membrane protein structure and dynamics. Nevertheless, the spectral resolution for uniformly labeled samples is often compromised by redundancy of the primary sequence and the presence of helical secondary structure that results in substantial resonance overlap. The ability to simplify the spectrum in order to obtain unambiguous site-specific assignments is a major bottleneck for structure determination. To address this problem, we used a combination of ¹⁵N reverse labeling, afterglow spectroscopic techniques, and frequency-selective dephasing experiments that dramatically improved the ability to resolve peaks in crowded spectra. This was demonstrated using the polytopic membrane protein EmrE, an efflux pump involved in multidrug resistance. Residues preceding the ¹⁵N reverse labeled amino acid were imaged using a 3D NCOCX afterglow experiment and those following were recorded using a frequency-selective dephasing experiment. Our approach reduced the spectral congestion and provided a sensitive way to obtain chemical shift assignments for a membrane protein where no high-resolution structure is available. This MAS methodology is widely applicable to the study of other polytopic membrane proteins in functional lipid bilayer environments.     

Structural effects of the binding of GTP to the wild-type and oncogenic forms of theras-gene-encoded p21 proteins

Molecular dynamics calculations have been performed to determine the average structures ofras-gene-encoded p21 proteins bound to GTP, i.e., the normal (wild-type) protein and two oncogenic forms of this protein, the Val 12- and Leu 61-p21 proteins. We find that the average structures for all of these proteins exhibit low coordinate fluctuations (which are highest for the normal protein), indicating convergence to specific structures. From previous dynamics calculations of the average structures of these proteins bound to GDP, major regional differences were found among these proteins (Monacoet al. (1995),J. Protein Chem., in press). We now find that the average structures of the oncogenic proteins are more similar to one another when the proteins are bound to GTP than when they are bound to GDP (Monacoet al. (1995),J. Protein Chem., in press). However, they still differ in structureat specific amino acid residues rather than in whole regions, in contradistinction to the results found for the p21-GDP complexes. Two exceptions are the regions 25–32, in anα-helical region, and 97–110. The two oncogenic (Val 12- and Leu 61-) proteins have similar structures which differ significantly in the region of residues 97–110. This region has recently been identified as being critical in the interaction of p21 with kinase target proteins. The differences in structure between the oncogenic proteins suggest the existence of more than one oncogenic form of the p21 protein that can activate different signaling pathways.     

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

The properties of disordered proteins are thought to depend on intrinsic conformational propensities for polyproline II (PP _II) structure. While intrinsic PP _II propensities have been measured for the common biological amino acids in short peptides, the ability of these experimentally determined propensities to quantitatively reproduce structural behavior in intrinsically disordered proteins (IDPs) has not been established. Presented here are results from molecular simulations of disordered proteins showing that the hydrodynamic radius (R _h) can be predicted from experimental PP _II propensities with good agreement, even when charge-based considerations are omitted. The simulations demonstrate that R _h and chain propensity for PP _II structure are linked via a simple power-law scaling relationship, which was tested using the experimental R _h of 22 IDPs covering a wide range of peptide lengths, net charge, and sequence composition. Charge effects on R _h were found to be generally weak when compared to PP _II effects on R _h. Results from this study indicate that the hydrodynamic dimensions of IDPs are evidence of considerable sequence-dependent backbone propensities for PP _II structure that qualitatively, if not quantitatively, match conformational propensities measured in peptides.