Genome-wide computational determination of the human metalloproteome.

Accurate prediction of protein function in humans is important for understanding biological processes at the molecular level in biomedicine and drug design. Over a third of proteins are commonly held to bind metal, and ∼10% of human proteins, to bind zinc. Therefore, an initial step in protein function prediction frequently involves predicting metal ion binding. In recent years, methods have been developed to predict a set of residues in 3D space forming the metal-ion binding site, often with a high degree of accuracy. Here, using extensions of these methods, we provide an extensive list of human proteins and their putative metal ion binding site residues, using translated gene sequences derived from the complete, resolved human genome. Under conditions of ∼90% selectivity, over 900 new human putative metal ion binding proteins are identified. A statistical analysis of resolved metal ion binding sites in the human metalloproteome is furnished and the importance of remote homology analysis is demonstrated. As an example, a novel metal-ion binding site involving a complex of a botulinum substrate with its inhibitor is presented. On the basis of the location of the predicted site and the interactions of the contacting residues at the complex interface, we postulate that metal ion binding in this region could influence complex formation and, consequently, the functioning of the protein. Thus, this work provides testable hypotheses about novel functions of known proteins.

Use of host-like peptide motifs in viral proteins is a prevalent strategy in host-virus interactions.

Viruses interact extensively with host proteins, but the mechanisms controlling these interactions are not well understood. We present a comprehensive analysis of eukaryotic linear motifs (ELMs) in 2,208 viral genomes and reveal that viruses exploit molecular mimicry of host-like ELMs to possibly assist in host-virus interactions. Using a statistical genomics approach, we identify a large number of potentially functional ELMs and observe that the occurrence of ELMs is often evolutionarily conserved but not uniform across virus families. Some viral proteins contain multiple types of ELMs, in striking similarity to complex regulatory modules in host proteins, suggesting that ELMs may act combinatorially to assist viral replication. Furthermore, a simple evolutionary model suggests that the inherent structural simplicity of ELMs often enables them to tolerate mutations and evolve quickly. Our findings suggest that ELMs may allow fast rewiring of host-virus interactions, which likely assists rapid viral evolution and adaptation to diverse environments.

Interactions of subunit CCT3 in the yeast chaperonin CCT/TRiC with Q/N-rich proteins revealed by high-throughput microscopy analysis.

The eukaryotic chaperonin containing t-complex polypeptide 1 (CCT/TRiC) is an ATP-fueled machine that assists protein folding. It consists of two back-to-back stacked rings formed by eight different subunits that are arranged in a fixed permutation. The different subunits of CCT are believed to possess unique substrate binding specificities that are still mostly unknown. Here, we used high-throughput microscopy analysis of yeast cells to determine changes in protein levels and localization as a result of a Glu to Asp mutation in the ATP binding site of subunits 3 (CCT3) or 6 (CCT6). The mutation in subunit CCT3 was found to induce cytoplasmic foci termed P-bodies where mRNAs, which are not translated, accumulate and can be degraded. Analysis of the changes in protein levels and structural modeling indicate that P-body formation in cells with the mutation in CCT3 is linked to the specific interaction of this subunit with Gln/Asn-rich segments that are enriched in many P-body proteins. An in vitro gel-shift analysis was used to show that the mutation in subunit CCT3 interferes with the ability of CCT to bind a Gln/Asn-rich protein aggregate. More generally, the strategy used in this work can be used to unravel the substrate specificities of other chaperone systems.

Modulation of folding kinetics of repeat proteins: interplay between intra- and interdomain interactions.

Repeat proteins have unique elongated structures that, unlike globular proteins, are quite modular. Despite their simple one-dimensional structure, repeat proteins exhibit intricate folding behavior with a complexity similar to that of globular proteins. Therefore, repeat proteins allow one to quantify fundamental aspects of the biophysics of protein folding. One important feature of repeat proteins is the interfaces between the repeating units. In particular, the distribution of stabilities within and between the repeats was previously suggested to affect their folding characteristics. In this study, we explore how the interface affects folding kinetics and cooperativity by investigating two families of repeat proteins, namely, the Ankyrin and tetratricopeptide repeat proteins, which differ in the number of interfacial contacts that are formed between their units as well as in their folding behavior. By using simple topology-based models, we show that modulating the energetic strength of the interface relative to that of the repeat itself can drastically change the protein stability, folding rate, and cooperativity. By further dissecting the interfacial contacts into several subsets, we isolated the effects of each of these groups on folding kinetics. Our study highlights the importance of interface connectivity in determining the folding behavior.

The origins and evolution of ubiquitination sites.

Protein ubiquitination is central to the regulation of various pathways in eukaryotes. The process of ubiquitination and its cellular outcome were investigated in hundreds of proteins to date. Despite this, the evolution of this regulatory mechanism has not yet been addressed comprehensively. Here, we quantify the rates of evolutionary changes of ubiquitination and SUMOylation (Small Ubiquitin-like MOdifier) sites. We estimate the time at which they first appeared, and compare them to acetylation and phosphorylation sites and to unmodified residues. We observe that the various modification sites studied exhibit similar rates. Mammalian ubiquitination sites are weakly more conserved than unmodified lysine residues, and a higher degree of relative conservation is observed when analyzing bona fide ubiquitination sites. Various reasons can be proposed for the limited level of excess conservation of ubiquitination, including shifts in locations of the sites, the presence of alternative sites, and changes in the regulatory pathways. We observe that disappearance of sites may be compensated by the presence of a lysine residue in close proximity, which is significant when compared to evolutionary patterns of unmodified lysine residues, especially in disordered regions. This emphasizes the importance of analyzing a window in the vicinity of functional residues, as well as the capability of the ubiquitination machinery to ubiquitinate residues in a certain region. Using prokaryotic orthologs of ubiquitinated proteins, we study how ubiquitination sites were formed, and observe that while sometimes sequence additions and rearrangements are involved, in many cases the ubiquitination machinery utilizes an already existing sequence without significantly changing it. Finally, we examine the evolution of ubiquitination, which is linked with other modifications, to infer how these complex regulatory modules have evolved. Our study gives initial insights into the formation of ubiquitination sites, their degree of conservation in various species, and their co-evolution with other posttranslational modifications.

What distinguishes GroEL substrates from other Escherichia coli proteins?

Experimental studies and theoretical considerations have shown that only a small subset of Escherichia coli proteins fold in vivo with the help of the GroE chaperone system. These proteins, termed GroE substrates, have been divided into three classes: (a) proteins that can fold independently, but are found to associate with GroEL; (b) proteins that require GroE when the cell is under stress; and (c) 'obligatory' proteins that require GroE assistance even under normal conditions. It remains unclear, however, why some proteins need GroE and others do not. Here, we review experimental and computational studies that addressed this question by comparing the sequences and structural, biophysical and evolutionary properties of GroE substrates with those of nonsubstrates. In general, obligatory substrates are found to have lower folding propensities and be more aggregation prone. GroE substrates are also more conserved than other proteins and tend to utilize more optimal codons, but this latter feature is less apparent for obligatory substrates. There is no evidence, however, for any specific sequence signatures although there is a tendency for sequence periodicity. Our review shows that reliable sequence- or structure-based predictions of GroE dependency remain a challenge. We suggest that the different classes of GroE substrates be studied separately and that proper control test sets (e.g. TIM barrel proteins that need GroE for folding versus TIM barrels that fold independently) be used more extensively in such studies.

Intrinsic disorder in ubiquitination substrates.

The ubiquitin-proteasome system is responsible for the degradation of numerous proteins in eukaryotes. Degradation is an essential process in many cellular pathways and involves the proteasome degrading a wide variety of unrelated substrates while retaining specificity in terms of its targets for destruction and avoiding unneeded proteolysis. How the proteasome achieves this task is the subject of intensive research. Many proteins are targeted for degradation by being covalently attached to a poly-ubiquitin chain. Several studies have indicated the importance of a disordered region for efficient degradation. Here, we analyze a data set of 482 in vivo ubiquitinated substrates and a subset in which ubiquitination is known to mediate degradation. We show that, in contrast to phosphorylation sites and other regulatory regions, ubiquitination sites do not tend to be located in disordered regions and that a large number of substrates are modified at structured regions. In degradation-mediated ubiquitination, there is a significant bias of ubiquitination sites to be in disordered regions; however, a significant number is still found in ordered regions. Moreover, in many cases, disordered regions are absent from ubiquitinated substrates or are located far away from the modified region. These surprising findings raise the question of how these proteins are successfully unfolded and ultimately degraded by the proteasome. They indicate that the folded domain must be perturbed by some additional factor, such as the p97 complex, or that ubiquitination may induce unfolding.

Searching DNA via a "Monkey Bar" mechanism: the significance of disordered tails.

The search through nonspecific DNA for a specific site by proteins is known to be facilitated by sliding, hopping, and intersegment transfer between separate DNA strands, yet the driving forces of these protein dynamics from the molecular perspective are unclear. In this study, molecular features of the DNA search mechanism were explored for three homologous proteins (the HoxD9, Antp, and NK-2 homeodomains) using a simple computational model in which protein-DNA interactions are represented solely by electrostatic forces. In particular, we studied the impact that disordered N-terminal tails (N-tails), which are more common in DNA-binding proteins than in other proteins, have on the efficiency of DNA search. While the three homeodomain proteins were found to use similar binding interfaces in specific and nonspecific interactions with DNAs, their different electrostatic potentials affect the nature of their sliding dynamics. The different lengths and net charges of the N-tails of the homeodomains affect their motion along the DNA. The presence of an N-tail increases sliding propensity but slows linear diffusion along the DNA. When the search is performed in the presence of two parallel DNA molecules, a direct transfer, which is facilitated by the protein tail, from one nonspecific DNA to another occurs. The tailed proteins jump between two DNA molecules through an intermediate in which the recognition helix of the protein is adsorbed to one DNA fragment and the N-tail is adsorbed to the second, suggesting a "monkey bar" mechanism. Our study illustrates how the molecular architecture of proteins controls the efficiency of DNA scanning.

Nonnative electrostatic interactions can modulate protein folding: molecular dynamics with a grain of salt.

In recent years, a growing number of protein folding studies have focused on the unfolded state, which is now recognized as playing a major role in the folding process. Some of these studies show that interactions occurring in the unfolded state can significantly affect the stability and kinetics of the protein folding reaction. In this study, we modeled the effect of electrostatic interactions, both native and nonnative, on the folding of three protein systems that underwent selective charge neutralization or reversal or complete charge suppression. In the case of the N-terminal L9 protein domain, our results directly attribute the increase in thermodynamic stability to destabilization of the unfolded ensemble, reaffirming the experimental observations. These results provide a deeper structural insight into the ensemble of the unfolded state and predict a new mutation site for increased protein stability. In the second case, charge reversal mutations of RNase Sa affected protein stability, with the destabilizing mutations being less destabilizing at higher salt concentrations, indicating the formation of charge-charge interactions in the unfolded state. In the N-terminal L9 and RNase Sa systems, changes in electrostatic interactions in the unfolded state that cause an increase in free energy had an overall compaction effect that suggests a decrease in entropy. In the third case, in which we compared the beta-lactalbumin and hen egg-white lysozyme protein homologues, we successfully eliminated differences between the folding kinetics of the two systems by suppressing electrostatic interactions, supporting previously reported findings. Our coarse-grained molecular dynamics study not only reproduces experimentally reported findings but also provides a detailed molecular understanding of the elusive unfolded-state ensemble and how charge-charge interactions can modulate the biophysical characteristics of folding.