A systematic search for protein signature sequences.

Signature sequences are contiguous patterns of amino acids 10-50 residues long that are associated with a particular structure or function in proteins. These may be of three types (by our nomenclature): superfamily signatures, remnant homologies, and motifs. We have performed a systematic search through a database of protein sequences to automatically and preferentially find remnant homologies and motifs. This was accomplished in three steps: 1. We generated a nonredundant sequence database. 2. We used BLAST3 (Altschul and Lipman, Proc. Natl. Acad. Sci. U.S.A. 87:5509-5513, 1990) to generate local pairwise and triplet sequence alignments for every protein in the database vs. every other. 3. We selected "interesting" alignments and grouped them into clusters. We find that most of the clusters contain segments from proteins which share a common structure or function. Many of them correspond to signatures previously noted in the literature. We discuss three previously recognized motifs in detail (FAD/NAD-binding, ATP/GTP-binding, and cytochrome b5-like domains) to demonstrate how the alignments generated by our procedure are consistent with previous work and make structural and functional sense. We also discuss two signatures (for N-acetyltransferases and glycerol-phosphate binding) which to our knowledge have not been previously recognized.     

Isolation of a cDNA clone encoding the RNase-superfamily-related gene highly expressed in chicken bone marrow cells.

The RNase gene superfamily combines functionally divergent proteins which share statistically significant sequence similarity. Known members assigned to this family include secretory and nonsecretory RNases; angiogenin; eosinophil cationic protein; eosinophil-derived neurotoxin; sialic-acid binding lectin and anti-tumor protein P-30. We report the cDNA cloning of the chicken RNase Super Family Related (RSFR) gene that is specifically overexpressed in normal bone marrow cells and bone marrow-derived AMV transformed monoblasts. It codes for a 139 amino acid protein with a putative signal peptide and remarkable conservation of active-site residues, other residues known to be important for substrate binding and catalytic activity and half-cystine residues common for all RNase family members. Phylogenetic tree analysis shows that RSFR defines a new group of genes within the family. We also conclude that an amino acid sequence block CKXXNTF(X) 11C is a "shortest RNase superfamily signature" which is both necessary and sufficient to identify all previously recognized family members as well as chicken RSFR.     

Prediction of amino acid sequence from structure

We have developed a method for the prediction of an amino acid sequence that is compatible with a three-dimensional backbone structure. Using only a backbone structure of a protein as input, the algorithm is capable of designing sequences that closely resemble natural members of the protein family to which the template structure belongs. In general, the predicted sequences are shown to have multiple sequence profile scores that are dramatically higher than those of random sequences, and sometimes better than some of the natural sequences that make up the superfamily. As anticipated, highly conserved but poorly predicted residues are often those that contribute to the functional rather than structural properties of the protein. Overall, our analysis suggests that statistical profile scores of designed sequences are a novel and valuable figure of merit for assessing and improving protein design algorithms.     

Subtilases: the superfamily of subtilisin-like serine proteases.

Subtilases are members of the clan (or superfamily) of subtilisin-like serine proteases. Over 200 subtilases are presently known, more than 170 of which with their complete amino acid sequence. In this update of our previous overview (Siezen RJ, de Vos WM, Leunissen JAM, Dijkstra BW, 1991, Protein Eng 4:719-731), details of more than 100 new subtilases discovered in the past five years are summarized, and amino acid sequences of their catalytic domains are compared in a multiple sequence alignment. Based on sequence homology, a subdivision into six families is proposed. Highly conserved residues of the catalytic domain are identified, as are large or unusual deletions and insertions. Predictions have been updated for Ca(2+)-binding sites, disulfide bonds, and substrate specificity, based on both sequence alignment and three-dimensional homology modeling.     

Combining structural and bioinformatics methods for the analysis of functionally important residues in DNA glycosylases

An essential function of DNA glycosylases is the recognition and excision of damaged bases in DNA, thereby preserving genomic integrity. Lesion recognition is a multistep process, which is only partially revealed by structural analysis of the catalytically competent complex. The functional role of additional residues can be predicted by combining structural data with analysis of amino acid conservation. The following postulate underlies this approach: if a family or superfamily can be broken into subgroups with different substrate specificities, residues highly conserved between these subgroups represent those important for enzyme catalysis and structure maintenance while residues highly conserved within a subgroup but not between subgroups represent residues important for substrate specificity. We review the bioinformatics approach used for this quantitative analysis and describe its application to the Nth superfamily and Fpg family of DNA glycosylases. These results serve as a starting point in planning site-directed mutagenesis experiments to elucidate the functional role of similar and dissimilar residues in DNA repair and other proteins.     

Combining protein sequence,structure, and dynamics: A novel approach for functional evolution analysis of PAS domain superfamily

PAS domains are widespread in archaea, bacteria, and eukaryota, and play important roles in various functions. In this study, we aim to explore functional evolutionary relationship among proteins in the PAS domain superfamily in view of the sequence‐structure‐dynamics‐function relationship. We collected protein sequences and crystal structure data from RCSB Protein Data Bank of the PAS domain superfamily belonging to three biological functions (nucleotide binding, photoreceptor activity, and transferase activity). Protein sequences were aligned and then used to select sequence‐conserved residues and build phylogenetic tree. Three‐dimensional structure alignment was also applied to obtain structure‐conserved residues. The protein dynamics were analyzed using elastic network model (ENM) and validated by molecular dynamics (MD) simulation. The result showed that the proteins with same function could be grouped by sequence similarity, and proteins in different functional groups displayed statistically significant difference in their vibrational patterns. Interestingly, in all three functional groups, conserved amino acid residues identified by sequence and structure conservation analysis generally have a lower fluctuation than other residues. In addition, the fluctuation of conserved residues in each biological function group was strongly correlated with the corresponding biological function. This research suggested a direct connection in which the protein sequences were related to various functions through structural dynamics. This is a new attempt to delineate functional evolution of proteins using the integrated information of sequence, structure, and dynamics.     

Integration of Evolutionary Features for the Identification of Functionally Important Residues in Major Facilitator Superfamily Transporters

The identification of functionally important residues is an important challenge for understanding the molecular mechanisms of proteins. Membrane protein transporters operate two-state allosteric conformational changes using functionally important cooperative residues that mediate long-range communication from the substrate binding site to the translocation pathway. In this study, we identified functionally important cooperative residues of membrane protein transporters by integrating sequence conservation and co-evolutionary information. A newly derived evolutionary feature, the co-evolutionary coupling number, was introduced to measure the connectivity of co-evolving residue pairs and was integrated with the sequence conservation score. We tested this method on three Major Facilitator Superfamily (MFS) transporters, LacY, GlpT, and EmrD. MFS transporters are an important family of membrane protein transporters, which utilize diverse substrates, catalyze different modes of transport using unique combinations of functional residues, and have enough characterized functional residues to validate the performance of our method. We found that the conserved cores of evolutionarily coupled residues are involved in specific substrate recognition and translocation of MFS transporters. Furthermore, a subset of the residues forms an interaction network connecting functional sites in the protein structure. We also confirmed that our method is effective on other membrane protein transporters. Our results provide insight into the location of functional residues important for the molecular mechanisms of membrane protein transporters.     

Crystal structure of an inactive duck delta II crystallin mutant with bound argininosuccinate

Delta-crystallin, the major soluble protein component of avian and reptilian eye lenses, is highly homologous to the urea cycle enzyme, argininosuccinate lyase (ASL). In duck lenses, there are two highly homologous delta crystallins, delta I and delta II, that are 94% identical in amino acid sequence. While delta II crystallin has been shown to exhibit ASL activity in vitro, delta I is enzymatically inactive. The X-ray structure of a His to Asn mutant of duck delta II crystallin (H162N) with bound argininosuccinate has been determined to 2.3 A resolution using the molecular replacement technique. The overall fold of the protein is similar to other members of the superfamily to which this protein belongs, with the active site located in a cleft formed by three different monomers in the tetramer. The active site of the H162N mutant structure reveals that the side chain of Glu 296 has a different orientation relative to the homologous residue in the H91N mutant structure [Abu-Abed et al. (1997) Biochemistry 36, 14012-14022]. This shift results in the loss of the hydrogen bond between His 162 and Glu 296 seen in the H91N and turkey delta I crystallin structures; this H-bond is believed to be crucial for the catalytic mechanism of ASL/delta II crystallin. Argininosuccinate was found to be bound to residues in each of the three monomers that form the active site. The fumarate moiety is oriented toward active site residues His 162 and Glu 296 and other residues that are part of two of the three highly conserved regions of amino acid sequence in the superfamily, while the arginine moiety of the substrate is oriented toward residues which belong to either domain 1 or domain 2. The analysis of the structure reveals that significant conformational changes occur on substrate binding. The comparison of this structure with the inactive turkey delta I crystallin reveals that the conformation of domain 1 is crucial for substrate affinity and that the delta I protein is almost certainly inactive because it can no longer bind the substrate.     

A Combinatorial Approach to Detect Coevolved Amino Acid Networks in Protein Families of Variable Divergence

Communication between distant sites often defines the biological role of a protein: amino acid long-range interactions are as important in binding specificity, allosteric regulation and conformational change as residues directly contacting the substrate. The maintaining of functional and structural coupling of long-range interacting residues requires coevolution of these residues. Networks of interaction between coevolved residues can be reconstructed, and from the networks, one can possibly derive insights into functional mechanisms for the protein family. We propose a combinatorial method for mapping conserved networks of amino acid interactions in a protein which is based on the analysis of a set of aligned sequences, the associated distance tree and the combinatorics of its subtrees. The degree of coevolution of all pairs of coevolved residues is identified numerically, and networks are reconstructed with a dedicated clustering algorithm. The method drops the constraints on high sequence divergence limiting the range of applicability of the statistical approaches previously proposed. We apply the method to four protein families where we show an accurate detection of functional networks and the possibility to treat sets of protein sequences of variable divergence.     

The 1.3-Angstrom-resolution crystal structure of beta-ketoacyl-acyl carrier protein synthase II from Streptococcus pneumoniae

The beta-ketoacyl-acyl carrier protein synthases are members of the thiolase superfamily and are key regulators of bacterial fatty acid synthesis. As essential components of the bacterial lipid metabolic pathway, they are an attractive target for antibacterial drug discovery. We have determined the 1.3 A resolution crystal structure of the beta-ketoacyl-acyl carrier protein synthase II (FabF) from the pathogenic organism Streptococcus pneumoniae. The protein adopts a duplicated betaalphabetaalphabetaalphabetabeta fold, which is characteristic of the thiolase superfamily. The two-fold pseudosymmetry is broken by the presence of distinct insertions in the two halves of the protein. These insertions have evolved to bind the specific substrates of this particular member of the thiolase superfamily. Docking of the pantetheine moiety of the substrate identifies the loop regions involved in substrate binding and indicates roles for specific, conserved residues in the substrate binding tunnel. The active site triad of this superfamily is present in spFabF as His 303, His 337, and Cys 164. Near the active site is an ion pair, Glu 346 and Lys 332, that is conserved in the condensing enzymes but is unusual in our structure in being stabilized by an Mg(2+) ion which interacts with Glu 346. The active site histidines interact asymmetrically with Lys 332, whose positive charge is closer to His 303, and we propose a specific role for the lysine in polarizing the imidazole ring of this histidine. This asymmetry suggests that the two histidines have unequal roles in catalysis and provides new insights into the catalytic mechanisms of these enzymes.     

Motif refinement of the peroxisomal targeting signal 1 and evaluation of taxon-specific differences

Eukaryote peroxisomes, plant glyoxysomes and trypanosomal glycosomes belong to the microbody family of organelles that compartmentalise a variety of biochemical processes. The interaction between the PTS1 signal and its cognate receptor Pex5 initiates the major import mechanism for proteins into the matrix of these organelles. Relying on the analysis of amino acid sequence variability of known PTS1-targeted proteins and PTS1-containing peptides that interact with Pex5 in the yeast two-hybrid assay, on binding site studies of the Pex5-ligand complex crystal structure, 3D models and sequences of Pex5 proteins from various taxa, we derived the requirements for a C-terminal amino acid sequence to interact productively with Pex5. We found evidence that, at least the 12 C-terminal residues of a given substrate protein are implicated in PTS1 signal recognition. This motif can be structurally and functionally divided into three regions: (i) the C-terminal tripeptide, (ii) a region interacting with the surface of Pex5 (about four residues further upstream), and (iii) a polar, solvent-accessible and unstructured region with linker function (the remaining five residues). Specificity differences are confined to taxonomic subgroups (metazoa and fungi) and are connected with amino acid type preferences in region 1 and deviating hydrophobicity patterns in region 2.     

Molecular cloning of a novel ras-like protein from chicken.

We have isolated a cDNA clone from a chicken DNA expression library which codes for a ras-like polypeptide of 216 amino acid residues. This polypeptide is closely related to the human protein TC4 and to the yeast protein Spil, two novel proteins that may be involved in the coordination of the cell cycle. In the amino-terminal region, the three polypeptides possess a P-loop motif characteristic of GTP-binding proteins. At the carboxy-terminal end, however, they lack the typical CAAX-box which is usually responsible for membrane anchorage of ras-like proteins. It is therefore likely that the three polypeptides define a new subclass of GTP-binding proteins within the ras-like superfamily.     

The LysE superfamily: topology of the lysine exporter LysE of Corynebacterium glutamicum, a paradyme for a novel superfamily of transmembrane solute translocators

In Corynebacterium glutamicum the LysE carrier protein exhibits the unique function of exporting L-lysine. We here analyze the membrane topology of LysE, a protein of 236 amino acyl residues, using PhoA- and LacZ-fusions. The amino-terminal end of LysE is located in the cytoplasm whereas the carboxy-terminal end is found in the periplasm. Although 6 hydrophobic domains were identified based on hydropathy analyses, only five transmembrane spanning helices appear to be present. The additional hydrophobic segment may dip into the membrane or be surface localized. We show that LysE is a member of a family of proteins found, for example, in Escherichia coil, Bacillus subtilis, Mycobacterium tuberculosis and Helicobacter pylori. This family, which we have designated the LysE family, is distantly related to two additional protein families which we have designated the YahN and CadD families. These three families, the members of which exhibit similar sizes, hydropathy profiles, and sequence motifs comprise the LysE superfamily. Functionally characterized members of the LysE superfamily export L-lysine, cadmium and possibly quarternary amines. We suggest that LysE superfamily members will prove to catalyze export of a variety of biologically important solutes.     

Estimation of amino acid residue substitution rates at local spatial regions and application in protein function inference: a Bayesian Monte Carlo approach

The amino acid sequences of proteins provide rich information for inferring distant phylogenetic relationships and for predicting protein functions. Estimating the rate matrix of residue substitutions from amino acid sequences is also important because the rate matrix can be used to develop scoring matrices for sequence alignment. Here we use a continuous time Markov process to model the substitution rates of residues and develop a Bayesian Markov chain Monte Carlo method for rate estimation. We validate our method using simulated artificial protein sequences. Because different local regions such as binding surfaces and the protein interior core experience different selection pressures due to functional or stability constraints, we use our method to estimate the substitution rates of local regions. Our results show that the substitution rates are very different for residues in the buried core and residues on the solvent-exposed surfaces. In addition, the rest of the proteins on the binding surfaces also have very different substitution rates from residues. Based on these findings, we further develop a method for protein function prediction by surface matching using scoring matrices derived from estimated substitution rates for residues located on the binding surfaces. We show with examples that our method is effective in identifying functionally related proteins that have overall low sequence identity, a task known to be very challenging.     

Docking and homology modeling explain inhibition of the human vesicular glutamate transporters

As membrane transporter proteins, VGLUT1-3 mediate the uptake of glutamate into synaptic vesicles at presynaptic nerve terminals of excitatory neural cells. This function is crucial for exocytosis and the role of glutamate as the major excitatory neurotransmitter in the central nervous system. The three transporters, sharing 76% amino acid sequence identity in humans, are highly homologous but differ in regional expression in the brain. Although little is known regarding their three-dimensional structures, hydropathy analysis on these proteins predicts 12 transmembrane segments connected by loops, a topology similar to other members in the major facilitator superfamily, where VGLUT1-3 have been phylogenetically classified. In this work, we present a three-dimensional model for the human VGLUT1 protein based on its distant bacterial homolog in the same superfamily, the glycerol-3-phosphate transporter from Escherichia coli. This structural model, stable during molecular dynamics simulations in phospholipid bilayers solvated by water, reveals amino acid residues that face its pore and are likely to affect substrate translocation. Docking of VGLUT1 substrates to this pore localizes two different binding sites, to which inhibitors also bind with an overall trend in binding affinity that is in agreement with previously published experimental data.     

Functional analysis of the glycero-manno-heptose 7-phosphate kinase domain from the bifunctional HldE protein, which is involved in ADP-L-glycero-D-manno-heptose biosynthesis

The core oligosaccharide component of the lipopolysaccharide can be subdivided into inner and outer core regions. In Escherichia coli, the inner core consists of two 3-deoxy-d-manno-octulosonic acid and three glycero-manno-heptose residues. The HldE protein participates in the biosynthesis of ADP-glycero-manno-heptose precursors used in the assembly of the inner core. HldE comprises two functional domains: an N-terminal region with homology to the ribokinase superfamily (HldE1 domain) and a C-terminal region with homology to the cytidylyltransferase superfamily (HldE2 domain). We have employed the structure of the E. coli ribokinase as a template to model the HldE1 domain and predict critical amino acids required for enzyme activity. Mutation of these residues renders the protein inactive as determined in vivo by functional complementation analysis. However, these mutations did not affect the secondary or tertiary structure of purified HldE1, as judged by fluorescence spectroscopy and circular dichroism. Furthermore, in vivo coexpression of wild-type, chromosomally encoded HldE and mutant HldE1 proteins with amino acid substitutions in the predicted ATP binding site caused a dominant negative phenotype as revealed by increased bacterial sensitivity to novobiocin. Copurification experiments demonstrated that HldE and HldE1 form a complex in vivo. Gel filtration chromatography resulted in the detection of a dimer as the predominant form of the native HldE1 protein. Altogether, our data support the notions that the HldE functional unit is a dimer and that structural components present in each HldE1 monomer are required for enzymatic activity.     

Crystal structures and proposed structural/functional classification of three protozoan proteins from the isochorismatase superfamily

We have determined the crystal structures of three homologous proteins from the pathogenic protozoans Leishmania donovani, Leishmania major, and Trypanosoma cruzi. We propose that these proteins represent a new subfamily within the isochorismatase superfamily (CDD classification cd004310). Their overall fold and key active site residues are structurally homologous both to the biochemically well-characterized N-carbamoylsarcosine-amidohydrolase, a cysteine hydrolase, and to the phenazine biosynthesis protein PHZD (isochorismase), an aspartyl hydrolase. All three proteins are annotated as mitochondrial-associated ribonuclease Mar1, based on a previous characterization of the homologous protein from L. tarentolae. This would constitute a new enzymatic activity for this structural superfamily, but this is not strongly supported by the observed structures. In these protozoan proteins, the extended active site is formed by inter-subunit association within a tetramer, which implies a distinct evolutionary history and substrate specificity from the previously characterized members of the isochorismatase superfamily. The characterization of the active site is supported crystallographically by the presence of an unidentified ligand bound at the active site cysteine of the T. cruzi structure.     

Inferring functional relationships of proteins from local sequence and spatial surface patterns

We describe a novel approach for inferring functional relationship of proteins by detecting sequence and spatial patterns of protein surfaces. Well-formed concave surface regions in the form of pockets and voids are examined to identify similarity relationship that might be directly related to protein function. We first exhaustively identify and measure analytically all 910,379 surface pockets and interior voids on 12,177 protein structures from the Protein Data Bank. The similarity of patterns of residues forming pockets and voids are then assessed in sequence, in spatial arrangement, and in orientational arrangement. Statistical significance in the form of E and p-values is then estimated for each of the three types of similarity measurements. Our method is fully automated without human intervention and can be used without input of query patterns. It does not assume any prior knowledge of functional residues of a protein, and can detect similarity based on surface patterns small and large. It also tolerates, to some extent, conformational flexibility of functional sites. We show with examples that this method can detect functional relationship with specificity for members of the same protein family and superfamily, as well as remotely related functional surfaces from proteins of different fold structures. We envision that this method can be used for discovering novel functional relationship of protein surfaces, for functional annotation of protein structures with unknown biological roles, and for further inquiries on evolutionary origins of structural elements important for protein function.     

Cleavage-site specificity of prolyl endopeptidase FAP investigated with a full-length protein substrate

Fibroblast activation protein (FAP) is a prolyl-cleaving endopeptidase proposed as an anti-cancer drug target. It is necessary to define its cleavage-site specificity to facilitate the identification of its in vivo substrates and to understand its biological functions. We found that the previously identified substrate of FAP, α(2)-anti-plasmin, is not a robust substrate in vitro. Instead, an intracellular protein, SPRY2, is cleavable by FAP and more suitable for investigation of its substrate specificity in the context of the full-length globular protein. FAP prefers uncharged residues, including small or bulky hydrophobic amino acids, but not charged amino acids, especially acidic residue at P1', P3 and P4 sites. Molecular modelling analysis shows that the substrate-binding site of FAP is surrounded by multiple tyrosine residues and some negatively charged residues, which may exert least preference for substrates with acidic residues. This provides an explanation why FAP cannot cleave interleukins, which have a glutamate at either P4 or P2', despite their P3-P2-P1 sites being identical to SPRY2 or α-AP. Our study provided new information on FAP cleavage-site specificity, which differs from the data obtained by profiling with a peptide library or with the denatured protein, gelatin, as the substrate. Furthermore, our study suggests that negatively charged residues should be avoided when designing FAP inhibitors.     

Structural insights into the lipase/esterase behavior in the Candida rugosa lipases family: crystal structure of the lipase 2 isoenzyme at 1.97A resolution

The yeast Candida rugosa produces several closely related extracellular lipases that differ in their substrate specificity. Here, we report the crystal structure of the isoenzyme lipase 2 at 1.97A resolution in its closed conformation. Lipase 2 shows a 79.4% amino acid sequence identity with lipase 1 and 82.2% with lipase 3, which makes it relevant to compare these three isoenzymes. Despite this high level of sequence identity, structural comparisons reveal several amino acid changes affecting the flap (residue 69), the substrate-binding pocket (residues 127, 132 and 450) and the mouth of the hydrophobic tunnel (residues 296 and 344), which may be responsible for the different substrate specificity and catalytic properties of this group of enzymes. Also, these comparisons reveal two distinct regions in the hydrophobic tunnel: a phenylalanyl-rich region and an aliphatic-rich region. Whereas this last region is essentially identical in the three isoenzymes, the phenylalanyl content in the first one is specific for each lipase, resulting in a different environment of the catalytic triad residues, which probably tunes finely their lipase/esterase character. The greater structural similarity observed between the monomeric form of lipase 3 and lipase 2 concerning the above-mentioned key residues led us to propose a significant esterase activity for this last protein. This enzymatic activity has been confirmed with biochemical experiments using cholesteryl [1-14C]oleate as substrate. Surprisingly, lipase 2 is a more efficient esterase than lipase 3, showing a twofold specific activity against cholesteryl [1-14C]oleate in our experimental conditions. These results show that subtle amino acid changes within a highly conserved protein fold may produce protein variants endowed with new enzymatic properties.