The importance of surface loops for stabilizing an eightfold beta alpha barrel protein.

An important step in understanding how a protein folds is to determine those regions of the sequence that are critical to both its stability and its folding pathway. We chose phosphoribosyl anthranilate isomerase from Escherichia coli, which is a monomeric representative of the (beta alpha)8 barrel family of proteins, to construct a variant that carries an internal tandem duplication of the fifth beta alpha module. This (beta alpha)9 variant was enzymically active and therefore must have a wild-type (beta alpha)8 core. It had a choice a priori to fold to three different folding frames, which are distinguished by carrying the duplicated segment as an insert into one out of three different loops. Steady-state kinetic constants, the fluorescence properties of a crucial tryptophan residue, and limited proteolysis showed that the stable (beta alpha)9 variant carries the insertion between beta-strand 5 and alpha-helix 5. This preference can be explained by the important role of loops between alpha helices and beta strands in stabilizing the structure of the enzyme.     

Analysis of interactive packing of secondary structural elements in alpha/beta units in proteins

An alpha-helix and a beta-strand are said to be interactively packed if at least one residue in each of the secondary structural elements loses 10% of its solvent accessible contact area on association with the other secondary structural element. An analysis of all such 5,975 nonidentical alpha/beta units in protein structures, defined at < or = 2.5 A resolution, shows that the interaxial distance between the alpha-helix and the beta-strand is linearly correlated with the residue-dependent function, log[(V/nda)/n-int], where V is the volume of amino acid residues in the packing interface, nda is the normalized difference in solvent accessible contact area of the residues in packed and unpacked secondary structural elements, and n-int is the number of residues in the packing interface. The beta-sheet unit (beta u), defined as a pair of adjacent parallel or antiparallel hydrogen-bonded beta-strands, packing with an alpha-helix shows a better correlation between the interaxial distance and log(V/nda) for the residues in the packing interface. This packing relationship is shown to be useful in the prediction of interaxial distances in alpha/beta units using the interacting residue information of equivalent alpha/beta units of homologous proteins. It is, therefore, of value in comparative modeling of protein structures.     

Sequence specificity, statistical potentials, and three-dimensional structure prediction with self-correcting distance geometry calculations of beta-sheet formation in proteins

A statistical analysis of a representative data set of 169 known protein structures was used to analyze the specificity of residue interactions between spatial neighboring strands in beta-sheets. Pairwise potentials were derived from the frequency of residue pairs in nearest contact, second nearest and third nearest contacts across neighboring beta-strands compared to the expected frequency of residue pairs in a random model. A pseudo-energy function based on these statistical pairwise potentials recognized native beta-sheets among possible alternative pairings. The native pairing was found within the three lowest energies in 73% of the cases in the training data set and in 63% of beta-sheets in a test data set of 67 proteins, which were not part of the training set. The energy function was also used to detect tripeptides, which occur frequently in beta-sheets of native proteins. The majority of native partners of tripeptides were distributed in a low energy range. Self-correcting distance geometry (SECODG) calculations using distance constraints sets derived from possible low energy pairing of beta-strands uniquely identified the native pairing of the beta-sheet in pancreatic trypsin inhibitor (BPTI). These results will be useful for predicting the structure of proteins from their amino acid sequence as well as for the design of proteins containing beta-sheets.     

Estimation of the number of alpha-helical and beta-strand segments in proteins using circular dichroism spectroscopy

A simple approach to estimate the number of alpha-helical and beta-strand segments from protein circular dichroism spectra is described. The alpha-helix and beta-sheet conformations in globular protein structures, assigned by DSSP and STRIDE algorithms, were divided into regular and distorted fractions by considering a certain number of terminal residues in a given alpha-helix or beta-strand segment to be distorted. The resulting secondary structure fractions for 29 reference proteins were used in the analyses of circular dichroism spectra by the SELCON method. From the performance indices of the analyses, we determined that, on an average, four residues per alpha-helix and two residues per beta-strand may be considered distorted in proteins. The number of alpha-helical and beta-strand segments and their average length in a given protein were estimated from the fraction of distorted alpha-helix and beta-strand conformations determined from the analysis of circular dichroism spectra. The statistical test for the reference protein set shows the high reliability of such a classification of protein secondary structure. The method was used to analyze the circular dichroism spectra of four additional proteins and the predicted structural characteristics agree with the crystal structure data.     

Two peptide fragments G55-I72 and K97-A109 from staphylococcal nuclease exhibit different behaviors in conformational preferences for helix formation

Two synthetic peptides, SNasealpha1 and SNasealpha2, corresponding to residues G55-I72 and K97-A109, respectively, of staphylococcal nuclease (SNase), are adopted for detecting the role of helix alpha1 (E57-A69) and helix alpha2 (M98-Q106) in the initiation of folding of SNase. The helix-forming tendencies of the two SNase peptide fragments are investigated using circular dichroism (CD) and two-dimensional (2D) nuclear magnetic resonance (NMR) methods in water and 40% trifluoroethanol (TFE) solutions. The coil-helix conformational transitions of the two peptides in the TFE-H2O mixture are different from each other. SNasealpha1 adopts a low population of localized helical conformation in water, and shows a gradual transition to helical conformation with increasing concentrations of TFE. SNasealpha2 is essentially unstructured in water, but undergoes a cooperative transition to a predominantly helical conformation at high TFE concentrations. Using the NMR data obtained in the presence of 40% TFE, an ensemble of alpha-helical structures has been calculated for both peptides in the absence of tertiary interactions. Analysis of all the experimental data available indicates that formation of ordered alpha-helical structures in the segments E57-A69 and M98-Q106 of SNase may require nonlocal interactions through transient contact with hydrophobic residues in other parts of the protein to stabilize the helical conformations in the folding. The folding of helix alpha1 is supposed to be effective in initiating protein folding. The formation of helix alpha2 depends strongly on the hydrophobic environment created in the protein folding, and is more important in the stabilization of the tertiary conformation of SNase.     

A helix-turn-strand structural motif common in alpha-beta proteins.

By exhaustive structural comparisons, we have found that about one-third of the alpha-helix-turn-beta-strand polypeptides in alpha-beta barrel domains share a common structural motif. The chief characteristics of this motif are that first, the geometry of the turn between the alpha-helix and the beta-strand is somewhat constrained, and second, the beta-strand contains a hydrophobic patch that fits into a hydrophobic pocket on the alpha-helix. The geometry of the turn does not seem to be a major determinant of the alpha-beta unit, because the turns vary in length from four to six residues. However, the motif does not occur when there are few constraints on the geometry of the turn-for instance, when the turns between the alpha-helix and the beta-strands are very long. It also occurs much less frequently in flat-sheet alpha-beta proteins, where the topology is much less regular and the amount of twist on the sheet varies considerably more than in the barrel proteins. The motif may be one of the basic building blocks from which alpha-beta barrels are constructed.     

Functional consequences of mutations at the allosteric interface in hetero- and homo-hemoglobin tetramers.

A seminal difference exists between the two types of chains that constitute the tetrameric hemoglobin in vertebrates. While alpha chains associate weakly into dimers, beta chains self-associate into tightly assembled tetramers. While heterotetramers bind ligands cooperatively with moderate affinity, homotetramers bind ligands with high affinity and without cooperativity. These characteristics lead to the conclusion that the beta 4 tetramer is frozen in a quaternary R-state resembling that of liganded HbA. X-ray diffraction studies of the liganded beta 4 tetramers and molecular modeling calculations revealed several differences relative to the native heterotetramer at the "allosteric" interface (alpha 1 beta 2 in HbA) and possibly at the origin of a large instability of the hypothetical deoxy T-state of the beta 4 tetramer. We have studied natural and artificial Hb mutants at different sites in the beta chains responsible for the T-state conformation in deoxy HbA with the view of restoring a low ligand affinity with heme-heme interaction in homotetramers. Functional studies have been performed for oxygen equilibrium binding and kinetics after flash photolysis of CO for both hetero- and homotetramers. Our conclusion is that the "allosteric" interface is so precisely tailored for maintaining the assembly between alpha beta dimers that any change in the side chains of beta 40 (C6), beta 99 (G1), and beta 101 (G3) involved in the interface results in increased R-state behavior. In the homotetramer, the mutations at these sites lead to the destabilization of the beta 4 hemoglobin and the formation of lower affinity noncooperative monomers.     

Dimerization of beta B2-crystallin: the role of the linker peptide and the N- and C-terminal extensions.

beta B2- and gamma B-crystallins of vertebrate eye lens are 2-domain proteins in which each domain consists of 2 Greek key motifs connected by a linker peptide. Although the folding topologies of beta B2- and gamma B-domains are very similar, gamma B-crystallin is always monomeric, whereas beta B2-crystallin associates to homodimers. It has been suggested that the linker or the protruding N- and C-terminal arms of beta B2-crystallin (not present in gamma B) are a necessary requirement for this association. In order to investigate the role of these segments for dimerization, we constructed two beta B2 mutants. In the first mutant, the linker peptide was replaced with the one from gamma B (beta B2 gamma L). In the second mutant, the N- and C-terminal arms of 15- and 12-residues length were deleted (beta B2 delta NC). The beta B2 gamma L mutant is monomeric, whereas the beta B2 delta NC mutant forms dimers and tetramers that cannot be interconverted without denaturation. The spectral properties of the beta B2 mutants, as well as their stabilities against denaturants, resemble those of wild-type beta B2-crystallin, thus indicating that the overall peptide fold of the subunits is not changed significantly. We conclude that the peptide linker in beta B2-crystallin is necessary for dimerization, whereas the N- and C-terminal arms appear to be involved in preventing the formation of higher homo-oligomers.     

Minimal surface as a model of beta-sheets

The purpose of this article is to present arguments based on experimental data that the beta-sheet structures in proteins are the result of the tendency to minimize surface areas. Thus, we propose the model that all beta-sheet structures are almost minimal surfaces, namely, their mean curvatures are nearly zero. To support this model, we chose 1740 disjoint beta-sheets with less than 10 strands from the all beta-protein class in a nonredundant 40% Structural Classification of Proteins (SCOP) database and applied the least-squares method to fit the minimal surface catenoid (and in some rare cases, the plane) to the beta-sheet structures. The fitting errors were extremely small: The error of 1729 beta-sheets with catenoid minimal surface is 0.90 +/- 0.55 A and the error of the remaining 11 flat sheets with the plane is 0.64 +/- 0.46 A. The fact that the commonly used models for some beta-sheet surfaces (i.e., the hyperboloid and strophoid) have very small mean curvatures (< 0.05) supports our model. Moreover, we showed that this model also includes the isotropically stressed configuration model proposed by Salemme, in which the intrastrand tendency of the individual chains to twist or coil is in equilibrium with the tendency of the interstrand hydrogen bonding to resist twisting of the sheet as a whole. As an application we used our model to quantify the two principal independent modes in the flexibility of beta-sheets, that is, the bending parameter of beta-sheets and the inclined angle of beta-strands in a sheet.     

Guanidine hydrochloride unfolding of a transmembrane beta-strand in FepA using site-directed spin labeling.

We have used the electron spin resonance (ESR) site-directed spin-labeling (SDSL) technique to examine the guanidine hydrochloride (Gdn-HCl) induced denaturation of several sites along a transmembrane beta-strand located in the ferric enterobactin receptor, FepA. In addition, we have continued the characterization of the beta-strand previously identified by our group (Klug CS et al., 1997, Biochemistry 36:13027-13033) to extend from the periplasm to the extracellular surface loop in FepA, an integral membrane protein containing a beta-barrel motif comprised of a series of antiparallel beta-strands that is responsible for transport of the iron chelate, ferric enterobactin (FeEnt), across the outer membrane of Escherichia coli and many related enteric bacteria. We have previously shown that a large surface loop in FepA containing the FeEnt binding site denatures independently of the beta-barrel domain (Klug CS et al., 1995, Biochemistry 34:14230-14236). The SDSL approach allows examination of the unfolding at individual residues independent of the global unfolding of the protein. This work shows that sites along the beta-strand that are exposed to the aqueous lumen of the channel denature more rapidly and with higher cooperativity than the surface loop, while sites on the hydrophobic side of the beta-strand undergo a limited degree of noncooperative unfolding and do not fully denature even at high (e.g., 4 M) Gdn-HCl concentrations. We conclude that, in a transmembrane beta-strand, the local environment of a given residue plays a significant role in the loss of structure at each site.     

Identification,classification, and analysis of beta-bulges in proteins.

A beta-bulge is a region of irregularity in a beta-sheet involving two beta-strands. It usually involves two or more residues in the bulged strand opposite to a single residue on the adjacent strand. These irregularities in beta-sheets were identified and classified automatically, extending the definition of beta-bulges given by Richardson et al. (Richardson, J.S., Getzoff, E.D., & Richardson, D.C., 1978, Proc. Natl. Acad. Sci. USA 75, 2574-2578). A set of 182 protein chains (170 proteins) was used, and a total of 362 bulges were extracted. Five types of beta-bulges were found: classic, G1, wide, bent, and special. Their characteristic amino acid preferences were found for most classes of bulges. Basically, bulges occur frequently in proteins; on average there are more than two bulges per protein. In general, beta-bulges produce two main changes in the structure of a beta-sheet: (1) disrupt the normal alternation of side-chain direction; (2) accentuate the twist of the sheet, altering the direction of the surrounding strands.     

Metal ion-binding ability of tetrapeptides containing alpha-aminoisobutyric acid.

alpha-Aminoisobutyric acid (Aib), one of the Calpha,alpha-disubstituted glycines, is a sterically hindered amino acid that acts as a conformational constraint in peptides. However, studies for the application of the ability of Aib to control conformation are quite few. The paper focuses on the molecular recognition ability of acyclic oligopeptides containing Aib. Liquid-liquid extraction of nine kinds of metal ions from aqueous layers to nonpolar organic layers with acyclic tetrapeptides, X-Trp-Xaa2-Gly-Xaa4-NH-Ar (X = H or C6H5CH2OCO (Z), Xaa2 = Aib or Gly, Xaa4 = Leu or Ala, Ar = phenyl or 3,5-dimethylphenyl) was examined using picrate as the anion of ion pairs. The extraction behaviour of the metal ions with the tetrapeptides was investigated in the pH range from 3 to 9. In the case of basic pH regions, Cu(II) and Ag(I) were effectively extracted with Trp-Aib-Gly-Leu-NH-Ar. Pd(II) was specifically extracted with Trp-Aib-Gly-Leu-NH-Ar in acidic pH regions. The extraction percent (%E) of the peptide host, which has a 3,5-dimethylphenyl group, was even larger than that of the host, which has a phenyl group. Moreover, Pd(II) was extracted with a peptide host which has Leu and a 3,5-dimethylphenyl group in the absence of picrate as the anion of ion pairs. The free alpha-amino group, the turn conformation and the hydrophobicity of peptide molecules were important factors for the extraction of the metals.     

Recurrent alpha beta loop structures in TIM barrel motifs show a distinct pattern of conserved structural features.

A systematic survey of seven parallel alpha/beta barrel protein domains, based on exhaustive structural comparisons, reveals that a sizable proportion of the alpha beta loops in these proteins--20 out of a total of 49--belong to either one of two loop types previously described by Thornton and co-workers. Six loops are of the alpha beta 1 type, with one residue between the alpha-helix and beta-strand, and 13 are of the alpha beta 3 type, with three residues between the helix and the strand. Protein fragments embedding the identified loops, and termed alpha beta connections since they contain parts of the flanking helix and strand, have been analyzed in detail revealing that each type of connection has a distinct set of conserved structural features. The orientation of the beta-strand relative to the helix and loop portions is different owing to a very localized difference in backbone conformation. In alpha beta 1 connections, the chain enters the beta-strand via a residue adopting an extended conformation, while in alpha beta 3 it does so via a residue in a near alpha-helical conformation. Other conserved structural features include distinct patterns of side chain orientation relative to the beta-sheet surface and of main chain H-bonds in the loop and the beta-strand moieties. Significant differences also occur in packing interactions of conserved hydrophobic residues situated in the last turn of the helix. Yet the alpha-helix surface of both types of connections adopts similar orientations relative to the barrel sheet surface. Our results suggest furthermore that conserved hydrophobic residues along the sequence of the connections, may be correlated more with specific patterns of interactions made with neighboring helices and sheet strands than with helix/strand packing within the connection itself. A number of intriguing observations are also made on the distribution of the identified alpha beta 1 and alpha beta 3 loops within the alpha/beta-barrel motifs. They often occur adjacent to each other; alpha beta 3 loops invariably involve even numbered beta-strands, while alpha beta 1 loops involve preferentially odd beta-strands; all the analyzed proteins contain at least one alpha beta 3 loop in the first half of the eightfold alpha/beta barrel. Possible origins of all these observations, and their relevance to the stability and folding of parallel alpha/beta barrel motifs are discussed.     

Constraints imposed by protease accessibility on the trans-membrane and surface topography of the colicin E1 ion channel.

The surface topography of a 190-residue COOH-terminal colicin E1 channel peptide (NH2-Met 333-Ile 522-COOH) bound to uniformly sized 0.2-micron liposomes was probed by accessibility of the peptide to proteases in order (1) to determine whether the channel structure contains trans-membrane segments in addition to the four alpha-helices previously identified and (2) to discriminate between different topographical possibilities for the surface-bound state. An unfolded surface-bound state is indicated by increased trypsin susceptibility of the bound peptide relative to that of the peptide in aqueous solution. The peptide is bound tightly to the membrane surface with Kd < 10(-7) M. The NH2-terminal 50 residues of the membrane-bound peptide are unbound or loosely bound as indicated by their accessibility to proteases, in contrast with the COOH-terminal 140 residues, which are almost protease inaccessible. The general protease accessibility of the NH2-terminal segment Ala 336-Lys 382 excludes any model for the closed channel state that would include trans-membrane helices on the NH2-terminal side of Lys 382. Lys 381-Lys 382 is a major site for protease cleavage of the surface-bound channel peptide. A site for proteinase K cleavage just upstream of the amphiphilic gating hairpin (K420-K461) implies the presence of a surface-exposed segment in this region. These protease accessibility data indicate that it is unlikely that there are any alpha-helices on the NH2-terminal side of the gating hairpin K420-K461 that are inserted into the membrane in the absence of a membrane potential. A model for the topography of an unfolded monomeric surface-bound intermediate of the colicin channel domain, including a trans-membrane hydrophobic helical hairpin and two or three long surface-bound helices, is proposed.     

Comparison of straight chain and cyclic unnatural amino acids embedded in the core of staphylococcal nuclease.

We have determined by X-ray crystallography the structures of several variants of staphylococcal nuclease with long flexible straight chain and equivalent length cyclic unnatural amino acid side chains embedded in the protein core. The terminal atoms in the straight side chains are not well defined by the observed electron density even though they remain buried within the protein interior. We have previously observed this behavior and have suggested that it may arise from the addition of side-chain vibrational and oscillational motions with each bond as a side chain grows away from the relatively rigid protein main chain and/or the population of multiple rotamers (Wynn R, Harkins P, Richards FM. Fox RO. 1996. Mobile unnatural amino acid side chains in the core of staphylococcal nuclease. Protein Sci 5:1026-1031). Reduction of the number of degrees of freedom by cyclization of a side chain would be expected to constrain these motions. These side chains are in fact well defined in the structures described here. Over-packing of the protein core results in a 1.0 A shift of helix 1 away from the site of mutation. Additionally, we have determined the structure of a side chain containing a single hydrogen to fluorine atom replacement on a methyl group. A fluorine atom is intermediate in size between methyl group and a hydrogen atom. The fluorine atom is observed in a single position indicating it does not rotate like methyl hydrogen atoms. This change also causes subtle differences in the packing interactions.     

Probing the alpha-complementing domain of E. coli beta-galactosidase with use of an insertional pentapeptide mutagenesis strategy based on Mu in vitro DNA transposition

Protein structure-function relationships can be studied by using linker insertion mutagenesis, which efficiently identifies essential regions in target proteins. Bacteriophage Mu in vitro DNA transposition was used to generate an extensive library of pentapeptide insertion mutants within the alpha-complementing domain 1 of Escherichia coli beta-galactosidase, yielding mutants at 100% efficiency. Each mutant contained an accurate 15-bp insertion that translated to five additional amino acids within the protein, and the insertions were distributed essentially randomly along the target sequence. Individual mutants (alpha-donors) were analyzed for their ability to restore (by alpha-complementation) beta-galactosidase activity of the M15 deletion mutant (alpha-acceptor), and the data were correlated to the structure of the beta-galactosidase tetramer. Most of the insertions were well tolerated, including many of those disrupting secondary structural elements even within the protein's interior. Nevertheless, certain sites were sensitive to mutations, indicating both known and previously unknown regions of functional importance. Inhibitory insertions within the N-terminus and loop regions most likely influenced protein tetramerization via direct local effects on protein-protein interactions. Within the domain 1 core, the insertions probably caused either lateral shifting of the polypeptide chain toward the protein's exterior or produced more pronounced structural distortions. Six percent of the mutant proteins exhibited temperature sensitivity, in general suggesting the method's usefulness for generation of conditional phenotypes. The method should be applicable to any cloned protein-encoding gene.     

Analysis of hydrophobicity in the alpha and beta chemokine families and its relevance to dimerization.

The chemokine family of chemotactic cytokines plays a key role in orchestrating the immune response. The family has been divided into 2 subfamilies, alpha and beta, based on the spacing of the first 2 cysteine residues, function, and chromosomal location. Members within each subfamily have 25-70% sequence identity, whereas the amino acid identity between members of the 2 subfamilies ranges from 20 to 40%. A quantitative analysis of the hydrophobic properties of 11 alpha and 9 beta chemokine sequences, based on the coordinates of the prototypic alpha and beta chemokines, interleukin-8 (IL-8), and human macrophage inflammatory protein-1 beta (hMIP-1 beta), respectively, is presented. The monomers of the alpha and beta chemokines have their strongest core hydrophobic cluster at equivalent positions, consistent with their similar tertiary structures. In contrast, the pattern of monomer surface hydrophobicity between the alpha and beta chemokines differs in a manner that is fully consistent with the observed differences in quaternary structure. The most hydrophobic surface clusters on the monomer subunits are located in very different regions of the alpha and beta chemokines and comprise in each case the amino acids that are buried at the interface of their respective dimers. The theoretical analysis of hydrophobicity strongly supports the hypothesis that the distinct dimers observed for IL-8 and hMIP-1 beta are preserved for all the alpha and beta chemokines, respectively. This provides a rational explanation for the lack of receptor crossbinding and reactivity between the alpha and beta chemokine subfamilies.     

Knowledge-based model of a glucosyltransferase from the oral bacterial group of mutans streptococci.

Mutans streptococci glucosyltransferases catalyze glucosyl transfer from sucrose to a glucan chain. We previously identified an aspartyl residue that participates in stabilizing the glucosyl transition state. The sequence surrounding the aspartate was found to have substantial sequence similarity with members of alpha-amylase family. Because little is known of the protein structure beyond the amino acid sequence, we used a knowledge-based interactive algorithm, MACAW, which provided significant level of homology with alpha-amylases and glucosyltransferase from Streptococcus downei gtfI (GTF). The significance of GTF similarity is underlined by GTF/alpha-amylase residues conserved in all but one alpha-amylase invariant residues. Site-directed mutagenesis of the three GTF catalytic residues are homologous with the alpha-amylase catalytic triad. The glucosyltransferases are members of the 4/7-superfamily that have a (beta/alpha)8-barrel structure and belong to family 13 of the glycohydralases.     

Catalytic function and affinity purification of site-directed mutant β-cyclodextrin glucanotransferase from alkalophilic Bacillus firmus var. alkalophilus

Subsites −3 and −7 in the active site of β-cyclodextrin glucanotransferase (β-CGTase) from alkalophilic Bacillus firmus var. alkalophilus were modified through site-directed mutagenesis to obtain novel mutant CGTases. Four mutant CGTases, H59Q, Y96M, 90-PPI-92, and Δ(154–160) were constructed and produced using a recombinant E. coli with a secretive expression system extracellularly. The secreted mutant β-CGTases were purified by one-step affinity adsorption chromatography using a β-cyclodextrin (CD) polymer as an adsorbent to nearly homogeneous purity. The catalytic activities were modified significantly compared to the wild-type. In particular, the Y96M and Δ(154–160) mutants increased cyclization activity around 1.5 times without any significant reduction of coupling and hydrolyzing activities. Meanwhile, the Y96M and Δ(154–160) mutants exhibited a much higher conversion yield into CDs from 28.6 to 39% without any recognizable change in the CD ratio. The conversion yield into linear maltooligosaccharides was also significantly reduced. The catalytic functions of subsites −3 and −7 in the active site of β-CGTase would appear to be related to the overall productivity of CDs rather than the product specificity.     

The effect of long-range interactions on the secondary structure formation of proteins

The influence of long-range residue interactions on defining secondary structure in a protein has long been discussed and is often cited as the current limitation to accurate secondary structure prediction. There are several experimental examples where a local sequence alone is not sufficient to determine its secondary structure, but a comprehensive survey on a large data set has not yet been done. Interestingly, some earlier studies denied the negative effect of long-range interactions on secondary structure prediction accuracy. Here, we have introduced the residue contact order (RCO), which directly indicates the separation of contacting residues in terms of the position in the sequence, and examined the relationship between the RCO and the prediction accuracy. A large data set of 2777 nonhomologous proteins was used in our analysis. Unlike previous studies, we do find that prediction accuracy drops as residues have contacts with more distant residues. Moreover, this negative correlation between the RCO and the prediction accuracy was found not only for beta-strands, but also for alpha-helices. The prediction accuracy of beta-strands is lower if residues have a high RCO or a low RCO, which corresponds to the situation that a beta-sheet is formed by beta-strands from different chains in a protein complex. The reason why the current study draws the opposite conclusion from the previous studies is examined. The implication for protein folding is also discussed.