Distinguishing structural and functional restraints in evolution in order to identify interaction sites

Structural genomics projects are producing many three-dimensional structures of proteins that have been identified only from their gene sequences. It is therefore important to develop computational methods that will predict sites involved in productive intermolecular interactions that might give clues about functions. Techniques based on evolutionary conservation of amino acids have the advantage over physiochemical methods in that they are more general. However, the majority of techniques neither use all available structural and sequence information, nor are able to distinguish between evolutionary restraints that arise from the need to maintain structure and those that arise from function. Three methods to identify evolutionary restraints on protein sequence and structure are described here. The first identifies those residues that have a higher degree of conservation than expected: this is achieved by comparing for each amino acid position the sequence conservation observed in the homologous family of proteins with the degree of conservation predicted on the basis of amino acid type and local environment. The second uses information theory to identify those positions where environment-specific substitution tables make poor predictions of the overall amino acid substitution pattern. The third method identifies those residues that have highly conserved positions when three-dimensional structures of proteins in a homologous family are superposed. The scores derived from these methods are mapped onto the protein three-dimensional structures and contoured, allowing identification clusters of residues with strong evolutionary restraints that are sites of interaction in proteins involved in a variety of functions. Our method differs from other published techniques by making use of structural information to identify restraints that arise from the structure of the protein and differentiating these restraints from others that derive from intermolecular interactions that mediate functions in the whole organism.     

SpaK/SpaR Two-component System Characterized by a Structure-driven Domain-fusion Method and in Vitro Phosphorylation Studies

Here we introduce a quantitative structure-driven computational domain-fusion method, which we used to predict the structures of proteins believed to be involved in regulation of the subtilin pathway in Bacillus subtilis, and used to predict a protein-protein complex formed by interaction between the proteins. Homology modeling of SpaK and SpaR yielded preliminary structural models based on a best template for SpaK comprising a dimer of a histidine kinase, and for SpaR a response regulator protein. Our LGA code was used to identify multi-domain proteins with structure homology to both modeled structures, yielding a set of domain-fusion templates then used to model a hypothetical SpaK/SpaR complex. The models were used to identify putative functional residues and residues at the protein-protein interface, and bioinformatics was used to compare functionally and structurally relevant residues in corresponding positions among proteins with structural homology to the templates. Models of the complex were evaluated in light of known properties of the functional residues within two-component systems involving His-Asp phosphorelays. Based on this analysis, a phosphotransferase complexed with a beryllofluoride was selected as the optimal template for modeling a SpaK/SpaR complex conformation. In vitro phosphorylation studies performed using wild type and site-directed SpaK mutant proteins validated the predictions derived from application of the structure-driven domain-fusion method: SpaK was phosphorylated in the presence of ³²P-ATP and the phosphate moiety was subsequently transferred to SpaR, supporting the hypothesis that SpaK and SpaR function as sensor and response regulator, respectively, in a two-component signal transduction system, and furthermore suggesting that the structure-driven domain-fusion approach correctly predicted a physical interaction between SpaK and SpaR. Our domain-fusion algorithm leverages quantitative structure information and provides a tool for generation of hypotheses regarding protein function, which can then be tested using empirical methods.     

Critical amino acids of p21 protein are located within beta-turns: further evaluation

It has been shown that malignant activation of ras proto-oncogenes was mediated by point mutations which resulted in the single amino acid conversions at positions 12, 13 or 61 of the ras gene products (p21 proteins). By analyzing randomly mutated ras genes, it has been demonstrated that amino acid substitutions at residues 12, 13, 59 and 63 activated p21. Furthermore, it has been shown that residues 16, 116 and 119 in p21 played critical roles in the guanine nucleotide binding and, consequently, the ability of the protein to induce changes characteristic of cellular transformation. By using the protein conformational prediction method of Chou and Fasman, the present work predicts that these critical amino acids, except glutamic acid at position 63, are located within beta-turns. The major "hot spots" for ras activation are codons 12 and 61. The author has predicted in an earlier paper that the single amino acid conversions at positions 12 and 61 would occur at beta-turn conformation consisting of residues 10-13 and 58-61, respectively. In the present study, probabilities of beta-turn occurrence at residues 10-13 or 58-61 of the p21 proteins encoded by various ras genes are compared. The probability for the normal p21 containing glycine as residue 12 is greatest, and the cancer-associated variants show less probabilities. The single amino acid substitutions at position 61 do not cause so decreased probabilities of beta-turn potential at residues 58-61, except the replacement by histidine. Histidine at position 61 is not predicted as occurring within a beta-turn.(ABSTRACT TRUNCATED AT 250 WORDS)     

Role of paired basic residues of protein C-termini in phospholipid binding

It is a well known phenomenon that the occurrence of several distinct amino acids at the C-terminus of proteins is non-random. We have analysed all Saccharomyces cerevisiae proteins predicted by computer databases and found lysine to be the most frequent residue both at the last (-1) and at the penultimate amino acid (-2) positions. To test the hypothesis that C-terminal basic residues efficiently bind to phospholipids we randomly expressed GST-fusion proteins from a yeast genomic library. Fifty-four different peptide fragments were found to bind phospholipids and 40% of them contained lysine/arginine residues at the (-1) or (-2) positions. One peptide showed high sequence similarity with the yeast protein Sip18p. Mutational analysis revealed that both C-terminal lysine residues of Sip18p are essential for phospholipid-binding in vitro. We assume that basic amino acid residues at the (-1) and (-2) positions in C-termini are suitable to attach the C-terminus of a given protein to membrane components such as phospholipids, thereby stabilizing the spatial structure of the protein or contributing to its subcellular localization. This mechanism could be an additional explanation for the C-terminal amino acid bias observed in proteins of several species.     

The protein structure code: what is its present status?

Current methods of prediction of protein conformation are reviewedand the algorithms on which they rely are presented. For non-homologousproteins and after cross-validation the reported methods exhibita probability index, i.e. the per cent of correctly predictedresidues per predicted residues, of 63–65% with a standarddeviation of the order of 7% for three conformational states—helix,ß-strand and coil. This present limitation in theaccuracy of predictions that use only the information of thelocal sequence can be related essentially to the effect of long-rangeinteractions specific for each protein family. The methods basedon sequence similarity can improve the accuracy of predictionby expressing explicitly the homology of the protein to be predictedwith proteins in the database. In these circumstances the probabilityindex can reach 87% with a standard deviation of 6.6%. Thisproperty can be used for modeling homologous proteins by aidingin amino acid sequence alignments. The prediction of the tertiarystructure of a protein is still limited to the case of modelinga structure based on the known three-dimensional structure ofa homologous protein.     

Determination of Amino Acid Residues Responsible for Specific Interaction of Protein Kinases with Small Molecule Inhibitors

Identifying amino acid positions that determine the specific interaction of proteins with small molecule ligands, is required for search of pharmaceutical targets, drug design, and solution of other biotechnology problems. We studied applicability of an original method SPrOS (specificity projection on sequence) developed to recognize functionally significant positions in amino acid sequences. The method allows residues specific to functional subgroups to be determined within the protein family based on their local surroundings in amino acid sequences. The efficiency of the method has been estimated on the protein kinase family. The residues associated with the protein specificity to inhibitors have been predicted. The results have been verified using 3D structures of protein–ligand complexes. Three small molecule inhibitors have been tested. Residues predicted with SPrOS either in contacted the inhibitor or influenced the conformation of the ligand–binding area. Excluding close homologues from the studied set makes it possible to decrease the number of difficult to interpret positions. The expediency of this procedure was determined by the relationship between an inhibitory spectrum and phylogenic partition. Thus, the method efficiency has been confirmed by matching the prediction results with the protein 3D structures.     

Statistical analyses of protein sequence alignments identify structures and mechanisms in signal activation of sensor histidine kinases

Statistical analyses of genome sequence‐derived protein sequence data can identify amino acid residues that interact between proteins or between domains of a protein. These statistical methods are based on evolution‐directed amino acid variation responding to structural and functional constraints in proteins. The identified residues form a basis for determining structure and folding of proteins as well as inferring mechanisms of protein function. When applied to two‐component systems, several research groups have shown they can be used to identify the amino acid interactions between response regulators and histidine kinases and the specificity therein. Recently, statistical studies between the HisKA and HATPase‐ATP‐binding domains of histidine kinases identified amino acid interactions for both the inactive and the active catalytic states of such kinases. The identified interactions generated a model structure for the domain conformation of the active state. This conformation requires an unwinding of a portion of the C‐terminal helix of the HisKA domain that destroys the inactive state residue contacts and suggests how signal‐binding determines the equilibrium between the inactive and active states of histidine kinases. The rapidly accumulating protein sequence databases from genome, metagenome and microbiome studies are an important resource for functional and structural understanding of proteins and protein complexes in microbes.     

Crystal structure of a recombinant form of the maltodextrin-binding protein carrying an inserted sequence of a B-cell epitope from the preS2 region of hepatitis B virus

We report the crystal structure of MalE-B133, a recombinant form of the maltodextrin-binding protein (MBP) of Escherichia coli carrying an inserted amino-acid sequence of a B-cell epitope from the preS2 region of the hepatitis B virus (HBV). The structure was determined by molecular replacement methods and refined to 2.7 Å resolution. MalE-B133 is an insertion/deletion mutant of MBP in which residues from positions 134 to 142, an external α helix in the wild-type structure, are replaced by a foreign peptide segment of 19 amino acids. The inserted residues correspond to the preS2 sequence from positions 132 to 145 and five flanking residues that arise from the creation of restriction sites. The conformation of the recombinant protein, excluding the inserted segment, closely resembles that of wild-type MBP in the closed maltose-bound form. MalE-B133 was shown by previous studies to display certain immunogenic and antigenic properties of the hepatitis B surface antigen (HBsAg), which contains the preS2 region. The crystal structure reveals the conformation of the first nine epitope residues (preS2 positions 132 to 140) exposed on the surface of the molecule. The remaining five epitope residues (preS2 positions 141 to 145) are not visible in electron density maps. The path of the polypeptide chain in the visible portion of the insert differs from that of the deleted segment in the structure of wild-type MBP, displaying a helical conformation at positions 134 to 140 (preS2 sequence numbering). A tripeptide (Asp-Pro-Arg) at the N terminus of the helix forms a stable structural motif that may be implicated in the cross-reactivity of anti-HBsAg antibodies with the hybrid protein. Proteins 27:1–8 © 1997 Wiley-Liss, Inc.     

A genetic screen to identify variants of bovine pancreatic trypsin inhibitor with altered folding energetics.

A genetic screening procedure has been developed to identify mutant forms of bovine pancreatic trypsin inhibitor (BPTI) that can fold to an active conformation but are inactivated more rapidly than the wild-type protein. Small cultures of Escherichia coli containing plasmids with mutagenized BPTI genes were grown in microtiter plates, lysed, and treated with dithiothreitol (DTT). Under these conditions, unfolding and inactivation of the wild-type protein has a half-time of about 10 hours. Variants of BPTI that are inactivated within 1 hour were identified by adding trypsin and a chromogenic substrate. Approximately 11,000 mutagenized clones were screened in this way and 75 clones that produce proteins that can fold but are inactivated by DTT were isolated. The genes coding for 68 "DTT-sensitive" mutant proteins were sequenced, and 25 different single amino acid substitutions at 15 of the 58 residues of the protein were identified. Most of the altered residues are largely buried in the core of the native wild-type structure and are highly conserved among proteins homologous to BPTI. These results indicate that a large fraction of the sequence of the protein contributes to the kinetic stability of the active conformation, but it also appears that substitutions can be tolerated at most sites without completely preventing folding. Because this genetic screen is based on changes in folding energetics, further studies of the isolated mutants are expected to provide information about the roles of the altered residues in folding and unfolding.     

Protein structure prediction

Current methods developed for predicting protein structure are reviewed. The most widely used algorithms of Chou and Fasman and Garnier et al for predicting secondary structure are compared to the most recent ones including sequence similarity methods, neural network, pattern recognition or joint prediction methods. The best of these methods correctly predict 63-65% of the residues in the database with cross-validation for 3 conformations, helix, beta strand and coli with a standard deviation of 6-8% per protein. However, when a homologous protein is already in the database, the accuracy of prediction by the similarity peptide method of Levin and Garnier reaches about 90%. Some conclusions can be drawn on the mechanism of protein folding. As all the prediction methods only use the local sequence for prediction (+/- 8 residues maximum) one can infer that 65% of the conformation of a residue is dictated on average by the local sequence, the rest is brought by the folding. The best predicted proteins or peptide segments are those for which the folding has less effect on the conformation. Presently, prediction of tertiary structure is only of practical use when the structure of a homologous protein is already known. Amino acid alignment to define residues of equivalent spatial position is critical for modelling of the protein. We showed for serine proteases that secondary structure prediction can help to define a better alignment. Non-homologous segments of the polypeptide chain, such as loops, libraries of known loops and/or energy minimization with various force fields, are used without yet giving satisfactory solutions. An example of modelling by homology, aided by secondary structure prediction on 2 regulatory proteins, Fnr and FixK is presented.     

Prediction of β-sheets in immunoglobulin chains. Comparison of various methods and an expanded 20 × 20 table for evaluation of the effects of nearest-neighbors on conformations of middle amino acids in proteins

The extraordinarily large number of immunoglobulins renders them an intriguing class of molecules for attempts to predict their conformations. The predictive method applied, using a 20 × 20 table of the observed effects of nearest-neighboring amino acids on the conformation (Φ,Ψ angles) of the middle residue in known proteins, indicates positions of tri-peptides that tend to break α-helices or regular β-sheets. This 20 × 20 table is derived from data on 19 proteins, as compared with the earlier version based on 12 proteins, and includes a separate listing of residues of β-turns that have helical Φ,Ψ values. Secondary conformations predicted by methods of Chou and Fasman, Lim and Burgess, Ponnuswamy, and Scheraga have also been compared; for all three methods, wrong predicitons of residues in β-sheet conformation exceed correct ones. Better predictions are obtained when there is agreement with two or three of the methods. If there is consistent overprediction of β-structure, as with the Chou and Fasman method, the use of the β-sheet-breaking tripeptides can improve pre-dictability somewhat.     

Comparison of the computed structures for the phosphate-binding loop of the p21 protein containing the oncogenic site Gly 12 with the X-ray crystallographic structures for this region in the p21 protein and EFtu. A model for the structure of the p21 protein in its oncogenic form

The GTP-binding p21 protein encoded by the ras-oncogene can be activated to cause malignant transformation of cells by substitution of a single amino acid at critical positions along the polypeptide chain. Substitution of any non-cyclic L-amino acid for Gly 12 in the normal protein results in a transforming protein. This substitution occurs in a hydrophobic sequence (residues 6-15) which is known to be involved in binding the phosphate moities of GTP (and GDP). We find, using conformational energy calculations, that the 6-15 segment of the normal protein (with Gly 12) adopts structures that contain a bend at residues 11 and 12 with the Gly in the D* conformation, not allowed energetically for L-amino acids. Substitution of non-cyclic L-amino acids for Gly 12 results in shifting this bend to residues 12 and 13. We show that many computed structures for the Gly 12-containing phosphate binding loop, segment 9-15, are superimposable on the corresponding segment of the recently determined X-ray crystallographic structure for residues 1-171 of the p21 protein. All such structures contain bends at residues 11 and 12 and most of these contain Gly 12 in the C* or D* conformational state. Other computed conformations for the 9-15 segment were superimposable on the structure of the corresponding 18-23 segment of EFtu, the bacterial chain elongation factor having structural similarities to the p21 protein in the phosphate-binding regions. This segment contains a Val residue where a Gly occurs in the p21 protein. As previously predicted, all of these superimposable conformations contain a bend at positions 12 and 13, not 11 and 12. If these structures that are superimposable on EFtu are introduced into the p21 protein structure, bad contacts occur between the sidechain of the residue (here Val) at position 12 and another phosphate binding loop region around position 61. These bad contacts between the two segments can be removed by changing the conformation of the 61 region in the p21 protein to the corresponding position of the homologous region in EFtu. In this new conformation, a large site becomes available for the binding of phosphate residues. In addition, such phenomena as autophosphorylation of the p21 protein by GTP can be explained with this new model structure for the activated protein which cannot be explained by the structure for the non-activated protein.     

Prediction of amyloid fibril-forming proteins

In Alzheimer's disease and spongiform encephalopathies proteins transform from their native states into fibrils. We find that several amyloid-forming proteins harbor an alpha-helix in a polypeptide segment that should form a beta-strand according to secondary structure predictions. In 1324 nonredundant protein structures, 37 beta-strands with > or =7 residues were predicted in segments where the experimentally determined structures show helices. These discordances include the prion protein (helix 2, positions 179-191), the Alzheimer amyloid beta-peptide (Abeta, positions 16-23), and lung surfactant protein C (SP-C, positions 12-27). In addition, human coagulation factor XIII (positions 258-266), triacylglycerol lipase from Candida antarctica (positions 256-266), and d-alanyl-d-alanine transpeptidase from Streptomyces R61 (positions 92-106) contain a discordant helix. These proteins have not been reported to form fibrils but in this study were found to form fibrils in buffered saline at pH 7.4. By replacing valines in the discordant helical part of SP-C with leucines, an alpha-helix is found experimentally and by secondary structure predictions. This analogue does not form fibrils under conditions where SP-C forms abundant fibrils. Likewise, when Abeta residues 14-23 are removed or changed to a nondiscordant sequence, fibrils are no longer formed. We propose that alpha-helix/beta-strand-discordant stretches are associated with amyloid fibril formation.     

Amino acid sequence homology applied to the prediction of protein secondary structures, and joint prediction with existing methods

The assumption that homologous segments in different proteins may share a similar conformation is applied to the prediction of secondary structures in proteins. Sequences homologous to a target protein are searched, without allowing any gap, and compared against a number of reference proteins of known three-dimensional structure, and then a conformational state (alpha, beta or coil) for each residue of the protein is predicted by looking at the secondary structure of corresponding homologous segments. This prediction is done in a statistical rather than 'deterministic' way, by assigning the most probable conformation state among homologous data to each residue site of a target protein. A test application for 22 sample proteins yields 60% correctness on the average, a better value in comparison with two other existing methods. Joint prediction combining three methods into one is shown to increase the reliability up to 70%, when only the regions identically predicted with the three methods are taken into account. Application of the present method to 10 proteins of unknown structure is demonstrated.     

A residue in helical conformation in the native state adopts a β-strand conformation in the folding transition state despite its high and canonical Φ-value

Delineating structures of the transition states in protein folding reactions has provided great insight into the mechanisms by which proteins fold. The most common method for obtaining this information is Φ-value analysis, which is carried out by measuring the changes in the folding and unfolding rates caused by single amino acid substitutions at various positions within a given protein. Canonical Φ-values range between 0 and 1, and residues displaying high values within this range are interpreted to be important in stabilizing the transition state structure, and to elicit this stabilization through native-like interactions. Although very successful in defining the general features of transition state structures, Φ-value analysis can be confounded when non-native interactions stabilize this state. In addition, direct information on backbone conformation within the transition state is not provided. In the work described here, we have investigated structure formation at a conserved β-bulge (with helical conformation) in the Fyn SH3 domain by characterizing the effects of substituting all natural amino acids at one position within this structural motif. By comparing the effects on folding rates of these substitutions with database-derived local structure propensity values, we have determined that this position adopts a non-native backbone conformation in the folding transition state. This result is surprising because this position displays a high and canonical Φ-value of 0.7. This work emphasizes the potential role of non-native conformations in folding pathways and demonstrates that even positions displaying high and canonical Φ-values may, nevertheless, adopt a non-native conformation in the transition state.     

Prediction of the three-dimensional structure of the rap-1A protein from its homology to theras-gene-encoded p21 protein

rap-1A, an anti-oncogene-encoded protein, is aras-p21-like protein whose sequence is over 80% homologous to p21 and which interacts with the same intracellular target proteins and is activated by the same mechanisms as p21, e.g., by binding GTP in place of GDP. Both interact with effector proteins in the same region, involving residues 32–47. However, activated rap-1A blocks the mitogenic signal transducing effects of p21. Optimal sequence alignment of p21 and rap-1A shows two insertions of rap-1A atras positions 120 and 138. We have constructed the three-dimensional structure of rap-1A bound to GTP by using the energy-minimized three-dimensional structure ofras-p21 as the basis for the modeling using a stepwise procedure in which identical and homologous amino acid residues in rap-1A are assumed to adopt the same conformation as the corresponding residues in p21. Side-chain conformations for homologous and nonhomologous residues are generated in conformations that are as close as possible to those of the corresponding side chains in p21. The entire structure has been subjected to a nested series of energy minimizations. The final predicted structure has an overall backbone deviation of 0.7 å from that ofras-p21. The effector binding domains from residues 32–47 are identical in both proteins (except for different side chains of different residues at position 45). A major difference occurs in the insertion region at residue 120. This region is in the middle of another effector loop of the p21 protein involving residues 115–126. Differences in sequence and structure in this region may contribute to the differences in cellular functions of these two proteins.     

Prediction of protein functional residues from sequence by probability density estimation

MOTIVATION: The prediction of ligand-binding residues or catalytically active residues of a protein may give important hints that can guide further genetic or biochemical studies. Existing sequence-based prediction methods mostly rank residue positions by evolutionary conservation calculated from a multiple sequence alignment of homologs. A problem hampering more wide-spread application of these methods is the low per-residue precision, which at 20% sensitivity is around 35% for ligand-binding residues and 20% for catalytic residues. RESULTS: We combine information from the conservation at each site, its amino acid distribution, as well as its predicted secondary structure (ss) and relative solvent accessibility (rsa). First, we measure conservation by how much the amino acid distribution at each site differs from the distribution expected for the predicted ss and rsa states. Second, we include the conservation of neighboring residues in a weighted linear score by analytically optimizing the signal-to-noise ratio of the total score. Third, we use conditional probability density estimation to calculate the probability of each site to be functional given its conservation, the observed amino acid distribution, and the predicted ss and rsa states. We have constructed two large data sets, one based on the Catalytic Site Atlas and the other on PDB SITE records, to benchmark methods for predicting functional residues. The new method FRcons predicts ligand-binding and catalytic residues with higher precision than alternative methods over the entire sensitivity range, reaching 50% and 40% precision at 20% sensitivity, respectively. AVAILABILITY: Server: http://frpred.tuebingen.mpg.de. Data sets: ftp://ftp.tuebingen.mpg.de/pub/protevo/FRpred/.     

Single-molecule nanopositioning: structural transitions of a helicase-DNA complex during ATP hydrolysis

The conformational states of Escherichia coli Rep helicase undergoing ATP hydrolysis while bound to a partial-duplex DNA (pdDNA) were studied using single-molecule FRET. Crystallographic studies showed that Rep bound to single-stranded DNA can exist in open and closed conformations that differ in the orientation of the 2B subdomain. FRET measurements between eight Rep mutants donor-labeled at different residues and pdDNA acceptor-labeled at the junction were conducted at each of the four nucleotide states. The positions of donor-labeled residues, based on crystal structure, and FRET measurements between these donor molecules and the acceptor fluorophore at the DNA junction were used to predict the most likely position for the DNA junction using a triangulation algorithm. These predicted junction positions are compared with the crystal structure to determine whether the open or closed conformation is more consistent with the FRET data. Our data revealed that there are two distinct Rep-pdDNA conformations in the ATPγS and ADP states, an unexpected finding. The primary conformation is similar to that observed in nucleotide-free and ADP.Pi states, and the secondary conformation is a novel conformation where the duplex DNA and 2B subdomain moved as a unit by 13 Å relative to the rest of the protein. The primary conformation found in all nucleotide states is consistent with the closed conformation of the crystal structure however; the secondary conformation is a new conformation that has not been observed before. We discuss the possible implications of this newly observed conformation.     

Processing of Escherichia coli alkaline phosphatase. Sequence requirements and possible conformations of the -6 to -4 region of the signal peptide

Analysis of the precursors of bacterial exported proteins revealed that those having bulky hydrophobic residues at position -5 have a high incidence of Pro residues at positions -6 and -4, Val at position -3, and Ser at positions -4 and -2. This led to a hypothesis that the previously observed inhibition of processing by bulky residues at position -5 can be suppressed by introduction of Pro, Ser, or Val in the corresponding nearby positions. Subsequent mutational analysis of Escherichia coli alkaline phosphatase showed that, as it was predicted, Pro on either side of bulky hydrophobic -5 Leu, Ile, or Tyr completely restores efficiency of the maturation. Introduction of Val at position -3 also partially suppresses the inhibition imposed by -5 Leu, while a Ser residue at position -4 or -2 does not restore processing. In addition, effective maturation of a mutant with Pro residues at positions from -6 throughout -4 proved that polyproline conformation of this region is permissive for processing. To understand the effects of the mutations, we modeled a peptide substrate into the active site of the signal peptidase using the known position of the beta-lactam inhibitor. The inhibitory effect of the -5 residue and its suppression by either Pro -6 or Pro -4 can be explained if we assume that Pro-containing -6 to -4 regions adopt a polyproline conformation whereas the region without Pro residues has a beta-conformation. These results permit us to specify sequence requirements at -6, -5, and -4 positions for efficient processing and to improve the prediction of yet unknown cleavage sites.