Dimensions and specificities of recognition sites on lectins and antibodies

A comparison is made of the specific combining sites of a number of lectins and of antibodies with emphasis on those reacting with blood group A, B, and H determinants. The ranges of site sizes and specificities of both groups are similar both from immunochemical studies and from the limited x-ray diffraction data available.     

Analysis and prediction of carbohydrate binding sites

An analysis of the characteristic properties of sugar binding sites was performed on a set of 19 sugar binding proteins. For each site six parameters were evaluated: solvation potential, residue propensity, hydrophobicity, planarity, protrusion and relative accessible surface area. Three of the parameters were found to distinguish the observed sugar binding sites from the other surface patches. These parameters were then used to calculate the probability for a surface patch to be a carbohydrate binding site. The prediction was optimized on a set of 19 non-homologous carbohydrate binding structures and a test prediction was carried out on a set of 40 protein-carbohydrate complexes. The overall accuracy of prediction achieved was 65%. Results were in general better for carbohydrate-binding enzymes than for the lectins, with a rate of success of 87%.     

Rice chitinases: sugar recognition specificities of the individual subsites

Sugar recognition specificities of class III (OsChib1a) and class I (OsChia1cDeltaChBD) chitinases from rice, Oryza sativa L., were investigated by analyzing (1)H- and (13)C-nuclear magnetic resonance spectra of the enzymatic products from partially N-acetylated chitosans. The reducing end residue of the enzymatic products obtained by the class III enzyme was found to be exclusively acetylated, whereas both acetylated and deacetylated units were found at the nearest neighbor to the reducing end residue. Both acetylated and deacetylated units were also found at the nonreducing end residue and its nearest neighbor of the class III enzyme products. Thus, only subsite (-1) among the contiguous subsites (-2) to (+2) of the class III enzyme was found to be specific to an acetylated residue. For the class I enzyme, the reducing end residue was preferentially acetylated, although the specificity was not absolute. The nearest neighbor to the acetylated reducing end residue was specifically acetylated. Moreover, the nonreducing end residue produced by the class I enzyme was exclusively acetylated, although there was a low but significant preference for deacetylated units at the nearest neighbor to the nonreducing end. These results suggest that the three contiguous subsites (-2), (-1), and (+1) of the class I enzyme are specific to three consecutive GlcNAc residues of the substrate. In rice plants, the target of the class I enzyme might be a consecutive GlcNAc sequence probably in the cell wall of fungal pathogen, whereas the class III enzyme might act toward an endogenous complex carbohydrate containing GlcNAc residue.     

A branch point consensus from Arabidopsis found by non-circular analysis allows for better prediction of acceptor sites.

Little knowledge exists about branch points in plants; it has even been claimed that plant introns lack conserved branch point sequences similar to those found in vertebrate introns. A putative branch point consensus sequence for Arabidopsis thaliana resembling the well known metazoan consensus sequence has been proposed, but this is based on search of sequences similar to those in yeast and metazoa. Here we present a novel consensus sequence found by a non-circular approach. A hidden Markov model with a fixed A nucleotide was trained on sequences upstream of the acceptor site. The consensus found by the Markov model shares features with the metazoan consensus, but differs in its details from the consensus proposed earlier. Despite the fact that branch point consensus sequences in plants are weak, we show that a prediction scheme incorporating them leads to a substantial improvement in the recognition of true acceptor sites; the false positive rate being reduced by a factor of 2. We take this as an indication that the consensus found here is the genuine one and that the branch point does play a role in the proper recognition of the acceptor site in plants.     

Analysis and prediction of functionally important sites in proteins

The rapidly increasing volume of sequence and structure information available for proteins poses the daunting task of determining their functional importance. Computational methods can prove to be very useful in understanding and characterizing the biochemical and evolutionary information contained in this wealth of data, particularly at functionally important sites. Therefore, we perform a detailed survey of compositional and evolutionary constraints at the molecular and biological function level for a large set of known functionally important sites extracted from a wide range of protein families. We compare the degree of conservation across different functional categories and provide detailed statistical insight to decipher the varying evolutionary constraints at functionally important sites. The compositional and evolutionary information at functionally important sites has been compiled into a library of functional templates. We developed a module that predicts functionally important columns (FIC) of an alignment based on the detection of a significant "template match score" to a library template. Our template match score measures an alignment column's similarity to a library template and combines a term explicitly representing a column's residue composition with various evolutionary conservation scores (information content and position-specific scoring matrix-derived statistics). Our benchmarking studies show good sensitivity/specificity for the prediction of functional sites and high accuracy in attributing correct molecular function type to the predicted sites. This prediction method is based on information derived from homologous sequences and no structural information is required. Therefore, this method could be extremely useful for large-scale functional annotation.     

Reassortment of DNA recognition domains and the evolution of new specificities

Type I restriction enzymes comprise three subunits only one of which, the S polypeptide, dictates the specificity of the DNA sequence recognized. Recombination between two different hsdS genes, SP and SB, led to the isolation of a system, SQ, which had a different specificity from that of either parent. The finding that the nucleotide sequence recognized by SQ is a hybrid containing components from both the SP and SB target sequences suggested that DNA recognition is carried out by two separable domains within each specificity polypeptide. To test this we have made the recombinant gene of reciprocal structure and demonstrate that it encodes a polypeptide whose recognition sequence, deduced in vivo, is as predicted by this model. We also report the sequence of the SB specificity gene, so that information is now available for the five known members of this family of enzymes. All show a similar organization of conserved and variable regions. Comparisons of the predicted amino acid sequences reveal large non-conserved areas which may not even be structurally similar. This is remarkable since these different S subunits are functionally identical, except for the specificity with respect to the DNA sequence with which they interact. We discuss the correlation of the variation in polypeptide sequence with recognition specificities.     

Mutational analysis of the engrailed homeodomain recognition helix by phage display.

The homeodomain (HD) is a ubiquitous protein fold that confers DNA binding function on a superfamily of eukaryotic gene regulatory proteins. Here, the DNA binding of recognition helix variants of the HD from the engrailed gene of Drosophila melanogaster was investigated by phage display. Nineteen different combinations of pairwise mutations at positions 50 and 54 were screened against a panel of four DNA sequences consisting of the engrailed consensus, a non-specific DNA control based on the lambda repressor operator OR1 and two model sequence targets con-taining imperfect versions of the 5'-TAAT-3' consensus. The resulting mutant proteins could be divided into four groups that varied with respect to their affinity for DNA and specificity for the engrailed consensus. The altered specificity phenotypes of several mutant proteins were confirmed by DNA mobility shift analysis. Lys50/Ala54 was the only mutant protein that exhibited preferential binding to a sequence other than the engrailed consensus. Arginine was also demonstrated to be a functional replacement for Ala54. The functional combinations at 50 and 54 identified by these experiments recapitulate the distribution of naturally occurring HD sequences and illustrate how the engrailed HD can be used as a framework to explore covariation among DNA binding residues.     

Reassortment of DNA recognition domains and the evolution of new specificities

Type I restriction enzymes comprise three subunits only one of which, the S polypeptide, dictates the specificity of the DNA sequence recognized. Recombination between two different hsdS genes, SP and SB, led to the isolation of a system, SQ, which had a different specificity from that of either parent. The finding that the nucleotide sequence recognized by SQ is a hybrid containing components from both the SP and SB target sequences suggested that DNA recognition is carried out by two separable domains within each specificity polypeptide. To test this we have made the recombinant gene of reciprocal structure and demonstrate that it encodes a polypeptide whose recognition sequence, deduced In vivo, is as predicted by this model. We also report the sequence of the SB specificity gene, so that information is now available for the five known members of this family of enzymes. Ali show a similar organization of conserved and variable regions. Comparisons of the predicted amino acid sequences reveal large non-conserved areas which may not even be structurally similar. This is remarkable since these different S subunits are functionally identical, except for the specificity with respect to the DNA sequence with which they interact. We discuss the correlation of the variation in polypeptide sequence with recognition specificities.     

Amino acid sequence specificities of an adhesive recognition signal

Synthetic peptides derived from the cell-binding domain of fibronectin have previously been found to inhibit fibronectin-mediated adhesion in vitro competitively and reversibly, as well as inhibiting cell migratory events in vivo. The amino acid sequence specificity required for this inhibitory activity has been examined further using variations of the originally identified active peptide sequences. The most active small peptide was found to be the pentapeptide Gly-Arg-Gly-Asp-Ser. Although the tetrapeptide Arg-Gly-Asp-Ser was found to retain substantial activity, it was approximately threefold less active. An "inverted" peptide sequence with these same four amino acids arranged in the mirror symmetrical sequence Ser-Asp-Gly-Arg was found to be nearly as active as the forward sequence. However, the same inverted tetrapeptide sequence embedded in a synthetic decapeptide derived from a sequence of histocompatibility antigens has minimal activity, suggesting the importance of adjacent sequences in modifying the activity of such peptides. Neither substitution of amino acids of the same charge nor reversal of the positions of the two charged amino acids retains biological activity. Decreasing the spacing between the charged residues also causes a loss of activity. Our results suggest the hypothesis that this adhesive recognition signal consists of a specific arrangement of one acidic and one basic charged group and additional information provided by adjacent amino acids.     

Using the recognition code to swap homeodomain target specificity in cell culture

The homeodomain (HD) is a 60 amino acid-long DNA-binding domain. A large fraction of HDs binds with high affinity sequences containing the 5′-TAAT-3′ core motif. However, NK-2 class HDs recognizes sequences containing the 5′-CAAG-3′ core motif. By using a cell transfection approach, here we show that modification of residues located in the N-terminal arm (at positions 6, 7 and 8) and in the recognition helix (at position 54) is enough to swap the “in vivo” binding specificity of TTF-1 HD (which is a member of the NK-2 class HD) from 5′-CAAG-3′ to 5′-TAAT-3′-containing targets. The role of residue at position 54 is also supported by data obtained with the HD of the Drosophila engrailed protein. These data support the notion that DNA-binding specificity “in vivo” is dictated by few critical residues.     

Identification of preferred actinomycin-DNA binding sites by the combinatorial method REPSA.

An important question in the study of ligand-DNA interactions is the determination of binding specificity. Here, we used the combinatorial method restriction endonuclease protection, selection, and amplification (REPSA) to identify the preferred duplex DNA-binding sites of the antineoplastic agent actinomycin D. After 10 rounds of REPSA, over 95% of the cloned DNAs exhibited significantly reduced FokI restriction endonuclease cleavage in the presence of 1 microM actinomycin. A chi(2) statistical analysis of their sequences found that 39 of the 45 clones contained one or more copies of the sequence 5'-(T/A)GC(A/T)-3', giving a p<0.001 for this consensus. A DNase I footprinting analysis of the cloned DNAs found that all possessed relatively high affinity actinomycin-binding sites with apparent dissociation constants between 12 and 258nM (average 98nM). The average footprint encompassed 7.6 bases and in most cases (90%) included one or more consensus sequences. Interestingly, several of the selected clones contained overlapping consensus sequences (e.g., 5'-TGCTGCT-3'), suggesting that such close proximity DNA-binding sites may actually be preferred by actinomycin under physiological conditions.