Sequence-similar, structure-dissimilar protein pairs in the PDB

It is often assumed that in the Protein Data Bank (PDB), two proteins with similar sequences will also have similar structures. Accordingly, it has proved useful to develop subsets of the PDB from which "redundant" structures have been removed, based on a sequence-based criterion for similarity. Similarly, when predicting protein structure using homology modeling, if a template structure for modeling a target sequence is selected by sequence alone, this implicitly assumes that all sequence-similar templates are equivalent. Here, we show that this assumption is often not correct and that standard approaches to create subsets of the PDB can lead to the loss of structurally and functionally important information. We have carried out sequence-based structural superpositions and geometry-based structural alignments of a large number of protein pairs to determine the extent to which sequence similarity ensures structural similarity. We find many examples where two proteins that are similar in sequence have structures that differ significantly from one another. The source of the structural differences usually has a functional basis. The number of such proteins pairs that are identified and the magnitude of the dissimilarity depend on the approach that is used to calculate the differences; in particular sequence-based structure superpositioning will identify a larger number of structurally dissimilar pairs than geometry-based structural alignments. When two sequences can be aligned in a statistically meaningful way, sequence-based structural superpositioning provides a meaningful measure of structural differences. This approach and geometry-based structure alignments reveal somewhat different information and one or the other might be preferable in a given application. Our results suggest that in some cases, notably homology modeling, the common use of nonredundant datasets, culled from the PDB based on sequence, may mask important structural and functional information. We have established a data base of sequence-similar, structurally dissimilar protein pairs that will help address this problem (http://luna.bioc.columbia.edu/rachel/seqsimstrdiff.htm).     

Efflux Pumps Represent Possible Evolutionary Convergence onto the β-Barrel Fold

1. Download high-res image (271KB)
2. Download full-size image

     

Protein structure comparison: implications for the nature of 'fold space', and structure and function prediction

The identification of geometric relationships between protein structures offers a powerful approach to predicting the structure and function of proteins. Methods to detect such relationships range from human pattern recognition to a variety of mathematical algorithms. A number of schemes for the classification of protein structure have found widespread use and these implicitly assume the organization of protein structure space into discrete categories. Recently, an alternative view has emerged in which protein fold space is seen as continuous and multidimensional. Significant relationships have been observed between proteins that belong to what have been termed different 'folds'. There has been progress in the use of these relationships in the prediction of protein structure and function.     

Characterization of the protein kinase activity of TRPM7/ChaK1, a protein kinase fused to the transient receptor potential ion channel

Channel-kinase TRPM7/ChaK1 is a member of a recently discovered family of protein kinases called alpha-kinases that display no sequence homology to conventional protein kinases. It is an unusual bifunctional protein that contains an alpha-kinase domain fused to an ion channel. The TRPM7/ChaK1 channel has been characterized using electrophysiological techniques, and recent evidence suggests that it may play a key role in the regulation of magnesium homeostasis. However, little is known about its protein kinase activity. To characterize the kinase activity of TRPM7/ChaK1, we expressed the kinase catalytic domain in bacteria. ChaK1-cat is able to undergo autophosphorylation and to phosphorylate myelin basic protein and histone H3 on serine and threonine residues. The kinase is specific for ATP and cannot use GTP as a substrate. ChaK1-cat is insensitive to staurosporine (up to 0.1 mM) but can be inhibited by rottlerin. Because the kinase domain is physically linked to an ion channel, we investigated the effect of ions on ChaK1-cat activity. The kinase requires Mg(2+) (optimum at 4-10 mM) or Mn(2+) (optimum at 3-5 mM), with activity in the presence of Mn(2+) being 2 orders of magnitude higher than in the presence of Mg(2+). Zn(2+) and Co(2+) inhibited ChaK1-cat kinase activity. Ca(2+) at concentrations up to 1 mM did not affect kinase activity. Considering intracellular ion concentrations, our results suggest that, among divalent metal ions, only Mg(2+) can directly modulate TRPM7/ChaK1 kinase activity in vivo.     

Genotype assignment in Gaucher disease by selective amplification of the active glucocerebrosidase gene.

Genomic DNA prepared from human cells in culture was amplified by the polymerase chain-reaction technique using two primers specific for the active human glucocerebrosidase gene. The 1,036-bp amplified fragment derived from the active gene was tested for the existence of three mutations--designated "370," "NciI," and "HhaI"--by allele-specific oligonucleotide hybridization. The results obtained from the cell lines examined permitted a clear distinction between homozygous affected, heterozygous, and normal genotypes. However, 28% of the possible affected loci were normal with respect to the three mutations, indicating the presence of additional mutations that remain to be elucidated. While the NciI mutation could be found in both Ashkenazi Jewish and non-Jewish type 1 patients, the only homozygotes with this mutation had the neurological (type 2 or type 3) form of the disease. The 370 mutation, on the other hand, was only present in type 1 patients and was not identified among any of the patients with neurologic forms of the disease.     

Molecular heterogeneity of late-onset forms of globoid-cell leukodystrophy.

Globoid-cell leukodystrophy (GLD) is an autosomal recessive inherited disorder caused by the deficiency of galactocerebrosidase, the lysosomal enzyme responsible for the degradation of the myelin glycolipid galactocerebroside. Although the most common form of the disease is the classical infantile form (Krabbe disease), later-onset forms also have been described. We have analyzed the galactocerebrosidase gene in 17 patients (nine families) with late-onset GLD and in 1 patient with classical Krabbe disease. Half of the patients were heterozygous for the large gene deletion associated with the 502C-->T polymorphism, the most common mutation in infantile patients. Several novel mutations that result in deficient galactocerebrosidase activity were also identified in these patients. They include the missense mutations R63H, G95S, M101L, G268S, Y298C, and I234T; the nonsense mutation S7X; a one-base deletion (805delG); a mutation that interferes with the splicing of intron 1; and a 34-nt insertion in the RNA, caused by the aberrant splicing of intron 6. All of these genetic defects are clustered in the first 10 exons of the galactocerebrosidase gene and therefore affect the 50-kD subunit of the mature enzyme. Studies on the distribution and enzymatic activity of the polymorphic alleles 1637T/C (I546/T546) provided support for previous data that had indicated the existence of two galactocerebrosidase forms with different catalytic activities in the general population. Our data also indicate that the mutations occur preferentially in the "low activity" 1637C allele.     

Metabolic activities in human skin fibroblasts preloaded with labeled GM2-ganglioside

Confluent cultures of human skin fibroblasts were maintained for 10 days with sphingosine labeled [3H]GM2. Labeled medium was then replaced with normal medium and the cells maintained for 42 days with weekly medium changes. Cells were harvested at regular intervals and cells, medium, and trypsin digest supernatant analyzed for [3H]GM2 and its metabolic products. The ganglioside can be membrane associated and removed by trypsin, or membrane incorporated and trypsin insensitive. The membrane incorporated material is apparently transported to the lysosomes slowly by membrane flow, where 80% of the cellular GM2 can be metabolized by day 42. [3H]GM2 as well as its metabolic products in control cells is continuously released into the medium, during which it can also become associated with the cell surface membrane. There is no detectable metabolism of the [3H]GM2 in GM2 gangliosidosis cell lines over the extended post-labeling period, indicating that there is no residual enzyme activity in these cells. Undegraded GM2 is continuously released into the medium and remains associated with the cell surface membrane as well.     

NMR study of the metabolic 15N isotopic enrichment of cyanophycin synthesized by the cyanobacterium Synechocystis sp. strain PCC 6308

1H, 13C and 15N nuclear magnetic resonance (NMR) spectroscopy has been used to characterize cyanophycin, a multi-l-arginyl-poly-[l-aspartic acid] polypeptide from the cyanobacterium Synechocystis sp. strain PCC 6308. 1H, 13C and 15N chemical shifts and 1JHN and 1JCN coupling constants were measured in isolated 15N-labeled cyanophycin, and showed chemical shift values and J-couplings consistent with the reported polypeptide structure. 15N enrichment levels were determined from the extent of 1H-15N J-coupling in 1H NMR spectra of cyanophycin. Similar experiments using 13C-15N coupling in 13C NMR spectra were not useful in determining enrichment levels.