Comparison of the hemocyanin beta-barrel with other Greek key beta-barrels: possible importance of the "beta-zipper" in protein structure and folding.

The Greek key beta-barrel topology is a folding motif observed in many proteins of widespread evolutionary origin. The arthropodan hemocyanins also have such a Greek key beta-barrel, which forms the core of the third domain of this protein. The hemocyanin beta-barrel was found to be structurally very similar to the beta-barrels of the immunoglobulin domains, Cu,Zn-superoxide dismutase and the chromophore carrying antitumor proteins. The structural similarity within this group of protein families is not accompanied by an evolutionary or functional relationship. It is therefore possible to study structure-sequence relations without bias from nonstructural constraints. The present study reports a conserved pattern of features in these Greek key beta-barrels that is strongly suggestive of a folding nucleation site. This proposed nucleation site, which we call a "beta-zipper," shows a pattern of well-conserved, large hydrophobic residues on two sequential beta-strands joined by a short loop. Each beta-zipper strand is near the center of one of the beta-sheets, so that the two strands face each other from opposite sides of the barrel and interact through their hydrophobic side chains, rather than forming a hydrogen-bonded beta-hairpin. Other protein families with Greek key beta-barrels that do not as strongly resemble the immunoglobulin fold--such as the azurins, plastocyanins, crystallins, and prealbumins--also contain the beta-zipper pattern, which might therefore be a universal feature of Greek key beta-barrel proteins.     

Topological diversity of artificial beta-barrels in water

Rigid-rod beta-barrels are composed of interdigitating, short, amphiphilic peptide strands flanked by stabilizing rigid-rod "staves". We here report studies on the topological diversity of these recently devised artificial beta-barrels with regard to their length. For this purpose, homologous p-octiphenyl, p-sexiphenyl, and p-quarterphenyl rods were equipped with complementary tripeptide strands based on the sequences Lys-Leu-Lys and Glu-Leu-Glu. The stability of rigid-rod beta-barrels of different length was determined by denaturation with guanidinium chloride. Free energies of delta GH2O = -5.2 kcalmol-1, delta GH2O = -2.9 kcalmol-1, and delta GH2O < -0.3 kcalmol-1 found for homologous p-octiphenyl, p-sexiphenyl, and p-quarterphenyl beta-barrels demonstrated strong dependence of beta-barrel stability on beta-barrel length. These results revealed a very qualitative minimal (approximately 23 A) and an "ideal" beta-barrel length (approximately 34 A), synergistic formation (alpha = 1.4) and remarkable stability for "ideal" p-octiphenyl beta-barrels exceeding that of several proteins and most synthetic models. Rigid-rod beta-barrels with p-oligophenyl "staves" longer than approximately 34 A will be very difficult to make and study because of rapidly decreasing rod solubilities. However, a strategy to bypass this apparent upper limitation of beta-barrel length is introduced: supramolecular matching of mismatched rods yielded elongated beta-barrels (61 A) of acceptable stability (delta GH2O = 2.2 - 3.1 kcalmol-1).     

Toward genomic identification of β-barrel membrane proteins: Composition and architecture of known structures

The amino acid composition and architecture of all beta-barrel membrane proteins of known three-dimensional structure have been examined to generate information that will be useful in identifying beta-barrels in genome databases. The database consists of 15 nonredundant structures, including several novel, recent structures. Known structures include monomeric, dimeric, and trimeric beta-barrels with between 8 and 22 membrane-spanning beta-strands each. For this analysis the membrane-interacting surfaces of the beta-barrels were identified with an experimentally derived, whole-residue hydrophobicity scale, and then the barrels were aligned normal to the bilayer and the position of the bilayer midplane was determined for each protein from the hydrophobicity profile. The abundance of each amino acid, relative to the genomic abundance, was calculated for the barrel exterior and interior. The architecture and diversity of known beta-barrels was also examined. For example, the distribution of rise-per-residue values perpendicular to the bilayer plane was found to be 2.7 +/- 0.25 A per residue, or about 10 +/- 1 residues across the membrane. Also, as noted by other authors, nearly every known membrane-spanning beta-barrel strand was found to have a short loop of seven residues or less connecting it to at least one adjacent strand. Using this information we have begun to generate rapid screening algorithms for the identification of beta-barrel membrane proteins in genomic databases. Application of one algorithm to the genomes of Escherichia coli and Pseudomonas aeruginosa confirms its ability to identify beta-barrels, and reveals dozens of unidentified open reading frames that potentially code for beta-barrel outer membrane proteins.     

The versatile beta-barrel membrane protein

The beta-barrel membrane protein is found in the outer membranes of bacteria, mitochondria and chloroplasts. Approximately 2-3% of the genes in Gram-negative bacterial genomes encode beta-barrels. Whereas there are fewer than 20 known three-dimensional beta-barrel structures, genomic databases currently contain thousands of beta-barrels belonging to dozens of families. New research is revealing the variety of beta-barrel structures and the variety of functions performed by these versatile proteins.     

The anatomy of protein beta-sheet topology

Here, we present a systematic analysis of the open-faced beta-sheet topologies in a set of non-redundant protein domain structures; in particular, we focus on the topological diversity of four-stranded beta-sheet motifs. Of the 96 topologies that are possible for a four-stranded beta-sheet, 42 were identified in known protein structures. Of these, four account for 50% of the structures that we have studied. Two sets of the topologies that were not observed may represent the section of the topological space that is not readily accessible to proteins on either thermodynamic or kinetic grounds. The first set contains topologies with alternating parallel and antiparallel beta-ladders. Their rare occurrence reflects the expectation that it is energetically unfavorable to match different hydrogen bonding patterns. The polypeptide chains in the second set of topologies go through convoluted paths and are expected to experience great kinetic frustrations during the folding processes. A knowledge of the potential causes for the topological preference of small beta-sheets also helps us to understand the topological properties of larger beta-sheet structures which frequently contain four-stranded motifs. The notion that protein topologies can only be taken from a confined and discrete space has important implications for structural genomics.     

Gene duplication of the eight-stranded beta-barrel OmpX produces a functional pore: a scenario for the evolution of transmembrane beta-barrels

The repeating unit of outer membrane beta-barrels from Gram-negative bacteria is the beta-hairpin, and representatives of this protein family always have an even strand number between eight and 22. Two dominant structural forms have eight and 16 strands, respectively, suggesting gene duplication as a possible mechanism for their evolution. We duplicated the sequence of OmpX, an eight-stranded beta-barrel protein of known structure, and obtained a beta-barrel, designated Omp2X, which can fold in vitro and in vivo. Using single-channel conductance measurements and PEG exclusion assays, we found that Omp2X has a pore size similar to that of OmpC, a natural 16-stranded barrel. Fusions of the homologous proteins OmpX, OmpA and OmpW were able to fold in vitro in all combinations tested, revealing that the general propensity to form a beta-barrel is sufficient to evolve larger barrels by simple genetic events.     

Transbilayer pores formed by beta-barrels: molecular modeling of pore structures and properties.

Transmembrane beta-barrels, first observed in bacterial porins, are possible models for a number of membrane channels. Restrained molecular dynamics simulations based on idealized C alpha beta templates have been used to generate models of such beta-barrels. Model beta-barrels have been analyzed in terms of their conformational, energetic, and pore properties. Model beta-barrels formed by N = 4, 8, 12 and 16 anti-parallel Ala10 strands have been developed. For each N, beta-barrels with shear numbers S = N to 2N have been modeled. In all beta-barrel models the constituent beta-strands adopt a pronounced right-handed twist. Interstrand interactions are of approximately equal stability for all models with N > or = 8, whereas such interactions are weaker for the N = 4 beta-barrels. In N = 4 beta-barrels the pore is too narrow (minimum radius approximately 0.6 A) to allow ion permeation. For N > or = 8, the pore radius depends on both N and S; for a given value of N an increase in S from N to 2N is predicted to result in an approximately threefold increase in pore conductance. Calculated maximal conductances for the beta-barrel models are compared with experimental values for porins and for K+ channels.     

The structure of bacterial outer membrane proteins

Integral membrane proteins come in two types, alpha-helical and beta-barrel proteins. In both types, all hydrogen bonding donors and acceptors of the polypeptide backbone are completely compensated and buried while nonpolar side chains point to the membrane. The alpha-helical type is more abundant and occurs in cytoplasmic (or inner) membranes, whereas the beta-barrels are known from outer membranes of bacteria. The beta-barrel construction is described by the number of strands and the shear number, which is a measure for the inclination angle of the beta-strands against the barrel axis. The common right-handed beta-twist requires shear numbers slightly larger than the number of strands. Membrane protein beta-barrels contain between 8 and 22 beta-strands and have a simple topology that is probably enforced by the folding process. The smallest barrels form inverse micelles and work as enzymes or they bind to other macromolecules. The medium-range barrels form more or less specific pores for nutrient uptake, whereas the largest barrels occur in active Fe(2+) transporters. The beta-barrels are suitable objects for channel engineering, because the structures are simple and because many of these proteins can be produced into inclusion bodies and recovered therefrom in the exact native conformation.     

Conformational and geometrical properties of idealized beta-barrels in proteins

An equation for calculating the distances between the atoms involved in forming an idealized hydrogen bond in a parallel or antiparallel beta-barrel has been derived by adjusting the corresponding data given by Pauling and Corey for a beta-sheet. Based on these distances, a geometrical optimization method was developed, by which one can generate various idealized beta-barrels: parallel or antiparallel, tilted or non-tilted, right-tilted or left-tilted. For each type of idealized beta-barrel thus obtained, the corresponding conformation and characteristic geometric parameters as well as their relationship are analyzed and discussed. Since the strand in a tilted beta-barrel traces a curve rather than a straight line on a cylinder-like surface, a regular chain in which the dihedral angles of each residue are the same cannot form a tilted beta-barrel but only a non-tilted beta-barrel. As observed, the strands of a right-tilted beta-barrel possess a very strong right-handed twist. The radii of the idealized tilted parallel and antiparallel beta-barrels are greater than those of the corresponding non-tilted ones by approximately 1 A and approximately 1.5 A, respectively. Consequently, there is relatively more room for a tilted beta-barrel to accommodate the internal side-chains, suggesting that a conformational change from a non-tilted beta-barrel to a tilted one would ease the repulsion among the crowded internal side-chains so as to make the structure more stable. The values of root-mean-square fits indicate that the idealized right-tilted beta-barrels coincide quite well with the observed beta-barrels in both parallel and antiparallel cases.     

Energetics of the structure and chain tilting of antiparallel beta-barrels in proteins

The preferred structural pattern of antiparallel beta-barrels in proteins, described as the right-handed tilting of the peptide strands with respect to the axis of the barrel, is accounted for in terms of intra- and interchain interaction energies. It is related to the preference of beta-sheets for right-handed twisting. Conformational energy computations have been carried out on three eight-stranded antiparallel beta-barrels composed of six-residue strands, in which L-Val and Gly alternate, and having a right-handed, a left-handed, or no tilt. After energy minimization, the relative energies of these structures were 0.0, 8.6, and 46.1 kcal/mol, respectively; i.e., the right-tilted beta-barrel is favored energetically, in agreement with anti-parallel beta-barrels observed in proteins. Tilting of the barrel is favored, relative to the nontilted structure, by both intra- and interstrand interactions, because tilting allows better packing of the bulky side chains. On the other hand, the energy difference between the left- and right-tilted barrels arises essentially from intrachain interactions. This is a consequence of the preference of beta-sheets for a right-handed twist. Space limitations inside the barrel are satisfied if there is an alternation of bulky residues and residues with small or no side chain (preferably Gly) in neighboring positions on adjacent strands. Such a pattern is seen frequently in antiparallel beta-barrels of globular proteins. The computations indicate that a structure with Val...Gly pairs can be accommodated in a beta-barrel with no distortion.     

Estimating the twist of beta-strands embedded within a regular parallel beta-barrel structure

The parallel beta-barrel is a recurrent structural motif found in a large variety of different enzymes belonging to the family of alpha/beta-proteins. It has been shown previously that the hyperboloid can be considered as a scaffold describing the parallel beta-barrel structure. To assess restraints on beta-strand twist imposed by a given scaffold geometry, the notion of scaffold twist, Ts, is introduced. From Ts, the beta-strand twist (Tw beta) expected for a given scaffold geometry can be derived and it is verified that this computed twist can be used to identify beta-barrels characterized by good hydrogen bonding. It is noted that Tw beta is only slightly affected for beta-barrels differing in the number (N) of beta-strands, suggesting that restraints on main-chain conformation of beta-strands are not likely to account for the N = 8 invariability observed in natural parallel beta-barrels thereby strengthening previous work rationalizing this constancy.     

OB(oligonucleotide/oligosaccharide binding)-fold: common structural and functional solution for non-homologous sequences.

A novel folding motif has been observed in four different proteins which bind oligonucleotides or oligosaccharides: staphylococcal nuclease, anticodon binding domain of asp-tRNA synthetase and B-subunits of heat-labile enterotoxin and verotoxin-1. The common fold of the four proteins, which we call the OB-fold, has a five-stranded beta-sheet coiled to form a closed beta-barrel. This barrel is capped by an alpha-helix located between the third and fourth strands. The barrel-helix frameworks can be superimposed with r.m.s. deviations of 1.4-2.2 A, but no similarities can be observed in the corresponding alignment of the four sequences. The nucleotide or sugar binding sites, known for three of the four proteins, are located in nearly the same position in each protein: on the side surface of the beta-barrel, where three loops come together. Here we describe the determinants of the OB-fold, based on an analysis of all four structures. These proposed determinants explain how very different sequences adopt the OB-fold. They also suggest a reinterpretation of the controversial structure of gene 5 ssDNA binding protein, which exhibits some topological and functional similarities with the OB-fold proteins.     

Tilt, twist, and coiling in beta-barrel membrane proteins: relation to infrared dichroism

The x-ray coordinates of beta-barrel transmembrane proteins from the porins superfamily and relatives are used to calculate the mean tilt of the beta-strands and their mean local twist and coiling angles. The 13 proteins examined correspond to beta-barrels with 8 to 22 strands, and shear numbers ranging from 8 to 24. The results are compared with predictions from the model of Murzin, Lesk, and Chothia for symmetrical regular barrels. Good agreement is found for the mean strand tilt, but the twist angles are smaller than those for open beta-sheets and beta-barrels with shorter strands. The model is reparameterised to account for the reduced twist characteristic of long-stranded transmembrane beta-barrels. This produces predictions of both twist and coiling angles that are in agreement with the mean values obtained from the x-ray structures. With the optimized parameters, the model can then be used to determine twist and coiling angles of transmembrane beta-barrels from measurements of the amide band infrared dichroism in oriented membranes. Satisfactory agreement is obtained for OmpF. The strand tilt obtained from the x-ray coordinates, or from the reparameterised model, can be combined with infrared dichroism measurements to obtain information on the orientation of the beta-barrel assembly in the membrane.     

Folding kinetics of the outer membrane proteins OmpA and FomA into phospholipid bilayers

The folding mechanism of outer membrane proteins (OMPs) of Gram-negative bacteria into lipid bilayers has been studied using OmpA of E. coli and FomA of F. nucleatum as examples. Both, OmpA and FomA are soluble in unfolded form in urea and insert and fold into phospholipid bilayers upon strong dilution of the denaturant urea. OmpA is a structural protein and forms a small ion channel, composed of an 8-stranded transmembrane beta-barrel domain. FomA is a voltage-dependent porin, predicted to form a 14 stranded beta-barrel. Both OMPs fold into a range of model membranes of very different phospholipid compositions. Three membrane-bound folding intermediates of OmpA were discovered in folding studies with dioleoylphosphatidylcholine bilayers that demonstrated a highly synchronized mechanism of secondary and tertiary structure formation of beta-barrel membrane proteins. A study on FomA folding into lipid bilayers indicated the presence of parallel folding pathways for OMPs with larger transmembrane beta-barrels.     

Origami in the outer membrane: the transmembrane arrangement of mitochondrial porins.

Voltage-dependent anion-selective channels (VDAC), also known as mitochondrial porins, are key regulators of metabolite flow across the mitochondrial outer membrane. Porins from a wide variety of organisms share remarkably similar electrophysiological properties, in spite of considerable sequence dissimilarity, indicating that they share a common structure. Based on primary sequence considerations, analogy with bacterial porins, and circular dichroism analysis, it is agreed that VDAC spans the outer membrane as a beta-barrel. However, the residues that form the antiparallel beta-strands comprising this barrel remain unknown. Various predictive methods, largely based on the known structures of bacterial beta-barrels, have been applied to the primary sequences of VDAC. Refinement and confirmation of these predictions have developed through numerous investigations of wild-type and variant porins, both in mitochondria and in artificial membranes. These experiments have involved VDAC from several sources, precluding the generation of a unified model. Herein, using the Neurospora VDAC sequence as a template, the published structural information and predictions have been reassessed to delineate a model that satisfies most of the available data.     

Calcium binding properties of gamma-crystallin: calcium ion binds at the Greek key beta gamma-crystallin fold

The beta- and gamma-crystallins are closely related lens proteins that are members of the betagamma-crystallin superfamily, which also include many non-lens members. Although beta-crystallin is known to be a calcium-binding protein, this property has not been reported in gamma-crystallin. We have studied the calcium binding properties of gamma-crystallin, and we show that it binds 4 mol eq of calcium with a dissociation constant of 90 microm. It also binds the calcium-mimic spectral probes, terbium and Stains-all. Calcium binding does not significantly influence protein secondary and tertiary structures. We present evidence that the Greek key crystallin fold is the site for calcium ion binding in gamma-crystallin. Peptides corresponding to Greek key motif of gamma-crystallin (42 residues) and their mutants were synthesized and studied for calcium binding. These peptides adopt beta-sheet conformation and form aggregates producing beta-sandwich. Our results with peptides show that, in Greek key motif, the amino acid adjacent to the conserved aromatic corner in the "a" strand and three amino acids of the "d" strand participate in calcium binding. We suggest that the betagamma superfamily represents a novel class of calcium-binding proteins with the Greek key betagamma-crystallin fold as potential calcium-binding sites. These results are of significance in understanding the mechanism of calcium homeostasis in the lens.     

Contributions of aromatic pairs to the folding and stability of long-lived human γD-crystallin

Human γD-crystallin (HγD-Crys) is a highly stable protein that remains folded in the eye lens for the majority of an individual's lifetime. HγD-Crys exhibits two homologous crystallin domains, each containing two Greek key motifs with eight β-strands. Six aromatic pairs (four Tyr/Tyr, one Tyr/Phe and one Phe/Phe) are present in the β-hairpin sequences of the Greek keys. Ultraviolet damage to the aromatic residues in lens crystallins may contribute to the genesis of cataract. Mutant proteins with these aromatic residues substituted with alanines were constructed and expressed in E. coli. All mutant proteins except F115A and F117A had lower thermal stability than the WT protein. In equilibrium experiments in guanidine hydrochloride (GuHCl), all mutant proteins had lower thermodynamic stability than the WT protein. N-terminal domain (N-td) substitutions shifted the N-td transition to lower GuHCl concentration, but the C-terminal domain (C-td) transition remained unaffected. C-td substitutions led to a more cooperative unfolding/refolding process, with both the N-td and C-td transitions shifted to lower GuHCl concentration. The aromatic pairs conserved for each Greek key motif (Greek key pairs) had larger contributions to both thermal stability and thermodynamic stability than the other pairs. Aromatic-aromatic interaction was estimated as 1.5-2.0 kcal/mol. In kinetic experiments, N-td substitutions accelerated the early phase of unfolding, while C-td substitutions accelerated the late phase, suggesting independent domain unfolding. Only substitutions of the second Greek key pair of each crystallin domain slowed refolding. The second Greek keys may provide nucleation sites during the folding of the double-Greek-key crystallin domains.     

The pore-lining region of shaker voltage-gated potassium channels: comparison of beta-barrel and alpha-helix bundle models.

Although there is a large body of site-directed mutagenesis data that identify the pore-lining sequence of the voltage-gated potassium channel, the structure of this region remains unknown. We have interpreted the available biochemical data as a set of topological and orientational restraints and employed these restraints to produce molecular models of the potassium channel pore region, H5. The H5 sequence has been modeled either as a tetramer of membrane-spanning beta-hairpins, thus producing an eight-stranded beta-barrel, or as a tetramer of incompletely membrane-spanning alpha-helical hairpins, thus producing an eight-staved alpha-helix bundle. In total, restraints-directed modeling has produced 40 different configurations of the beta-barrel model, each configuration comprising an ensemble of 20 structures, and 24 different configurations of the alpha-helix bundle model, each comprising an ensemble of 24 structures. Thus, over 1300 model structures for H5 have been generated. Configurations have been ranked on the basis of their predicted pore properties and on the extent of their agreement with the biochemical data. This ranking is employed to identify particular configurations of H5 that may be explored further as models of the pore-lining region of the voltage-gated potassium channel pore.     

The structure of the translational initiation factor IF1 from E.coli contains an oligomer-binding motif.

The structure of the translational initiation factor IF1 from Escherichia coli has been determined with multidimensional NMR spectroscopy. Using 1041 distance and 78 dihedral constraints, 40 distance geometry structures were calculated, which were refined by restrained molecular dynamics. From this set, 19 structures were selected, having low constraint energy and few constraint violations. The ensemble of 19 structures displays a root-mean-square deviation versus the average of 0.49 A for the backbone atoms and 1.12 A for all atoms for residues 6-36 and 46-67. The structure of IF1 is characterized by a five-stranded beta-barrel. The loop connecting strands three and four contains a short 3(10) helix but this region shows considerably higher flexibility than the beta-barrel. The fold of IF1 is very similar to that found in the bacterial cold shock proteins CspA and CspB, the N-terminal domain of aspartyl-tRNA synthetase and the staphylococcal nuclease, and can be identified as the oligomer-binding motif. Several proteins of this family are nucleic acid-binding proteins. This suggests that IF1 plays its role in the initiation of protein synthesis by nucleic acid interactions. Specific changes of NMR signals of IF1 upon titration with 30S ribosomal subunit identifies several residues that are involved in the interaction with ribosomes.