Combinations of turns in proteins.

We observed that beta- and gamma-turns in protein structure may be associated as peptides representing combinations of turns that span between nine and 26 amino acid residues along the polypeptide backbone chain and often correspond to loops in the protein structure. Around 475 peptides resulted from the analysis of a non-redundant data set corresponding to 248 protein crystal structures selected from the Protein Data Bank. Nearly 40% protein chains are associated with two or more peptides and the peptides with nine and 10 amino acid residues are more frequent. A maximum of four distinct peptides varying in number of amino acid residues were observed in at least 10 proteins along the same protein chain. Nearly 80% peptides comprise type IV beta-turns that are associated with irregular dihedral angle values suggesting this may be important for the conformational diversity associated with the loops in proteins. In general, predominant interactions that possibly stabilize these peptides involve main-chain and side-chain interactions with solvent, in addition to hydrogen bond, salt-bridge and non-bonded interactions. Majority of the peptides were observed in hydrolase, oxidoreductase, transferase, serine proteinase/inhibitor complex, electron transport/electron transfer and lyase proteins.     

A detailed consideration of a principal domain of vertebrate fibrinogen and its relatives.

Vertebrate fibrinogen is a complex multidomained protein, the structure of which has been inferred mainly from electron microscopy and amino acid sequence studies. Among its most prominent features are two terminal globules, moieties that are mostly composed of the carboxyl-terminal two-thirds of the beta and gamma chains. Sequences homologous to the latter segments are found in several other animal proteins, always as the carboxyl-terminal contributions. An alignment of 15 amino acid sequences from various fibrinogens and related proteins has been used to make judgments about secondary structure. The nature of amino acids at each position in the alignment was used to distinguish alpha helices and beta structure on the one hand from loops and turns on the other, and the resulting assignments compared with predictions of secondary structure by other methods. Additionally, constraints imposed by the locations of cystines, carbohydrate attachment residues, and proteinase-sensitive points provided further insights into the general organization of the postulated secondary structures. Other ancillary data, including the effects of bound calcium and the locations of labeled or variant residues, were also considered. An intriguing similarity to a portion of the recently reported structure of a calcium-dependent lectin is noted.     

The conformation of glucagon: predictions and consequences.

It is proposed that glucagon, a polypeptide hormone, is delicately balanced between two major conformational states. Utilizing a new predictive model [Chou, P.Y., and Fasman, G.D. (1974), Biochemistry 13, 222] which considers all the conformational states in proteins (helix, beta sheet, random coil, and beta turns), the secondary structural regions of glucagon are computed herein. The conformational sensitivity of glucagon may be due to residues 19-27 which have both alpha-helical potential (mean value of Palpha = 1.19) as well as beta-sheet potential (mean value of Pbeta = 1.25). Two conformational states are predicted for glucagon. In predicted form (a), residues 5-10 form a beta-sheet region while residues 19-27 form an alpha-helical region (31% alpha, 21% beta) agreeing well with the circular dichroism (CD) spectra of glucagon. The similarity in the CD spectra of glucagon and insulin further suggests the presence of beta structure in glucagon, since X-ray analysis of insulin showed 24% beta sheet. In predicted form (b), both regions, residues 5-10 and residues 19-27, are beta sheets sheets (0% alpha, 52% beta) in agreement with the infrared spectral evidence that glucagon gels and fibrils have a predominant beta-sheet conformation. Since three reverse beta turns are predicted at residues 2-5, 10-13, and 15-18, glucagon may possess tertiary structure in agreement with viscosity and tritium-hydrogen exchange experiments. A proposal is offered concerning an induced alpha yields beta transition at residues 22-27 in glucagon during receptor site binding. Amino acid substitutions are proposed which should disrupt the beta sheets of glucagon with concomitant loss of biological activity. The experimental findings that glucagon aggregates to form dimers, trimers, and hexamers can be explained in terms of beta-sheet interactions as outlined in the present predictive model. Thus the conflicting conclusions of previous workers, concerning the conformation of glucagon in different environments, can be rationalized by the suggested conformational transition occurring within the molecule.     

Tight Turns of Outer Membrane Proteins: An Analysis of Sequence,Structure, and Hydrogen Bonding

As a structural class, tight turns can control molecular recognition, enzymatic activity, and nucleation of folding. They have been extensively characterized in soluble proteins but have not been characterized in outer membrane proteins (OMPs), where they also support critical functions. We clustered the 4 to 6 residue tight turns of 110 OMPs to characterize the phi/psi angles, sequence, and hydrogen bonding of these structures. We find significant differences between reports of soluble protein tight turns and OMP tight turns. Since OMP strands are less twisted than soluble strands, they favor different turn structures types. Moreover, the membrane localization of OMPs yields different sequence hallmarks for their tight turns relative to soluble protein turns. We also characterize the differences in phi/psi angles, sequence, and hydrogen bonding between OMP extracellular loops and OMP periplasmic turns. As previously noted, the extracellular loops tend to be much longer than the periplasmic turns. We find that this difference in length is due to the broader distribution of lengths of the extracellular loops not a large difference in the median length. Extracellular loops also tend to have more charged residues as predicted by the charge-out rule. Finally, in all OMP tight turns, hydrogen bonding between the side chain and backbone 2 to 4 residues away from that side chain plays an important role. These bonds preferentially use an Asp, Asn, Ser, or Thr residue in a beta or pro phi/psi conformation. We anticipate that this study will be applicable to future design and structure prediction of OMPs.     

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.     

Structural analysis of the PsbQ protein of photosystem II by Fourier transform infrared and circular dichroic spectroscopy and by bioinformatic methods

The structure of PsbQ, one of the three main extrinsic proteins associated with the oxygen-evolving complex (OEC) of higher plants and green algae, is examined by Fourier transform infrared (FTIR) and circular dichroic (CD) spectroscopy and by computational structural prediction methods. This protein, together with two other lumenally bound extrinsic proteins, PsbO and PsbP, is essential for the stability and full activity of the OEC in plants. The FTIR spectra obtained in both H(2)O and D(2)O suggest a mainly alpha-helix structure on the basis of the relative areas of the constituents of the amide I and I' bands. The FTIR quantitative analyses indicate that PsbQ contains about 53% alpha-helix, 7% turns, 14% nonordered structure, and 24% beta-strand plus other beta-type extended structures. CD analyses indicate that PsbQ is a mainly alpha-helix protein (about 64%), presenting a small percentage assigned to beta-strand ( approximately 7%) and a larger amount assigned to turns and nonregular structures ( approximately 29%). Independent of the spectroscopic analyses, computational methods for protein structure prediction of PsbQ were utilized. First, a multiple alignment of 12 sequences of PsbQ was obtained after an extensive search in the public databases for protein and EST sequences. Based on this alignment, computational prediction of the secondary structure and the solvent accessibility suggest the presence of two different structural domains in PsbQ: a major C-terminal domain containing four alpha-helices and a minor N-terminal domain with a poorly defined secondary structure enriched in proline and glycine residues. The search for PsbQ analogues by fold recognition methods, not based on the secondary structure, also indicates that PsbQ is a four alpha-helix protein, most probably folding as an up-down bundle. The results obtained by both the spectroscopic and computational methods are in agreement, all indicating that PsbQ is mainly an alpha protein, and show the value of using both methodologies for protein structure investigation.     

Electrostatic interactions at the C-terminal domain of nucleoplasmin modulate its chromatin decondensation activity

The chromatin decondensation activity, thermal stability, and secondary structure of recombinant nucleoplasmin, of two deletion mutants, and of the protein isolated from Xenopus oocytes have been characterized. As previously reported, the chromatin decondensation activity of recombinant, unphosphorylated nucleoplasmin is almost negligible. Our data show that deletion of 50 residues at the C-terminal domain of the protein, containing the positively charged nuclear localization sequence, activates its chromatin decondensation ability and decreases its stability. Interestingly, both the decondensation activity and thermal stability of this deletion mutant resemble those of the phosphorylated protein isolated from Xenopus oocytes. Deletion of 80 residues at the C-terminal domain, containing the above-mentioned positively charged region and a poly(Glu) tract, inactivates the protein and increases its thermal stability. These findings, along with the effect of salt on the thermal stability of these proteins, suggest that electrostatic interactions between the positive nuclear localization sequence and the poly(Glu) tract, at the C-terminal domain, modulate protein activity and stability.     

Use of secondary structural information and C alpha-C alpha distance restraints to model protein structures with MODELLER

Protein secondary structure predictions and amino acid long range contact map predictions from primary sequence of proteins have been explored to aid in modelling protein tertiary structures. In order to evaluate the usefulness of secondary structure and 3D-residue contact prediction methods to model protein structures we have used the known Q3 (alpha-helix,beta-strands and irregular turns/loops) secondary structure information, along with residue-residue contact information as restraints for MODELLER. We present here results of our modelling studies on 30 best resolved single domain protein structures of varied lengths. The results shows that it is very difficult to obtain useful models even with 100% accurate secondary structure predictions and accurate residue contact predictions for up to 30% of residues in a sequence. The best models that we obtained for proteins of lengths 37, 70, 118, 136 and 193 amino acid residues are of RMSDs 4.17, 5.27, 9.12, 7.89 and 9.69,respectively. The results show that one can obtain better models for the proteins which have high percent of alpha-helix content. This analysis further shows that MODELLER restrain optimization program can be useful only if we have truly homologous structure(s) as a template where it derives numerous restraints, almost identical to the templates used. This analysis also clearly indicates that even if we satisfy several true residue-residue contact distances, up to 30%of their sequence length with fully known secondary structural information, we end up predicting model structures much distant from their corresponding native structures.     

Logical analysis of the mechanism of protein folding II. The nucleation process

Regions of secondary structure are predicted, without using information about the conformation of the protein itself, and compared with crystallographic assignments for seven proteins of recently published sequence and conformation (Table 1). It is observed in Table 3 that the prediction of helices is good (78.7% for %cor.ass.3), except for proteins having large antiparallel pleated sheets, and the prediction of β-structure is quite good (51.2% for %cor.ass.3) except for helix-rich proteins.The prediction of secondary structure from sequence, and a survey of all protein structures analysed so far by X-ray crystallography, suggest that nuceleation starts in almost all cases from interactions in the medium range between the regions having helical potential (α-candidate) and β-structural potential (β-candidate), which are very close to each other but separated by at least three hydrophilic or neutral residues in four consecutive residues on the polypeptide chain. Predictability of loops or turns is enhanced to 71.3% (%cor.ass.2) from 64.4% by taking into account the contiguous α-β interactions. Such a medium-range interaction is called here a probable nucleus. There are a lot of nuclei in large proteins such as carboxypeptidase Aα, while there exists at least one in small proteins like the trypsin inhibitor, Moreover, such an interaction could be a transitionary state towards a helix-rich protein, and towards a helix-deficient protein having a large antiparallel pleated sheet β-structure as well.The analysis of the relation between probable nuclei with regard to their mutual spatial proximity strongly suggests that the topological pathway of the polypeptide chain in three-dimensional space might be decided by the long-range interactions between an α-candidate and a β-candidate. An empirical rule is observed that almost all parallel pleated sheets are accompanied by helices in their neighbourhood. An accumulation of chemical facts, such as complementation experiments, combinations of disulphide bonds, etc., seems also to be elucidated by the proposed mechanism of protein folding.     

A machine learning approach for the prediction of protein surface loop flexibility

Proteins often undergo conformational changes when binding to each other. A major fraction of backbone conformational changes involves motion on the protein surface, particularly in loops. Accounting for the motion of protein surface loops represents a challenge for protein-protein docking algorithms. A first step in addressing this challenge is to distinguish protein surface loops that are likely to undergo backbone conformational changes upon protein-protein binding (mobile loops) from those that are not (stationary loops). In this study, we developed a machine learning strategy based on support vector machines (SVMs). Our SVM uses three features of loop residues in the unbound protein structures-Ramachandran angles, crystallographic B-factors, and relative accessible surface area-to distinguish mobile loops from stationary ones. This method yields an average prediction accuracy of 75.3% compared with a random prediction accuracy of 50%, and an average of 0.79 area under the receiver operating characteristic (ROC) curve using cross-validation. Testing the method on an independent dataset, we obtained a prediction accuracy of 70.5%. Finally, we applied the method to 11 complexes that involve members from the Ras superfamily and achieved prediction accuracy of 92.8% for the Ras superfamily proteins and 74.4% for their binding partners.     

Spectroscopic methods for analysis of protein secondary structure

Several methods for determination of the secondary structure of proteins by spectroscopic measurements are reviewed. Circular dichroism (CD) spectroscopy provides rapid determinations of protein secondary structure with dilute solutions and a way to rapidly assess conformational changes resulting from addition of ligands. Both CD and Raman spectroscopies are particularly useful for measurements over a range of temperatures. Infrared (IR) and Raman spectroscopy require only small volumes of protein solution. The frequencies of amide bands are analyzed to determine the distribution of secondary structures in proteins. NMR chemical shifts may also be used to determine the positions of secondary structure within the primary sequence of a protein. However, the chemical shifts must first be assigned to particular residues, making the technique considerably slower than the optical methods. These data, together with sophisticated molecular modeling techniques, allow for refinement of protein structural models as well as rapid assessment of conformational changes resulting from ligand binding or macromolecular interactions. A selected number of examples are given to illustrate the power of the techniques in applications of biological interest.     

Prediction of secondary structure by evolutionary comparison: application to the alpha subunit of tryptophan synthase

The amino acid sequences of the a subunits of tryptophan synthase from ten different microorganisms were aligned by standard procedures. The alpha helices, beta strands and turns of each sequence were predicted separately by two standard prediction algorithms and averaged at homologous sequence positions. Additional evidence for conserved secondary structure was derived from profiles of average hydropathy and chain flexibility values, leading to a joint prediction. There is good agreement between (1) predicted beta strands, maximal hydropathy and minimal flexibility, and (2) predicted loops, great chain flexibility, and protein segments that accept insertions of various lengths in individual sequences. The a subunit is predicted to have eight repeated beta-loop-alpha-loop motifs with an extra N-terminal alpha helix and an intercalated segment of highly conserved residues. This pattern suggests that the territory structure of the a subunit is an eightfold alpha/beta barrel. The distribution of conserved amino acid residues and published data on limited proteolysis, chemical modification, and mutagenesis are consistent with the alpha/beta barrel structure. Both the active site of the a subunit and the combining site for the beta 2 subunit are at the end of the barrel formed by the carboxyl-termini of the beta strands.     

Conformational analysis and clustering of short and medium size loops connecting regular secondary structures: a database for modeling and prediction.

Loops are regions of nonrepetitive conformation connecting regular secondary structures. We identified 2,024 loops of one to eight residues in length, with acceptable main-chain bond lengths and peptide bond angles, from a database of 223 protein and protein-domain structures. Each loop is characterized by its sequence, main-chain conformation, and relative disposition of its bounding secondary structures as described by the separation between the tips of their axes and the angle between them. Loops, grouped according to their length and type of their bounding secondary structures, were superposed and clustered into 161 conformational classes, corresponding to 63% of all loops. Of these, 109 (51% of the loops) were populated by at least four nonhomologous loops or four loops sharing a low sequence identity. Another 52 classes, including 12% of the loops, were populated by at least three loops of low sequence similarity from three or fewer nonhomologous groups. Loop class suprafamilies resulting from variations in the termini of secondary structures are discussed in this article. Most previously described loop conformations were found among the classes. New classes included a 2:4 type IV hairpin, a helix-capping loop, and a loop that mediates dinucleotide-binding. The relative disposition of bounding secondary structures varies among loop classes, with some classes such as beta-hairpins being very restrictive. For each class, sequence preferences as key residues were identified; those most frequently at these conserved positions than in proteins were Gly, Asp, Pro, Phe, and Cys. Most of these residues are involved in stabilizing loop conformation, often through a positive phi conformation or secondary structure capping. Identification of helix-capping residues and beta-breakers among the highly conserved positions supported our decision to group loops according to their bounding secondary structures. Several of the identified loop classes were associated with specific functions, and all of the member loops had the same function; key residues were conserved for this purpose, as is the case for the parvalbumin-like calcium-binding loops. A significant number, but not all, of the member loops of other loop classes had the same function, as is the case for the helix-turn-helix DNA-binding loops. This article provides a systematic and coherent conformational classification of loops, covering a broad range of lengths and all four combinations of bounding secondary structure types, and supplies a useful basis for modelling of loop conformations where the bounding secondary structures are known or reliably predicted.     

Use of secondary structural information and C<Emphasis Type="Italic">α</Emphasis>-C<Emphasis Type="Italic">α</Emphasis> distance restraints to model protein structures with MODELLER

Protein secondary structure predictions and amino acid long range contact map predictions from primary sequence of proteins have been explored to aid in modelling protein tertiary structures. In order to evaluate the usefulness of secondary structure and 3D-residue contact prediction methods to model protein structures we have used the known Q3 (alpha-helix, beta-strands and irregular turns/loops) secondary structure information, along with residue-residue contact information as restraints for MODELLER. We present here results of our modelling studies on 30 best resolved single domain protein structures of varied lengths. The results shows that it is very difficult to obtain useful models even with 100% accurate secondary structure predictions and accurate residue contact predictions for up to 30% of residues in a sequence. The best models that we obtained for proteins of lengths 37, 70, 118, 136 and 193 amino acid residues are of RMSDs 4.17, 5.27, 9.12, 7.89 and 9.69, respectively. The results show that one can obtain better models for the proteins which have high percent of alpha-helix content. This analysis further shows that MODELLER restrain optimization program can be useful only if we have truly homologous structure(s) as a template where it derives numerous restraints, almost identical to the templates used. This analysis also clearly indicates that even if we satisfy several true residue-residue contact distances, up to 30% of their sequence length with fully known secondary structural information, we end up predicting model structures much distant from their corresponding native structures.     

Analysis of the stability and function of nucleoplasmin through cysteine mutants

Xenopus laevis nucleoplasmin is a pentameric nuclear chaperone. The relation between the structure and the multifunctional aspects of the molecule has not yet been clearly established. In the course of analysing a C-terminally His-tagged recombinant version of the region equivalent to the trypsin resistant core (r-NP142) of the molecule, we found that this domain exhibited a substantially decreased oligomerization potential. To better understand the role of the three cysteines of nucleoplasmin on its pentameric functional structure, we have selectively mutated these residues to serine and generated three mutants (C15S, C35S, and C45S) both for the complete recombinant nucleoplasmin (r-NP) and the truncated r-NP142 non-tagged forms. We demonstrate that there are no disulphide bridges stabilizing either the monomer or the pentamer. Neither C15S nor C35S has any structural effects, while the mutation C45S abolishes the ability of r-NP142 to pentamerize. This structural impairment suggests that hydrophobic interactions of Cys 45 are critical for the stability of the protein. Our studies allow to analyse for the first time the structural and functional properties of nucleoplasmin in its monomeric form.     

Protein loop structure prediction with flexible stem geometries

The structure prediction of loops with flexible stem residues is addressed in this article. While the secondary structure of the stem residues is assumed to be known, the geometry of the protein into which the loop must fit is considered to be unknown in our methodology. As a consequence, the compatibility of the loop with the remainder of the protein is not used as a criterion to reject loop decoys. The loop structure prediction with flexible stems is more difficult than fitting loops into a known protein structure in that a larger conformational space has to be covered. The main focus of the study is to assess the precision of loop structure prediction if no information on the protein geometry is available. The proposed approach is based on (1) dihedral angle sampling, (2) structure optimization by energy minimization with a physically based energy function, (3) clustering, and (4) a comparison of strategies for the selection of loops identified in (3). Steps (1) and (2) have similarities to previous approaches to loop structure prediction with fixed stems. Step (3) is based on a new iterative approach to clustering that is tailored for the loop structure prediction problem with flexible stems. In this new approach, clustering is not only used to identify conformers that are likely to be close to the native structure, but clustering is also employed to identify far-from-native decoys. By discarding these decoys iteratively, the overall quality of the ensemble and the loop structure prediction is improved. Step (4) provides a comparative study of criteria for loop selection based on energy, colony energy, cluster density, and a hybrid criterion introduced here. The proposed method is tested on a large set of 3215 loops from proteins in the Pdb-Select25 set and to 179 loops from proteins from the Casp6 experiment.     

Protein structure prediction

Current methods developed for predicting protein structure are reviewed. The most widely used algorithms of Chou and Fasman and Garnier et al for predicting secondary structure are compared to the most recent ones including sequence similarity methods, neural network, pattern recognition or joint prediction methods. The best of these methods correctly predict 63-65% of the residues in the database with cross-validation for 3 conformations, helix, beta strand and coli with a standard deviation of 6-8% per protein. However, when a homologous protein is already in the database, the accuracy of prediction by the similarity peptide method of Levin and Garnier reaches about 90%. Some conclusions can be drawn on the mechanism of protein folding. As all the prediction methods only use the local sequence for prediction (+/- 8 residues maximum) one can infer that 65% of the conformation of a residue is dictated on average by the local sequence, the rest is brought by the folding. The best predicted proteins or peptide segments are those for which the folding has less effect on the conformation. Presently, prediction of tertiary structure is only of practical use when the structure of a homologous protein is already known. Amino acid alignment to define residues of equivalent spatial position is critical for modelling of the protein. We showed for serine proteases that secondary structure prediction can help to define a better alignment. Non-homologous segments of the polypeptide chain, such as loops, libraries of known loops and/or energy minimization with various force fields, are used without yet giving satisfactory solutions. An example of modelling by homology, aided by secondary structure prediction on 2 regulatory proteins, Fnr and FixK is presented.     

Prediction of the secondary and tertiary structure of plastocyanin.

The amino acid sequences of nine plastocyanins were examined using four published methods for the prediction of secondary structure in proteins. The results of the four methods were combined in such a way as to maximize agreement, and the position of alpha helices, beta sheets, and beta turns in plastocyanin was predicted. From this result and other information, such as the position of conserved residues and the requirements for coordination of copper, a preliminary model for the mainchain folding of the molecule was presented.