Characterization and comparison of protein structures. Part I-characterization.

The quantitative criteria characterizing the regularity of Calpha-backbones in the protein structures are presented. A technique is based on the Fourier remapping of the Cartesian coordinates for the Calpha-chain. The Fourier spectra identify the hidden periodicities and symmetries in protein structures, while the integral regularity is assessed via the spectral structural entropies. The formal unification of digitizing and the similarities in statistics for the random counterparts allow study of the direct correlations between the distribution of physico-chemical characteristics along the amino acid sequence and the spatial conformation of the polypeptide chain. The significant correlations are found for both hydrophobicity and side-chain volumes, though, as expected, the effects for hydrophobicity turn out essentially stronger. A scheme is illustrated by the set of 120 protein structures comprising the representatives from the main superfamilies and superfolds.     

Size-independent comparison of protein three-dimensional structures

Protein structures are routinely compared by their root-mean-square deviation (RMSD) in atomic coordinates after optimal rigid body superposition. What is not so clear is the significance of different RMSD values, particularly above the customary arbitrary cutoff for obvious similarity of 2–3 Å. Our earlier work argued for an intrinsic cutoff for protein similarity that varied with the number of residues in the polypeptide chains being compared. Here we introduce a new measure, ρ, of structural similarity based on RMSD that is independent of the sizes of the molecules involved, or of any other special properties of molecules. When ρ is less than 0.4–0.5, protein structures are visually recognized to be obviously similar, but the mathematically pleasing intrinsic cutoff of ρ>1.0 corresponds to overall similarity in folding motif at a level not usually recognized until smoothing of the polypeptide chain path makes it striking. When the structures are scaled to unit radius of gyration and equal principle moments of inertia, the comparisons are even more universal, since they are no longer obscured by differences in overall size and ellipticity. With increasing chain length, the distribution of ρ for pairs of random structures is skewed to higher values, but the value for the best 1% of the comparisons rises only slowly with the number of residues. This level is close to an intrinsic cutoff between similar and dissimilar comparisons, namely the maximal scaled ρ possible for the two structures to be more similar to each other than one is to the other's mirror image. The intrinsic cutoff is independent of the number of residues or points being compared. For proteins having fewer than 100 residues, the 1% ρ falls below the intrinsic cutoff, so that for very small proteins, geometrically significant similarity can often occur by chance. We believe these ideas will be helpful in judging success in NMR structure determination and protein folding modeling. © 1995 Wiley-Liss, Inc.     

Robust probabilistic superposition and comparison of protein structures

Background  

Protein structure comparison is a central issue in structural bioinformatics. The standard dissimilarity measure for protein structures is the root mean square deviation (RMSD) of representative atom positions such as α-carbons. To evaluate the RMSD the structures under comparison must be superimposed optimally so as to minimize the RMSD. How to evaluate optimal fits becomes a matter of debate, if the structures contain regions which differ largely - a situation encountered in NMR ensembles and proteins undergoing large-scale conformational transitions.     

Alternative alignments from comparison of protein structures

Comparison of two protein structures often results in not only a global alignment but also a number of distinct local alignments; the latter, referred to as alternative alignments, are however usually ignored in existing protein structure comparison analyses. Here, we used a novel method of protein structure comparison to extensively identify and characterize the alternative alignments obtained for structure pairs of a fold classification database. We showed that all alternative alignments can be classified into one of just a few types, and with which illustrated the potential of using alternative alignments to identify recurring protein substructures, including the internal structural repeats of a protein. Furthermore, we showed that among the alternative alignments obtained, permuted alignments, which included both circular and scrambled permutations, are as prevalent as topological alignments. These results demonstrated that the so far largely unattended alternative alignments of protein structures have implications and applications for research of protein classification and evolution.     

Characterization and comparison of membrane-associated and cytosolic cAMP-dependent protein kinases. Studies on human erythrocyte protein kinases.

Cyclic AMP-dependent protein kinase from human erythrocyte plasma membranes was solubilized with Triton X-100, partially purified, and systematically characterized by a series of physicochemical studies. Sedimentation and gel filtration experiments showed that the 6.6 S holoenzyme had a Stokes radius (a) of 5.7 nm and was dissociated into native 4.8 S cAMP-binding (a = 4.5 nm) and 3.2 S catalytic (a = 2.6 nm) subunits. A minimum subunit molecular weight of 48,000 was established for the regulatory subunit by photoaffinity labeling with 8-azido[32P]cAMP, sodium dodecyl sulfate-polyacrylamide gel electrophoresis, and autoradiography. These data suggest an asymmetric tetrameric (R2C2) structure (Mr approximately equal to 160,000) for the membrane-derived enzyme. Membrane-derived protein kinase was characterized as a type I enzyme on the basis of its R subunit molecular weight, pI values (R, 4.9; holoenzyme, 5.75 and 5.95), dissociation by 0.5 M NaCl and 50 microgram/ml of protamine, 20-fold reduced affinity for cAMP in the presence of 0.3 mM MgATP, elution from DEAE-cellulose at low ionic strength, and kinetic and cAMP-binding properties. The physicochemical properties of the membrane protein kinase closely parallel the characteristics of erythrocyte cytosolic protein kinase I but are clearly dissimilar from those of the soluble type II enzyme. Moreover, regulatory subunits of the membrane-associated and cytosolic type I kinases were indistinguishable in size, shape, subunit molecular weight, charge, binding and reassociation properties, and peptide maps of the photoaffinity-labeled cAMP-binding site, suggesting a high degree of structural and functional homology in this pair of enzymes. In view of the predominant occurrence of particulate type II protein kinases in rabbit heart and bovine cerebral cortex, the present results suggest that the distribution of membrane-associated protein kinases may be tissue- or species-specific, but not isoenzyme-specific.     

Characterization of novel proteins based on known protein structures

The genome sciences face the challenge to characterize structure and function of a vast number of novel genes. Sequence search techniques are used to infer functional and structural information from similarities to experimentally characterized genes or proteins. The persistent goal is to refine these techniques and to develop alternative and complementary methods to increase the range of reliable inference.Here, we focus on the structural and functional assignments that can be inferred from the known three-dimensional structures of proteins. The study uses all structures in the Protein Data Bank that were known by the end of 1997. The protein structures released in 1998 were then characterized in terms of functional and structural similarity to the previously known structures, yielding an estimate of the maximum amount of information on novel protein sequences that can be obtained from inference techniques.The 147 globular proteins corresponding to 196 domains released in 1998 have no clear sequence similarity to previously known structures. However, 75 % of the domains have extensive structure similarity to previously known folds, and most importantly, in two out of three cases similarity in structure coincides with related function. In view of this analysis, full utilization of existing structure data bases would provide information for many new targets even if the relationship is not accessible from sequence information alone. Currently, the most sophisticated techniques detect of the order of one-third of these relationships.     

Characterization of non-planar peptide groups in protein crystal structures

Peptide groups are generally assumed to be planar in protein structure, due to 'rigid' partial double bond character of peptide bonds, thus the value of peptide torsion angle omega should be restricted to 180 degrees for the usual trans form of peptide unit. However, on analyzing the ultra-high resolution protein crystal database, we find that in some cases, omega deviates >10 degrees from its usual value of 180 degrees, indicating significant non-planarity of peptide groups. Moreover, the non-planarity for most of the amino acids is found to be 'biased' towards values of omega smaller than 180 degrees. Similar trend for to is confirmed by the neutron diffraction data for proteins. The neutron diffraction database also reveals that non-planar peptide groups are generally correlated to 'pyramidal' structure of the peptide-nitrogen bonds. Consequently, the hydrogen atom of peptide group deviates from its planar position, as measured by the 'improper' torsion angle theta. Thus, we find that both the angles omega and theta point towards a significant amount of non-planarity of peptide groups, which cannot be ignored. The role of peptide nonplanarity in protein function is, however, not yet clear.     

A systematic approach to the comparison of protein structures

A systematic method has been developed for comparing the backbone conformations of proteins (Remington & Matthews, 1978). Two proteins are compared by successively optimizing the agreement between all possible segments of a chosen length from one protein, and all possible segments of the same length from the other protein. The method reveals any similarities between the two proteins, and provides an estimate of the statistical significance of any given structure agreement that is obtained.The method has been tested in a number of cases, including comparisons of the dehydrogenases and of the pancreatic and bacterial serine proteases. These examples were chosen to test the ability of the comparison method to detect structural similarities in the presence of large insertions and deletions. The results suggest that the detection of the “nucleotide binding fold” in the dehydrogenases is at the limit of the capability of the comparison technique in its original form, although it may be possible to generalize the method to allow for insertions and deletions in proteins.The results of many protein comparisons, made with different probe lengths, are summarized. For medium and long probe lengths, the average value of the structural agreement does not depend very much on the type of protein being compared. The average value of the structure agreement increases with the square root of the probe length, but for probe lengths above about 40 residues, the standard deviation is independent of probe length. From these observations it is possible to construct a generalized probability diagram to evaluate the significance of any structure agreement that might be obtained in comparing two proteins.     

Chemically accurate protein structures: Validation of protein NMR structures by comparison of measured and predicted pK a values

A new method is presented for evaluating the quality of protein structures obtained by NMR. This method exploits the dependence between measurable chemical properties of a protein, namely pK _a values of acidic residues, and protein structure. The accurate and fast empirical computational method employed by the PROPKA program () allows the user to test the ability of a given structure to reproduce known pK _a values, which in turn can be used as a criterion for the selection of more accurate structures. We demonstrate the feasibility of this novel idea for a series of proteins for which both␣NMR and X-ray structures, as well as pK _a values of all ionizable residues, have been determined. For the 17 NMR ensembles used in this study, this criterion is shown effective in the elimination of a large number of NMR structure ensemble members.     

Characterization of the backbone geometry of protein native state structures

We present a simple method for analyzing the geometry of noncovalent residue-residue interactions stabilizing the protein structure, which takes into account the constraints on the local backbone geometry. We find that the principal geometrical constraints are amino acid aspecific and are associated with hydrogen bond formation in helices and sheets. In contrast, amino acid residues in nonhelical and nonextended conformations, which make noncovalent interactions stabilizing the protein tertiary structure, display greater flexibility. We apply the method to an analysis of the packing of helices in helical bundle proteins requiring an efficient packing of amino acid side-chains of the interacting helices.     

Characterization and comparison of arcelin seed protein variants from common bean

Four variants of arcelin, an insecticidal seed storage protein of bean, Phaseolus vulgaris L., were investigated. Each variant (arcelin-1, -2, -3, and -4) was purified, and solubilities and M_rs were determined. For arcelins-1, -2, and -4, the isoelectric points, hemagglutinating activities, immunological cross-reactivities, and N-terminal amino acid sequences were determined. On the basis of native and denatured M_rs, the variants were classified as being composed of dimer protein (arcelin-2), tetramer protein (arcelins-3 and -4), or both dimer and tetramer proteins (arcelin-1). Although the dimer proteins (arcelins-1d and -2) could be distinguished by M_rs and isoelectric points, they were identical for their first 37 N-terminal amino acids and had similar immunological cross-reactions, and bean lines containing these variants had a DNA restriction fragment in common. The tetramer proteins arcelin-1t and arcelin-4 also could be distinguished from each other based on M_rs and isoelectric points; however, they had similar immunological cross-reactions and they were 77 to 93% identical for N-terminal amino acid composition. The similarities among arcelin variants, phytohemagglutinin, and a bean α-amylase inhibitor suggest that they are all encoded by related members of a lectin gene family.     

Evolution of protein sequences and structures.

The relationship between sequence similarity and structural similarity has been examined in 36 protein families with five or more diverse members whose structures are known. The structural similarity within a family (as determined with the DALI structure comparison program) is linearly related to sequence similarity (as determined by a Smith-Waterman search of the protein sequences in the structure database). The correlation between structural similarity and sequence similarity is very high; 18 of the 36 families had linear correlation coefficients r>/=0.878, and only nine had correlation coefficients r相似文献   

LOPAL and SCAMP: techniques for the comparison and display of protein structures

This paper describes two computer programs designed to assist in the comparison of protein structures. LOPAL (LOoP ALignment) applies a dynamic programming algorithm to the comparison of regions of protein three-dimensional (3D) structure and gives a similarity score and suggested sequence alignment with that score. SCAMP (Structure Comparison and Alignment of Multiple Proteins) is an interactive graphics program for the Evans and Sutherland PS300 graphics terminal that allows the simultaneous display, manipulation and pairwise least-squares fitting of up to nine independent structures. Together, LOPAL and SCAMP provide an integrated system for characterizing structural similarities in proteins with the aim of improving the accuracy of predicted protein structures. An application of these programs to loop regions in the immunoglobulin constant domains is illustrated.     

Characterization and comparison of chloramphenicol acetyltransferase variants.

1. Variants of chloramphenicol acetyltransferase from a variety of bacterial species have been isolated and purified to homogeneity. They constitute a heterogeneous group of proteins as judged by analytical affinity and hydrophobic ('detergent') chromatography, native and sodium dodecyl sulfate electrophoresis, sensitivity to sulfhydryl specific reagents, steady state kinetic analysis, and reaction with antisera. 2. The most striking observation is that three variants of chloramphenicol acetyltransferase (R factor type III, Streptomyces acrimycini, and Agrobacterium tumefaciens) possess an apparent subunit molecular weight (24,500) which is significantly greater than that of all other variants examined (22,500). The three atypical variants are not identical since they show marked differences in a number of important parameters. 3. Although the fundamental mechanism of catalysis may prove to be identical for all chloramphenicol acetyltransferase variants, there is a wide range of sensitivity to thiol-directed inhibitors among the enzymes studied. 4. Amino acid sequence analysis of the N-termini of selected variants suggests that the qualitative differences among chloramphenicol acetyltransferase variants is a reflection of structural heterogeneity which is most marked in comparisons between variants from Gram-positive and Gram-negative species.     

Biological insights from topology independent comparison of protein 3D structures

Comparing and classifying the three-dimensional (3D) structures of proteins is of crucial importance to molecular biology, from helping to determine the function of a protein to determining its evolutionary relationships. Traditionally, 3D structures are classified into groups of families that closely resemble the grouping according to their primary sequence. However, significant structural similarities exist at multiple levels between proteins that belong to these different structural families. In this study, we propose a new algorithm, CLICK, to capture such similarities. The method optimally superimposes a pair of protein structures independent of topology. Amino acid residues are represented by the Cartesian coordinates of a representative point (usually the C(α) atom), side chain solvent accessibility, and secondary structure. Structural comparison is effected by matching cliques of points. CLICK was extensively benchmarked for alignment accuracy on four different sets: (i) 9537 pair-wise alignments between two structures with the same topology; (ii) 64 alignments from set (i) that were considered to constitute difficult alignment cases; (iii) 199 pair-wise alignments between proteins with similar structure but different topology; and (iv) 1275 pair-wise alignments of RNA structures. The accuracy of CLICK alignments was measured by the average structure overlap score and compared with other alignment methods, including HOMSTRAD, MUSTANG, Geometric Hashing, SALIGN, DALI, GANGSTA(+), FATCAT, ARTS and SARA. On average, CLICK produces pair-wise alignments that are either comparable or statistically significantly more accurate than all of these other methods. We have used CLICK to uncover relationships between (previously) unrelated proteins. These new biological insights include: (i) detecting hinge regions in proteins where domain or sub-domains show flexibility; (ii) discovering similar small molecule binding sites from proteins of different folds and (iii) discovering topological variants of known structural/sequence motifs. Our method can generally be applied to compare any pair of molecular structures represented in Cartesian coordinates as exemplified by the RNA structure superimposition benchmark.     

Proteasomes: protein and gene structures.

Proteasomes are ring- or cylinder-shaped particles that have a sedimentation coefficient of 20S and are composed of a characteristic set of small polypeptides. These particles have a latent multicatalytic proteinase activity. Recently, proteasomes were found to combine reversibly with multiple protein components to form 26S proteolytic complexes that catalyze ATP-dependent, selective breakdown of proteins ligated with ubiquitin. This suggests that the 26S complexes are a new type of ATP-requiring protease in eukaryotic cells. We have studied the structures of various eukaryotic proteasomes at the molecular level by physicochemical and recombinant DNA techniques and have proposed that the gross structures of proteasomes, such as their size and shape, have been highly conserved during evolution. Proteasome subunits appear to be encoded by a family of homologous genes named the "proteasome gene family," which may have evolved from a common ancestral gene. Evidence obtained by genetic analyses in yeast and studies on the levels of proteasome expression in various eukaryotic cells indicates that proteasomes have essential roles in the cell. In this review, we summarize available information on the protein and gene structures of proteasomes and discuss the biological functions of proteasomes.     

Characterization of cysteine thiol modifications based on protein microenvironments and local secondary structures

We have demonstrated earlier that protein microenvironments were conserved around disulfide‐bridged cystine motifs with similar functions, irrespective of diversity in protein sequences. Here, cysteine thiol modifications were characterized based on protein microenvironments, secondary structures and specific protein functions. Protein microenvironment around an amino acid was defined as the summation of hydrophobic contributions from the surrounding protein fragments and the solvent molecules present within its first contact shell. Cysteine functions (modifications) were grouped into enzymatic and non‐enzymatic classes. Modifications studied were—disulfide formation, thio‐ether formation, metal‐binding, nitrosylation, acylation, selenylation, glutathionylation, sulfenylation, and ribosylation. 1079 enzymatic proteins were reported from high‐resolution crystal structures. Protein microenvironments around cysteine thiol, derived from above crystal structures, were clustered into 3 groups—buried‐hydrophobic, intermediate and exposed‐hydrophilic clusters. Characterization of cysteine functions were statistically meaningful for 4 modifications (disulfide formation, thioether formation, sulfenylation, and iron/zinc binding) those have sufficient amount of data in the current dataset. Results showed that protein microenvironment, secondary structure and protein functions were conserved for enzymatic cysteine functions, in contrast to the same function from non‐enzymatic cysteines. Disulfide forming enzymatic cysteines were tightly packed within intermediate protein microenvironment cluster, have alpha‐helical conformation and mostly belonged to CxxC motif of electron transport proteins. Disulfide forming non‐enzymatic cysteines did not belong to conserved motif and have variable secondary structures. Similarly, enzymatic thioether forming cysteines have conserved microenvironment compared to non‐enzymatic cystienes. Based on the compatibility between protein microenvironment and cysteine modifications, more efficient drug molecules could be designed against cysteine‐related diseases.     

Characterization of a hyaluronic acid-binding protein from sheep brain comparison with human brain hyaluronectin.

1. A hyaluronic acid (HA)-binding glycoprotein from sheep brain was characterized. 2. The specific affinity for HA was shown in vitro by high performance liquid chromatography, polyacrylamide gel electrophoresis and ELISA methods. 3. The KD for high molecular weight HA was 5.4 10(-9) M at 37 degrees C and lower than 10(-10) M at 4 degrees C. 4. No link protein was found and HA molecules could bind up to 10 times their weight of the glycoprotein. 5. The specific site for interaction was the HA-derived decasaccharide HA10. 6. The protein is composed of one polypeptidic chain. Tryptophan and lysine play a prominent role in the conformation of the binding site to HA. 7. Enzyme analysis indicated that the protein different forms are due to differences in glycosylation and that N- and O-linkages coexist in the molecules. 8. Immunohistochemistry localized the glycoprotein at the nodes of Ranvier and at the periphery of neurons. The perineuronal labeling was seen around all neurons studied in the cerebellum whereas it was almost undetectable in the cerebral hemispheres. 9. HA is not saturated by hyaluronectin (HN) in the sheep nervous system. 10. The glycoprotein is largely similar to human brain HN, and different from the hyaluronate-binding protein characterized in the cartilage.