A predicted consensus structure for the protein kinase C2 homology (C2H) domain, the repeating unit of synaptotagmin

Alignments grow, secondary structure prediction improves

Using information from sequence alignments significantly improves protein secondary structure prediction. Typically, more divergent profiles yield better predictions. Recently, various groups have shown that accuracy can be improved significantly by using PSI-BLAST profiles to develop new prediction methods. Here, we focused on the influences of various alignment strategies on two 8-year-old PHD methods. The following results stood out. (i) PHD using pairwise alignments predicts about 72% of all residues correctly in one of the three states: helix, strand, and other. Using larger databases and PSI-BLAST raised accuracy to 75%. (ii) More than 60% of the improvement originated from the growth of current sequence databases; about 20% resulted from detailed changes in the alignment procedure (substitution matrix, thresholds, and gap penalties). Another 20% of the improvement resulted from carefully using iterated PSI-BLAST searches. (iii) It is of interest that we failed to improve prediction accuracy further when attempting to refine the alignment by dynamic programming (MaxHom and ClustalW). (iv) Improvement through family growth appears to saturate at some point. However, most families have not reached this saturation. Hence, we anticipate that prediction accuracy will continue to rise with database growth.

Predicted structure and fold recognition for the glutamyl tRNA reductase family of proteins

The conserved residues of glutamyl tRNA reductase (GTR) from Hordeum vulgare (GTRhorvu) were found from an alignment/pile-up of 24 homologous sequences found using BLAST searches. A multiple alignment of sequences was used to obtain a prediction of the secondary structure of the GTR's. This secondary structure was submitted to the THREADER program to find possible homologous 3D structures. To help select the template for predicting the fold for GTRhorvu, we employed both molecular-biological and biochemical information about GTRhorvu. After fitting the secondary structure of GTRhorvu to the selected template, the MODELLER program was used to determine the fold for GTRhorvu. This model was built using the B subunit of succinyl CoA synthetase, 1scuB, as a template for the 3D structure of GTRhorvu. From the predicted structure, possible regions were identified for the binding of glutamyl-tRNA, NADPH and a heme inhibitor. The predicted structure was used to propose a detailed biochemical mechanism for the GTR, involving Mg catalyzed thioester formation and reduction by NADPH to glutamate-1-semialdehyde. Sites for these reactions are identified. The predicted structure has been deposited in the Brookhaven database as ID 1b61.

Evolutionary, mechanistic, and predictive analyses of the hydroxymethyldihydropterin pyrophosphokinase family of proteins

A prediction has been prepared ab initio for the secondary structure of the hydroxymethyldihydropterin pyrophosphokinase (HPPK) family of proteins starting from a set of aligned homologous protein sequences. Attempts to identify a fold by threading failed, judging by the inability to find a threading "hit" that had a secondary structure that was plausibly congruent to the predicted secondary structure for the HPPK family. Therefore, a set of tertiary structure models was assembled ab initio, where alternative models were built and used to select between alternative secondary structure models. This prediction report illustrates the importance of non-computational approaches to structure prediction at its present frontier, which is to obtain medium resolution models of tertiary structure.

Structure prediction in a post-genomic environment: a secondary and tertiary structural model for the initiation factor 5A family

Two predictions have been prepared for the fold of initiation factor 5A (IF5A) starting from a set of homologous sequences. In the first, a secondary structural model was predicted for the protein in 1994, when only eleven homologs (and no eubacterial homologs) had been sequenced. The second was made recently, after genome projects had generated a total of 33 sequences for the protein family from species of all three kingdoms of life. With the second set of sequences, but not with the first, it was possible to predict that the N-terminal domain of the protein folds in a possibly open beta-barrel/sandwich core structure, with a short helix capping one side of the barrel. We place the pair of predictions in the public domain before an experimental structure is known. This example illustrates the impact of genome sequencing projects on structure prediction from sequence alignments.

Structural aspects of the functional modules in human protein kinase-C alpha deduced from comparative analyses

Three-dimensional models of the five functional modules in human protein kinase C alpha (PKC alpha) have been generated on the basis of known related structures. The catalytic region at the C-terminus of the sequence and the N-terminal auto-inhibitory pseudo-substrate have been modeled using the crystal structure complex of cAMP-dependent protein kinase (cAPK) and PKI peptide. While the N-terminal helix of the catalytic region of PKC alpha is predicted to be in a different location compared with cAPK, the C-terminal extension is modeled like that in the cAPK. The predicted permissive phosphorylation site of PKC alpha, Thr 497, is found to be entirely consistent with the mutagenesis studies. Basic Lys and Arg residues in the pseudo-substrate make several specific interactions with acidic residues in the catalytic region and may interact with the permissive phosphorylation site. Models of the two zinc-binding modules of PKC alpha are based on nuclear magnetic resonance and crystal structures of such modules in other PKC isoforms while the calcium phospholipid binding module (C2) is based on the crystal structure of a repeating unit in synaptotagmin I. Phorbol ester binding regions in zinc-binding modules and the calcium binding region in the C2 domain are similar to those in the basis structures. A hypothetical model of the relative positions of all five modules has the putative lipid binding ends of the C2 and the two zinc-binding domains pointing in the same direction and may serve as a basis for further experiments.

The phospho-beta-galactosidase and synaptotagmin predictions

Two bona fide consensus predictions of secondary and tertiary structure in a protein family, made and announced before experimental structures were known, are evaluated in light of the subsequently determined experimental structures. The first, for phospho-beta-galactosidase, identified the core strands of an 8-fold alpha-beta barrel, and identified the 8-fold alpha-beta barrel itself, which was found in the subsequently determined experimental structure to be the core folding domain. The second, for synaptotagmin, identified seven out of eight beta-strands in the structure correctly, missing only a noncore strand. Three preferred "topologies" were selected from several hundred thousand possible topologies of these seven predicted strands using a rule-based analysis. The subsequently determined experimental structure showed that these seven strands in synaptotagmin adopt one of the three preferred topologies. We were unable, however, to identify the correct topology from among these three topologies.