A systematic comparative and structural analysis of protein phosphorylation sites based on the mtcPTM database
© Jiménez et al. licensee BioMed Central Ltd. 2007
Received: 3 January 2007
Accepted: 23 May 2007
Published: 23 May 2007
mtcPTM is an online repository of human and mouse phosphosites in which data are hierarchically organized to preserve biologically relevant experimental information, thus allowing straightforward comparisons of phosphorylation patterns found under different conditions. The database also contains the largest available collection of atomic models of phosphorylatable proteins. Detailed analysis of this structural dataset reveals that phosphorylation sites are found in a heterogeneous range of structural and sequence contexts. mtcPTM is available on the web http://www.mitocheck.org/cgi-bin/mtcPTM/search.
In recent years, several sequencing projects have revealed the complete transcriptomes and proteomes for a number of organisms, including human [1, 2]. The current challenge is to place this information within the dynamic context of the cell in order to elucidate how individual molecules interact to achieve the complex behavior of cellular processes, which translates into the ability of living organisms to adapt and thrive in a myriad of environments and conditions. Thus, much effort has been invested in identifying, for example, the transcription patterns of genes and the interacting partners of proteins in order to determine the connections that establish the intricate cellular pathways [3, 4]. To understand these networks fully, however, we must also comprehend how their connections are regulated when the states of individual components are altered, for example by means of post-translational modifications (PTMs). It is therefore crucial to identify which proteins can be modified as well as the effect and lifetime of the PTMs.
Among PTMs, reversible protein phosphorylation is known to play a key role in regulating a variety of processes in eukaryotes, from the cell division cycle to neuronal plasticity [5, 6]. The most commonly observed phosphorylations affect serine, threonine, and tyrosine residues [7, 8], although phosphorylation of histidines and aspartates has also been reported (for review ). Protein phosphorylation is catalyzed by enzymes called protein kinases, which are usually specific for either tyrosine or serine/threonine, with few of them being able to modify all three residues indistinguishably [10–12]. The human genome encodes 518 protein kinases [13, 14], and recent estimates suggest that around one-third of cellular proteins could undergo phosphorylation . Despite the progress made during the past few decades, our knowledge about regulation of protein function by phosphorylation and the basis of kinase specificity remains incomplete, mainly because of lack of data. High-throughput proteomic approaches are expected to help fill this gap because they can identify large amounts of in vivo modified peptides (for review [16, 17]).
Protein kinases catalyze the formation of a covalent bond between a phosphate group and a hydroxyl moiety of an amino-acid side chain. Because of the size and charge of the phosphate groups, their introduction could have a local, and potentially global, effect on the modified proteins. This effect may translate into modulation of protein activity, subcellular localization, half-life, and ability to interact with other molecules [8, 11]. Undoubtedly, the best characterized examples of the molecular effects of phosphorylation on proteins are from high-resolution structural studies (for review [18–20]). For example, some modifications that affect residues that are part of or in the vicinity of catalytic sites and protein docking interfaces may promote or disrupt substrate binding by a combination of steric and electrostatic effects, without apparent major local structural rearrangements Histidine-containing phosphocarrier protein (HPr) , isocitrate dehydrogenase , signal transducer and activator of transcription [STAT]3B , STAT-1 , and Stage II sporulation protein (SpoII)AA/SpoIIAB ). On the other hand, the modifications could cause conformational changes that result either disorder-to-order transitions (glycogen phosphorylase [26, 27]) or increased local flexibility if the native amino-acid packing is disrupted (protein kinase A [28, 29], mitogen-activated protein kinase , ubiquitin-protein ligase E3 , and potassium channel inactivation domain ). However, because of technical challenges, few atomic structures of proteins are available in their phosphorylated state.
Although atomic models of the proteins in their nonphosphorylated form can provide invaluable clues that may enhance understanding of the molecular impact of modifications on proteins or allow us to predict them , no public resource is available that routinely stores and provides this information. Furthermore, current phosphosite databases only address the storage and display of phosphosites [33–35], disregarding the experimental context of the phosphorylation. We have developed the mtcPTM (MitoCheck's post-translational modifications) database to address these needs. The mtcPTM database is a repository of PTMs in human and mouse proteins that aims to preserve and present the experimental evidence that led to the identification of each modification. We show that the graphical display of these data allows intuitive comparisons between phosphorylation patterns from different sources or experiments. The database also contains structural information on those modified protein domains for which the actual structure, or the structure of a close homolog, is available. In addition, we have analyzed in detail this large structural collection to investigate the molecular characteristics of phosphorylatable sites in terms of solvent accessibility, secondary structure preference, and degree of conservation. We report that, in general, modified residues are in flexible/exposed regions and, although they are no more conserved than expected, they present highly variable degrees of conservation. Finally, we elaborate on those cases of phosphorylatable residues that were found buried in the structures, predicting the structural/functional effect of their modification on these proteins. As part of the MitoCheck programme, a European Union-funded project whose overall aim is to study the regulation of mitosis by phosphorylation , mtcPTM was originally developed for the study of differential phosphorylation in mitosis. However, its general design is readily applicable to any data, regardless of experimental source. The database is publicly available online, and experimentalists are encouraged to submit their data for storage and display.
Handling and storage of phosphosite data
The mtcPTM database contains data retrieved from literature, protein annotations, and other databases. In the future, the database will also display phosphorylation sites that have been mapped as part of the MitoCheck project. The mtcPTM database therefore handles quite different datasets, for which the available information varies. For example, modifications retrieved from literature and protein annotation are usually recorded as individual residues, in which experimental information can only be recovered by reading the original report. By contrast, high-throughput mass spectrometry (MS) data take the form of phosphorylated positions within peptide sequences. In this case, mtcPTM preserves the experimental context of the phosphosites by grouping the MS peptides into sets according to individual experiments and assigning to each group a hierarchical data structure that summarizes the experimental information. This simple hierarchy comprises data source (for instance, a research group or programme), experimental category (for example, label describing a set of experiments that are undertaken with a combined aim), and individual experiments (data obtained from the same sample). Thus, two experiments undertaken, for example, by MitoCheck to determine the differential phosphorylation state of a protein along the cell cycle would receive the following common labels: 'MitoCheck', 'timing', and a specific label, for example interphase or mitosis.
At present, the mtcPTM database stores 13,051 and 7,930 peptides from human and mouse, respectively, corresponding to 13,116 (serine: 9839; threonine: 2067; tyrosine: 1210) and 8,889 (serine: 6942; threonine: 1470; tyrosine: 477) phosphorylations. The human-related data comprise 3842 genes and 7753 proteins, whereas for mouse they represent 2721 genes and 3866 proteins.
Display of protein phosphorylation data
The database also contains structural models for proteins and protein domains that contain modified residues. These models have been automatically built by homology modeling to empirically determined atomic co-ordinates. A conservative criterion for assignment of sequences to structures was used in order to minimize errors in the modeled domains (see Materials and methods, below). The coordinates of the models are provided as RasMol scripts , including the pair-wise alignments between the modeled Ensembl sequences and its structural templates. The mtcPTM database currently contains 2,599 structural models, 658 for mouse proteins (comprising 529 genes), and 1,191 for human (686 genes). On comparing the phosphosite dataset with these models, only a small proportion (10% in both human and mouse) of phosphosites were found in structurally defined regions. This finding is not expected to be caused by bias resulting from the type of structural data currently available, because similar proportions were observed when counting modified positions within the far more diverse Pfam domains (85% for both human and mouse proteins fell outside defined Pfam domains) . This suggested that phosphorylated sites tend to be found in flexible, unstructured segments and linkers between domains, which is in agreement with previous observations .
Interestingly, the distribution of residues between linkers and (structured) domains was not even. Phosphorylated threonine and serine residues were mainly located outside domains (structures). In mouse, 86% (91%) of serines and 83% (87%) of threonines were found in linkers between domains (structures), and similar numbers were obtained in human, specifically 87% (92%) serines and 83% (90%) threonines. However, this distribution was less biased for tyrosines, in which 37% (34%) in human and 31% (31%) in mouse were found within domains (structures). At present, it is unknown whether these differences between tyrosine and serine/threonine residues correlate with their propensity to appear in structured and flexible regions, respectively, or whether it actually reflects a biologically distinct feature of their regulation, such as specific properties in kinase recognition. Of note, the existence of different structural rules for substrate binding between serine/threonine and tyrosine protein kinases has previously been suggested .
As mentioned previously, atomic information from modified and unmodified forms of the proteins is invaluable in rationalizing the molecular effect and functional impact of phosphorylations. Therefore, even though a considerable proportion of phosphosites is situated away from structured regions, we wished to take advantage of the large structural dataset collected here to undertake a detailed study of the properties of these residues, as well as the potential effect of their modification on the domains.
Compiling a nonredundant set of structural models
For this analysis, we first defined a nonredundant (NR) set from all of the structural models stored in the database in order to preclude potential biases arising from the comparisons of highly similar structures (see Materials and methods, below). The NR set comprised 324 structural models, representing a wide range of Pfam domains, and contained 264 modified serine/threonine and 219 tyrosine residues. Half of the models were less than 150 amino acids long, indicating that half of the models represented single domains and the other half multidomain structures. Regardless of their length, the majority of the models (72%) contained only one phosphorylated residue. Furthermore, 70% of the models shared at least 80% sequence identity with their templates and only 15% less than 40%; therefore, the overall quality of the models is expected to be high.
For the structural analyses, the phosphorylated sites were clustered into two groups: one composed of serine and threonine residues, the other of tyrosines. This grouping is based on the similar characteristics of serine and threonine, and the fact that they are usually targeted by the same protein kinases. The study focused on the following structural features of the phosphosites: relative position within structured domains, solvent accessibility, secondary structure preference, and degree of conservation.
Phosphosites can accumulate at the flanks of structured domains
The second group contained residues that were in short linkers joining adjacent domains. Examples of these are the human Zinc finger protein 174 (ENSP00000268655) and the mouse discs large homolog 4 (ENSMUSP00000018700). In the first example, the phosphorylation  can take place between two zinc-finger motifs (Figure 4c). Modifications targeting the short linkers joining zinc-finger domains were also found in other proteins (data not shown), and they may regulate oligonucleotide binding because the phosphosites are part of the putative DNA binding interface. In the second example, a number of phosphosites  accumulated between the PDZ and SH3 domains of the mouse discs large homolog 4 (Figure 4d), and it is tempting to speculate that the phosphorylated state of the residues may affect the relative positioning or allosteric communication between the domains.
The last group corresponded to those sites located in long and unstructured termini relatively far away from the domains. These models were mainly built from NMR structures. For these cases, it is difficult to predict the effect that the phosphorylations could have. However, by analogy to the effect observed in other examples and considering that disordered regions appear to play important roles in protein-protein recognition events , the phosphorylation state of these sites may regulate the interaction of additional effectors to these regions, which may be especially important for those in closer proximity to the structured domains.
Phosphorylatable residues are not always accessible to solvent
Phosphorylated serine/threonines show a marginal preference for loops, whereas tyrosines do not
Phosphosites are not more conserved than expected
Because PTMs can play functional roles, phosphosites would be expected to be under purifying selection, and thus conserved through evolution. To investigate this, multiple sequence alignments were calculated from homologs to the modeled structures , and the conservation of each position corresponding to the phosphosites was assessed. The initial alignments, which can be retrieved via the mtcPTM web interface, contained nonredundant sequences sharing at least 30% sequence identity with the model. Although the inclusion of sequences that were up to 30% identical to the query domain ensured that they would adopt nearly identical structural arrangements to it , the alignments could present not only orthologous but also paralogous domains . For the latter, the phosphorylation patterns may be different or absent because of functional divergence. Furthermore, the alignments may also contain sequences from distantly related organisms in which the phosphorylation patterns may have evolved differently. To account for these potential sources of variability, the degree of conservation of each phosphosite was assessed for several subdivisions of the initial alignments. Briefly, conservation scores were calculated for the full alignments (all sequences at least 30% identical to the query) and for three subsets containing only sequences that were at least 40%, 50%, or 60% identical to the query. In alignments obtained from sequence identity cut-offs equal to or higher than 40%, most sequences are expected to be orthologous .
What happens when phosphorylatable sites are buried
As mentioned above, most phosphorylatable sites were considerably exposed to solvent and thus potentially accessible by protein kinases. However, for a few phosphosites, their side chains were found to present not only low solvent accessibility but to be actually packed into the domain core. Modification of these buried residues is likely to have structural implications because the intramolecular packing between the two states may be different. Depending on the amount of atomic interactions involved, the conformational changes could have local or global effects, from rigid body displacements to partial or total unfolding. In fact, our dataset contained some examples of proteins that have already been shown to undergo conformational changes upon phosphorylation (mitogen-activated protein kinase  and ubiquitin-protein ligase E3 ). In both cases the structural rearrangements are critical for activation of the proteins.
Given the intriguing nature of the buried phosphorylatable residues, we studied them systematically to elucidate how the phosphorylation could take place and what its potential structural impact could be. During the analysis, in order to ensure that the conformation of the residues of interest was likely to be native, only models in which the phosphorylatable side chains had been built based on the same or similar residues from the templates were considered. We also checked the consistency of poor solvent accessibility for those residues in which there existed other available models, with similar sequence identity to the templates, in the redundant set. We found 13 examples of this in which ten exhibited similar low accessibility (at a 10% cut-off) and three examples in which both the exposed and buried versions could exist, depending on the conformational states of the proteins. The latter included the active and auto-inhibitory conformations of human tyrosine-protein kinase c-Src [52, 53]. The other two examples are discussed below.
The analysis of phosphorylatable buried residues revealed that their modifications could have three major structural/functional effects on the structures: regulation of function by affecting functional sites directly or indirectly; spatial rearrangements, presumably by rigid body movements, of domains within a protein; and opening of the structure, leading to local flexibility.
Phosphorylation of buried residues found at or close to functional sites
It is tempting to hypothesize that phosphorylation of the carboxyl-terminal C2 domain of human tricalbin (ENSP00000267113) could also regulate its substrate-binding capability. This domain is phosphorylated at position Y1009 , which is located at the concave face of the β-sandwich, a region rich in positively charged residues (Figure 8c). In homologous domains, including the template used for the modeling , this poly-basic region can have phospholipid-binding capabilities . Therefore, covalent attachment of a phosphate to this tyrosine would alter the net charge of the region and thus could affect its putative ability to interact with lipids. Of note, the equivalent tyrosine (Y822) in another C2 domain within this protein can also be phosphorylated . Interestingly, these tyrosines are very conserved across C2 domains, including distantly related paralogs. However, phosphorylations at these positions have not been reported for other well characterized C2 domain subfamilies. Therefore, phosphorylation of this tyrosine may be an exclusive feature of some tricalbin C2 domains.
In addition to phosphosites in binding pockets for small substrates, modifications were also found in areas involved in binding to larger molecules. This was the case for histones (ENSP00000350159) and the zinc-finger domain of the human estrogen receptor (ENSP00000343925). In the latter, the residue that can be phosphorylated, namely S165, is found at a short loop, which is in direct contact with DNA (Figure 8e). Phosphorylation at S165 is critical for activation of the transcription factor, perhaps by precluding interactions with inhibitors that may occlude the DNA-binding interface . In the histone, the phosphorylatable Y52 residue  is packed into one of the two small cores of the histone domain (Figure 8d). This region is involved in both DNA binding and interactions with another histone monomer. In this case, the introduction of phosphate groups may also regulate the DNA-binding or protein-binding capabilities of the domains.
Phosphorylation of residues buried between domains within the same protein
The phosphorylation state of buried residues could influence the structural conformation of the protein
When buried residues that participate in the packing of the protein core or of isolated secondary structure elements to the globular domain can be phosphorylated, their modifications could lead to structural instability. In the case of isolated elements of secondary structure, for example those at the domain termini, this instability may lead to their detachment, perhaps without compromising the structural integrity of the core domain. On the other hand, when the modified residues are part of the hydrophobic core of the domain, the modification could cause considerable structural rearrangements, including local or even total unfolding, if the packing of the unmodified residue were critical for the maintenance of the overall structure and no alternative stable packing for its phosphorylated form could be established.
Additional examples of phosphorylatable buried residues at terminal regions were found in the regulator of G-protein signaling 16 and the serine/threonine protein phosphatase PP1-β catalytic subunit. In the former (ENSP00000356529), the phosphorylation of Y177 , which is found in the carboxyl-terminal helix packing tightly against the amino-terminal helix (Figure 10d), may disrupt the interaction between these two helices. This effect could be reinforced by additional phosphorylation of Y168, which also packs against the carboxyl-terminus of the amino-terminal helix (Figure 10d). Both tyrosines are well conserved and only occasionally replaced by phenylalanines, which is indicative of their important structural roles. In the human serine/threonine protein phosphatase PP1-β catalytic subunit (ENSP00000351298), S41  is found in the second amino-terminal helix (Figure 10e), whose packing to the domain appears critical in maintaining an optimal arrangement of the two terminal helices. In the crystal structure, the loop joining these two helices participates in the binding to a PEG molecule and, to a lesser extent, to the 130 kDa myosin-binding subunit of smooth muscle myosin phosphatase . Of note, S41 is only conserved in closely related sequences and is otherwise replaced by small hydrophobic residues that could still allow the helix to pack against the domain.
There were also examples of buried residues in nonterminal regions for which their phosphorylation may induce local structural rearrangements. The first example was the structure of human DJ-1, a protein that is related to male fertility and Parkinson's disease (ENSP00000340278) . Y67 is found in a loop that contributes to the packing of a short helix against the domain (Figure 10f). This helix is the most disordered region in the electron density of the crystal structure . Its flexibility could increase upon phosphorylation of Y67, because the phosphorylated amino acid may preclude its packing to the domain. This could ultimately affect the preceding loop, which participates in the dimer interface in the crystal structure. Interestingly, Y67 is poorly conserved, although the important role of this position is evident because the preferred replacement is phenylalanine. The second example was the human 60S ribosomal protein L7-A (ENSP00000339795). Here, Y139 packs against a helical motif on the surface of the protein (Figure 10g). Its phosphorylation may affect the conformation of this motif, which is involved in RNA-protein interactions . Finally, the structure of the human elongation factor EEF1A was an excellent example of both potential local conformational rearrangement and domain re-orientation (ENSP00000339053). The first involves Y29 (Figure 10h), which is found packing in the core of a helical subdomain formed at the amino-terminal elongation factor Tu GTP-binding domain. The second involves T432 (Figure 10h), which is packed between two domains, namely elongation factor Tu GTP-binding domain and elongation factor Tu carboxyl-terminal domain, of this multi-domain protein. Phosphorylation of T432 may affect the relative orientations between domains within the protein, whereas Y29 may modulate the conformation of the helical subdomain. The orientations of the domains and of the helical bundle are likely to be critical for an effective interaction of the protein with the catalytical carboxyl-terminal domain of EEF1BA (Figure 10h).
We have presented a database, mtcPTM, that stores human and mouse phosphosites . The database integrates data from low-throughput and high-throughput screenings and, in contrast to other, similar databases, it explicitly preserves the experimental context for each phosphosite by means of a three-level hierarchy annotation. Furthermore, the phosphosites are stored as relative positions within experimentally determined peptides, which allows automatic updates without loss of information, against new genome assemblies or gene builds. The data are publicly accessible via a web interface, in which the user can retrieve the mapping of the peptides onto the Ensembl genome as well as extensive information about the phosphorylated proteins including graphical comparisons of phosphorylation patterns under different conditions. At present, only mass spectrometric data are explicitly referenced to their experimental sources. In the future, the low-throughput data could be integrated likewise by, for example, using the relevant literature references to manually assign the data hierarchy.
Another major asset of the mtcPTM database is that it contains atomic models for a considerable number of phosphosites. These models have been automatically built by homology to experimentally determined structures using a conservative procedure in order to minimize modeling errors. Although a higher number of structural models could have been obtained by the use of more sensitive fold recognition techniques, the use of templates sharing low sequence similarity with the targets may decrease the overall quality of the models. To our knowledge, the structural set provided here is the largest freely available collection of phosphorylatable proteins.
From the study of this large structural collection, some general trends have been observed. For example, we find that some termini or regions close to the termini of structured domains may be preferentially phosphorylated. These phosphosites would be readily accessible by protein kinases, and their behavior could be similar to those found in the most frequent case of unstructured regions linking domains. Also, the degree of conservation of phosphosites varies considerably from protein to protein not only with respect to distantly related species and paralogs but also between closely related organisms. Few cases of highly conserved sites were found across alignments containing very diverse sequences, and some sites were also poorly conserved even in alignments of highly similar sequences. Overall, this indicates that the evolution of phosphosites could be less constrained than that of other functional motifs, such as catalytic sites. This variation could be organism-specific for a number of cellular processes or it may reflect the different importance of phosphorylation as a regulatory strategy between organisms. Alternatively, the precise position of phosphosites may in some cases not be critical for its action, especially in proteins regulated by multisite phosphorylation . Furthermore, we have also noticed that high conservation does not necessarily indicate that a given site will be phosphorylated in all homologs (for example, the C2 domains of tricalbin). Further studies will be needed to clarify and elaborate on these observations.
We have also reported that the side chains of phosphorylatable amino acids are usually exposed to solvent. However, we found a significant number of sites that were buried. It could be possible that some of these cases were artifactual, caused by experimental in vitro phosphorylation of non-full-length constructs, ambiguous peptide-protein assignment, or modeling artefacts. However, in most cases the models were built from structures that share high similarity to the target sequences and thus, except in regions with intrinsic high flexibility or poor steric constraints, the conformation of the modeled side chains are probably nearly native, especially for conserved buried residues. Also, all instances considered here in which there was available structural information for a phosphopeptide matching to different genes (mainly genes with multiple, highly similar copies such as histones, elongation factors and ribosomal proteins), the sequences of the modeled domains were nearly, if not completely, identical and therefore the equivalent phosphosites presented similar low solvent accessibility because they were built based on the same structural templates. Furthermore, some of these buried phosphosites were either identified under in vivo conditions or have already been empirically characterized in some detail, suggesting that most buried phosphorylatable residues reported here are likely to be bona fide phosphosites.
By exploring systematically the examples of phosphorylatable buried sites, we have been able to predict the functional impact of the modifications to some proteins and to classify the residues into three categories according to their predicted effect. These categories are residues in or close to binding pockets or intermolecular interfaces; residues at interfaces or in hinge regions between structural units of multidomain proteins; and residues whose phosphorylated state is likely to result in conformational rearrangements, including detachment of secondary structure elements. In the first two cases, the disturbance of binding sites or interfaces would introduce steric and electrostatic constraints for intermolecular binding leading to loss of activity and rigid body rearrangements between the structural units of the multidomain proteins, respectively. For the third category, it is tempting to suggest that isolated elements of secondary structure joined to the main domain by long disordered, usually loosely packed, linkers could be released to solvent upon modification of residues that pack them against the main domain. This release could serve to re-orient domains within the protein. Alternatively, the now flexible fragments or newly exposed areas on the core domain could play a role in the downstream recruitment of additional effectors. We have also observed cases in which the modifications are likely to result in local unfolding/refolding. However, we have not observed cases in which the modification could lead to major unfolding, perhaps because of the limited size of our dataset or because these modifications are rare. In any case, it is likely that these modifications could regulate the half-life of the protein.
How protein kinases could have access to the buried residues remains to be addressed. However, several possible mechanisms could facilitate kinase accessibility. If residues to be phosphorylated were in relatively loose or in disordered structural elements, then this intrinsic flexibility may result in spontaneous opening of the regions, temporarily exposing the residues to protein kinases. We have seen examples of phosphosites sitting in naturally flexible regions as identified by examining the temperature factor of the crystal structure templates (also see ) or the variability in the ensembles of various atomic models solved by NMR. However, it must be remembered that in some cases additional co-factors/effectors may be required to actively extract and present the residues to the phosphorylating enzymes.
Finally, although this structural study has revealed some general features, it has also made it clear that the impact of phosphorylations will depend on the atomic environment of the modified residue and thus is likely to be case specific. Therefore, we believe that making these structural data publicly available via the mtcPTM database will be of great importance to experimentalists who wish to examine in detail their particular proteins of interest. This will also provide access to the numerous examples that have not been dealt with here, including exposed functionally relevant phosphosites and those that could be silent (nonfunctional) phosphorylations.
We have implemented a database of phosphorylated residues that allows straightforward comparison of phosphorylation patterns obtained from different experiments. In addition, clues about the molecular effect of phosphorylation on some proteins domains have been obtained by examining their structural models stored on the database. It is hoped that the mtcPTM database will serve as a working tool for experimentalists, but it should also be a useful resource to theoretical biologists because it houses a large collection of protein phosphorylation data. Ultimately, the knowledge obtained from mining this database will be important to further our understanding of the regulation of protein function by PTMs. Researchers are encouraged to submit their data for storage and public display.
Materials and methods
Genomic peptide matching
The genomic matching was performed separately for each set of peptides from individual experiments. First, the peptides were compared as strings for perfect matches to Ensembl proteins. Unmatched peptides were then compared against all Uniprot entries from a given organism . When perfect matches to Uniprot were found, the matched Uniprot entry was used to extend the peptide sequences by 10 residues at both termini. The extended peptides were then BLASTed against Ensembl proteins . This sequence extension was important to improve the signal-to-noise ratio of the BLAST searches, especially for short peptides, and to increase the chance of identifying the variation that precluded a perfect match, especially when the variation took place at the ends or involved insertion/deletions that could have truncated the reported alignment. Only matches covering the initial, unextended peptides were considered further.
The adequacy of the matches was then assessed by their length coverage, sequence identity and the existence of previous perfect matches to the same proteins. The philosophy of the assessment was initially to weight more favorably high sequence similarity between genomic regions and experimentally determined peptides, and only to consider low similarity cases that involved proteins previously matched perfectly by another peptide from the same experiment. Thus, this assessment had three hierarchical levels of stringency against which all of the matches for a given peptide were examined until a match was considered to be adequate. First, if the number of mismatches was zero and the length covered was at least 85% of the peptide, or the full-length peptide was matched with a sequence identity of at least 85%, then the peptide-protein match was accepted independently on whether there existed previous matches for that protein. These 85% cut-offs minimized spurious matches without losing many real ones, as deduced from a test dataset (data not shown). If no adequate peptide matches were found, then the second assessment would admit matches with values of at least 85% for both length coverage and sequence identity but only involving proteins previously identified. In the final and least strict assignment, protein-peptide matches were accepted if the protein had previously been identified and the matches had either 85% length coverage or sequence identity. At present, the peptide matching procedure does not penalize disagreements between the flanking sequences from the proteins matched by the peptides with respect to those expected by the characteristic digestion pattern of the protease employed in the experiment.
Phosphorylated sites obtained from literature and other resources were handled in a similar way because they were actually stored in the database as artificially long peptide sequences, extracted from their original source entries, with the modified residues situated at their central position. This facilitates the complete reassignment of all protein-peptide matches every time that a new version of Ensembl is released.
Clustering peptide-protein matches from each experiment into minimal protein lists
Following the matching of peptides to genes, it was necessary to analyze further those peptides with multiple matches, either to more than one protein from a single gene or to proteins from different genes, in order to determine whether the assignment could be further restricted to a single peptide-protein or peptide-gene pair, respectively. Thus, after the initial assignment of every peptide to one or more proteins, the peptides from an experiment were grouped according to the matched proteins in order to find a minimal list of proteins that could explain all of the matches observed . A protein was regarded as part of the minimal list if at least one of its peptides was unique or if it represented a unique combination of peptides that was not a subgroup of the peptide set of another protein. When the experiment provides high sequence coverage for each protein, this assignment can be unambiguous. However, when the sequence coverage is low, as is sometimes the case for high-throughput data, it can be difficult to distinguish between splicing variants or cross-matches. Because of this, even though every matched protein will eventually be described as being part or not of a minimal list, the mtcPTM database still keeps all the matches originally found, which can be visualized online for further critical inspection (Figure 1).
Structural modeling of protein domains containing phosphorylatable residues
Only proteins belonging to the minimal lists were considered for homology modeling. The feasibility of building models for these proteins, or for fragments of them, was determined according to their significant sequence similarity to available high-resolution structures. The sequence set of determined structures was obtained from Protein Data Bank (PDB) co-ordinates . Positions for which no electron density was available were given the amino acid letter code 'X'. The final set included a set of redundant sequences containing all the full-length sequences from PDB structures as well as those of their individual domains as defined by Structural Classification of Proteins (SCOP) [77, 78]. This redundancy was necessary because our fold recognition procedure mandated that the regions matched onto the protein queries should cover the entire length of the structural units in order to avoid missing structural elements in the three-dimensional models.
The structural assignment was performed as follows. Each protein sequence was first BLASTed against the structural database. Matches were initially considered adequate if they covered 80% of the length of the structural unit sharing no less than 30% sequence identity. Subsequently, the regions matched on the protein were extended in both directions to account for the full-length of the structure, and both sequences were realigned using ClustalX . The extension/alignment procedure was iterated until the protein covered exactly the entire length of the structure. After convergence, all pair-wise alignments with less than 30% sequence identity (including gaps) or more than 5% gapped regions in both sequences were discarded. If two possible models from the same protein covered exactly the same set of modified residues, only that with the highest sequence identity to their structural template was considered. The whole procedure was automatic and therefore the pair-wise alignments were not inspected or adjusted manually. However, given the strict sequence similarity and length coverage cut-offs imposed, it is likely that most alignments will be correct. The homology models were built from the pair-wise alignments using MODELLER with default values .
Creating a nonredundant set of structural models
To make a nonredundant set of structural models the sequences of all the modeled structures were clustered using BLASTclust with cut-offs of 40% sequence similarity and 80% length coverage. After the clustering, the member with the highest number of phosphorylated sites (or highest sequence identity to the template if their number of phosphosites was identical) from each group was taken. The nonredundant set used in the structural analysis (and provided in the Additional data files) was based on the mtcPTM release corresponding to Ensembl version 40. Updated versions of the structural models can be found in the live site.
Additional data files
JLJ, BH, JH, JMP are funded by the 6th Framework Programme of the European Union via the Integrated Project MitoCheck; and RD and the Sanger Institute by the Wellcome Trust.
- Lander ES, Linton LM, Birren B, Nusbaum C, Zody MC, Baldwin J, Devon K, Dewar K, Doyle M, FitzHugh W, et al: Initial sequencing and analysis of the human genome. Nature. 2001, 409: 860-921. 10.1038/35057062.PubMedView ArticleGoogle Scholar
- Venter JC, Adams MD, Myers EW, Li PW, Mural RJ, Sutton GG, Smith HO, Yandell M, Evans CA, Holt RA, et al: The sequence of the human genome. Science. 2001, 291: 1304-1351. 10.1126/science.1058040.PubMedView ArticleGoogle Scholar
- Goll J, Uetz P: The elusive yeast interactome. Genome Biol. 2006, 7: 223-PubMedPubMed CentralView ArticleGoogle Scholar
- Lee NH: Genomic approaches for reconstructing gene networks. Pharmacogenomics. 2005, 6: 245-258. 10.1517/146224188.8.131.52.PubMedView ArticleGoogle Scholar
- Mansuy IM, Shenolikar S: Protein serine/threonine phosphatases in neuronal plasticity and disorders of learning and memory. Trends Neurosci. 2006, 29: 679-686. 10.1016/j.tins.2006.10.004.PubMedView ArticleGoogle Scholar
- Pawson T, Scott JD: Protein phosphorylation in signaling: 50 years and counting. Trends Biochem Sci. 2005, 30: 286-290. 10.1016/j.tibs.2005.04.013.PubMedView ArticleGoogle Scholar
- Cohen P: Dissection of the protein phosphorylation cascades involved in insulin and growth factor action. Biochem Soc Trans. 1993, 21: 555-567.PubMedView ArticleGoogle Scholar
- Cohen P: The regulation of protein function by multisite phosphorylation - a 25 year update. Trends Biochem Sci. 2000, 25: 596-601. 10.1016/S0968-0004(00)01712-6.PubMedView ArticleGoogle Scholar
- Attwood PV, Piggott MJ, Zu XL, Besant PG: Focus on phosphohistidine. Amino Acids. 2007, 32: 145-156. 10.1007/s00726-006-0443-6.PubMedView ArticleGoogle Scholar
- Lindberg RA, Quinn AM, Hunter T: Dual-specificity protein kinases: will any hydroxyl do?. Trends Biochem Sci. 1992, 17: 114-119. 10.1016/0968-0004(92)90248-8.PubMedView ArticleGoogle Scholar
- Stroud RM: Mechanisms of biological control by phosphorylation. Curr Opin Struct Biol. 1991, 1: 826-835. 10.1016/0959-440X(91)90186-W.View ArticleGoogle Scholar
- Kobe B, Kampmann T, Forwood JK, Listwan P, Brinkworth RI: Substrate specificity of protein kinases and computational prediction of substrates. Biochim Biophys Acta. 2005, 1754: 200-209.PubMedView ArticleGoogle Scholar
- Manning G, Whyte DB, Martinez R, Hunter T, Sudarsanam S: The protein kinase complement of the human genome. Science. 2002, 298: 1912-1934. 10.1126/science.1075762.PubMedView ArticleGoogle Scholar
- Park J, Hu Y, Murthy TV, Vannberg F, Shen B, Rolfs A, Hutti JE, Cantley LC, Labaer J, Harlow E, et al: Building a human kinase gene repository: bioinformatics, molecular cloning, and functional validation. Proc Natl Acad Sci USA. 2005, 102: 8114-8119. 10.1073/pnas.0503141102.PubMedPubMed CentralView ArticleGoogle Scholar
- Kreegipuu A, Blom N, Brunak S: PhosphoBase, a database of phosphorylation sites: release 2.0. Nucleic Acids Res. 1999, 27: 237-239. 10.1093/nar/27.1.237.PubMedPubMed CentralView ArticleGoogle Scholar
- Johnson SA, Hunter T: Kinomics: methods for deciphering the kinome. Nat Methods. 2005, 2: 17-25. 10.1038/nmeth731.PubMedView ArticleGoogle Scholar
- Mumby M, Brekken D: Phosphoproteomics: new insights into cellular signaling. Genome Biol. 2005, 6: 230-10.1186/gb-2005-6-9-230.PubMedPubMed CentralView ArticleGoogle Scholar
- Johnson LN, Barford D: The effects of phosphorylation on the structure and function of proteins. Annu Rev Biophys Biomol Struct. 1993, 22: 199-232. 10.1146/annurev.bb.22.060193.001215.PubMedView ArticleGoogle Scholar
- Johnson LN, Lewis RJ: Structural basis for control by phosphorylation. Chem Rev. 2001, 101: 2209-2242. 10.1021/cr000225s.PubMedView ArticleGoogle Scholar
- Krupa A, Preethi G, Srinivasan N: Structural modes of stabilization of permissive phosphorylation sites in protein kinases: distinct strategies in Ser/Thr and Tyr kinases. J Mol Biol. 2004, 339: 1025-1039. 10.1016/j.jmb.2004.04.043.PubMedView ArticleGoogle Scholar
- Herzberg O, Reddy P, Sutrina S, Saier MH, Reizer J, Kapadia G: Structure of the histidine-containing phosphocarrier protein HPr from Bacillus subtilis at 2.0-A resolution. Proc Natl Acad Sci USA. 1992, 89: 2499-2503. 10.1073/pnas.89.6.2499.PubMedPubMed CentralView ArticleGoogle Scholar
- Hurley JH, Dean AM, Thorsness PE, Koshland DE, Stroud RM: Regulation of isocitrate dehydrogenase by phosphorylation involves no long-range conformational change in the free enzyme. J Biol Chem. 1990, 265: 3599-3602.PubMedGoogle Scholar
- Becker S, Groner B, Muller CW: Three-dimensional structure of the Stat3beta homodimer bound to DNA. Nature. 1998, 394: 145-151. 10.1038/28980.PubMedView ArticleGoogle Scholar
- Chen X, Vinkemeier U, Zhao Y, Jeruzalmi D, Darnell JE, Kuriyan J: Crystal structure of a tyrosine phosphorylated STAT-1 dimer bound to DNA. Cell. 1998, 93: 827-839. 10.1016/S0092-8674(00)81443-9.PubMedView ArticleGoogle Scholar
- Kovacs H, Comfort D, Lord M, Campbell ID, Yudkin MD: Solution structure of SpoIIAA, a phosphorylatable component of the system that regulates transcription factor sigmaF of Bacillus subtilis. Proc Natl Acad Sci USA. 1998, 95: 5067-5071. 10.1073/pnas.95.9.5067.PubMedPubMed CentralView ArticleGoogle Scholar
- Barford D, Hu SH, Johnson LN: Structural mechanism for glycogen phosphorylase control by phosphorylation and AMP. J Mol Biol. 1991, 218: 233-260. 10.1016/0022-2836(91)90887-C.PubMedView ArticleGoogle Scholar
- Barford D, Johnson LN: The allosteric transition of glycogen phosphorylase. Nature. 1989, 340: 609-616. 10.1038/340609a0.PubMedView ArticleGoogle Scholar
- Knighton DR, Zheng JH, Ten Eyck LF, Ashford VA, Xuong NH, Taylor SS, Sowadski JM: Crystal structure of the catalytic subunit of cyclic adenosine monophosphate-dependent protein kinase. Science. 1991, 253: 407-414. 10.1126/science.1862342.PubMedView ArticleGoogle Scholar
- Knighton DR, Zheng JH, Ten Eyck LF, Xuong NH, Taylor SS, Sowadski JM: Structure of a peptide inhibitor bound to the catalytic subunit of cyclic adenosine monophosphate-dependent protein kinase. Science. 1991, 253: 414-420. 10.1126/science.1862343.PubMedView ArticleGoogle Scholar
- Canagarajah BJ, Khokhlatchev A, Cobb MH, Goldsmith EJ: Activation mechanism of the MAP kinase ERK2 by dual phosphorylation. Cell. 1997, 90: 859-869. 10.1016/S0092-8674(00)80351-7.PubMedView ArticleGoogle Scholar
- Zheng N, Wang P, Jeffrey PD, Pavletich NP: Structure of a c-Cbl-UbcH7 complex: RING domain function in ubiquitin-protein ligases. Cell. 2000, 102: 533-539. 10.1016/S0092-8674(00)00057-X.PubMedView ArticleGoogle Scholar
- Antz C, Bauer T, Kalbacher H, Frank R, Covarrubias M, Kalbitzer HR, Ruppersberg JP, Baukrowitz T, Fakler B: Control of K+ channel gating by protein phosphorylation: structural switches of the inactivation gate. Nat Struct Biol. 1999, 6: 146-150. 10.1038/5833.PubMedView ArticleGoogle Scholar
- Diella F, Cameron S, Gemund C, Linding R, Via A, Kuster B, Sicheritz-Ponten T, Blom N, Gibson TJ: Phospho.ELM: a database of experimentally verified phosphorylation sites in eukaryotic proteins. BMC Bioinformatics. 2004, 5: 79-10.1186/1471-2105-5-79.PubMedPubMed CentralView ArticleGoogle Scholar
- Hornbeck PV, Chabra I, Kornhauser JM, Skrzypek E, Zhang B: PhosphoSite: a bioinformatics resource dedicated to physiological protein phosphorylation. Proteomics. 2004, 4: 1551-1561. 10.1002/pmic.200300772.PubMedView ArticleGoogle Scholar
- Lee TY, Huang HD, Hung JH, Huang HY, Yang YS, Wang TH: dbPTM: an information repository of protein post-translational modification. Nucleic Acids Res. 2006, D622-627. 10.1093/nar/gkj083. 34 Database
- MitoCheck web site. [http://www.mitocheck.org]
- Birney E, Andrews D, Caccamo M, Chen Y, Clarke L, Coates G, Cox T, Cunningham F, Curwen V, Cutts T, et al: Ensembl 2006. Nucleic Acids Res. 2006, D556-D561. 10.1093/nar/gkj133. 34 Database
- Schwartz D, Gygi SP: An iterative statistical approach to the identification of protein phosphorylation motifs from large-scale data sets. Nat Biotechnol. 2005, 23: 1391-1398. 10.1038/nbt1146.PubMedView ArticleGoogle Scholar
- Sayle RA, Milner-White EJ: RASMOL: biomolecular graphics for all. Trends Biochem Sci. 1995, 20: 374-10.1016/S0968-0004(00)89080-5.PubMedView ArticleGoogle Scholar
- Finn RD, Mistry J, Schuster-Bockler B, Griffiths-Jones S, Hollich V, Lassmann T, Moxon S, Marshall M, Khanna A, Durbin R, et al: Pfam: clans, web tools and services. Nucleic Acids Res. 2006, D247-D251. 10.1093/nar/gkj149. 34 Database
- Iakoucheva LM, Radivojac P, Brown CJ, O'Connor TR, Sikes JG, Obradovic Z, Dunker AK: The importance of intrinsic disorder for protein phosphorylation. Nucleic Acids Res. 2004, 32: 1037-1049. 10.1093/nar/gkh253.PubMedPubMed CentralView ArticleGoogle Scholar
- Brinkworth RI, Breinl RA, Kobe B: Structural basis and prediction of substrate specificity in protein serine/threonine kinases. Proc Natl Acad Sci USA. 2003, 100: 74-79. 10.1073/pnas.0134224100.PubMedPubMed CentralView ArticleGoogle Scholar
- Kabsch W, Sander C: Dictionary of protein secondary structure: pattern recognition of hydrogen-bonded and geometrical features. Biopolymers. 1983, 22: 2577-2637. 10.1002/bip.360221211.PubMedView ArticleGoogle Scholar
- Hubbard SJ, Campbell SF, Thornton JM: Molecular recognition. Conformational analysis of limited proteolytic sites and serine proteinase protein inhibitors. J Mol Biol. 1991, 220: 507-530. 10.1016/0022-2836(91)90027-4.PubMedView ArticleGoogle Scholar
- Rush J, Moritz A, Lee KA, Guo A, Goss VL, Spek EJ, Zhang H, Zha XM, Polakiewicz RD, Comb MJ: Immunoaffinity profiling of tyrosine phosphorylation in cancer cells. Nat Biotechnol. 2005, 23: 94-101. 10.1038/nbt1046.PubMedView ArticleGoogle Scholar
- Pekarsky Y, Hallas C, Palamarchuk A, Koval A, Bullrich F, Hirata Y, Bichi R, Letofsky J, Croce CM: Akt phosphorylates and regulates the orphan nuclear receptor Nur77. Proc Natl Acad Sci USA. 2001, 98: 3690-3694. 10.1073/pnas.051003198.PubMedPubMed CentralView ArticleGoogle Scholar
- Collins MO, Yu L, Coba MP, Husi H, Campuzano I, Blackstock WP, Choudhary JS, Grant SG: Proteomic analysis of in vivo phosphorylated synaptic proteins. J Biol Chem. 2005, 280: 5972-5982. 10.1074/jbc.M411220200.PubMedView ArticleGoogle Scholar
- Dosztanyi Z, Chen J, Dunker AK, Simon I, Tompa P: Disorder and sequence repeats in hub proteins and their implications for network evolution. J Proteome Res. 2006, 5: 2985-2995. 10.1021/pr060171o.PubMedView ArticleGoogle Scholar
- Jimenez JL: Does structural and chemical divergence play a role in precluding undesirable protein interactions?. Proteins. 2005, 59: 757-764. 10.1002/prot.20448.PubMedView ArticleGoogle Scholar
- Chothia C, Lesk AM: The relation between the divergence of sequence and structure in proteins. EMBO J. 1986, 5: 823-826.PubMedPubMed CentralGoogle Scholar
- Wilson CA, Kreychman J, Gerstein M: Assessing annotation transfer for genomics: quantifying the relations between protein sequence, structure and function through traditional and probabilistic scores. J Mol Biol. 2000, 297: 233-249. 10.1006/jmbi.2000.3550.PubMedView ArticleGoogle Scholar
- Cowan-Jacob SW, Fendrich G, Manley PW, Jahnke W, Fabbro D, Liebetanz J, Meyer T: The crystal structure of a c-Src complex in an active conformation suggests possible steps in c-Src activation. Structure. 2005, 13: 861-871. 10.1016/j.str.2005.03.012.PubMedView ArticleGoogle Scholar
- Xu W, Doshi A, Lei M, Eck MJ, Harrison SC: Crystal structures of c-Src reveal features of its autoinhibitory mechanism. Mol Cell. 1999, 3: 629-638. 10.1016/S1097-2765(00)80356-1.PubMedView ArticleGoogle Scholar
- Beausoleil SA, Jedrychowski M, Schwartz D, Elias JE, Villen J, Li J, Cohn MA, Cantley LC, Gygi SP: Large-scale characterization of HeLa cell nuclear phosphoproteins. Proc Natl Acad Sci USA. 2004, 101: 12130-12135. 10.1073/pnas.0404720101.PubMedPubMed CentralView ArticleGoogle Scholar
- Liu L, Song X, He D, Komma C, Kita A, Virbasius JV, Huang G, Bellamy HD, Miki K, Czech MP, et al: Crystal structure of the C2 domain of class II phosphatidylinositide 3-kinase C2alpha. J Biol Chem. 2006, 281: 4254-4260. 10.1074/jbc.M510791200.PubMedView ArticleGoogle Scholar
- Jimenez JL, Smith GR, Contreras-Moreira B, Sgouros JG, Meunier FA, Bates PA, Schiavo G: Functional recycling of C2 domains throughout evolution: a comparative study of synaptotagmin, protein kinase C and phospholipase C by sequence, structural and modelling approaches. J Mol Biol. 2003, 333: 621-639. 10.1016/j.jmb.2003.08.052.PubMedView ArticleGoogle Scholar
- Joel PB, Smith J, Sturgill TW, Fisher TL, Blenis J, Lannigan DA: pp90rsk1 regulates estrogen receptor-mediated transcription through phosphorylation of Ser-167. Mol Cell Biol. 1998, 18: 1978-1984.PubMedPubMed CentralView ArticleGoogle Scholar
- Misura KM, Scheller RH, Weis WI: Three-dimensional structure of the neuronal-Sec1-syntaxin 1a complex. Nature. 2000, 404: 355-362. 10.1038/35006120.PubMedView ArticleGoogle Scholar
- Picciotto MR, Cohn JA, Bertuzzi G, Greengard P, Nairn AC: Phosphorylation of the cystic fibrosis transmembrane conductance regulator. J Biol Chem. 1992, 267: 12742-12752.PubMedGoogle Scholar
- Lewis HA, Zhao X, Wang C, Sauder JM, Rooney I, Noland BW, Lorimer D, Kearins MC, Conners K, Condon B, et al: Impact of the deltaF508 mutation in first nucleotide-binding domain of human cystic fibrosis transmembrane conductance regulator on domain folding and structure. J Biol Chem. 2005, 280: 1346-1353. 10.1074/jbc.M410968200.PubMedView ArticleGoogle Scholar
- el Benna J, Faust LP, Babior BM: The phosphorylation of the respiratory burst oxidase component p47phox during neutrophil activation. Phosphorylation of sites recognized by protein kinase C and by proline-directed kinases. J Biol Chem. 1994, 269: 23431-23436.PubMedGoogle Scholar
- Groemping Y, Lapouge K, Smerdon SJ, Rittinger K: Molecular basis of phosphorylation-induced activation of the NADPH oxidase. Cell. 2003, 113: 343-355. 10.1016/S0092-8674(03)00314-3.PubMedView ArticleGoogle Scholar
- Ogura K, Nobuhisa I, Yuzawa S, Takeya R, Torikai S, Saikawa K, Sumimoto H, Inagaki F: NMR solution structure of the tandem Src homology 3 domains of p47phox complexed with a p22phox-derived proline-rich peptide. J Biol Chem. 2006, 281: 3660-3668. 10.1074/jbc.M505193200.PubMedView ArticleGoogle Scholar
- Rosengarth A, Gerke V, Luecke H: X-ray structure of full-length annexin 1 and implications for membrane aggregation. J Mol Biol. 2001, 306: 489-498. 10.1006/jmbi.2000.4423.PubMedView ArticleGoogle Scholar
- Weng X, Luecke H, Song IS, Kang DS, Kim SH, Huber R: Crystal structure of human annexin I at 2.5 A resolution. Protein Sci. 1993, 2: 448-458.PubMedPubMed CentralView ArticleGoogle Scholar
- Derrien A, Druey KM: RGS16 function is regulated by epidermal growth factor receptor-mediated tyrosine phosphorylation. J Biol Chem. 2001, 276: 48532-48538.PubMedGoogle Scholar
- Terrak M, Kerff F, Langsetmo K, Tao T, Dominguez R: Structural basis of protein phosphatase 1 regulation. Nature. 2004, 429: 780-784. 10.1038/nature02582.PubMedView ArticleGoogle Scholar
- Honbou K, Suzuki NN, Horiuchi M, Niki T, Taira T, Ariga H, Inagaki F: The crystal structure of DJ-1, a protein related to male fertility and Parkinson's disease. J Biol Chem. 2003, 278: 31380-31384. 10.1074/jbc.M305878200.PubMedView ArticleGoogle Scholar
- Huai Q, Sun Y, Wang H, Chin LS, Li L, Robinson H, Ke H: Crystal structure of DJ-1/RS and implication on familial Parkinson's disease. FEBS Lett. 2003, 549: 171-175. 10.1016/S0014-5793(03)00764-6.PubMedView ArticleGoogle Scholar
- Spahn CM, Gomez-Lorenzo MG, Grassucci RA, Jorgensen R, Andersen GR, Beckmann R, Penczek PA, Ballesta JP, Frank J: Domain movements of elongation factor eEF2 and the eukaryotic 80S ribosome facilitate tRNA translocation. EMBO J. 2004, 23: 1008-1019. 10.1038/sj.emboj.7600102.PubMedPubMed CentralView ArticleGoogle Scholar
- mtcPTM database. [http://www.mtocheck.org/cgi-bin/mtcPTM/search]
- Blom N, Gammeltoft S, Brunak S: Sequence and structure-based prediction of eukaryotic protein phosphorylation sites. J Mol Biol. 1999, 294: 1351-1362. 10.1006/jmbi.1999.3310.PubMedView ArticleGoogle Scholar
- Wu CH, Apweiler R, Bairoch A, Natale DA, Barker WC, Boeckmann B, Ferro S, Gasteiger E, Huang H, Lopez R, et al: The Universal Protein Resource (UniProt): an expanding universe of protein information. Nucleic Acids Res. 2006, D187-D191. 10.1093/nar/gkj161. 34 Database
- Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.PubMedPubMed CentralView ArticleGoogle Scholar
- Nesvizhskii AI, Aebersold R: Interpretation of shotgun proteomic data: the protein inference problem. Mol Cell Proteomics. 2005, 4: 419-1440. 10.1074/mcp.R500005-MCP200.View ArticleGoogle Scholar
- Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H, Shindyalov IN, Bourne PE: The Protein Data Bank. Nucleic Acids Res. 2000, 28: 235-242. 10.1093/nar/28.1.235.PubMedPubMed CentralView ArticleGoogle Scholar
- Andreeva A, Howorth D, Brenner SE, Hubbard TJ, Chothia C, Murzin AG: SCOP database in 2004: refinements integrate structure and sequence family data. Nucleic Acids Res. 2004, D226-D229. 10.1093/nar/gkh039. 32 Database
- Chandonia JM, Hon G, Walker NS, Lo Conte L, Koehl P, Levitt M, Brenner SE: The ASTRAL Compendium in 2004. Nucleic Acids Res. 2004, D189-D192. 10.1093/nar/gkh034. 32 Database
- Aiyar A: The use of CLUSTAL W and CLUSTAL X for multiple sequence alignment. Methods Mol Biol. 2000, 132: 221-241.PubMedGoogle Scholar
- Sali A, Blundell TL: Comparative protein modelling by satisfaction of spatial restraints. J Mol Biol. 1993, 234: 779-815. 10.1006/jmbi.1993.1626.PubMedView ArticleGoogle Scholar
- Scheffzek K, Stephan I, Jensen ON, Illenberger D, Gierschik P: The Rac-RhoGDI complex and the structural basis for the regulation of Rho proteins by RhoGDI. Nat Struct Biol. 2000, 7: 122-126. 10.1038/72392.PubMedView ArticleGoogle Scholar
- Meinke G, Sigler PB: DNA-binding mechanism of the monomeric orphan nuclear receptor NGFI-B. Nat Struct Biol. 1999, 6: 471-477. 10.1038/8276.PubMedView ArticleGoogle Scholar
- Appleby TC, Erion MD, Ealick SE: The structure of human 5'-deoxy-5'-methylthioadenosine phosphorylase at 1.7 A resolution provides insights into substrate binding and catalysis. Structure. 1999, 7: 629-641. 10.1016/S0969-2126(99)80084-7.PubMedView ArticleGoogle Scholar
- Davey CA, Sargent DF, Luger K, Maeder AW, Richmond TJ: Solvent mediated interactions in the structure of the nucleosome core particle at 1.9 a resolution. J Mol Biol. 2002, 319: 1097-1113. 10.1016/S0022-2836(02)00386-8.PubMedView ArticleGoogle Scholar
- Schwabe JW, Chapman L, Finch JT, Rhodes D: The crystal structure of the estrogen receptor DNA-binding domain bound to DNA: how receptors discriminate between their response elements. Cell. 1993, 75: 567-578. 10.1016/0092-8674(93)90390-C.PubMedView ArticleGoogle Scholar
- Schalk I, Zeng K, Wu SK, Stura EA, Matteson J, Huang M, Tandon A, Wilson IA, Balch WE: Structure and mutational analysis of Rab GDP-dissociation inhibitor. Nature. 1996, 381: 42-48. 10.1038/381042a0.PubMedView ArticleGoogle Scholar
- Andersen GR, Pedersen L, Valente L, Chatterjee I, Kinzy TG, Kjeldgaard M, Nyborg J: Structural basis for nucleotide exchange and competition with tRNA in the yeast elongation factor complex eEF1A:eEF1Balpha. Mol Cell. 2000, 6: 1261-1266. 10.1016/S1097-2765(00)00122-2.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.