Aplysia attractin: biophysical characterization and modeling of a water-borne pheromone.

Attractin, a 58-residue protein secreted by the mollusk Aplysia californica, stimulates sexually mature animals to approach egg cordons. Attractin from five different Aplysia species are approximately 40% identical in sequence. Recombinant attractin, expressed in insect cells and purified by reverse-phase high-performance liquid chromatography (RP-HPLC), is active in a bioassay using A. brasiliana; its circular dichroism (CD) spectrum indicates a predominantly alpha-helical structure. Matrix-assisted laser desorption/ionization mass spectrometry (MALDI-MS) characterization of proteolytic fragments identified disulfide bonds between the six conserved cysteines (I-VI, II-V, III-IV, where the Roman numeral indicates the order of occurrence in the primary sequence). Attractin has no significant similarity to any other sequence in the database. The protozoan Euplotes pheromones were selected by fold recognition as possible templates. These diverse proteins have three alpha-helices, with six cysteine residues disulfide-bonded in a different pattern from attractin. Model structures with good stereochemical parameters were prepared using the EXDIS/DIAMOD/FANTOM program suite and constraints based on sequence alignments with the Euplotes templates and the attractin disulfide bonds. A potential receptor-binding site is suggested based on these data. Future structural characterization of attractin will be needed to confirm these models.     

Membrane proteins: the 'Wild West' of structural biology

Historically, the task of determining the structure of membrane proteins has been hindered by experimental difficulties associated with their lipid-embedded domains. Here, we provide an overview of recently developed experimental and predictive tools that are changing our view of this largely unexplored territory - the 'Wild West' of structural biology. Crystallography, single-particle methods and atomic force microscopy are being used to study huge membrane proteins with increasing detail. Solid-state nuclear magnetic resonance strategies provide orientational constraints for structure determination of transmembrane (TM) alpha-helices and accurate measurements of intramolecular distances, even in very complex systems. Longer distance constraints are determined by site-directed spin-labelling electron paramagnetic resonance, but current labelling strategies still constitute some limitation. Other methods, such as site-specific infrared dichroism, enable orientational analysis of TM alpha-helices in aligned bilayers and, combined with novel computational and predictive tools that use evolutionary conservation data, are being used to analyze TM alpha-helical bundles.     

GXXXG and AXXXA: common alpha-helical interaction motifs in proteins, particularly in extremophiles

The GXXXG motif is a frequently occurring sequence of residues that is known to favor helix-helix interactions in membrane proteins. Here we show that the GXXXG motif is also prevalent in soluble proteins whose structures have been determined. Some 152 proteins from a non-redundant PDB set contain at least one alpha-helix with the GXXXG motif, 41 +/- 9% more than expected if glycine residues were uniformly distributed in those alpha-helices. More than 50% of the GXXXG-containing alpha-helices participate in helix-helix interactions. In fact, 26 of those helix-helix interactions are structurally similar to the helix-helix interaction of the glycophorin A dimer, where two transmembrane helices associate to form a dimer stabilized by the GXXXG motif. As for the glycophorin A structure, we find backbone-to-backbone atomic contacts of the C alpha-H...O type in each of these 26 helix-helix interactions that display the stereochemical hallmarks of hydrogen bond formation. These glycophorin A-like helix-helix interactions are enriched in the general set of helix-helix interactions containing the GXXXG motif, suggesting that the inferred C alpha-H...O hydrogen bonds stabilize the helix-helix interactions. In addition to the GXXXG motif, some 808 proteins from the non-redundant PDB set contain at least one alpha-helix with the AXXXA motif (30 +/- 3% greater than expected). Both the GXXXG and AXXXA motifs occur frequently in predicted alpha-helices from 24 fully sequenced genomes. Occurrence of the AXXXA motif is enhanced to a greater extent in thermophiles than in mesophiles, suggesting that helical interaction based on the AXXXA motif may be a common mechanism of thermostability in protein structures. We conclude that the GXXXG sequence motif stabilizes helix-helix interactions in proteins, and that the AXXXA sequence motif also stabilizes the folded state of proteins.     

Left handed alpha-helix formation by a bacterial peptide

The alpha-helix is a common element of secondary structure in proteins and peptides. In eukaryotic organisms, which exclusively incorporate L-amino acids into such molecules, stereochemical interactions make such alpha-helices, invariably right-handed. Pseudomonas tolaasii Paine is the causal organism of the economically significant brown blotch disease of the cultivated mushroom Agaricus bisporus (Lange) Imbach. P. Tolaasii proceduces an extracellular lipodepsipeptide toxin, tolaasin, which causes the brown pitted lesions on the mushroom cap. Circular dichroism studies on tolaasin in a membrane-like environment indicate the presence of a left-handed alpha-helix, probably formed by a sequence of 7 D-amino acids in the peptide. P. tolaasii represents the first reported example of an organism which has evolved the ability to biosynthesize a left-handed alpha-helix.     

Comparison of the structures of cro and lambda repressor proteins from bacteriophage lambda

The three-dimensional structures of cro repressor protein and of the amino-terminal domain of lambda repressor protein, both from bacteriophage lambda, are compared. The second and third alpha-helices, alpha 2 and alpha 3, are shown to have essentially identical conformations in the two proteins, confirming the significance of the amino acid sequence homology previously noted between these and other DNA binding proteins in the region corresponding to these helices. The correspondence between the two-helical units in cro and lambda repressor protein is better than the striking agreement noted previously between two-helical units in cro and catabolite gene-activator protein. Parts of the first alpha-helices of repressor and cro show a structural correspondence that suggests a revised sequence homology between the two proteins in their extreme amino-terminal regions. In particular, there is a short loop between the alpha 1 and alpha 2 helices of lambda repressor that is missing from cro. This structural difference may account for the observed differences found with different cros and repressors in the pattern of phosphates whose ethylation prevents the binding of these proteins to their specific recognition sites. Although the two proteins have strikingly similar alpha 2-alpha 3 helical units that are presumed to bind to DNA in an essentially similar manner, stereochemical restrictions prevent the alpha 2-alpha 3 units of the respective proteins aligning on the DNA in exactly the same way.     

Towards protein tertiary fold prediction using distance and motif constraints

Based on a simplified model of the all-alpha class of protein, all packing arrangements of alpha-helices were generated and assessed by both general and specific structural rules. The method was applied to myoglobin and parvalbumin, which were both ranked in the top 4% of folds under the general packing constraints. Incorporation of the restrictions implied by the EF-hand motifs of parvalbumin were sufficient to select the correct fold as one of two (equal scoring) possibilities. Myoglobin scored well under the general packing constraints and the addition of a single distance constraint, implied by haem binding, was sufficient to select the correct fold as one of several candidates. Incorporation of a score for complementary hydrophobic packing between helices further selected myoglobin as a unique fold but did not improve the ranking of parvalbumin. For both proteins, the alpha-helices were predicted from multiply aligned sequences using pattern-matching methods and no specific aspect of the known X-ray structures influenced this or the prediction of the correct folds. Although the method is currently of limited generality, its further applications and extension to a more detailed structural level are discussed.     

A test of proposed rules for helix capping: implications for protein design

alpha-helices within proteins are often terminated (capped) by distinctive configurations of the polypeptide chain. Two common arrangements are the Schellman motif and the alternative alpha(L) motif. Rose and coworkers developed stereochemical rules to identify the locations of such motifs in proteins of unknown structure based only on their amino acid sequences. To check the effectiveness of these rules, they made specific predictions regarding the structural and thermodynamic consequences of certain mutations in T4 lysozyme. We have constructed these mutants and show here that they have neither the structure nor the stability that was predicted. The results show the complexity of the protein-folding problem. Comparison of known protein structures may show that a characteristic sequence of amino acids (a sequence motif) corresponds to a conserved structural motif. In any particular protein, however, changes in other parts of the sequence may result in a different conformation. The structure is determined by sequence as a whole, not by parts considered in isolation.     

Stereochemical theory of the 3-dimensional structure of globular proteins. I. Highly helical intermediate structures

The authors analyze the physical prerequisites on which the proposed stereochemical theory of the three-dimensional structure of globular proteins is based. The theory represents a stereochemical modelling of the mechanism of protein self-organization suggested earlier by one of the authors. According to this mechanism, a highly helical intermediate structure(s) is formed at first and then it passes into the native one. In the highly-helical intermediate structure the arrangement of the polypeptide chain in space is the same as in the native structure. These two structures differ mainly by the secondary structure of the chain. The transition into the native structure proceeds under the effect of long-range interactions which transform the excess alpha-helices into beta-structural and irregular conformations. The so-called s-helices are considered (the alpha-helix, whose hydrophobic groups form a separate cluster on its surface). s-Helices can be obtained on the greater part of the polypeptide chain of any globular protein. In the unfolded protein chain they are the most stable and rapidly formed structures. It has been shown that namely s-helices are the initial blocks for the formation of the highly-helical intermediate structure. Stereochemical principles of the s-helix packing that permit to predict the three-dimensional structure of highly helical proteins have been found. According to these principles the highly helical structure represents the packing of hydrophobic surfaces and s-helices. In their turn, hydrophobic surfaces are formed as a result of complementary interaction of borders of hydrophobic clusters of two s-helices according to the "knob-hole" principle.     

Limits on alpha-helix prediction with neural network models.

Using a backpropagation neural network model we have found a limit for secondary structure prediction from local sequence. By including only sequences from whole alpha-helix and non-alpha-helix structures in our training and test sets--sequences spanning boundaries between these two structures were excluded--it was possible to investigate directly the relationship between sequence and structure for alpha-helix. A group of non-alpha-helix sequences, that was disrupting overall prediction success, was indistinguishable to the network from alpha-helix sequences. These sequences were found to occur at regions adjacent to the termini of alpha-helices with statistical significance, suggesting that potentially longer alpha-helices are disrupted by global constraints. Some of these regions spanned more than 20 residues. On these whole structure sequences, 10 residues in length, a comparatively high prediction success of 78% with a correlation coefficient of 0.52 was achieved. In addition, the structure of the input space, the distribution of beta-sheet in this space, and the effect of segment length were also investigated.     

Bundles of amphipathic transmembrane alpha-helices as a structural motif for ion-conducting channel proteins: studies on sodium channels and acetylcholine receptors

Channel proteins are transmembrane symmetric (or pseudosymmetric) oligomers organized around a central ionic pore. We present here a molecular model of the pore forming structures of two channel proteins with different primary structures and oligomeric size: the voltage-sensitive sodium channel and the nicotinic cholinergic receptor. We report low-energy arrangements of alpha-helical bundles calculated by semiempiricial potential energy functions and optimization routines and further refined using molecular dynamics. The ion-conducting pore is considered to be a symmetric or pseudosymmetric homooligomer of 3-5 amphipathic alpha-helices arranged such that the polar residues line a central hydrophilic pathway and the apolar residues face the hydrophobic bilayer interior. The channel lining exposes either charged (Asp, Glu, Arg, Lys) or polar-neutral (Ser, Thr) residues. A bundle of four parallel helices constrained to C4 symmetry, the helix axis aligned with the symmetry axis, and the helices constrained to idealized dihedral angles, produces a structure with a pore of the size inferred for the sodium channel protein (area approximately 16 A2). Similarly, a pentameric array optimized with constraints to maintain C5 symmetry and backbone torsions characteristic of alpha-helices adopts a structure that appears well suited to form the lining of the nicotinic cholinergic receptor (pore area approximately 46 A2). Thus, bundles of amphipathic alpha-helices satisfy the structural, energetic, and dynamic requirements to be the molecular structural motif underlying the function of ionic channels.     

Helix geometry in proteins

In this report we describe a general survey of all helices found in 57 of the known protein crystal structures, together with a detailed analysis of 48 alpha-helices found in 16 of the structures that are determined to high resolution. The survey of all helices reveals a total of 291 alpha-helices, 71 3(10)-helices and no examples of pi-helices. The conformations of the observed helices are significantly different from the "ideal" linear structures. The mean phi, psi angles for the alpha- and 3(10)-helices found in proteins are, respectively, (-62 degrees, -41 degrees) and (-71 degrees, -18 degrees). A computer program, HBEND, is used to characterize and to quantify the different types of helix distortion. alpha-Helices are classified as regular or irregular, linear, curved or kinked. Of the 48 alpha-helices analysed, only 15% are considered to be linear; 17% are kinked, and 58% are curved. The curvature of helices is caused by differences in the peptide hydrogen bonding on opposite faces of the helix, reflecting carbonyl-solvent/side-chain interactions for the exposed residues, and packing constraints for residues involved in the hydrophobic core. Kinked helices arise either as a result of included proline residues, or because of conflicting requirements for the optimal packing of the helix side-chains. In alpha-helices where there are kinks caused by proline residues, we show that the angle of kink is relatively constant (approximately 26 degrees), and that there is minimal disruption of the helix hydrogen bonding. The proline residues responsible for the kinks are highly conserved, suggesting that these distortions may be structurally/functionally important.     

Single proline substitutions in predicted alpha-helices of murine granulocyte-macrophage colony-stimulating factor result in a loss in bioactivity and altered glycosylation

Contributions of alpha-helices to biological activity in murine granulocyte-macrophage colony-stimulating factor were analyzed using site-directed mutagenesis and protein expression in COS-1 cells. A series of single proline substitutions were made for residues within the four predicted alpha-helices as a means of disrupting local helical secondary structure. Mutations in three of the four helices resulted in marked reductions in bioactivity. Five mutants E21P, L56P, E60P, L63P, and L107P showed 10(2)-10(4)-fold reduction in bioactivity as well as hyperglycosylation. The same Pro substitutions made on non-N-glycosylated molecules had a similar loss in bioactivity implying that a Pro-induced structural change and not hyperglycosylation was responsible for the major decrease in bioactivity. Additional amino acid substitutions at these residues which conserved charge or hydrophobicity, or replaced the original residue with an Ala, verified that conformational changes in the protein structure were specifically due to steric constraints imposed by the Pro residue rather than loss of important side chain functions.     

Discovery of a significant,nontopological preference for antiparallel alignment of helices with parallel regions in sheets

To help elucidate the role of secondary structure packing preferences in protein folding, here we present an analysis of the packing geometry observed between alpha-helices and between alpha-helices and beta-sheets in 1316 diverse, nonredundant protein structures. Finite-length vectors were fit to the alpha-carbon atoms in each of the helices and strands, and the packing angle between the vectors, Omega, was determined at the closest point of approach within each helix-helix or helix-sheet pair. Helix-sheet interactions were found in 391 of the proteins, and the distributions of Omega values were calculated for all the helix-sheet and helix-helix interactions. The packing angle preferences for helix-helix interactions are similar to those previously observed. However, analysis of helix-strand packing preferences uncovered a remarkable tendency for helices to align antiparallel to parallel regions of beta-sheets, independent of the topological constraints or prevalence of beta-alpha-beta motifs in the proteins. This packing angle preference is significantly diminished in helix interactions involving mixed and antiparallel beta-sheets, suggesting a role for helix-sheet dipole alignment in guiding supersecondary structure formation in protein folding. This knowledge of preferred packing angles can be used to guide the engineering of stable protein modules.     

Recurring loop motif in proteins that occurs in right-handed and left-handed forms. Its relationship with alpha-helices and beta-bulge loops

A common feature of alpha-helices in proteins is a loop at the C-terminal end, with a characteristic hydrogen bond pattern. It is noted that several loops with the same structural features occur independently of alpha-helices; two are even situated at the loop ends of beta-hairpins. The name paperclip is suggested for loops possessing the appropriate hydrogen bonds. A number of features of paperclips are described: they exist in two classes, depending on the number of residues at the loop end; one class is very much commoner than the other. Two paperclips are found that belong to the common class, except that the main-chain conformation of each is the mirror image of that normally found. The majority of paperclips are shown to have tightly clustered sets of main-chain dihedral angles. These are somewhat similar to, but distinct from, a subgroup of another common family of loops that have been called beta-bulge loops; in the latter, the dihedral angles are also tightly clustered. The high degree of clustering in both cases is likely to be a result of steric constraints associated with hydrogen bond patterns at the ends of loops.     

New Insights into the Interdependence between Amino Acid Stereochemistry and Protein Structure

To successfully design new proteins and understand the effects of mutations in natural proteins, we must understand the geometric and physicochemical principles underlying protein structure. The side chains of amino acids in peptides and proteins adopt specific dihedral angle combinations; however, we still do not have a fundamental quantitative understanding of why some side-chain dihedral angle combinations are highly populated and others are not. Here we employ a hard-sphere plus stereochemical constraint model of dipeptide mimetics to enumerate the side-chain dihedral angles of leucine (Leu) and isoleucine (Ile), and identify those conformations that are sterically allowed versus those that are not as a function of the backbone dihedral angles ? and ψ. We compare our results with the observed distributions of side-chain dihedral angles in proteins of known structure. With the hard-sphere plus stereochemical constraint model, we obtain agreement between the model predictions and the observed side-chain dihedral angle distributions for Leu and Ile. These results quantify the extent to which local, geometrical constraints determine protein side-chain conformations.     

Properties of polyproline II, a secondary structure element implicated in protein-protein interactions

The polyproline II (PPII) conformation of protein backbone is an important secondary structure type. It is unusual in that, due to steric constraints, its main-chain hydrogen-bond donors and acceptors cannot easily be satisfied. It is unable to make local hydrogen bonds, in a manner similar to that of alpha-helices, and it cannot easily satisfy the hydrogen-bonding potential of neighboring residues in polyproline conformation in a manner analogous to beta-strands. Here we describe an analysis of polyproline conformations using the HOMSTRAD database of structurally aligned proteins. This allows us not only to determine amino acid propensities from a much larger database than previously but also to investigate conservation of amino acids in polyproline conformations, and the conservation of the conformation itself. Although proline is common in polyproline helices, helices without proline represent 46% of the total. No other amino acid appears to be greatly preferred; glycine and aromatic amino acids have low propensities for PPII. Accordingly, the hydrogen-bonding potential of PPII main-chain is mainly satisfied by water molecules and by other parts of the main-chain. Side-chain to main-chain interactions are mostly nonlocal. Interestingly, the increased number of nonsatisfied H-bond donors and acceptors (as compared with alpha-helices and beta-strands) makes PPII conformers well suited to take part in protein-protein interactions.     

New Insights into the Interdependence between Amino Acid Stereochemistry and Protein Structure

To successfully design new proteins and understand the effects of mutations in natural proteins, we must understand the geometric and physicochemical principles underlying protein structure. The side chains of amino acids in peptides and proteins adopt specific dihedral angle combinations; however, we still do not have a fundamental quantitative understanding of why some side-chain dihedral angle combinations are highly populated and others are not. Here we employ a hard-sphere plus stereochemical constraint model of dipeptide mimetics to enumerate the side-chain dihedral angles of leucine (Leu) and isoleucine (Ile), and identify those conformations that are sterically allowed versus those that are not as a function of the backbone dihedral angles ϕ and ψ. We compare our results with the observed distributions of side-chain dihedral angles in proteins of known structure. With the hard-sphere plus stereochemical constraint model, we obtain agreement between the model predictions and the observed side-chain dihedral angle distributions for Leu and Ile. These results quantify the extent to which local, geometrical constraints determine protein side-chain conformations.     

How many fold types of protein are there in nature?

Many protein structures have now been determined and reveal that protein molecules can adopt the same fold despite having very different sequences. It has been suggested that, owing to different stereochemical constraints, the number of ways that a sequence can fold may be limited. Therefore, it is reasonable to ask how many fold types exist in nature. Several groups have tackled this problem with very different results. In the present study, a novel statistical sampling approach is used to reestimate this number. The results suggest that the number of protein folds in nature is probably several hundreds. © 1996 Wiley-Liss, Inc.     

A molecular dynamics study of the C-terminal fragment of the L7/L12 ribosomal protein. Secondary structure motion in a 150 picosecond trajectory

A 150 picosecond molecular dynamics computer simulation of the C-terminal fragment of the L7/L12 ribosomal protein from Escherichia coli is reported. The molecular dynamics results are compared with the available high-resolution X-ray data in terms of atomic positions, distances and positional fluctuations. Good agreement is found between the molecular dynamics results and the X-ray data. The form and parameters of the interaction potential energy function and the procedures for deriving it are discussed. Some current misunderstandings concerning the ways of evaluating the efficiency of molecular dynamics algorithms and of application of bond-length constraints in protein simulations are cleared up. The 150 picosecond trajectory has been scanned in a search for correlated motions within and between secondary structure elements. The beta-strands have diffusional stretching modes, and uncorrelated transversal displacements. The dynamic analysis of alpha-helices shows a variety of features. The atomic fluctuations differ between the helix ends; this effect reflects long time-scale motions. Two alpha-helices, alpha A and alpha C, show diffusive longitudinal stretching modes. The third helix, alpha B, has a correlated asymmetric longitudinal stretching; the N-terminal part dominates this behaviour. Furthermore, alpha B presents a librational motion with respect to the other parts of the molecule with a frequency of approximately 5 cm-1. This motion is coupled to helix stretching. Interestingly, the regions of highly conserved residues contain the most mobile parts of the molecule.