Consensus sequence for processing of peptide precursors at monobasic sites

Many regulatory peptide precursors undergo post-translational processing at mono- and/or dibasic residues. Comparison of amino acids around the monobasic cleavage sites suggests that these cleavages follow certain sequence motifs and can be described as the rules that govern monobasic cleavages: (i) a basic amino acid it present at either 3, 5, or 7 amino acids N-terminal to the cleavage site, (ii) hydrophobic aliphatic amino acids (leucine, isoleucine, valine, or methionine) are never present in the position C-terminal to the monobasic amino acid at the cleavage site, (iii) a cysteine is never present in the vicinity of the cleavage site, and (iv) an aromatic amino acid is never present at the position N-terminal to the monobasic amino acid at the cleavage site. In addition to these rules, the monobasic cleavages follow certain tendencies: (i) the amino acid at the cleavage site tends to be predominantly arginine, (ii) the amino acid at the position C-terminal to the cleavage site tends to be serine, alanine or glycine in more than 60% of the cases, (iii) the amino acid at either 3, 5, or 7 position N-terminal to the cleavage site tends to be arginine, (iv) aromatic amino acids are rare at the position C-terminal to the monobasic amino acid at the cleavage site, and (v) aliphatic amino acids tend to be in the two positions N-terminal to and the two positions C-terminal to the cleavage site, except as noted above. When compared with a large number of sequence containing single basic amino acids, these rules and tendencies are capable of not only correctly predicting the processing sites, but also are capable of excluding most of the single basic sequences that are known to be uncleaved. Many or these rules can also be applied to correctly predict the dibasic and multibasic cleavage sites suggesting that the rules and tendencies could govern endoproteolytic processing at the monobasic, dibasic and multibasic sites.     

Detecting coevolution in and among protein domains

Correlated changes of nucleic or amino acids have provided strong information about the structures and interactions of molecules. Despite the rich literature in coevolutionary sequence analysis, previous methods often have to trade off between generality, simplicity, phylogenetic information, and specific knowledge about interactions. Furthermore, despite the evidence of coevolution in selected protein families, a comprehensive screening of coevolution among all protein domains is still lacking. We propose an augmented continuous-time Markov process model for sequence coevolution. The model can handle different types of interactions, incorporate phylogenetic information and sequence substitution, has only one extra free parameter, and requires no knowledge about interaction rules. We employ this model to large-scale screenings on the entire protein domain database (Pfam). Strikingly, with 0.1 trillion tests executed, the majority of the inferred coevolving protein domains are functionally related, and the coevolving amino acid residues are spatially coupled. Moreover, many of the coevolving positions are located at functionally important sites of proteins/protein complexes, such as the subunit linkers of superoxide dismutase, the tRNA binding sites of ribosomes, the DNA binding region of RNA polymerase, and the active and ligand binding sites of various enzymes. The results suggest sequence coevolution manifests structural and functional constraints of proteins. The intricate relations between sequence coevolution and various selective constraints are worth pursuing at a deeper level.     

Protein-DNA interactions: amino acid conservation and the effects of mutations on binding specificity

We investigate the conservation of amino acid residue sequences in 21 DNA-binding protein families and study the effects that mutations have on DNA-sequence recognition. The observations are best understood by assigning each protein family to one of three classes: (i) non-specific, where binding is independent of DNA sequence; (ii) highly specific, where binding is specific and all members of the family target the same DNA sequence; and (iii) multi-specific, where binding is also specific, but individual family members target different DNA sequences. Overall, protein residues in contact with the DNA are better conserved than the rest of the protein surface, but there is a complex underlying trend of conservation for individual residue positions. Amino acid residues that interact with the DNA backbone are well conserved across all protein families and provide a core of stabilising contacts for homologous protein-DNA complexes. In contrast, amino acid residues that interact with DNA bases have variable levels of conservation depending on the family classification. In non-specific families, base-contacting residues are well conserved and interactions are always found in the minor groove where there is little discrimination between base types. In highly specific families, base-contacting residues are highly conserved and allow member proteins to recognise the same target sequence. In multi-specific families, base-contacting residues undergo frequent mutations and enable different proteins to recognise distinct target sequences. Finally, we report that interactions with bases in the target sequence often follow (though not always) a universal code of amino acid-base recognition and the effects of amino acid mutations can be most easily understood for these interactions.     

DNA binding specificity of mutant glucocorticoid receptor DNA-binding domains

Mutation of a small number of amino acids in the DNA-binding domain of the estrogen receptor to the corresponding sequence of the glucocorticoid receptor switches the specificity of the receptor in transactivation assays (Mader, S., Kumar, V., de Verneuil, H., and Chambon, P. (1989) Nature 338, 271-274). We have made the corresponding reciprocal mutations in the context of the glucocorticoid receptor DNA-binding domain and studied the binding of wild type and mutant purified proteins to palindromic glucocorticoid and estrogen response elements as well as to elements of intermediate sequence, using gel mobility shift assays. We show here that a protein with two altered amino acids binds glucocorticoid and estrogen response elements with a low but equal affinity, whereas a protein with an additional changed residue has a high affinity for estrogen response elements but still retains a considerable affinity for glucocorticoid response elements. Using binding sites of intermediate sequence we have further characterized the interaction with DNA. The in vitro DNA binding results are confirmed by in vivo transactivation assays in yeast. Finally we suggest a testable model for amino acid/base pair interactions involved in recognition by the glucocorticoid receptor DNA-binding domain of its target sequence.     

Identification of low density lipoprotein receptor binding domains of human apolipoprotein B-100: a proposed consensus LDL receptor binding sequence of apoB-100

The human liver apoB-100 gene cloned in the lambda gt-11 expression vector expresses fusion proteins reacting with apoB antibodies. A fusion protein induced from a apoB-lambda gt-11 clone reacted with apoB-100 monoclonal antibodies known to block the binding of LDL to the LDL receptor. The fusion protein contains an amino acid sequence domain enriched in positively charged residues which is complementary to the negatively charged amino acids present in the consensus LDL receptor binding domain. This sequence of apoB-100 is proposed as a binding domain for the interaction with the LDL receptor. Comparison of derived amino acid sequences from the entire structure of apoB-100 molecule revealed several similar domains enriched in positively charged amino acids. A consensus sequence of the potential LDL binding domain was identified which contained positively charged amino acids at positions 1, 5 and 8 and a loop of 8-11 amino acids followed by two adjacent positively charged amino acids. These results are interpreted as indicating that there are several potential LDL receptor binding domains in apoB-100.     

Position‐coded multivalent peptide–peptide interactions revealed by tryptophan‐scanning mutagenesis

We demonstrate in this contribution the evidence that significant cooperative binding effect can be identified for the amino acid sites that are determinant to the binding characteristics in peptide–peptide interactions. The analysis of tryptophan‐scanning mutagenesis of the 14‐mer peptide consisting only of glycine provides a mapping of position‐dependent contributions to the binding energy. The pronounced tryptophan‐associated peptide–peptide interactions are originated from the indole moieties with the main chains of 14‐mer glycines containing N–H and C?O moieties. Specifically, with the presence of two tryptophans as determinant amino acids, cooperative binding can be observed, which are dependent on relative positions of the two tryptophans with a “volcano”‐like characteristics. An optimal separation of 6–10 amino acids between two adjacent binding sites can be identified to achieve maximal binding interactions.     

Mutational analysis of a blood coagulation factor VIII-binding peptide.

A peptide screened from a combinatorial peptide library with the sequence EYKSWEYC performed best as a ligand for affinity chromatography of human blood coagulation factor VIII (FVIII). With this peptide immobilized on monolithic CIM columns via epoxy groups we were able to capture FVIII from diluted plasma. Rational substitution of amino acids by spot synthesis revealed that lysine and cysteine can be exchanged for almost all other proteinogenic amino acids without loss of affinity to FVIII. This offers the possibility of site-specific attachment via either one of these residues or the N- or C-terminus. The aliphatic positions O5 (tryptophan) and O7 (tyrosine), together with the charged position O6 (glutamic acid), seem to form the core of the binding unit. In the positions with aliphatic amino acids, substitution by tyrosine or phenylalanine, and in the positions with charged amino acids, substitution by aspartic acid or lysine, preserved the affinity to FVIII. The functionality of the selected peptides was confirmed by affinity chromatography. Selective binding and elution could be achieved.     

Amino acid propensities for the collagen triple-helix

Determination of the tendencies of amino acids to form alpha-helical and beta-sheet structures has been important in clarifying stabilizing interactions, protein design, and the protein folding problem. In this study, we have determined for the first time a complete scale of amino acid propensities for another important protein motif: the collagen triple-helix conformation with its Gly-X-Y repeating sequence. Guest triplets of the form Gly-X-Hyp and Gly-Pro-Y are used to quantitate the conformational propensities of all 20 amino acids for the X and Y positions in the context of a (Gly-Pro-Hyp)(8) host peptide. The rankings for both the X and Y positions show the highly stabilizing nature of imino acids and the destabilizing effects of Gly and aromatic residues. Many residues show differing propensities in the X versus Y position, related to the nonequivalence of these positions in terms of interchain interactions and solvent exposure. The propensity of amino acids to adopt a polyproline II-like conformation plays a role in their triple-helix rankings, as shown by a moderate correlation of triple-helix propensity with frequency of occurrence in polyproline II-like regions. The high propensity of ionizable residues in the X position suggests the importance of interchain hydrogen bonding directly or through water to backbone carbonyls or hydroxyprolines. The low propensity of side chains with branching at the C(delta) in the Y position supports models suggesting these groups block solvent access to backbone C=O groups. These data provide a first step in defining sequence-dependent variations in local triple-helix stability and binding, and are important for a general understanding of side chain interactions in all proteins.     

Yeast Puf3 mutants reveal the complexity of Puf-RNA binding and identify a loop required for regulation of mRNA decay

The eukaryotic Puf proteins regulate mRNA translation and degradation by binding the 3' untranslated regions of target mRNAs. Crystal structure analysis of a human Puf bound to RNA suggested a modular mode of binding, with specific amino acids within each of eight repeat domains contacting a single nucleotide of the target RNA. Here we study the mechanism by which the yeast Puf3p binds and stimulates the degradation of COX17 mRNA. Mutation of the predicted RNA-binding positions of Puf3p to those found in Puf5p demonstrated that a single amino acid change in Puf3p abolished detectable binding to COX17. Since this amino acid position in both Puf3p and Puf5p is predicted to contact an adenine in the respective target RNAs, the amino acid in Puf3p must play a more critical role in promoting COX17 interaction. In contrast, an amino acid change in the third repeat of Puf3p, which interacts with the only divergent nucleotide between the Puf3p and Puf5p targets, had no effect on binding COX17. These results argue that a simple set of rules cannot reliably link specific amino acid positions with target specificity. Each of these amino acid changes in Puf3p enhanced binding to the Puf5p target HO RNA, suggesting a different mode of binding to this target. Finally, we identified an outer surface loop that was dispensable for binding but was required to promote both rapid deadenylation and subsequent decapping of the COX17 mRNA, most likely as a point of protein-protein interactions.     

Probing the amino acids critical for protein oligomerisation and protein–nucleotide interaction in Mycobacterium tuberculosis PII protein through integration of computational and experimental approaches

We investigated the interacting amino acids critical for the stability and ATP binding of Mycobacterium tuberculosis PII protein through a series of site specific mutagenesis experiments. We assessed the effect of mutants using glutaraldehyde crosslinking and size exclusion chromatography and isothermal titration calorimetry. Mutations in the amino acid pair R60–E62 affecting central electrostatic interaction resulted in insoluble proteins. Multiple sequence alignment of PII orthologs displayed a conserved pattern of charged residues at these positions. Mutation of amino acid D97 to a neutral residue was tolerated whereas positive charge was not acceptable. Mutation of R107 alone had no effect on trimer formation. However, the combination of neutral residues both at positions 97 and 107 was not acceptable even with the pair at 60–62 intact. Reversal of charge polarity could partially restore the interaction. The residues including K90, R101 and R103 with potential to form H-bonds to ATP are conserved throughout across numerous orthologs of PII but when mutated to Alanine, they did not show significant differences in the total free energy change of the interaction as examined through isothermal titration calorimetry. The ATP binding pattern showed anti-cooperativity using three-site binding model. We observed compensatory effect in enthalpy and entropy changes and these may represent structural adjustments to accommodate ATP in the cavity even in absence of some interactions to perform the requisite function. In this respect these small differences between the PII orthologs may have evolved to suite species specific physiological niches.     

Characterization of the DNA binding protein encoded by the N-specific filamentous Escherichia coli phage IKe. Binding properties of the protein and nucleotide sequence of the gene

A DNA binding protein encoded by the filamentous single-stranded DNA phage IKe has been isolated from IKe-infected Escherichia coli cells. Fluorescence and in vitro binding studies have shown that the protein binds co-operatively and with a high specificity to single-stranded but not to double-stranded DNA. From titration of the protein to poly(dA) it has been calculated that approximately four bases of the DNA are covered by one monomer of protein. These binding characteristics closely resemble those of gene V protein encoded by the F-specific filamentous phages M13 and fd. The nucleotide sequence of the gene specifying the IKe DNA binding protein has been established. When compared to the nucleotide sequence of gene V of phage M13 it shows an homology of 58%, indicating that these two phages are evolutionarily related. The IKe DNA binding protein is 88 amino acids long which is one amino acid residue larger than the gene V protein sequence. When the IKe DNA binding protein sequence is compared with that of gene V protein it was found that 39 amino acid residues have identical positions in both proteins. The positions of all five tyrosine residues, a number of which are known to be involved in DNA binding, are conserved. Secondary structure predictions indicate that the two proteins contain similar structural domains. It is proposed that the tyrosine residues which are involved in DNA binding are the ones in or next to a beta-turn, at positions 26, 41 and 56 in gene V protein and at positions 27, 42 and 57 in the IKe DNA binding protein.     

Synthesis of alanyl nucleobase amino acids and their incorporation into proteins

Proteins which bind to nucleic acids and regulate their structure and functions are numerous and exceptionally important. Such proteins employ a variety of strategies for recognition of the relevant structural elements in their nucleic acid substrates, some of which have been shown to involve rather subtle interactions which might have been difficult to design from first principles. In the present study, we have explored the preparation of proteins containing unnatural amino acids having nucleobase side chains. In principle, the introduction of multiple nucleobase amino acids into the nucleic acid binding domain of a protein should enable these modified proteins to interact with their nucleic acid substrates using Watson-Crick and other base pairing interactions. We describe the synthesis of five alanyl nucleobase amino acids protected in a fashion which enabled their attachment to a suppressor tRNA, and their incorporation into each of two proteins with acceptable efficiencies. The nucleobases studied included cytosine, uracil, thymine, adenine and guanine, i.e. the major nucleobase constituents of DNA and RNA. Dihydrofolate reductase was chosen as one model protein to enable direct comparison of the facility of incorporation of the nucleobase amino acids with numerous other unnatural amino acids studied previously. The Klenow fragment of DNA polymerase I was chosen as a representative DNA binding protein whose mode of action has been studied in detail.     

Amino acid substitutions in the V domain of nectin-1 (HveC) that impair entry activity for herpes simplex virus types 1 and 2 but not for Pseudorabies virus or bovine herpesvirus 1

The entry of herpes simplex virus (HSV) into cells requires the interaction of viral glycoprotein D (gD) with a cellular gD receptor to trigger the fusion of viral and cellular membranes. Nectin-1, a member of the immunoglobulin superfamily, can serve as a gD receptor for HSV types 1 and 2 (HSV-1 and HSV-2, respectively) as well as for the animal herpesviruses porcine pseudorabies virus (PRV) and bovine herpesvirus 1 (BHV-1). The HSV-1 gD binding domain of nectin-1 is hypothesized to overlap amino acids 64 to 104 of the N-terminal variable domain-like immunoglobulin domain. Moreover, the HSV-1 and PRV gDs compete for binding to nectin-1. Here we report that two amino acids within this region, at positions 77 and 85, are critical for HSV-1 and HSV-2 entry but not for the entry of PRV or BHV-1. Replacement of either amino acid 77 or amino acid 85 reduced HSV-1 and HSV-2 gD binding but had a lesser effect on HSV entry activity, suggesting that weak interactions between gD and nectin-1 are sufficient to trigger the mechanism of HSV entry. Substitution of both amino acid 77 and amino acid 85 in nectin-1 significantly impaired entry activity for HSV-1 and HSV-2 and eliminated binding to soluble forms of HSV-1 and HSV-2 gDs but did not impair the entry of PRV and BHV-1. Thus, amino acids 77 and 85 of nectin-1 form part of the interface with HSV gD or influence the conformation of that interface. Moreover, the binding sites for HSV and PRV or BHV-1 gDs on nectin-1 may overlap but are not identical.     

Rat brain calbindin D28: six domain structure and extensive amino acid homology with chicken calbindin D28

Calbindin D28 cDNA clones were isolated from a rat brain library using a chicken intestinal Calbindin D28 cDNA probe. Nucleotide sequence analysis of these clones shows an open reading frame of 78 nucleotide coding for a 261 amino acid 29,994 dalton protein. The predicted amino acid sequence contains six repeats of a domain with the feature of an EF-hand calcium binding site. In domains II and VI, two of the five oxygen-containing amino acids important for the coordination of calcium are absent, suggesting that these two sites have lost their calcium-binding capability. Comparing the amino acid sequence to that recently reported for the chicken Calbindin D28 there is 79% homology. Tolerating conservative differences, the homology increases to 93%. Interestingly, domains II and VI which have presumably lost their calcium binding ability are very conserved among the two species (81% and 78%, respectively). Since an EF hand calcium binding site requires only certain types of amino acids at certain positions, rather than a specific amino acid sequence, maintaining a calcium binding site is a weak conservation pressure. To explain the high degree of homology of rat and chicken Calbindin D28, and in particular the conservation of the two degenerated domains over the 300 million years since divergence of birds and mammals, additional function(s) of the Calbindin D28 are postulated.     

The ORFa protein, the putative transposase of maize transposable element Ac, has a basic DNA binding domain.

Ac encodes the 807 amino acid ORFa protein which binds specifically to multiple AAACGG motifs that are subterminally located in both ends of Ac. The wild-type ORFa protein and a number of deletion and amino acid exchange mutants were expressed in Escherichia coli, renatured and used for mobility shift assays. At least 136 amino acids from the N-terminus and 537 C-terminal amino acids may be removed from the ORFa protein without destroying the DNA binding domain, whereas a protein starting at amino acid 189 is DNA binding deficient. Certain basic amino acids between positions 190 and 200 are essential for DNA binding, as their substitution with uncharged amino acids leads to the loss of this function. The DNA binding domain of ORFa protein has an overall basic character, but no obvious sequence homology to any other known DNA binding protein. The homologies to the major open reading frames of transposable elements Tam3 from Antirrhinum majus and Hobo from Drosophila are found between the C-terminal two thirds of the three proteins. The ORFa protein forms discrete complexes with target DNA that appear, depending on the protein concentration, as a 'ladder' of bands on gels, indicating the occupation of target DNA by multiple ORFa protein molecules.     

The amino acid sequence of satyr tragopan lysozyme and its activity

The amino acid sequence of satyr tragopan lysozyme and its activity was analyzed. Carboxymethylated lysozyme was digested with trypsin and the resulting peptides were sequenced. The established amino acid sequence had three amino acid substitutions at positions 103 (Asn to Ser), 106 (Ser to Asn), and 121 (His to Gln) comparing with Temminck's tragopan lysozyme and five amino acid substitutions at positions 3 (Phe to Tyr), 15 (His to Leu), 41 (Gln to His), 101 (Asp to Gly) and 103 (Asn to Ser) with chicken lysozyme. The time course analysis using N-acetylglucosamine pentamer as a substrate showed a decrease of binding free energy change, 1.1 kcal/mol at subsite A and 0.2 kcal/mol at subsite B, between satyr tragopan and chicken lysozymes. This was assumed to be responsible for the amino acid substitutions at subsite A-B at position 101 (Asp to Gly), however another substitution at position 103 (Asn to Ser) considered not to affect the change of the substrate binding affinity by the observation of identical time course of satyr tragopan lysozyme with turkey and Temminck's tragopan lysozymes that carried the identical amino acids with chicken lysozyme at this position. These results indicate that the observed decrease of binding free energy change at subsites A-B of satyr tragopan lysozyme was responsible for the amino acid substitution at position 101 (Asp to Gly).     

Kinetic and modeling studies of S3-S3' subsites of HIV proteinases.

Kinetic analysis and modeling studies of HIV-1 and HIV-2 proteinases were carried out using the oligopeptide substrate [formula: see text] and its analogs containing single amino acid substitutions in P3-P3' positions. The two proteinases acted similarly on the substrates except those having certain hydrophobic amino acids at P2, P1, P2', and P3' positions (Ala, Leu, Met, Phe). Various amino acids seemed to be acceptable at P3 and P3' positions, while the P2 and P2' positions seemed to be more restrictive. Polar uncharged residues resulted in relatively good binding at P3 and P2 positions, while at P2' and P3' positions they gave very high Km values, indicating substantial differences in the respective S and S' subsites of the enzyme. Lys prevented substrate hydrolysis at any of the P2-P2' positions. The large differences for subsite preference at P2 and P2' positions seem to be at least partially due to the different internal interactions of P2 residue with P1', and P2' residue with P1. As expected on the basis of amino acid frequency in the naturally occurring cleavage sites, hydrophobic residues at P1 position resulted in cleavable peptides, while polar and beta-branched amino acids prevented hydrolysis. On the other hand, changing the P1' Pro to other amino acids prevented substrate hydrolysis, even if the substituted amino acid had produced a good substrate in other oligopeptides representing naturally occurring cleavage sites. The results suggest that the subsite specificity of the HIV proteinases may strongly depend on the sequence context of the substrate.     

Signal peptide of Bacillus subtilis alpha-amylase

Mature alpha-amylase of Bacillus subtilis is known to be formed from its precursor by the removal of the NH2-terminal 41 amino acid sequence (41 amino acid leader sequence). DNA fragments coding for short sequences consisting of 28 (Pro as the COOH terminus) 29 (Ala), 31 (Ala), and 33 (Ala) amino acids from the translation initiator, Met, in the leader sequence were prepared and fused in frame to the DNA encoding the mature alpha-amylase. The secretion activity of the 33 amino acid sequence was nearly twice as high as that of the parental 41 amino acid sequence, whereas the activity of the 31 amino acid sequence was 75% of that of the parent. In contrast, almost no secretion activity was observed with the 28 and 29 amino acid sequences. The signal peptide cleavage site of the precursor expressed from the plasmid encoding the 33 amino acid sequence was located between Ala and Leu at positions 33 and 34 and that from the 31 amino acid sequence between Thr and Ala at positions 33 and 34. The NH2-terminal amino acid from the latter corresponded to the 3rd amino acid of the mature enzyme. These results indicated that the functional signal peptide of the B. subtilis beta-amylase consists of the first 33 amino acids from the initiator, Met.     

Segmentally variable genes: a new perspective on adaptation

Genomic sequence variation is the hallmark of life and is key to understanding diversity and adaptation among the numerous microorganisms on earth. Analysis of the sequenced microbial genomes suggests that genes are evolving at many different rates. We have attempted to derive a new classification of genes into three broad categories: lineage-specific genes that evolve rapidly and appear unique to individual species or strains; highly conserved genes that frequently perform housekeeping functions; and partially variable genes that contain highly variable regions, at least 70 amino acids long, interspersed among well-conserved regions. The latter we term segmentally variable genes (SVGs), and we suggest that they are especially interesting targets for biochemical studies. Among these genes are ones necessary to deal with the environment, including genes involved in host-pathogen interactions, defense mechanisms, and intracellular responses to internal and environmental changes. For the most part, the detailed function of these variable regions remains unknown. We propose that they are likely to perform important binding functions responsible for protein-protein, protein-nucleic acid, or protein-small molecule interactions. Discerning their function and identifying their binding partners may offer biologists new insights into the basic mechanisms of adaptation, context-dependent evolution, and the interaction between microbes and their environment.