Identification of members of the Bex gene family as olfactory marker protein (OMP) binding partners

Olfactory marker protein (OMP) expression is a hallmark of mature vertebrate olfactory receptor neurons (ORNs). Evidence for OMP function derives from altered behavioral and electrophysiological activities of OMP-KO mice. The molecular basis for the altered phenotype following the deletion of OMP is still unclear. Recent structural studies predict the involvement of OMP in protein-protein interaction. Here we report the identification of an OMP partner, Bex2, by phage-display screening of an olfactory mucosal cDNA-library. In situ hybridization demonstrates cellular co-localization of OMP mRNA with mRNAs for Bex1, Bex2, and Bex3 in ORNs of olfactory tissue of the mouse. The OMP/Bex interaction has been confirmed by demonstrating the chemical cross-linking of recombinant rat OMP with a synthetic peptide derived from the Bex amino acid sequence. The subcellular localization of Bex and OMP proteins was evaluated in transfected HEK293 cells. Bex is visualized in the nucleus and cytoplasm. Following co-transfection we observed the unexpected presence of some OMP in the nucleus along with Bex. Together, these data argue convincingly that we have identified Bex as an OMP partner whose further characterization will provide insight to the role of OMP and to the mechanism of the OMP/Bex interaction in ORN differentiation and function.     

The interaction of Bex and OMP reveals a dimer of OMP with a short half-life

Olfactory marker protein (OMP) participates in the olfactory signal transduction pathway. This is evident from the behavioral and electrophysiological deficits of OMP-null mice, which can be reversed by intranasal infection of olfactory sensory neurons with an OMP-expressing adenovirus. Bex, brain expressed X-linked protein, has been identified as a protein that interacts with OMP. We have now further characterized the interaction of OMP and Bex1/2 by in vitro binding assays and by immuno-coprecipitation experiments. OMP is a 19 kDa protein but these immunoprecipitation studies have revealed the unexpected presence of a 38 kDa band in addition to the expected 19 kDa band. Furthermore, the 38 kDa form was preferentially co-immunoprecipitated with Bex from cell extracts. In-gel tryptic digestion, mass spectrometry, and two-dimensional gel electrophoresis indicate that the 38 kDa protein behaves as a covalently cross-linked OMP-homodimer. The 38 kDa band was also identified in western blots of olfactory epithelium demonstrating its presence in vivo. The stabilities and subcellular localizations of the OMP-monomer and -dimer were studied in transfected cells. These results demonstrated that the OMP-dimer is much less stable than the monomer, and that while the monomer is present both in the nuclear and cytosolic compartments, the dimer is preferentially located in a Triton X-100 insoluble cytoskeletal fraction. These novel observations led us to hypothesize that regulation of the level of the rapidly turning-over OMP-dimer and its interaction with Bex1/2 is critical for OMP function in sensory transduction.     

Dependence of alpha-helical and beta-sheet amino acid propensities on the overall protein fold type

ABSTRACT: BACKGROUND: A large number of studies have been carried out to obtain amino acid propensities for alpha- helices and beta-sheets. The obtained propensities for alpha-helices are consistent with each other, and the pair-wise correlation coefficient is frequently high. On the other hand, the beta-sheet propensities obtained by several studies differed significantly, indicating that the context significantly affects beta-sheet propensity. RESULTS: We calculated amino acid propensities for alpha-helices and beta-sheets for 39 and 24 protein folds, respectively, and addressed whether they correlate with the fold. The propensities were also calculated for exposed and buried sites, respectively. Results showed that alpha-helix propensities do not differ significantly by fold, but beta-sheet propensities are diverse and depend on the fold. The propensities calculated for exposed sites and buried sites are similar for alpha-helix, but such is not the case for the beta-sheet propensities. We also found some fold dependence on amino acid frequency in beta-strands. Folds with a high Ser, Thr and Asn content at exposed sites in beta-strands tend to have a low Leu, Ile, Glu, Lys and Arg content (correlation coefficient = 0.90) and to have flat beta-sheets. At buried sites in beta-strands, the content of Tyr, Trp, Gln and Ser correlates negatively with the content of Val, Ile and Leu (correlation coefficient = 0.93). "All-beta" proteins tend to have a higher content of Tyr, Trp, Gln and Ser, whereas alpha/beta proteins tend to have a higher content of Val, Ile and Leu. CONCLUSIONS: The alpha-helix propensities are similar for all folds and for exposed and buried residues. However, beta-sheet propensities calculated for exposed residues differ from those for buried residues, indicating that the exposed-residue fraction is one of the major factors governing amino acid composition in beta-strands. Furthermore, the correlations we detected suggest that amino acid composition is related to folding properties such as the twist of a beta-strand or association between two beta sheets.     

Solution structure of the coxsackievirus and adenovirus receptor domain 2

The coxsackievirus and adenovirus receptor (CAR) mediates entry of coxsackievirus and adenovirus. CAR possesses an extracellular region that is comprised of 2 immunoglobulin domains termed CAR-D1 and CAR-D2. In the present work, the solution structure of CAR-D2, consisting of residues 142-235 of human CAR, has been determined by NMR spectroscopy. CAR-D2 is shown to be a beta-sandwich motif comprised of two beta-sheets, which are stabilized by two disulfide bonds. The first beta-sheet is comprised of beta-strands A, B, and E, and the second beta-sheet is comprised of beta-strands C, F, and G. A relatively hydrophobic helix is found between beta-strands C and E, which replaces beta-strand D of the classical c-type immunoglobulin fold.     

Calorimetric dissection of thermal unfolding of OspA,a predominantly beta-sheet protein containing a single-layer beta-sheet

Outer surface protein A (OspA) from Borrelia burgdorferi is a predominantly beta-sheet protein comprised of beta-strands beta1-beta21 and a short C-terminal alpha-helix. It contains two globular domains (N and C-terminal domains) and a unique single-layer beta-sheet (central beta-sheet) that connects the two domains. OspA contains an unusually large number of charged amino acid residues. To understand the mechanism of stabilization of this unique beta-sheet protein, thorough thermodynamic investigations of OspA and its truncated mutant lacking a part of the C-terminal domain were conducted using calorimetry and circular dichroism. The stability of OspA was found to be sensitive to pH and salt concentration. The heat capacity curve clearly consisted of two components, and all the thermodynamic parameters were obtained for each step. The thermodynamic parameters associated with the two transitions are consistent with a previously proposed model, in which the first transition corresponds to the unfolding of the C-terminal domain and the last two beta-strands of the central beta-sheet, and the second transition corresponds to that of the N-terminal domain and the first beta-strand of the central beta-sheet in the second peak. The ratio of calorimetric and van't Hoff enthalpies indicates that the first peak includes another thermodynamic intermediate state. Large heat capacity changes were observed for both transitions, indicative of large changes in the exposure of hydrophobic surfaces associated with the transitions. This observation demonstrates that hydrophobic parts are buried efficiently in the native structure in spite of the low content of hydrophobic residues in OspA. By decomposing the enthalpy, entropy, and Gibbs free energy into contributions from different interactions, we found that the enthalpy changes for hydrogen bonding and polar interactions are exceptionally large, indicating that OspA maintains its stability by making full use of its unique beta-sheet and high content of polar residues. These thermodynamic analyses demonstrated that it is possible to maintain protein tertiary structure by making effective use of an unusual amino acid composition.     

Dichroic ratios in polarized Fourier transform infrared for nonaxial symmetry of beta-sheet structures

The transition moments for the amide bands from beta-sheet peptide structures generally do not exhibit axial symmetry about the director in linearly polarized Fourier transform infrared (FTIR) measurements on oriented systems. The angular dependences of the dichroic ratios of the amide bands are derived for beta-sheet structures in attenuated total reflection (ATR) and polarized transmission experiments on samples that are oriented with respect to the normal to the substrate and are randomly distributed with respect to the azimuthal angle in the plane of the orienting substrate. The orientational distributions of both the beta-strands and the beta-sheets are considered, and explicit expressions are given for the dichroic ratios of the amide I and amide II bands. The dichroic ratio of the amide II band, which is parallel polarized, can yield the orientation of the beta-strands directly, but to specify the orientations of the beta-sheets completely requires measurement of the dichroic ratios of both the amide I and amide II bands, or generally two bands with parallel and perpendicular polarizations. A random distribution in tilt of the planes of the beta-sheets does not give rise to equal dichroic ratios for bands with perpendicular and parallel polarizations, such as the amide I and amide II bands. The results are applied to previous ATR and polarized transmission FTIR measurements on a potassium channel-associated peptide, the Escherichia coli outer membrane protein OmpA, and the E. coli OmpF porin protein in oriented membranes.     

Solution structure and antibody binding studies of the envelope protein domain III from the New York strain of West Nile virus

The solution structure of domain III from the New York West Nile virus strain 385-99 (WN-rED3) has been determined by NMR methods. The West Nile domain III structure is a beta-barrel structure formed from seven anti-parallel beta-strands in two beta-sheets. One anti-parallel beta-sheet consists of beta-strands beta1 (Phe(299)-Asp(307)), beta2 (Val(313)-Tyr(319)), beta4 (Arg(354)-Leu(355)), and beta5 (Lys(370)-Glu(376)) arranged so that beta2 is flanked on either side by beta1 and beta5. The short beta4 flanks the end of the remaining side of beta5. The remaining anti-parallel beta-sheet is formed from strands beta3 (Ile(340)-Val(343)), beta6 (Gly(380)-Arg(388)), and beta7 (Gln(391)-Lys(399)) arranged with beta6 at the center. Residues implicated in antigenic differences between different West Nile virus strains (and other flaviviruses) and neutralization are located on the outer surface of the protein. Characterization of the binding of monoclonal antibodies to WN-rED3 mutants, which were identified through neutralization escape experiments, indicate that antibody neutralization directly correlates with binding affinities. These studies provide an insight into theoretical virus-receptor interaction points, structure of immunogenic determinants, and potential targets for antiviral agents against West Nile virus and highlight differences between West Nile virus and other flavivirus structures that may represent critical determinants of virulence.     

Structural plasticity associated with the beta-propeller architecture

The beta-propeller architecture observed in protein tertiary structure and classified into the five different types according to number of 'blades' (or beta-sheets) and a sixth type classified according to the secondary structure composition of the blades (the beta beta alpha beta-molecular unit) is characterized by variations (or plasticity) in the structure. These correspond to the number of beta-strands associated with the blade, the number of amino acid residues associated with equivalent beta-strands in the different blades and the presence of alpha-helices and twisted beta-strands. We have generated a beta-sheet associated beta-strand pattern that may be important for protein structure prediction and modeling. Analysis of the beta-propellers extracted primarily from the SCOP database revealed there are 179 beta-propellers. The examination of the secondary structure corresponding to the beta-propeller using PDBsum that was useful to define the beta-sheet associated beta-strand pattern, combined with visualization on graphics display revealed structural plasticity associated with the beta-propeller architecture. Particularly, the type 6- and 7-bladed beta-propellers known to be associated with sequence and functional diversity are more common and associated with relatively more structural variations compared to the other beta-propeller types.     

Microscopic stability of cold shock protein A examined by NMR native state hydrogen exchange as a function of urea and trimethylamine N-oxide

Native state hydrogen exchange of cold shock protein A (CspA) has been characterized as a function of the denaturant urea and of the stabilizing agent trimethylamine N-oxide (TMAO). The structure of CspA has five strands of beta-sheet. Strands beta1-beta4 have strongly protected amide protons that, based on experiments as a function of urea, exchange through a simple all-or-none global unfolding mechanism. By contrast, the protection of amide protons from strand beta5 is too weak to measure in water. Strand beta5 is hydrogen bonded to strands beta3 and beta4, both of which afford strong protection from solvent exchange. Gaussian network model (GNM) simulations, which assume that the degree of protection depends on tertiary contact density in the native structure, accurately predict the strong protection observed in strands beta1-beta4 but fail to account for the weak protection in strand beta5. The most conspicuous feature of strand beta5 is its low sequence hydrophobicity. In the presence of TMAO, there is an increase in the protection of strands beta1-beta4, and protection extends to amide protons in more hydrophilic segments of the protein, including strand beta5 and the loops connecting the beta-strands. TMAO stabilizes proteins by raising the free energy of the denatured state, due to highly unfavorable interactions between TMAO and the exposed peptide backbone. As such, the stabilizing effects of TMAO are expected to be relatively independent of sequence hydrophobicity. The present results suggest that the magnitude of solvent exchange protection depends more on solvent accessibility in the ensemble of exchange susceptible conformations than on the strength of hydrogen-bonding interactions in the native structure.     

On the interaction of ribosomal protein L5 with 5S rRNA

Ribosomal protein L5, a 5S rRNA binding protein in the large subunit, is composed of a five-stranded antiparallel beta-sheet and four alpha-helices, and folds in a way that is topologically similar to the ribonucleprotein (RNP) domain [Nakashima et al., RNA 7, 692-701, 20011. The crystal structure of ribosomal protein L5 (BstL5) from Bacillus stearothermophilus suggests that a concave surface formed by an anti-parallel beta-sheet and long loop structures are strongly involved in 5S rRNA binding. To identify amino acid residues responsible for 5S rRNA binding, we made use of Ala-scanning mutagenesis of evolutionarily conserved amino acids occurred at beta-strands and loop structures in BstL5. The mutation of Lys33 at the beta 1-strand caused a significant reduction in 5S rRNA binding. In addition, the Arg92, Phe122, and Glu134 mutations on the beta2-strand, the alpha3-beta4 loop, and the beta4-beta5 loop, respectively, resulted in a moderate decrease in the 5S rRNA binding affinity. In contrast, mutation of the conserved residue Pro65 at the beta2-strand had little effect on the 5S rRNA binding activity. These results, taken together with previous results, identified Lys33, Asn37, Gln63, and Thr90 on the beta-sheet structure, and Phe77 at the beta2-beta3 loop as critical residues for the 5S rRNA binding. The contribution of these amino acids to 5S rRNA binding was further quantitatively evaluated by surface plasmon resonance (SPR) analysis by the use of BIAcore. The results showed that the amino acids on the beta-sheet structure are required to decrease the dissociation rate constant for the BstL5-5S rRNA complex, while those on the loops are to increase the association rate constant for the BstL5-5S rRNA interaction.     

Ca2+ Extrusion by NCX Is Compromised in Olfactory Sensory Neurons of OMP?/? Mice

Background

The role of olfactory marker protein (OMP), a hallmark of mature olfactory sensory neurons (OSNs), has been poorly understood since its discovery. The electrophysiological and behavioral phenotypes of OMP knockout mice indicated that OMP influences olfactory signal transduction. However, the mechanism by which this occurs remained unknown.

Principal Findings

We used intact olfactory epithelium obtained from WT and OMP^−/− mice to monitor the Ca²⁺ dynamics induced by the activation of cyclic nucleotide-gated channels, voltage-operated Ca²⁺ channels, or Ca²⁺ stores in single dendritic knobs of OSNs. Our data suggested that OMP could act to modulate the Ca²⁺-homeostasis in these neurons by influencing the activity of the plasma membrane Na⁺/Ca²⁺-exchanger (NCX). Immunohistochemistry verifies colocalization of NCX1 and OMP in the cilia and knobs of OSNs. To test the role of NCX activity, we compared the kinetics of Ca²⁺ elevation by stimulating the reverse mode of NCX in both WT and OMP^−/− mice. The resulting Ca²⁺ responses indicate that OMP facilitates NCX activity and allows rapid Ca²⁺ extrusion from OSN knobs. To address the mechanism by which OMP influences NCX activity in OSNs we studied protein-peptide interactions in real-time using surface plasmon resonance technology. We demonstrate the direct interaction of the XIP regulatory-peptide of NCX with calmodulin (CaM).

Conclusions

Since CaM also binds to the Bex protein, an interacting protein partner of OMP, these observations strongly suggest that OMP can influence CaM efficacy and thus alters NCX activity by a series of protein-protein interactions.     

NMR structure of the conserved hypothetical protein TM0487 from Thermotoga maritima: implications for 216 homologous DUF59 proteins

The NMR structure of the conserved hypothetical protein TM0487 from Thermotoga maritima represents an alpha/beta-topology formed by the regular secondary structures alpha1-beta1-beta2-alpha2-beta3-beta4-alpha3- beta5-3(10)-alpha4, with a small anti-parallel beta-sheet of beta-strands 1 and 2, and a mixed parallel/anti-parallel beta-sheet of beta-strands 3-5. Similar folds have previously been observed in other proteins, with amino acid sequence identity as low as 3% and a variety of different functions. There are also 216 sequence homologs of TM0487, which all have the signature sequence of domains of unknown function 59 (DUF59), for which no three-dimensional structures have as yet been reported. The TM0487 structure thus presents a platform for homology modeling of this large group of DUF59 proteins. Conserved among most of the DUF59s are 13 hydrophobic residues, which are clustered in the core of TM0487. A putative active site of TM0487 consisting of residues D20, E22, L23, T51, T52, and C55 is conserved in 98 of the 216 DUF59 sequences. Asp20 is buried within the proposed active site without any compensating positive charge, which suggests that its pK(a) value may be perturbed. Furthermore, the DUF59 family includes ORFs that are part of a conserved chromosomal group of proteins predicted to be involved in Fe-S cluster metabolism.     

A general model of the P2 protein of peripheral nervous system myelin based on secondary structure predictions, tertiary folding principles, and experimental observations

The amino acid sequence of the P2 protein of peripheral myelin was analyzed with regard to regions of probable alpha-helix, beta-structure, beta-turn, and unordered conformation by means of several algorithms commonly used to predict secondary structure in proteins. Because of the high beta-sheet content and virtual absence of alpha-helix shown by the circular dichroic spectra of the protein, a bias was introduced into the algorithms to favor the beta-structure over the alpha-helical conformation. In order to define those beta-sheet residues that could lie on the external hydrophilic surface of the protein and those that could lie in its hydrophobic interior, the predicted beta-strands were examined for charged and uncharged amino acids located at alternating positions in the sequence. The sequential beta-strands in the predicted secondary structure were then ordered into beta-sheets and aligned according to generally accepted tertiary folding principles and certain chemical properties peculiar to the P2 protein. The general model of the P2 protein that emerged was a "Greek key" beta-barrel, consisting of eight antiparallel beta-strands with a two-stranded ribbon of antiparallel beta-structure emerging from one end. The model has an uncharged, hydrophobic core and a highly hydrophilic surface. The two Cys residues, which form a disulfide, occur in a loop connecting two adjacent antiparallel strands. Two hydrophilic loops, each containing a cluster of acidic residues and a single Phe, protrude from one end of the molecule. The general model is consistent with many of the properties of the actual protein, including the relatively weak nature of its association with myelin lipids and the positions of amino acid substitutions. Alternative beta-strand orderings yield three specific models having different interstrand connections across the barrel ends.     

Twist and shear in beta-sheets and beta-ribbons

The structures of the beta-sheets and the beta-ribbons have been analysed using high-resolution protein structure data. Systematic asymmetries measured in both parallel and antiparallel beta-structures include the sheet twist and the strand shear. In order to determine the origin of these asymmetries, numerous interactions and correlations were examined. The strongest correlations are observed for residues in antiparallel beta-sheets and beta-ribbons that form non-H-bonded pairs. For these residues, the sheet twist is correlated to the backbone phi angle but not to the psi angle. Our analysis supports the existence of an inter-strand C(alpha)H(alpha)...O weak H-bond, which, together with the CO...HN H-bond, constitutes a bifurcated H-bond that links neighbouring beta-strands. Residues of beta-sheets and beta-ribbons in high-resolution protein structures form a distinct region of the Ramachandran plot, which is determined by the formation of the bifurcated H-bond, the formation of an intra-strand O...H(alpha) non-bonded polar interaction, and an intra-strand O...C(beta) steric clash. Using beta-strands parameterised by phi-psi values from the allowed beta-sheet region of the Ramachandran plot, the shear and the right-hand twist can be reproduced in a simple model of the antiparallel and parallel beta-ribbon that models the bifurcated H-bonds specifically. The conformations of interior residues of beta-sheets are shown to be subsets of the conformations of residues of beta-ribbons.     

A cross-strand Trp Trp pair stabilizes the hPin1 WW domain at the expense of function

Using the human Pin1 WW domain (hPin1 WW), we show that replacement of two nearest neighbor non-hydrogen-bonded residues on adjacent beta-strands with tryptophan (Trp) residues increases beta-sheet thermodynamic stability by 4.8 kJ mol(-1) at physiological temperature. One-dimensional NMR studies confirmed that introduction of the Trp-Trp pair does not globally perturb the structure of the triple-stranded beta-sheet, while circular dichroism studies suggest that the engineered cross-strand Trp-Trp pair adopts a side-chain conformation similar to that first reported for a designed "Trp-zipper" beta-hairpin peptide, wherein the indole side chains stack perpendicular to each other. Even though the mutated side chains in wild-type hPin1 WW are not conserved among WW domains and compose the beta-sheet surface opposite to that responsible for ligand binding, introduction of the cross-strand Trp-Trp pair effectively eliminates hPin1 WW function as assessed by the loss of binding affinity toward a natural peptide ligand. Maximizing both thermodynamic stability and the domain function of hPin1 WW by the above mentioned approach appears to be difficult, analogous to the situation with loop 1 optimization explored previously. That introduction of a non-hydrogen-bonded cross-strand Trp-Trp pair within the hPin1 WW domain eliminates function may provide a rationale for why this energetically favorable pairwise interaction has not yet been identified in WW domains or any other biologically evolved protein with known three-dimensional structure.     

Mean curvature as a major determinant of beta-sheet propensity

MOTIVATION: Despite the importance of beta-sheets as building blocks in proteins and also toxic elements in the pathological disorders, ranging from Alzheimer's disease to mad cow disease, the principles underlying their stability are not well understood. Non-random beta-sheet propensities of amino acids have been revealed both by their distinct statistical preferences within known protein structures and by the relative thermodynamic scales through the experimental host-guest systems. However, recent fitting analysis has proved that a native beta-sheet conforms to a minimal surface with zero mean curvature, like the physical model of soap films. RESULTS: We here suggest that the stability of a residue in the all beta-sheet proteins can be measured with its mean curvature parameter, using discrete differential geometry. The sharply decreasing mean curvature with increasing number of beta-strands identifies a significant cooperative effect whereby the interstrand interaction increases in strength with the number of beta-strands. Furthermore, strong correlations of mean curvatures with previous beta-sheet propensities of amino acids show that their intrinsic differences in adopting the ideal beta-sheet structure are affected by the water-accessible area of side-chains, and result in the distinct statistical and thermodynamic beta-sheet propensities. Therefore, we conclude that mean curvature should be considered as the significant stability index of a beta-sheet structure.     

Structural similarity of YbeD protein from Escherichia coli to allosteric regulatory domains

Lipoic acid is an essential prosthetic group in several metabolic pathways. The biosynthetic pathway of protein lipoylation in Escherichia coli involves gene products of the lip operon. YbeD is a conserved bacterial protein located in the dacA-lipB intergenic region. Here, we report the nuclear magnetic resonance structure of YbeD from E. coli. The structure includes a beta alpha beta beta alpha beta fold with two alpha-helices on one side of a four-strand antiparallel beta-sheet. The beta 2-beta 3 loop shows the highest sequence conservation and is likely functionally important. The beta-sheet surface contains a patch of conserved hydrophobic residues, suggesting a role in protein-protein interactions. YbeD shows striking structural homology to the regulatory domain from d-3-phosphoglycerate dehydrogenase, hinting at a role in the allosteric regulation of lipoic acid biosynthesis or the glycine cleavage system.     

Correlations among morphology, beta-sheet stability, and molecular structure in prion peptide aggregates

The misfolding of proteins into beta-sheets and the subsequent aggregation of these sheets into fibrous networks underlies many diseases. In this paper, the role of peptide structure in determining the ordering of beta-sheet aggregates and the morphology of fibrils and protofibrils is dissected. Using a series of peptides based on residues 109-122 of the Syrian hamster prion protein (H1) with a range of substitutions at position 117, the link between side chain interactions and beta-sheet thermal stability has been investigated. The thermal stability of beta-sheets is associated with the peptides' ability to adopt the same alignment as wild-type H1, with residue 117 in register across all beta-strands [Silva, R. A. G. D., Barber-Armstrong, W., and Decatur, S. M. (2003) J. Am. Chem. Soc. 125, 13674-13675]. These aligned strands are capable of forming long, rigid, and twisted fibrils (as visualized by atomic force microscopy) which are thermostable. Peptides which do not adopt this strand alignment aggregate to form thin, flexible, and smooth protofibrils. The ability to form ordered aggregates, and thus to form twisted fibrils, is modulated by the structure of the side chain of residue 117.     

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.