Cooperative effects in hydrogen-bonding of protein secondary structure elements: a systematic analysis of crystal data using Secbase

A systematic analysis of the hydrogen-bonding geometry in helices and beta sheets has been performed. The distances and angles between the backbone carbonyl O and amide N atoms were correlated considering more than 1500 protein chains in crystal structures determined to a resolution better than 1.5 A. They reveal statistically significant trends in the H-bond geometry across the different secondary structural elements. The analysis has been performed using Secbase, a modular extension of Relibase (Receptor Ligand Database) which integrates information about secondary structural elements assigned to individual protein structures with the various search facilities implemented into Relibase. A comparison of the mean hydrogen-bond distances in alpha helices and 3(10) helices of increasing length shows opposing trends. Whereas in alpha helices the mean H-bond distance shrinks with increasing helix length and turn number, the corresponding mean dimension in 3(10) helices expands in a comparable series. Comparing similarly the hydrogen-bond lengths in beta sheets there is no difference to be found between the mean H-bond length in antiparallel and parallel beta sheets along the strand direction. In contrast, an interesting systematic trend appears to be given for the hydrogen bonds perpendicular to the strands bridging across an extended sheet. With increasing number of accumulated strands, which results in a growing number of back-to-back piling hydrogen bonds across the strands, a slight decrease of the mean H-bond distance is apparent in parallel beta sheets whereas such trends are obviously not given in antiparallel beta sheets. This observation suggests that cooperative effects mutually polarizing spatially well-aligned hydrogen bonds are present either in alpha helices and parallel beta sheets whereas such influences seem to be lacking in 3(10) helices and antiparallel beta sheets.     

Distributions of beta sheets in proteins with application to structure prediction

We recently developed the Rosetta algorithm for ab initio protein structure prediction, which generates protein structures from fragment libraries using simulated annealing. The scoring function in this algorithm favors the assembly of strands into sheets. However, it does not discriminate between different sheet motifs. After generating many structures using Rosetta, we found that the folding algorithm predominantly generates very local structures. We surveyed the distribution of beta-sheet motifs with two edge strands (open sheets) in a large set of non-homologous proteins. We investigated how much of that distribution can be accounted for by rules previously published in the literature, and developed a filter and a scoring method that enables us to improve protein structure prediction for beta-sheet proteins. Proteins 2002;48:85-97.     

Topological structure analysis of the protein-protein interaction network in budding yeast

Interaction detection methods have led to the discovery of thousands of interactions between proteins, and discerning relevance within large-scale data sets is important to present-day biology. Here, a spectral method derived from graph theory was introduced to uncover hidden topological structures (i.e. quasi-cliques and quasi-bipartites) of complicated protein-protein interaction networks. Our analyses suggest that these hidden topological structures consist of biologically relevant functional groups. This result motivates a new method to predict the function of uncharacterized proteins based on the classification of known proteins within topological structures. Using this spectral analysis method, 48 quasi-cliques and six quasi-bipartites were isolated from a network involving 11,855 interactions among 2617 proteins in budding yeast, and 76 uncharacterized proteins were assigned functions.     

Graph identification of proteins in tomograms (GRIP-Tomo)

In this study, we present a method of pattern mining based on network theory that enables the identification of protein structures or complexes from synthetic volume densities, without the knowledge of predefined templates or human biases for refinement. We hypothesized that the topological connectivity of protein structures is invariant, and they are distinctive for the purpose of protein identification from distorted data presented in volume densities. Three-dimensional densities of a protein or a complex from simulated tomographic volumes were transformed into mathematical graphs as observables. We systematically introduced data distortion or defects such as missing fullness of data, the tumbling effect, and the missing wedge effect into the simulated volumes, and varied the distance cutoffs in pixels to capture the varying connectivity between the density cluster centroids in the presence of defects. A similarity score between the graphs from the simulated volumes and the graphs transformed from the physical protein structures in point data was calculated by comparing their network theory order parameters including node degrees, betweenness centrality, and graph densities. By capturing the essential topological features defining the heterogeneous morphologies of a network, we were able to accurately identify proteins and homo-multimeric complexes from 10 topologically distinctive samples without realistic noise added. Our approach empowers future developments of tomogram processing by providing pattern mining with interpretability, to enable the classification of single-domain protein native topologies as well as distinct single-domain proteins from multimeric complexes within noisy volumes.     

Computer simulations aimed at structure prediction of supersecondary motifs in proteins

It is well established that protein structures are more conserved than protein sequences. One-third of all known protein structures can be classified into ten protein folds, which themselves are composed mainly of alpha-helical hairpin, beta hairpin, and betaalphabeta supersecondary structural elements. In this study, we explore the ability of a recent Monte Carlo-based procedure to generate the 3D structures of eight polypeptides that correspond to units of supersecondary structure and three-stranded antiparallel beta sheet. Starting from extended or misfolded compact conformations, all Monte Carlo simulations show significant success in predicting the native topology using a simplified chain representation and an energy model optimized on other structures. Preliminary results on model peptides from nucleotide binding proteins suggest that this simple protein folding model can help clarify the relation between sequence and topology.     

Prediction of protein coarse contact maps

Prediction of topological representations of proteins that are geometrically invariants can contribute towards the solution of fundamental open problems in structural genomics like folding. In this paper we focus on coarse grained protein contact maps, a representation that describes the spatial neighborhood relation between secondary structure elements such as helices, beta sheets, and random coils. Our methodology is based on searching the graph space. The search algorithm is guided by an adaptive evaluation function computed by a specialized noncausal recursive connectionist architecture. The neural network is trained using candidate graphs generated during examples of successful searches. Our results demonstrate the viability of the approach for predicting coarse contact maps.     

An amino acid code for β‐sheet packing structure

To understand the relationship between protein sequence and structure, this work extends the knob‐socket model in an investigation of β‐sheet packing. Over a comprehensive set of β‐sheet folds, the contacts between residues were used to identify packing cliques: sets of residues that all contact each other. These packing cliques were then classified based on size and contact order. From this analysis, the two types of four‐residue packing cliques necessary to describe β‐sheet packing were characterized. Both occur between two adjacent hydrogen bonded β‐strands. First, defining the secondary structure packing within β‐sheets, the combined socket or XY:HG pocket consists of four residues i, i+2 on one strand and j, j+2 on the other. Second, characterizing the tertiary packing between β‐sheets, the knob‐socket XY:H+B consists of a three‐residue XY:H socket (i, i+2 on one strand and j on the other) packed against a knob B residue (residue k distant in sequence). Depending on the packing depth of the knob B residue, two types of knob‐sockets are found: side‐chain and main‐chain sockets. The amino acid composition of the pockets and knob‐sockets reveal the sequence specificity of β‐sheet packing. For β‐sheet formation, the XY:HG pocket clearly shows sequence specificity of amino acids. For tertiary packing, the XY:H+B side‐chain and main‐chain sockets exhibit distinct amino acid preferences at each position. These relationships define an amino acid code for β‐sheet structure and provide an intuitive topological mapping of β‐sheet packing. Proteins 2014; 82:2128–2140. © 2014 Wiley Periodicals, Inc.     

Protein-protein recognition and the association of immunoglobulin constant domains

Four immunoglobulin constant domain interfaces were analysed to determine how the light and heavy chains of the domains recognize each other and secondarily, why the packing of the beta sheets is different from that normally observed in the interior of domains of beta sheet proteins. It was found that two bulky groups in the centre of each of the domain interfaces prevent the alignment of sheet axes by intercalating with residues in the sheet opposite. This intercalation forms the basis of constant domain light chain-heavy chain recognition.     

BuildBeta—A system for automatically constructing beta sheets

We describe a method that can thoroughly sample a protein conformational space given the protein primary sequence of amino acids and secondary structure predictions. Specifically, we target proteins with β‐sheets because they are particularly challenging for ab initio protein structure prediction because of the complexity of sampling long‐range strand pairings. Using some basic packing principles, inverse kinematics (IK), and β‐pairing scores, this method creates all possible β‐sheet arrangements including those that have the correct packing of β‐strands. It uses the IK algorithms of ProteinShop to move α‐helices and β‐strands as rigid bodies by rotating the dihedral angles in the coil regions. Our results show that our approach produces structures that are within 4–6 Å RMSD of the native one regardless of the protein size and β‐sheet topology although this number may increase if the protein has long loops or complex α‐helical regions. Proteins 2010. © Published 2009 Wiley‐Liss, Inc.     

Structures of single‐layer β‐sheet proteins evolved from β‐hairpin repeats

Free‐standing single‐layer β‐sheets are extremely rare in naturally occurring proteins, even though β‐sheet motifs are ubiquitous. Here we report the crystal structures of three homologous, single‐layer, anti‐parallel β‐sheet proteins, comprised of three or four twisted β‐hairpin repeats. The structures reveal that, in addition to the hydrogen bond network characteristic of β‐sheets, additional hydrophobic interactions mediated by small clusters of residues adjacent to the turns likely play a significant role in the structural stability and compensate for the lack of a compact hydrophobic core. These structures enabled identification of a family of secreted proteins that are broadly distributed in bacteria from the human gut microbiome and are putatively involved in the metabolism of complex carbohydrates. A conserved surface patch, rich in solvent‐exposed tyrosine residues, was identified on the concave surface of the β‐sheet. These new modular single‐layer β‐sheet proteins may serve as a new model system for studying folding and design of β‐rich proteins.     

A new protein folding algorithm based on hydrophobic compactness: Rigid Unconnected Secondary Structure Iterative Assembly (RUSSIA). I: Methodology

We present an algorithm that is able to propose compact models of protein 3D structures, only starting from the prediction of the nature and length of regular secondary structures. Helices are modeled by cylinders and sheets by helicoid surfaces, all strands of a sheet being considered as a single block. It means that relative topology of the strands inside one sheet is a prerequisite. Loops are only considered as constraints, given by the maximal distance between their Calpha extremities according to their sequence length. Unconnected regular secondary structures are reduced to a single point, the center of their hydrophobic faces. These centers are then repeatedly moved in order to obtain a compact hydrophobic core. To prevent secondary structures from interpenetrating, a repulsive term is introduced in the function whose minimization leads to the compact structure. This RUSSIA (Rigid Unconnected Secondary Structure Assembly) algorithm has the advantage of relying on a small number of variables and therefore many initial conformations can be tested. Flexibility is produced in the following way: helices or sheets are allowed to rotate around the direction leading to the center of the model; residues in a sheet can slide along the main direction of the strand where they are embedded. RUSSIA is fast and simple and it produces on a test set several neighbor good models with an r.m.s. to the native structures in the range 1.4-3.7 A. These models can be further treated by statistical potentials used in threading approaches in order to detect the best candidate. The limits of the present method are the following: small proteins with few secondary structures are excluded; multi domain proteins must be split into several compact globular domains from their sequences; sheets of more than five strands and completely buried helices are not treated. In this first paper the algorithm is developed and in Part II, which follows, some applications are presented and the program is evaluated.     

Conformation of alpha zeins in solid state by Fourier transform IR

The major maize storage proteins (alpha zeins) are deposited as an insoluble mass in the protein bodies of the endosperm. Because they are insoluble in water, most structural studies are performed in alcohol solutions. To solve the question raised by several authors about denaturation of the alpha zein structure by alcohol, we analyze the secondary structure of alpha zeins prepared with and without solubilization in alcohol (corn gluten meal and protein bodies with high concentrations of alpha zeins and traces of beta zeins). The secondary structures of alpha zeins are analyzed in the solid state by Fourier transform IR spectroscopy (FTIR) in KBr pellets and solid-state 13C-NMR spectroscopy. The proportion of secondary structures obtained by FTIR of alpha zeins prepared with and without solubilization in alcohol yield almost identical proportions of alpha helices and beta sheets. The proportion of alpha helices (43%) agrees with that measured by circular dichroism in an alcohol solution. However, the proportion of beta sheets (28%) is higher than the one measured by the same technique. Gluten and protein body samples with high beta zein content showed higher beta sheet and lower alpha helix proportions than that obtained for alpha zein preparations. The solid-state 13C-NMR spectra show the carbonyl peak for the alpha zeins at delta 176 and for the sample rich in beta zeins at delta 172, which demonstrates the presence of a high content of alpha helices and beta sheets, respectively. These results indicate that alcohol solubilization does not affect the conformation of alpha zeins, validating the secondary structure measurements in solution.     

Fourier transform infared spectroscopy investigation of protein conformation in spray-dried protein/trehalose powders

Spray drying is a way to generate protein solids (powders), which is also true for lyophilization. Sugars are used to protect proteins from conformational changes and chemical degradations arising from drying processes and storage conditions such as the humidity. The influence of trehalose and humidity on the conformation and hydration of spray-dried recombinant human granolucyte colony stimulating factor (rhG-CSF) and recombinant consensus interferon-alpha (rConIFN) was investigated using Fourier transform IR spectroscopy. The spectral analysis of spray-dried powders in the amide I region demonstrated that trehalose stabilized the alpha-helical conformation of both rhG-CSF and rConIFN proteins. Exposure of the pure protein powders to 33% relative humidity (RH) resulted in the formation of beta sheets and loss of turns but no change in alpha-helical structure. Trehalose reduced the magnitude of the changes in beta sheets and turns. Exposure of the pure protein powders to 75% RH resulted in the loss of alpha-helical conformation with a corresponding increase in beta structures (beta sheets and turns). Trehalose did not protect proteins from the loss of alpha-helical structures, but it reduced the formation of antiparallel beta sheets. Hydrogen-deuterium exchange (H-D exchange) was used to further characterize these hydration-induced conformational changes. At 33% RH the percent exchange of the protein decreased with increasing trehalose content, indicating a greater protection of the protein from H-D exchange by a higher concentration of trehalose. Such protection correlates with decreased conformational changes of the protein by trehalose at this humidity. At 75% RH the degree of H-D exchange of the protein was insensitive to the powder composition in all powders. Surprisingly, the H-D exchange of trehalose was low at about 20-25%, which was nearly independent of the protein/trehalose ratio and humidity, indicating that the exchangeable protons on trehalose molecules are highly protected in protein-containing powders. The observed protein hydration is related to the effect of trehalose on the conformational changes of the protein under humidity.     

Relative stability of protein structures determined by X-ray crystallography or NMR spectroscopy: a molecular dynamics simulation study

The relative stability of protein structures determined by either X-ray crystallography or nuclear magnetic resonance (NMR) spectroscopy has been investigated by using molecular dynamics simulation techniques. Published structures of 34 proteins containing between 50 and 100 residues have been evaluated. The proteins selected represent a mixture of secondary structure types including all alpha, all beta, and alpha/beta. The proteins selected do not contain cysteine-cysteine bridges. In addition, any crystallographic waters, metal ions, cofactors, or bound ligands were removed before the systems were simulated. The stability of the structures was evaluated by simulating, under identical conditions, each of the proteins for at least 5 ns in explicit solvent. It is found that not only do NMR-derived structures have, on average, higher internal strain than structures determined by X-ray crystallography but that a significant proportion of the structures are unstable and rapidly diverge in simulations.     

Layered structures of hydrophobically modified chitosan derivatives

Two series of chitosan derivatives, N-aliphatic acyl chitosans (C_n-CS) and N-aliphatic-O-dicinnamoyl-chitosans (C_n-CinCS) with n=2, 4, 8, 12 and 18, were prepared through regioselective reactions. The solid state structures were studied by FTIR and X-ray diffraction techniques. Two different types of layered structures were found to exist in the powder samples of C_n-CS and C_n-CinCS. The C_n-CS series crystallized into a sheet-type structure, in which all the flexible side chains lied down in the sheet in a direction normal to the backbones and partially interdigitated with each other. For the C_n-CinCS series, the polar backbones were stacked into sheets, and the flexible side chains occupied the space between sheets in a direction inclined to the backbones. The relationship between the structures and the solubilities are discussed.     

Looking at proteins: representations, folding, packing, and design. Biophysical Society National Lecture, 1992.

Looking at proteins is an active process of interpretation and selection, emphasizing some features and deleting others. Multiple representations are needed, for such purposes as showing motions or conveying both the chain connectivity and the three-dimensional shape simultaneously. In studying and comparing protein structures, ideas are suggested about the determinants of tertiary structure and of folding (e.g., that Greek key beta barrels may fold up two strands at a time). The design and synthesis of new proteins "from scratch" provides a route toward the experimental testing of such ideas. It has also been a fruitful new perspective from which to look at structures, requiring such things as statistics on very narrowly defined structural categories and explicit attention to "negative design" criteria that actively block unwanted alternatives (e.g., reverse topology of a helix bundle, or edge-to-edge aggregation of beta sheets). Recently, the field of protein design has produced a rather unexpected general result: apparently we do indeed know enough to successfully design proteins that fold into approximately correct structures, but not enough to design unique, native-like structures. The degree of order varies considerably, but even the best designed material shows multiple conformations by NMR, more similar to a "molten globule" folding intermediate than to a well ordered native tertiary structure. In response to this conclusion, we are now working on systems that test useful questions with approximate structures (such as determining which factors most influence the choice of helix-bundle topology) and also analyzing how natural proteins achieve unique core conformations (e.g., for side chains on the interior side of a beta sheet, illustrated in the kinemages).     

Mechanisms determining the morphology of the peripheral ER

The endoplasmic reticulum (ER) consists of the nuclear envelope and a peripheral network of tubules and membrane sheets. The tubules are shaped by the curvature-stabilizing proteins reticulons and DP1/Yop1p, but how the sheets are formed is unclear. Here, we identify several sheet-enriched membrane proteins in the mammalian ER, including proteins that translocate and modify newly synthesized polypeptides, as well as coiled-coil membrane proteins that are highly upregulated in cells with proliferated ER sheets, all of which are localized by membrane-bound polysomes. These results indicate that sheets and tubules correspond to rough and smooth ER, respectively. One of the coiled-coil proteins, Climp63, serves as a "luminal ER spacer" and forms sheets when overexpressed. More universally, however, sheet formation appears to involve the reticulons and DP1/Yop1p, which localize to sheet edges and whose abundance determines the ratio of sheets to tubules. These proteins may generate sheets by stabilizing the high curvature of edges.     

Systematic application of two-dimensional 1H nuclear-magnetic-resonance techniques for studies of proteins. 2. Combined use of correlated spectroscopy and nuclear Overhauser spectroscopy for sequential assignments of backbone resonances and elucidation of polypeptide secondary structures

This paper describes a new nuclear magnetic resonance approach for the determination of secondary structure in globular proteins. To illustrate the practical application of the new procedure, two-dimensional correlated spectroscopy and two-dimensional nuclear Overhauser enhancement spectroscopy were used to obtain individual assignments for all the backbone protons of the beta-sheet secondary structures in the basic pancreatic trypsin inhibitor. First, combined connectivity diagrams of these two methods recorded in both 2H2O solution and H2O solution of the inhibitor were employed to obtain sequential, individual resonance assignments for the separate strands in the beta sheet. Second, a 2D nuclear Overhauser enhancement spectrum recorded with a long mixing time was used to determine how the separate, extended polypeptide strands are linked by hydrogen bonds in the sheet structures. By combination of these results with the identifications of the amino acid side-chain resonances described in the preceding paper, the beta-sheet structures can, without reference to data on the spatial structure obtained with other techniques, be localized in the amino acid sequence. This investigation confirms results on limited regions of the beta sheet in the inhibitor obtained previously with one-dimensional nuclear magnetic resonance experiments and demonstrates that the entire beta-sheet structure seen in single crystals of the inhibitor is preserved in aqueous solution.     

Flexibility of beta-sheets: principal component analysis of database protein structures

Protein folds are built primarily from the packing together of two types of structures: alpha-helices and beta-sheets. Neither structure is rigid, and the flexibility of helices and sheets is often important in determining the final fold (e.g., coiled coils and beta-barrels). Recent work has quantified the flexibility of alpha-helices using a principal component analysis (PCA) of database helical structures (J. Mol. Bio. 2003, 327, pp. 229-237). Here, we extend the analysis to beta-sheet flexibility using PCA on a database of beta-sheet structures. For sheets of varying dimension and geometry, we find two dominant modes of flexibility: twist and bend. The distributions of amplitudes for these modes are found to be Gaussian and independent, suggesting that the PCA twist and bend modes can be identified as the soft elastic normal modes of sheets. We consider the scaling of mode eigenvalues with sheet size and find that parallel beta-sheets are more rigid than antiparallel sheets over the entire range studied. Finally, we discuss the application of our PCA results to modeling and design of beta-sheet proteins.     

Comparison of solvent-inaccessible cores of homologous proteins: definitions useful for protein modelling

The three-dimensional structures of 41 homologous proteins (belonging to eight families) were compared by pairwise superposition. A subset of 'core' residues was defined as those whose side chains have less than 7% of their surface exposed to solvent. This subset has significantly higher sequence identity and lower root mean square (RMS) alpha carbon separation than for all topologically equivalent residues in the structure, when members of a protein family are superposed. For such superpositions the relationship between RMS distance and percentage sequence identity of this subset of residues is similar to that for all equivalent residues, although some variation is observed between families of proteins which are predominantly beta sheet and those which are mainly alpha helix. The definition of a structurally more conserved core may be useful in model building proteins from an homologous family. The RMS differences of coordinates of structures of proteins with identical sequences are found to be related to the resolutions of the structures.