Exploring the sequence determinants of amyloid structure using position-specific scoring matrices.

Protein aggregation results in beta-sheet-like assemblies that adopt either a variety of amorphous morphologies or ordered amyloid-like structures. These differences in structure also reflect biological differences; amyloid and amorphous beta-sheet aggregates have different chaperone affinities, accumulate in different cellular locations and are degraded by different mechanisms. Further, amyloid function depends entirely on a high intrinsic degree of order. Here we experimentally explored the sequence space of amyloid hexapeptides and used the derived data to build Waltz, a web-based tool that uses a position-specific scoring matrix to determine amyloid-forming sequences. Waltz allows users to identify and better distinguish between amyloid sequences and amorphous beta-sheet aggregates and allowed us to identify amyloid-forming regions in functional amyloids.

Quantifying information transfer by protein domains: analysis of the Fyn SH2 domain structure.

BACKGROUND: Efficient communication between distant sites within a protein is essential for cooperative biological response. Although often associated with large allosteric movements, more subtle changes in protein dynamics can also induce long-range correlations. However, an appropriate formalism that directly relates protein structural dynamics to information exchange between functional sites is still lacking. RESULTS: Here we introduce a method to analyze protein dynamics within the framework of information theory and show that signal transduction within proteins can be considered as a particular instance of communication over a noisy channel. In particular, we analyze the conformational correlations between protein residues and apply the concept of mutual information to quantify information exchange. Mapping out changes of mutual information on the protein structure then allows visualizing how distal communication is achieved. We illustrate the approach by analyzing information transfer by the SH2 domain of Fyn tyrosine kinase, obtained from Monte Carlo dynamics simulations. Our analysis reveals that the Fyn SH2 domain forms a noisy communication channel that couples residues located in the phosphopeptide and specificity binding sites and a number of residues at the other side of the domain near the linkers that connect the SH2 domain to the SH3 and kinase domains. We find that for this particular domain, communication is affected by a series of contiguous residues that connect distal sites by crossing the core of the SH2 domain. CONCLUSION: As a result, our method provides a means to directly map the exchange of biological information on the structure of protein domains, making it clear how binding triggers conformational changes in the protein structure. As such it provides a structural road, next to the existing attempts at sequence level, to predict long-range interactions within protein structures.

How evolutionary pressure against protein aggregation shaped chaperone specificity.

As protein aggregation is potentially lethal, control of protein conformation by molecular chaperones is essential for cellular organisms. This is especially important during protein expression and translocation, since proteins are then unfolded and therefore most susceptible to form non-native interactions. Using TANGO, a statistical mechanics algorithm to predict protein aggregation, we here analyse the aggregation propensities of 28 complete proteomes. Our results show that between 10% and 20% of the residues in these proteomes are within aggregating protein segments and that this represents a lower limit for the aggregation tendency of globular proteins. Further, we show that not only evolution strongly pressurizes aggregation downwards by minimizing the amount of strongly aggregating sequences but also by selectively capping strongly aggregating hydrophobic protein sequences with arginine, lysine and proline. These residues are strongly favoured at these positions as they function as gatekeepers that are most efficient at opposing aggregation. Finally, we demonstrate that the substrate specificity of different unrelated chaperone families is geared by these gatekeepers. Chaperones face the difficulty of having to combine substrate affinity for a broad range of hydrophobic sequences and selectivity for those hydrophobic sequences that aggregate most strongly. We show that chaperones achieve these requirements by preferentially binding hydrophobic sequences that are capped by positively charged gatekeeper residues. In other words, targeting evolutionarily selected gatekeepers allows chaperones to prioritize substrate recognition according to aggregation propensity.

Domain swapping in p13suc1 results in formation of native-like, cytotoxic aggregates.

The field of protein aggregation has been occupied mainly with the study of beta-strand self-association that occurs as a result of misfolding and leads to the formation of toxic protein aggregates and amyloid fibers. However, some of these aggregates retain native-like structural and enzymatic properties suggesting mechanisms other than beta-strand assembly. p13suc1 is a small protein that can exist as a monomer or a domain-swapped dimer. Here, we show that, under native conditions, p13suc1 forms three-dimensional domain-swapped aggregates, and that these aggregates are cytotoxic. Thus, toxicity of protein aggregates is not only associated with beta-rich assemblies and amyloid fibers, involving non-native interactions, but it can be induced by oligomeric misassembly that maintains predominantly native-like interactions.

The unfolding story of three-dimensional domain swapping.

Three-dimensional domain swapping is the event by which a monomer exchanges part of its structure with identical monomers to form an oligomer where each subunit has a similar structure to the monomer. The accumulating number of observations of this phenomenon in crystal structures has prompted speculation as to its biological relevance. Domain swapping was originally proposed to be a mechanism for the emergence of oligomeric proteins and as a means for functional regulation, but also to be a potentially harmful process leading to misfolding and aggregation. We highlight experimental studies carried out within the last few years that have led to a much greater understanding of the mechanism of domain swapping and of the residue- and structure-specific features that facilitate the process. We discuss the potential biological implications of domain swapping in light of these findings.

The structure of the transition state for folding of domain-swapped dimeric p13suc1.

suc1 has two native states, a monomer and a domain-swapped dimer, in which one molecule exchanges a beta strand with an identical partner. Thus, monomer and dimer have the same structures but are topologically distinct. Importantly, residues that exchange are part of the folding nucleus of the monomer and therefore forming these interactions in the dimer would be expected to incur a large entropic cost. Here we present the transition state for folding/unfolding of domain-swapped dimeric suc1 and compare it with its monomeric counterpart. The same overall structure is observed in the two transition states but the phi values are consistently higher for the domain-swapped dimer. Thus, a greater entropic penalty for bringing together the key interactions in the dimer is overcome by mobilizing more contacts in the transition state, thereby achieving a greater enthalpic gain.

Prediction of water and metal binding sites and their affinities by using the Fold-X force field.

The empirical force field Fold-X was developed previously to allow rapid free energy calculations in proteins. Here, we present an enhanced version of the force field allowing prediction of the position of structural water molecules and metal ions, together called single atom ligands. Fold-X picks up 76% of water molecules found to interact with two or more polar atoms of proteins in high-resolution crystal structures and predicts their position to within 0.8 A on average. The prediction of metal ion-binding sites have success rates between 90% and 97% depending on the metal, with an overall standard deviation on the position of binding of 0.3-0.6 A. The following metals were included in the force field: Mg2+, Ca2+, Zn2+, Mn2+, and Cu2+. As a result, the current version of Fold-X can accurately decorate a protein structure with biologically important ions and water molecules. Additionally, the free energy of binding of Ca2+ and Zn2+ (i.e., the natural logarithm of the dissociation constant) and its dependence on ionic strength correlate reasonably well with the experimental data available in the literature, allowing one to discriminate between high- and low-affinity binding sites. Importantly, the accuracy of the energy prediction presented here is sufficient to efficiently discriminate between Mg2+, Ca2+, and Zn2+ binding.

Intermediates control domain swapping during folding of p13suc1.

The 13-kDa protein p13(suc1) has two folded states, a monomer and a structurally similar domain-swapped dimer formed by exchange of a beta-strand. The refolding reaction of p13(suc1) is multiphasic, and in this paper we analyze the kinetics as a function of denaturant and protein concentration and compare the behavior of wild type and a set of mutants previously designed with dimerization propensities that span 9 orders of magnitude. We show that the folding reactions of wild type and all mutants produce the monomer predominantly despite their very different equilibrium behavior. However, the addition of low concentrations of denaturant in the refolding buffer leads to thermodynamic control of the folding reaction with products that correspond to the wild type and mutant equilibrium dimerization propensities. We present evidence that the kinetic control in the absence of urea arises because of the population of the folding intermediates. Intermediates are usually considered to be detrimental to folding because they slow down the reaction; however, our work shows that intermediates buffer the monomeric folding pathway against the effect of mutations that favor the nonfunctional, dimeric state at equilibrium.

Folding and association of the human cell cycle regulatory proteins ckshs1 and ckshs2.

The two human proteins ckshs1 and ckshs2 are each 79 amino acids in length and consist of a four-stranded beta-sheet capped at one end by two alpha-helices. They are members of the cks family of essential cell cycle regulatory proteins that can adopt two native states, a monomer and a domain-swapped dimer formed by exchange of a C-terminal beta-strand. ckshs1 and ckshs2 both have marginal thermodynamic stability (the free energies of unfolding at 25 degrees C are 3.0 and 2.5 kcal/mol, respectively) and low kinetic stability (the rates of unfolding in water are approximately 1 s(-1)). Refolding of their denatured states to the monomeric forms of the proteins is slowed by transient oligomerization that is likely to occur via domain swapping. The folding behavior of ckshs1 and ckshs2 is markedly different from that of suc1, the cks protein from Schizosaccharomyces pombe, but the domain swapping propensities are similar. The greater thermodynamic and kinetic stability of suc1 and the population of a folding intermediate are most likely a consequence of its larger size (113 residues). The similarity in the domain swapping propensities, despite the contrast in other biophysical properties, may be attributable to the common double-proline motif in the hinge loop that connects the swapped domain to the rest of the protein. The motif was shown previously for suc1 to control the equilibrium between the monomer and the domain-swapped dimer. Finally, according to our model, the kinetic barrier separating the monomer and the domain-swapped dimer arises because the protein must unfold for beta-strand exchange to occur. Consistent with this, interconversion between the two states is much faster in the human proteins than it is for suc1, reflecting the faster unfolding rates of the former.