A consensus prediction of the secondary structure for the 6-phospho-beta-D-galactosidase superfamily

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.

Crystal structure of the beta-glycosidase from the hyperthermophilic archeon Sulfolobus solfataricus: resilience as a key factor in thermostability

Enzymes from hyperthermophilic organisms must operate at temperatures which rapidly denature proteins from mesophiles. The structural basis of this thermostability is still poorly understood. Towards a further understanding of hyperthermostability, we have determined the crystal structure of the beta-glycosidase (clan GH-1A, family 1) from the hyperthermophilic archaeon Sulfolobus solfataricus at 2.6 A resolution. The enzyme is a tetramer with subunit molecular mass at 60 kDa, and crystallises with half of the tetramer in the asymmetric unit. The structure is a (betaalpha)8 barrel, but with substantial elaborations between the beta-strands and alpha-helices in each repeat. The active site occurs at the centre of the top face of the barrel and is connected to the surface by a radial channel which becomes a blind-ended tunnel in the tetramer, and probably acts as the binding site for extended oligosaccharide substrates. Analysis of the structure reveals two features which differ significantly from mesophile proteins; (1) an unusually large proportion of surface ion-pairs involved in networks that cross-link sequentially separate structures on the protein surface, and (2) an unusually large number of solvent molecules buried in hydrophilic cavities between sequentially separate structures in the protein core. These factors suggest a model for hyperthermostability via resilience rather than rigidity.

Improving protein secondary structure prediction with aligned homologous sequences

Most recent protein secondary structure prediction methods use sequence alignments to improve the prediction quality. We investigate the relationship between the location of secondary structural elements, gaps, and variable residue positions in multiple sequence alignments. We further investigate how these relationships compare with those found in structurally aligned protein families. We show how such associations may be used to improve the quality of prediction of the secondary structure elements, using the Quadratic-Logistic method with profiles. Furthermore, we analyze the extent to which the number of homologous sequences influences the quality of prediction. The analysis of variable residue positions shows that surprisingly, helical regions exhibit greater variability than do coil regions, which are generally thought to be the most common secondary structure elements in loops. However, the correlation between variability and the presence of helices does not significantly improve prediction quality. Gaps are a distinct signal for coil regions. Increasing the coil propensity for those residues occurring in gap regions enhances the overall prediction quality. Prediction accuracy increases initially with the number of homologues, but changes negligibly as the number of homologues exceeds about 14. The alignment quality affects the prediction more than other factors, hence a careful selection and alignment of even a small number of homologues can lead to significant improvements in prediction accuracy.

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.

The three-dimensional structure of 6-phospho-beta-galactosidase from Lactococcus lactis

BACKGROUND

The enzyme 6-phospho-beta-galactosidase hydrolyzes phospholactose, the product of a phosphor-enolpyruvate-dependent phosphotransferase system. It belongs to glycosidase family 1 and no structure has yet been published for a member of this family.

RESULTS

The crystal structure of 6-phospho-beta-galactosidase was determined at 2.3 A resolution by multiple isomorphous replacement, using the wild-type enzyme and a designed cysteine mutant. A second crystal form, found with the mutant enzyme, was solved by molecular replacement, yielding the conformation of two chain loops that are invisible in the first crystal form. The active center, located through catalytic residues identified in previous studies, cannot be accessed by the substrate if the two loops are in their defined conformation. The enzyme contains a (beta alpha)8 barrel and the relationship of its chain fold to that of other glycosidases has been quantified. As a side issue, we observed that a cysteine point mutant designed for X-ray analysis crystallized mainly as a symmetric dimer around an intermolecular disulfide bridge formed by the newly introduced cysteine.

CONCLUSIONS

The presented analysis provides a basis on which to model all other family 1 members and thereby will help in elucidating the catalytic mechanisms of these sequence-related enzymes. Moreover, this enzyme belongs to a superfamily of glycosidases sharing a (beta alpha)8 barrel with catalytic glutamates/aspartates at the ends of the fourth and the seventh strands of the beta barrel.

A preference-based free-energy parameterization of enzyme-inhibitor binding. Applications to HIV-1-protease inhibitor design

The interface between protein receptor-ligand complexes has been studied with respect to their binary interatomic interactions. Crystal structure data have been used to catalogue surfaces buried by atoms from each member of a bound complex and determine a statistical preference for pairs of amino-acid atoms. A simple free energy model of the receptor-ligand system is constructed from these atom-atom preferences and used to assess the energetic importance of interfacial interactions. The free energy approximation of binding strength in this model has a reliability of about +/- 1.5 kcal/mol, despite limited knowledge of the unbound states. The main utility of such a scheme lies in the identification of important stabilizing atomic interactions across the receptor-ligand interface. Thus, apart from an overall hydrophobic attraction (Young L, Jernigan RL, Covell DG, 1994, Protein Sci 3:717-729), a rich variety of specific interactions is observed. An analysis of 10 HIV-1 protease inhibitor complexes is presented that reveals a common binding motif comprised of energetically important contacts with a rather limited set of atoms. Design improvements to existing HIV-1 protease inhibitors are explored based on a detailed analysis of this binding motif.