Conformational contagion in a protein: structural properties of a chameleon sequence

Dimer dissociation is a key energetic event in the fold-switch pathway of KaiB

Cyanobacteria possesses the simplest circadian clock, composed of three proteins that act as a phosphorylation oscillator: KaiA, KaiB, and KaiC. The timing of this oscillator is determined by the fold-switch of KaiB, a structural rearrangement of its C-terminal half that is accompanied by a change in the oligomerization state. During the day, KaiB forms a stable tetramer (gsKaiB), whereas it adopts a monomeric thioredoxin-like fold during the night (fsKaiB). Although the structures and functions of both native states are well studied, little is known about the sequence and structure determinants that control their structural interconversion. Here, we used confinement molecular dynamics (CCR-MD) and folding simulations using structure-based models to show that the dissociation of the gsKaiB dimer is a key energetic event for the fold-switch. Hydrogen-deuterium exchange mass spectrometry (HDXMS) recapitulates the local stability of protein regions reported by CCR-MD, with both approaches consistently indicating that the energy and backbone flexibility changes are solely associated with the region that fold-switches between gsKaiB and fsKaiB and that the localized regions that differentially stabilize gsKaiB also involve regions outside the dimer interface. Moreover, two mutants (R23C and R75C) previously reported to be relevant for altering the rhythmicity of the Kai clock were also studied by HDXMS. Particularly, R75C populates dimeric and monomeric states with a deuterium incorporation profile comparable to the one observed for fsKaiB, emphasizing the importance of the oligomerization state of KaiB for the fold-switch. These findings suggest that the information necessary to control the rhythmicity of the cyanobacterial biological clock is, to a great extent, encoded within the KaiB sequence.

Turncoat Polypeptides: We Adapt to Our Environment

A common interpretation of Anfinsen's hypothesis states that one amino acid sequence should fold into a single, native, ordered state, or a highly similar set thereof, coinciding with the global minimum in the folding-energy landscape, which, in turn, is responsible for the function of the protein. However, this classical view is challenged by many proteins and peptide sequences, which can adopt exchangeable, significantly dissimilar conformations that even fulfill different biological roles. The similarities and differences of concepts related to these proteins, mainly chameleon sequences, metamorphic proteins, and switch peptides, which are all denoted herein "turncoat" polypeptides, are reviewed. As well as adding a twist to the conventional view of protein folding, the lack of structural definition adds clear versatility to the activity of proteins and can be used as a tool for protein design and further application in biotechnology and biomedicine.

An overview of 25 years of research on Thermococcus kodakarensis, a genetically versatile model organism for archaeal research

Almost 25 years have passed since the discovery of a planktonic, heterotrophic, hyperthermophilic archaeon named Thermococcus kodakarensis KOD1, previously known as Pyrococcus sp. KOD1, by Imanaka and coworkers. T. kodakarensis is one of the most studied archaeon in terms of metabolic pathways, available genomic resources, established genetic engineering techniques, reporter constructs, in vitro transcription/translation machinery, and gene expression/gene knockout systems. In addition to all these, ease of growth using various carbon sources makes it a facile archaeal model organism. Here, in this review, an attempt is made to reflect what we have learnt from this hyperthermophilic archaeon.

The CDR1 and Other Regions of Immunoglobulin Light Chains are Hot Spots for Amyloid Aggregation

Immunoglobulin light chain-derived (AL) amyloidosis is a debilitating disease without known cure. Almost nothing is known about the structural factors driving the amyloidogenesis of the light chains. This study aimed to identify the fibrillogenic hotspots of the model protein 6aJL2 and in pursuing this goal, two complementary approaches were applied. One of them was based on several web-based computational tools optimized to predict fibrillogenic/aggregation-prone sequences based on different structural and biophysical properties of the polypeptide chain. Then, the predictions were confirmed with an ad-hoc synthetic peptide library. In the second approach, 6aJL2 protein was proteolyzed with trypsin, and the products incubated in aggregation-promoting conditions. Then, the aggregation-prone fragments were identified by combining standard proteomic methods, and the results validated with a set of synthetic peptides with the sequence of the tryptic fragments. Both strategies coincided to identify a fibrillogenic hotspot located at the CDR1 and β-strand C of the protein, which was confirmed by scanning proline mutagenesis analysis. However, only the proteolysis-based strategy revealed additional fibrillogenic hotspots in two other regions of the protein. It was shown that a fibrillogenic hotspot associated to the CDR1 is also encoded by several κ and λ germline variable domain gene segments. Some parts of this study have been included in the chapter “The Structural Determinants of the Immunoglobulin Light Chain Amyloid Aggregation”, published in Physical Biology of Proteins and Peptides, Springer 2015 (ISBN 978-3-319-21687-4).

Chameleon sequences in neurodegenerative diseases

Chameleon sequences can adopt either alpha helix sheet or a coil conformation. Defining chameleon sequences in PDB (Protein Data Bank) may yield to an insight on defining peptides and proteins responsible in neurodegeneration. In this research, we benefitted from the large PDB and performed a sequence analysis on Chameleons, where we developed an algorithm to extract peptide segments with identical sequences, but different structures. In order to find new chameleon sequences, we extracted a set of 8315 non-redundant protein sequences from the PDB with an identity less than 25%. Our data was classified to "helix to strand (HE)", "helix to coil (HC)" and "strand to coil (CE)" alterations. We also analyzed the occurrence of singlet and doublet amino acids and the solvent accessibility in the chameleon sequences; we then sorted out the proteins with the most number of chameleon sequences and named them Chameleon Flexible Proteins (CFPs) in our dataset. Our data revealed that Gly, Val, Ile, Tyr and Phe, are the major amino acids in Chameleons. We also found that there are proteins such as Insulin Degrading Enzyme IDE and GTP-binding nuclear protein Ran (RAN) with the most number of chameleons (640 and 405 respectively). These proteins have known roles in neurodegenerative diseases. Therefore it can be inferred that other CFP's can serve as key proteins in neurodegeneration, and a study on them can shed light on curing and preventing neurodegenerative diseases.

Virtual screening on an α-helix to β-strand switchable region of the FGFR2 extracellular domain revealed positive and negative modulators

The secondary structure of some protein segments may vary between α-helix and β-strand. To predict these switchable segments, we have developed an algorithm, Switch-P, based solely on the protein sequence. This algorithm was used on the extracellular parts of FGF receptors. For FGFR2, it predicted that β4 and β5 strands of the third Ig-like domain were highly switchable. These two strands possess a high number of somatic mutations associated with cancer. Analysis of PDB structures of FGF receptors confirmed the switchability prediction for β5. We thus evaluated if compound-driven α-helix/β-strand switching of β5 could modulate FGFR2 signaling. We performed the virtual screening of a library containing 1.4 million of chemical compounds with two models of the third Ig-like domain of FGFR2 showing different secondary structures for β5, and we selected 32 compounds. Experimental testing using proliferation assays with FGF7-stimulated SNU-16 cells and a FGFR2-dependent Erk1/2 phosphorylation assay with FGFR2-transfected L6 cells, revealed activators and inhibitors of FGFR2. Our method for the identification of switchable proteinic regions, associated with our virtual screening approach, provides an opportunity to discover new generation of drugs with under-explored mechanism of action.

A survey of machine learning methods for secondary and supersecondary protein structure prediction

In this chapter we provide a survey of protein secondary and supersecondary structure prediction using methods from machine learning. Our focus is on machine learning methods applicable to β-hairpin and β-sheet prediction, but we also discuss methods for more general supersecondary structure prediction. We provide background on the secondary and supersecondary structures that we discuss, the features used to describe them, and the basic theory behind the machine learning methods used. We survey the machine learning methods available for secondary and supersecondary structure prediction and compare them where possible.

Certain heptapeptide and large sequences representing an entire helix, strand or coil conformation in proteins are associated as chameleon sequences

Helices, strands and coils in proteins of known three-dimensional structure, corresponding to heptapeptide and large sequences (‘probe’ peptides), were scanned against peptide sequences of variable length, comprising seven or more residues that correspond to a different conformation (‘target’ peptides) in protein crystal structures available from the Protein Data Bank (PDB). Where the ‘probe’ and ‘target’ peptide sequences exactly match, they correspond to ‘chameleon’ sequences in protein structures. We observed ∼548 heptapeptide and large chameleon sequences that included peptides in the coil conformation from 53,794 PDB files that were analyzed. However, after excluding several chameleon peptides based on the quality of protein structure data, redundancy and peptides associated with cloning artifacts, such as, histidine-tags, we observed only ten chameleon peptides in structurally different proteins and the maximum length comprised seven amino acid residues. Our analysis suggests that the quality of protein structure data is important for identifying possibly, the ‘true chameleons’ in PDB. Majority of the chameleon sequences correspond to an entire strand in one protein that is observed as part of helix sequence in another protein. The heptapeptide chameleons are characterized with a high propensity of alanine, leucine and valine amino acid residues. The total hydropathy values range between −11.2 and 22.9, the difference in solvent accessibility between 2.0 Å² and 373 Å² units and the difference in total number of residue neighbor contacts between 0 and 7 residues. Our work identifies for the first time heptapeptide and large sequences that correspond to a single complete helix, strand or coil, which adopt entirely different secondary structures in another protein.

Intrinsically disordered proteins from A to Z

The ideas that proteins might possess specific functions without being uniquely folded into rigid 3D-structures and that these floppy polypeptides might constitute a noticeable part of any given proteome would have been considered as a preposterous fiction 15 or even 10 years ago. The situation has changed recently, and the existence of functional yet intrinsically disordered proteins and regions has become accepted by a significant number of protein scientists. These fuzzy objects with fuzzy structures and fuzzy functions are among the most interesting and attractive targets for modern protein research. This review summarizes some of the major discoveries and breakthroughs in the field of intrinsic disorder by representing related concepts and definitions.

Stabilization by fusion to the C-terminus of hyperthermophile Sulfolobus tokodaii RNase HI: a possibility of protein stabilization tag

RNase HI from the hyperthermophile Sulfolobus tokodaii (Sto-RNase HI) is stabilized by its C-terminal residues. In this work, the stabilization effect of the Sto-RNase HI C-terminal residues was investigated in detail by thermodynamic measurements of the stability of variants lacking the disulfide bond (C58/145A), or the six C-terminal residues (ΔC6) and by structural analysis of ΔC6. The results showed that the C-terminal does not affect overall structure and stabilization is caused by local interactions of the C-terminal, suggesting that the C-terminal residues could be used as a “stabilization tag.” The Sto-RNase HI C-terminal residues (-IGCIILT) were introduced as a tag on three proteins. Each chimeric protein was more stable than its wild-type protein. These results suggested the possibility of a simple stabilization technique using a stabilization tag such as Sto-RNase HI C-terminal residues.

Short-lived alpha-helical intermediates in the folding of beta-sheet proteins

Several lines of evidence point strongly toward the importance of highly alpha-helical intermediates in the folding of all globular proteins, regardless of their native structure. However, experimental refolding studies demonstrate no observable alpha-helical intermediate during refolding of some beta-sheet proteins and have dampened enthusiasm for this model of protein folding. In this study, beta-sheet proteins were hypothesized to have potential to form amphiphilic helices at a period of <3.6 residues/turn that matches or exceeds the potential at 3.6 residues/turn. Hypothetically, such potential is the basis for an effective and unidirectional mechanism by which highly alpha-helical intermediates might be rapidly disassembled during folding and potentially accounts for the difficulty in detecting highly alpha-helical intermediates during the folding of some proteins. The presence of this potential was confirmed, indicating that a model entailing ubiquitous formation of alpha-helical intermediates during the folding of globular proteins predicts previously unrecognized features of primary structure. Further, the folding of fatty acid binding protein, a predominantly beta-sheet protein that exhibits no apparent highly alpha-helical intermediate during folding, was dramatically accelerated by 2,2,2-trifluoroethanol, a solvent that stabilizes alpha-helical structure. This observation suggests that formation of an alpha-helix can be a rate-limiting step during folding of a predominantly beta-sheet protein and further supports the role of highly alpha-helical intermediates in the folding of all globular proteins.

Analysis of protein chameleon sequence characteristics

Conversion of local structural state of a protein from an α-helix to a β-strand is usually associated with a major change in the tertiary structure. Similar changes were observed during the self assembly of amyloidogenic proteins to form fibrils, which are implicated in severe diseases conditions, e.g., Alzheimer disease. Studies have emphasized that certain protein sequence fragments known as chameleon sequences do not have a strong preference for either helical or the extended conformations. Surprisingly, the information on the local sequence neighborhood can be used to predict their secondary at a high accuracy level. Here we report a large scale-analysis of chameleon sequences to estimate their propensities to be associated with different local structural states such as α -helices, β-strands and coils. With the help of the propensity information derived from the amino acid composition, we underline their complexity, as more than one quarter of them prefers coil state over to the regular secondary structures. About half of them show preference for both α-helix and β-sheet conformations and either of these two states is favored by the rest.

Slow unfolding of monomeric proteins from hyperthermophiles with reversible unfolding

Based on the differences in their optimal growth temperatures microorganisms can be classified into psychrophiles, mesophiles, thermophiles, and hyperthermophiles. Proteins from hyperthermophiles generally exhibit greater stability than those from other organisms. In this review, we collect data about the stability and folding of monomeric proteins from hyperthermophilies with reversible unfolding, from the equilibrium and kinetic aspects. The results indicate that slow unfolding is a general strategy by which proteins from hyperthermophiles adapt to higher temperatures. Hydrophobic interaction is one of the factors in the molecular mechanism of the slow unfolding of proteins from hyperthermophiles.

Proline effect on the thermostability and slow unfolding of a hyperthermophilic protein

Ribonuclease HII from hyperthermophile Thermococcus kodakaraensis (Tk-RNase HII) is a robust monomeric protein under kinetic control, which possesses some proline residues at the N-terminal of alpha-helices. Proline residue at the N-terminal of an alpha-helix is thought to stabilize a protein. In this work, the thermostability and folding kinetics of Tk-RNase HII were measured for mutant proteins in which a proline residue is introduced (Xaa to Pro) or removed (Pro to Ala) at the N-terminal of alpha-helices. In the folding experiments, the mutant proteins examined exhibit little influence on the remarkably slow unfolding of Tk-RNase HII. In contrast, E111P and K199P exhibit some thermostabilization, whereas P46A, P70A and P174A have some thermodestabilization. E111P/K199P and P46A/P70A double mutations cause cumulative changes in stability. We conclude that the proline effect on protein thermostability is observed in a hyperthermophilic protein, but each proline residue at the N-terminal of an alpha-helix slightly contributes to the thermostability. The present results also mean that even a natural hyperthermophilic protein can acquire improved thermostability.

Flexible nets: disorder and induced fit in the associations of p53 and 14-3-3 with their partners

Background

Proteins are involved in many interactions with other proteins leading to networks that regulate and control a wide variety of physiological processes. Some of these proteins, called hub proteins or hubs, bind to many different protein partners. Protein intrinsic disorder, via diversity arising from structural plasticity or flexibility, provide a means for hubs to associate with many partners (Dunker AK, Cortese MS, Romero P, Iakoucheva LM, Uversky VN: Flexible Nets: The roles of intrinsic disorder in protein interaction networks. FEBS J 2005, 272:5129-5148).

Results

Here we present a detailed examination of two divergent examples: 1) p53, which uses different disordered regions to bind to different partners and which also has several individual disordered regions that each bind to multiple partners, and 2) 14-3-3, which is a structured protein that associates with many different intrinsically disordered partners. For both examples, three-dimensional structures of multiple complexes reveal that the flexibility and plasticity of intrinsically disordered protein regions as well as induced-fit changes in the structured regions are both important for binding diversity.

Conclusions

These data support the conjecture that hub proteins often utilize intrinsic disorder to bind to multiple partners and provide detailed information about induced fit in structured regions.