Predicting transmembrane beta-barrels and interstrand residue interactions from sequence

Transmembrane beta-barrel (TMB) proteins are embedded in the outer membrane of Gram-negative bacteria, mitochondria, and chloroplasts. The cellular location and functional diversity of beta-barrel outer membrane proteins (omps) makes them an important protein class. At the present time, very few nonhomologous TMB structures have been determined by X-ray diffraction because of the experimental difficulty encountered in crystallizing transmembrane proteins. A novel method using pairwise interstrand residue statistical potentials derived from globular (nonouter membrane) proteins is introduced to predict the supersecondary structure of transmembrane beta-barrel proteins. The algorithm transFold employs a generalized hidden Markov model (i.e., multitape S-attribute grammar) to describe potential beta-barrel supersecondary structures and then computes by dynamic programming the minimum free energy beta-barrel structure. Hence, the approach can be viewed as a "wrapping" component that may capture folding processes with an initiation stage followed by progressive interaction of the sequence with the already-formed motifs. This approach differs significantly from others, which use traditional machine learning to solve this problem, because it does not require a training phase on known TMB structures and is the first to explicitly capture and predict long-range interactions. TransFold outperforms previous programs for predicting TMBs on smaller (相似文献   

Sequence motifs and antimotifs in beta-barrel membrane proteins from a genome-wide analysis: the Ala-Tyr dichotomy and chaperone binding motifs

Beta-barrel membrane proteins are found in the outer membrane of gram-negative bacteria, mitochondria, and chloroplasts. Although sequence motifs have been studied in alpha-helical membrane proteins and have been shown to play important roles in their assembly, it is not clear whether over-represented motifs and under-represented anti-motifs exist in beta-barrel membrane proteins. We have developed probabilistic models to identify sequence motifs of residue pairs on the same strand separated by an arbitrary number of residues. A rigorous statistical model is essential for this study because of the difficulty associated with the short length of the strands and the small amount of structural data. By comparing to the null model of exhaustive permutation of residues within the same beta-strand, propensity values of sequence patterns of two residues and p-values measuring statistical significance are calculated exactly by several analytical formulae we have developed or by enumeration. We find that there are characteristic sequence motifs and antimotifs in transmembrane (TM) beta-strands. The amino acid Tyr plays an important role in several such motifs. We find a general dichotomy consisting of favorable Aliphatic-Tyr sequence motifs and unfavorable Tyr-Aliphatic antimotifs. Tyr is also part of a terminal motif, YxF, which is likely to be important for chaperone binding. Our results also suggest several experiments that can help to elucidate the mechanisms of in vitro and in vivo folding of beta-barrel membrane proteins.     

Interplay between hydrophobic cluster and loop propensity in beta-hairpin formation

Autonomously folding beta-hairpins have recently emerged as powerful tools for elucidating the origins of antiparallel beta-sheet folding preferences. Analysis of such model systems has suggested four potential sources of beta-sheet stability: (1) the conformational propensity of the loop segment that connects adjacent strands; (2) favorable contacts between side-chains on adjacent strands; (3) interstrand hydrogen bonds; and (4) the intrinsic beta-sheet propensities of the strand residues. We describe the design and analysis of a series of isomeric 20 residue peptides in which factors (1)-(4) are identical. Differences in beta-hairpin formation within this series demonstrate that these four factors, individually, are not sufficient to explain beta-sheet stability. In agreement with the prediction of a simple statistical mechanical model for beta-hairpin formation, our results show that the separation between the loop segment and an interstrand cluster of hydrophobic side-chains strongly influences beta-hairpin size and stability, with a smaller separation leading to greater stability.     

Folding type-specific secondary structure propensities of amino acids,derived from α-helical, β-sheet, α/β, and α+β proteins of known structures

Folding type-specific secondary structure propensities of 20 naturally occurring amino acids have been derived from α-helical, β-sheet, α/β, and α+β proteins of known structures. These data show that each residue type of amino acids has intrinsic propensities in different regions of secondary structures for different folding types of proteins. Each of the folding types shows markedly different rank ordering, indicating folding type-specific effects on the secondary structure propensities of amino acids. Rigorous statistical tests have been made to validate the folding type-specific effects. It should be noted that α and β proteins have relatively small α-helices and β-strands forming propensities respectively compared with those of α+β and α/β proteins. This may suggest that, with more complex architectures than α and β proteins, α+β and α/β proteins require larger propensities to distinguish from interacting α-helices and β-strands. Our finding of folding type-specific secondary structure propensities suggests that sequence space accessible to each folding type may have differing features. Differing sequence space features might be constrained by topological requirement for each of the folding types. Almost all strong β-sheet forming residues are hydrophobic in character regardless of folding types, thus suggesting the hydrophobicities of side chains as a key determinant of β-sheet structures. In contrast, conformational entropy of side chains is a major determinant of the helical propensities of amino acids, although other interactions such as hydrophobicities and charged interactions cannot be neglected. These results will be helpful to protein design, class-based secondary structure prediction, and protein folding. © 1998 John Wiley & Sons, Inc. Biopoly 45: 35–49, 1998     

Correct folding of the beta-barrel of the human membrane protein VDAC requires a lipid bilayer

Spontaneous membrane insertion and folding of beta-barrel membrane proteins from an unfolded state into lipid bilayers has been shown previously only for few outer membrane proteins of Gram-negative bacteria. Here we investigated membrane insertion and folding of a human membrane protein, the isoform 1 of the voltage-dependent anion-selective channel (hVDAC1) of mitochondrial outer membranes. Two classes of transmembrane proteins with either alpha-helical or beta-barrel membrane domains are known from the solved high-resolution structures. VDAC forms a transmembrane beta-barrel with an additional N-terminal alpha-helix. We demonstrate that similar to bacterial OmpA, urea-unfolded hVDAC1 spontaneously inserts and folds into lipid bilayers upon denaturant dilution in the absence of folding assistants or energy sources like ATP. Recordings of the voltage-dependence of the single channel conductance confirmed folding of hVDAC1 to its active form. hVDAC1 developed first beta-sheet secondary structure in aqueous solution, while the alpha-helical structure was formed in the presence of lipid or detergent. In stark contrast to bacterial beta-barrel membrane proteins, hVDAC1 formed different structures in detergent micelles and phospholipid bilayers, with higher content of beta-sheet and lower content of alpha-helix when inserted and folded into lipid bilayers. Experiments with mixtures of lipid and detergent indicated that the content of beta-sheet secondary structure in hVDAC1 decreased at increased detergent content. Unlike bacterial beta-barrel membrane proteins, hVDAC1 was not stable even in mild detergents such as LDAO or dodecylmaltoside. Spontaneous folding of outer membrane proteins into lipid bilayers indicates that in cells, the main purpose of membrane-inserted or associated assembly factors may be to select and target beta-barrel membrane proteins towards the outer membrane instead of actively assembling them under consumption of energy as described for the translocons of cytoplasmic membranes.     

PA63 channel of anthrax toxin: an extended beta-barrel

Anthrax toxin consists of three protein components: protective antigen (PA), lethal factor (LF), and edema factor (EF). PA(63), generated by protease "nicking" of whole PA, is responsible for delivering the toxin's catalytic fragments (LF and EF) to the target cell's cytosol. In planar bilayer membranes, trypsin-nicked PA makes cation-selective voltage-gated channels with a pore diameter of > or =12 A. The channels are presumed to be heptameric "mushrooms", with an extracellular "cap" region and a membrane-inserted, beta-barrel "stem". Although the crystal structure of the water-soluble monomeric form has been resolved to 2.1 A and that of the heptameric "prepore" to 4.5 A, the structure for the membrane-bound channel (pore) has not been determined. We have engineered mutant channels that are cysteine-substituted in residues in the putative beta-barrel, and identified the residues lining the channel lumen by their accessibility to a water-soluble sulfhydryl-specific reagent. The reaction with lumen-exposed cysteinyl side chains causes a drop in channel conductance, which we used to map the residues that line the pore. Our results indicate that the beta-barrel structure extends beyond the bilayer and involves residues that are buried in the monomer. The implication is that major rearrangement of domains in the prepore cap region is required for membrane insertion of the beta-barrel stem.     

Empirical lipid propensities of amino acid residues in multispan alpha helical membrane proteins

Characterizing the interactions between amino acid residues and lipid molecules is important for understanding the assembly of transmembrane helices and for studying membrane protein folding. In this study we develop TMLIP (TransMembrane helix-LIPid), an empirically derived propensity of individual residue types to face lipid membrane based on statistical analysis of high-resolution structures of membrane proteins. Lipid accessibilities of amino acid residues within the transmembrane (TM) region of 29 structures of helical membrane proteins are studied with a spherical probe of radius of 1.9 A. Our results show that there are characteristic preferences for residues to face the headgroup region and the hydrocarbon core region of lipid membrane. Amino acid residues Lys, Arg, Trp, Phe, and Leu are often found exposed at the headgroup regions of the membrane, where they have high propensity to face phospholipid headgroups and glycerol backbones. In the hydrocarbon core region, the strongest preference for interacting with lipids is observed for Ile, Leu, Phe and Val. Small and polar amino acid residues are usually buried inside helical bundles and are strongly lipophobic. There is a strong correlation between various hydrophobicity scales and the propensity of a given residue to face the lipids in the hydrocarbon region of the bilayer. Our data suggest a possibly significant contribution of the lipophobic effect to the folding of membrane proteins. This study shows that membrane proteins have exceedingly apolar exteriors rather than highly polar interiors. Prediction of lipid-facing surfaces of boundary helices using TMLIP1 results in a 54% accuracy, which is significantly better than random (25% accuracy). We also compare performance of TMLIP with another lipid propensity scale, kPROT, and with several hydrophobicity scales using hydrophobic moment analysis.     

Folding and aggregation are selectively influenced by the conformational preferences of the alpha-helices of muscle acylphosphatase

The native state of human muscle acylphosphatase (AcP) presents two alpha-helices. In this study we have investigated folding and aggregation of a number of protein variants having mutations aimed at changing the propensity of these helical regions. Equilibrium and kinetic measurements of folding indicate that only helix-2, spanning residues 55-67, is largely stabilized in the transition state for folding therefore playing a relevant role in this process. On the contrary, the aggregation rate appears to vary only for the variants in which the propensity of the region corresponding to helix-1, spanning residues 22-32, is changed. Mutations that stabilize the first helix slow down the aggregation process while those that destabilize it increase the aggregation rate. AcP variants with the first helix destabilized aggregate with rates increased to different extents depending on whether the introduced mutations also alter the propensity to form beta-sheet structure. The fact that the first alpha-helix is important for aggregation and the second helix is important for folding indicates that these processes are highly specific. This partitioning does not reflect the difference in intrinsic alpha-helical propensities of the two helices, because helix-1 is the one presenting the highest propensity. Both processes of folding and aggregation do not therefore initiate from regions that have simply secondary structure propensities favorable for such processes. The identification of the regions involved in aggregation and the understanding of the factors that promote such a process are of fundamental importance to elucidate the principles by which proteins have evolved and for successful protein design.     

Analysis and prediction of helix-helix interactions in membrane channels and transporters

Membrane proteins span a large variety of different functions such as cell-surface receptors, redox proteins, ion channels, and transporters. Proteins with functional pores show different characteristics of helix-helix packing as other helical membrane proteins. We found that the helix-helix contacts of 13 nonhomologous high-resolution structures of membrane channels and transporters are mainly accomplished by weakly polar amino acids (G > S > T > F) that preferably create contacts every fourth residue, typical for right-handed helix crossings. There is a strong correlation between the now available biological hydrophobicity scale and the propensities of the weakly polar and hydrophobic residues to be buried at helix-helix interfaces or to be exposed to the lipids in membrane channels and transporters. The polar residues, however, make no major contribution towards the packing of their transmembrane helices, and are therefore subsumed to be primarily exposed to the polar milieu during the folding process. The contact formation of membrane channels and transporters is therefore ruled by the solubility of the residues, which we suppose to be the driving force for the assembly of their transmembrane helices. By contrast, in 14 nonhomologous high-resolution structures of other membrane protein coils, also large and polar amino acids (D > S > M > Q) create characteristic contacts every 3.5th residues, which is a signature for left-handed helix crossings. Accordingly, it seems that dependent on the function, different concepts of folding and stabilization are realized for helical membrane proteins. Using a sequence-based matrix prediction method these differences are exploited to improve the prediction of buried and exposed residues of transmembrane helices significantly. When the sequence motifs typical for membrane channels and transporters were applied for the prediction of helix-helix contacts the quality of prediction rises by 16% to an average value of 76%, compared to the same approach when only single amino acid positions are taken into account.     

Folding and assembly of beta-barrel membrane proteins

Beta-barrel membrane proteins occur in the outer membranes of Gram-negative bacteria, mitochondria and chloroplasts. The membrane-spanning sequences of beta-barrel membrane proteins are less hydrophobic than those of alpha-helical membrane proteins, which is probably the main reason why completely different folding and membrane assembly pathways have evolved for these two classes of membrane proteins. Some beta-barrel membrane proteins can be spontaneously refolded into lipid bilayer model membranes in vitro. They may also have this ability in vivo although lipid and protein chaperones likely assist with their assembly in appropriate target membranes. This review summarizes recent work on the thermodynamic stability and the mechanism of membrane insertion of beta-barrel membrane proteins in lipid model and biological membranes. How lipid compositions affect folding and assembly of beta-barrel membrane proteins is also reviewed. The stability of these proteins in membranes is not as large as previously thought (<10 kcal/mol) and is modulated by elastic forces of the lipid bilayer. Detailed kinetic studies indicate that beta-barrel membrane proteins fold in distinct steps with several intermediates that can be characterized in vitro. Formation of the barrel is synchronized with membrane insertion and all beta-hairpins insert simultaneously in a concerted pathway.     

Sequence periodicity and secondary structure propensity in model proteins

We explore the question of whether local effects (originating from the amino acids intrinsic secondary structure propensities) or nonlocal effects (reflecting the sequence of amino acids as a whole) play a larger role in determining the fold of globular proteins. Earlier circular dichroism studies have shown that the pattern of polar, non polar amino acids (nonlocal effect) dominates over the amino acid intrinsic propensity (local effect) in determining the secondary structure of oligomeric peptides. In this article, we present a coarse grained computational model that allows us to quantitatively estimate the role of local and nonlocal factors in determining both the secondary and tertiary structure of small, globular proteins. The amino acid intrinsic secondary structure propensity is modeled by a dihedral potential term. This dihedral potential is parametrized to match with experimental measurements of secondary structure propensity. Similarly, the magnitude of the attraction between hydrophobic residues is parametrized to match the experimental transfer free energies of hydrophobic amino acids. Under these parametrization conditions, we systematically explore the degree of frustration a given polar, non polar pattern can tolerate when the secondary structure intrinsic propensities are in opposition to it. When the parameters are in the biophysically relevant range, we observe that the fold of small, globular proteins is determined by the pattern of polar, non polar amino acids regardless of their instrinsic secondary structure propensities. Our simulations shed new light on previous observations that tertiary interactions are more influential in determining protein structure than secondary structure propensity. The fact that this can be inferred using a simple polymer model that lacks most of the biochemical details points to the fundamental importance of binary patterning in governing folding.     

A propensity scale for type II polyproline helices (PPII): aromatic amino acids in proline-rich sequences strongly disfavor PPII due to proline-aromatic interactions

Type II polyproline helices (PPII) are a fundamental secondary structure of proteins, common in globular and nonglobular regions and important in cellular signaling. We developed a propensity scale for PPII using a host-guest system with sequence Ac-GPPXPPGY-NH(2), where X represents any amino acid. We found that proline has the highest PPII propensity, but most other amino acids display significant PPII propensities. The PPII propensity of leucine was the highest of all propensities of non-proline residues. Alanine and residues with linear side chains displayed the next highest PPII propensities. Three classes of residues displayed lower PPII propensities: β-branched amino acids (Thr, Val, and Ile), short amino acids with polar side chains (Asn, protonated Asp, Ser, Thr, and Cys), and aromatic amino acids (Phe, Tyr, and Trp). tert-Leucine particularly disfavored PPII. The basis of the low PPII propensities of aromatic amino acids in this context was significant cis-trans isomerism, with proline-rich peptides containing aromatic residues exhibiting 45-60% cis amide bonds, due to Pro-cis-Pro-aromatic and aromatic-cis-Pro amide bonds.     

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.     

Folding of beta-sheets in membranes: specificity and promiscuity in peptide model systems

The interactions that drive the folding of beta-barrel membrane proteins have not been well studied because there have been few available model systems for membrane beta-sheets. In this work, we expand on a recently described model system to explore the contributions of interstrand hydrogen bonds, side-chain/side-chain interactions and side-chain/membrane interactions to beta-sheet formation in membranes. These experiments are based on the observation that the hydrophobic hexapeptide acetyl-Trp-Leu-Leu-Leu-Leu-Leu-OH (AcWLLLLL) folds, cooperatively and reversibly, into oligomeric, antiparallel beta-sheets in phosphatidylcholine membranes. To systematically characterize the important interactions that drive beta-sheet formation in membranes, we have used circular dichroism spectroscopy to determine the membrane secondary structure of each member of a complete host-guest family of related peptides of the form AcWLL-X-LL, where X is one of the natural amino acids. Peptides with hydrophobic X-residues of any size or character (X=Ala, Val, Ile, Leu, Cys, Met, Phe and Trp) form similar beta-sheets in membranes, while peptides with any polar X-residue or Gly or Pro at the X-position are random-coils, even when bound to membranes at high concentrations. The observed membrane sheet preferences correlate poorly with intrinsic sheet propensity scales measured in soluble proteins, but they correlate well with several membrane hydrophobicity scales. These results support the idea that the predominant interactions of the side-chains in membrane-bound beta-sheets are with the membrane lipids, and that backbone hydrogen bonding is the major driving force for the stabilization of beta-sheets in membranes.     

Peptide Partitioning Properties from Direct Insertion Studies

Partitioning properties of polypeptides are at the heart of biological membrane phenomena and their precise quantification is vital for ab-initio structure prediction and the accurate simulation of membrane protein folding and function. Recently the cellular translocon machinery has been employed to determine membrane insertion propensities and transfer energetics for a series of polyleucine segments embedded in a carrier sequence. We show here that the insertion propensity, pathway, and transfer energetics into synthetic POPC bilayers can be fully described by direct atomistic peptide partitioning simulations. The insertion probability as a function of peptide length follows two-state Boltzmann statistics, in agreement with the experiments. The simulations expose a systematic offset between translocon-mediated and direct insertion free energies. Compared to the experiment the insertion threshold is shifted toward shorter peptides by ∼2 leucine residues. The simulations reveal many hitherto unknown atomic-resolution details about the partitioning process and promise to provide a powerful tool for urgently needed calibration of lipid parameters to match experimentally observed peptide transfer energies.     

Structural and mechanistic investigations on Salmonella typhimurium acetate kinase (AckA): identification of a putative ligand binding pocket at the dimeric interface

Background

A large number of studies have been carried out to obtain amino acid propensities for ??-helices and ??-sheets. The obtained propensities for ??-helices are consistent with each other, and the pair-wise correlation coefficient is frequently high. On the other hand, the ??-sheet propensities obtained by several studies differed significantly, indicating that the context significantly affects ??-sheet propensity.

Results

We calculated amino acid propensities for ??-helices and ??-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 ??-helix propensities do not differ significantly by fold, but ??-sheet propensities are diverse and depend on the fold. The propensities calculated for exposed sites and buried sites are similar for ??-helix, but such is not the case for the ??-sheet propensities. We also found some fold dependence on amino acid frequency in ??-strands. Folds with a high Ser, Thr and Asn content at exposed sites in ??-strands tend to have a low Leu, Ile, Glu, Lys and Arg content (correlation coefficient = ?0.90) and to have flat ??-sheets. At buried sites in ??-strands, the content of Tyr, Trp, Gln and Ser correlates negatively with the content of Val, Ile and Leu (correlation coefficient = ?0.93). "All-??" proteins tend to have a higher content of Tyr, Trp, Gln and Ser, whereas "??/??" proteins tend to have a higher content of Val, Ile and Leu.

Conclusions

The ??-helix propensities are similar for all folds and for exposed and buried residues. However, ??-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 ??-strands. Furthermore, the correlations we detected suggest that amino acid composition is related to folding properties such as the twist of a ??-strand or association between two ?? sheets.     

Residue-specific α-helix propensities from molecular simulation

Formation of α-helices is a fundamental process in protein folding and assembly. By studying helix formation in molecular simulations of a series of alanine-based peptides, we obtain the temperature-dependent α-helix propensities of all 20 naturally occurring residues with two recent additive force fields, Amber ff03w and Amber ff99SB¹. Encouragingly, we find that the overall helix propensity of many residues is captured well by both energy functions, with Amber ff99SB¹ being more accurate. Nonetheless, there are some residues that deviate considerably from experiment, which can be attributed to two aspects of the energy function: i), variations of the charge model used to determine the atomic partial charges, with residues whose backbone charges differ most from alanine tending to have the largest error; ii), side-chain torsion potentials, as illustrated by the effect of modifications to the torsion angles of I, L, D, N. We find that constrained refitting of residue charges for charged residues in Amber ff99SB¹ significantly improves their helix propensity. The resulting parameters should more faithfully reproduce helix propensities in simulations of protein folding and disordered proteins.     

Protein folding kinetics provides a context-independent assessment of beta-strand propensity in the Fyn SH3 domain

Structural database-derived propensities for amino acids to adopt particular local protein structures, such as alpha-helix and beta-strand, have long been recognized and effectively exploited for the prediction of protein secondary structure. However, the experimental verification of database-derived propensities using mutagenesis studies has been problematic, especially for beta-strand propensities, because local structural preferences are often confounded by non-local interactions arising from formation of the native tertiary structure. Thus, the overall thermodynamic stability of a protein is not always altered in a predictable manner by changes in local structural propensity at a single position. In this study, we have undertaken an investigation of the relationship between beta-strand propensity and protein folding kinetics. By characterizing the effects of a wide variety of amino acid substitutions at two different beta-strand positions in an SH3 domain, we have found that the observed changes in protein folding rates are very well correlated to beta-strand propensities for almost all of the substitutions examined. In contrast, there is little correlation between propensities and unfolding rates. These data indicate that beta-strand conformation is well formed in the structured portion of the SH3 domain transition state, and that local structure propensity strongly influences the stability of the transition state. Since the transition state is known to be packed more loosely than the native state and likely lacks many of the non-local stabilizing interactions seen in the native state, we suggest that folding kinetics studies may generally provide an effective means for the experimental validation of database-derived local structural propensities.     

Conformational properties of a peptide model for unfolded alpha-helices

Models of protein folding often hypothesize that the first step is local secondary structure formation. The assumption is that unfolded polypeptide chains possess an intrinsic propensity to form these local secondary structures. On the basis of this idea, it is tempting to model the local conformational properties of unfolded proteins using well-established residue secondary structure propensities, in particular, alpha-helix forming propensities. We have used spectroscopic methods to investigate the conformational behavior of a host-guest series of peptides designed to model unfolded alpha-helices. A suitable peptide model for unfolded alpha-helices was determined from studies of the length dependence of the conformational properties of alanine-based peptides. The chosen host peptide possessed a small, detectable, alpha-helix content. Substituting various representative guest residues into the central position of the host peptide at times changed the conformational behavior dramatically, and often in ways that could not be predicted from known alpha-helix forming propensities. The data presented can be used to rationalize some of these propensities. However, it is clear that secondary structure propensities cannot be used to predict the local conformational properties of unfolded proteins.     

The beta-barrel finder (BBF) program,allowing identification of outer membrane beta-barrel proteins encoded within prokaryotic genomes

Many outer membrane proteins (OMPs) in Gram-negative bacteria possess known beta-barrel three-dimensional (3D) structures. These proteins, including channel-forming transmembrane porins, are diverse in sequence but exhibit common structural features. We here report computational analyses of six outer membrane proteins of known 3D structures with respect to (1) secondary structure, (2) hydropathy, and (3) amphipathicity. Using these characteristics, as well as the presence of an N-terminal targeting sequence, a program was developed allowing prediction of integral membrane beta-barrel proteins encoded within any completely sequenced prokaryotic genome. This program, termed the beta-barrel finder (BBF) program, was used to analyze the proteins encoded within the Escherichia coli genome. Out of 4290 sequences examined, 118 (2.8%) were retrieved. Of these, almost all known outer membrane proteins with established beta-barrel structures as well as many probable outer membrane proteins were identified. This program should be useful for predicting the occurrence of outer membrane proteins in bacteria with completely sequenced genomes.