Statistical analyses of the vibrational circular dichroism of selected proteins and relationship to secondary structures

The vibrational circular dichroism (VCD) spectra of 20 proteins dissolved in D2O are presented in the amide I' region. These data are decomposed into a linear combination of orthogonal subspectra generated by the principal component method of factor analysis, and the results for 13 of them are compared to their secondary structures as determined from X-ray crystallography. Factor analysis of the VCD yields six statistically significant subspectra that can be used to reproduce the spectra. Their coefficients can then be used to characterize a given protein. Comparison of cluster analyses of these VCD coefficients and of the secondary structure fractional coefficients from X-ray crystallography showed that proteins clustered in the VCD analysis were also clustered in the X-ray analysis. The relative fractions of alpha-helix and beta-sheet in the protein dominate the clustering in both data sets. Qualitative characterization of the secondary structure of a given protein is obtained from its clustering on the basis of spectral characteristics. A strong linear correlation was found between the coefficient of the second subspectrum and the alpha-helical fraction for the proteins studied. The second coefficient also correlated to the beta-sheet fraction, and the first coefficient weakly correlated to the fraction for "other". Subsequent multiple-parameter regression analyses of the VCD factor analysis coefficients, constrained to include only significant dependencies, yielded reliable determination of the alpha-helix fraction and somewhat less confident determination of beta-sheet, bend, and "other" components. Predictive capability for proteins not in the regression was good. Varimax rotation of the coefficients transformed the subspectra and gave simple correlations to secondary structure components but had less reliability and more restrictions than the multiple regression on the original coefficients. The partial least-squares analysis method was also used to predict fractional secondary structures for the training set proteins but resulted in somewhat higher average error, particularly for beta-sheet, than the multiple regression. The turn fraction was effectively undetermined in both the regression and partial least-squares analyses. These statistical analyses represent the first determination of a quantitative relationship between VCD spectra and secondary structure in proteins.     

Estimation of the number of alpha-helical and beta-strand segments in proteins using circular dichroism spectroscopy

A simple approach to estimate the number of alpha-helical and beta-strand segments from protein circular dichroism spectra is described. The alpha-helix and beta-sheet conformations in globular protein structures, assigned by DSSP and STRIDE algorithms, were divided into regular and distorted fractions by considering a certain number of terminal residues in a given alpha-helix or beta-strand segment to be distorted. The resulting secondary structure fractions for 29 reference proteins were used in the analyses of circular dichroism spectra by the SELCON method. From the performance indices of the analyses, we determined that, on an average, four residues per alpha-helix and two residues per beta-strand may be considered distorted in proteins. The number of alpha-helical and beta-strand segments and their average length in a given protein were estimated from the fraction of distorted alpha-helix and beta-strand conformations determined from the analysis of circular dichroism spectra. The statistical test for the reference protein set shows the high reliability of such a classification of protein secondary structure. The method was used to analyze the circular dichroism spectra of four additional proteins and the predicted structural characteristics agree with the crystal structure data.     

Vibrational circular dichroism studies of epidermal growth factor and basic fibroblast growth factor.

Vibrational circular dichroism (VCD) studies are reported for two unrelated recombinant growth factor proteins: epidermal growth factor and basic fibroblast growth factor (bFGF). NMR, electronic CD, and bFGF X-ray studies indicate that these two proteins are primarily composed of beta-sheet and loop secondary structure elements with no detectable alpha-helices. Two reports on solution conformation of these proteins using FTIR absorption spectroscopy with subsequent resolution enhancement confirmed the presence of a large fraction of a beta-sheet conformation but in addition indicated the presence of large absorption bands in the 1650-1656 cm-1 region, which are typically assigned to alpha-helices. The VCD spectra of both proteins have band shapes that strongly resemble those of other high beta-sheet fraction proteins, such as the trypsin family of proteins. Quantitative analysis of the VCD spectra also indicates that these proteins are predominantly in beta-sheet and extended ("other") conformations with very little alpha-helix fraction. These results agree with the CD interpretation and affirm that the FTIR peaks in the region 1650-1656 cm-1 can be assigned to loops. This study provides an example of the limitations of using FTIR frequencies alone for examination of protein secondary structure.     

Quantitative analysis of protein far UV circular dichroism spectra by neural networks.

A new method based on neural network theory is presented to analyze and quantify the information content of far UV circular dichroism spectra. Using a backpropagation network model with a single hidden layer between input and output, it was possible to deduce five different secondary structure fractions (helix, parallel and antiparallel beta-sheet, beta-turn and random coil) with satisfactory correlations between calculated and measured secondary structure data. We demonstrate that for each wavelength interval a specific network is suitable. The remaining discrepancy between the secondary structure data from neural network prediction and crystallography may be attributed to errors in the determination of protein concentration and random noise in the CD signal, as indicated by simulations.     

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.     

Biochemical analysis of the modular enzyme inosine 5'-monophosphate dehydrogenase.

Two prominent domains have been identified in the X-ray crystal structure of inosine-5'-monophosphate dehydrogenase (IMPDH), a core domain consisting of an alpha/beta barrel which contains the active site and an inserted subdomain whose structure is less well defined. The core domain encompassing amino acids 1-108 and 244-514 of wild-type human IMPDH (II) connected by the tetrapeptide linker Ile-Arg-Thr-Gly was expressed. The subdomain including amino acids 99-244 of human wild-type IMPDH (II) was expressed as a His-tagged fusion protein, where the His-tag was removable by enterokinase cleavage. These two proteins as well as wild-type human IMPDH (II), all proteins expressed in Escherichia coli, have been purified to apparent homogeneity. Both the wild-type and core domain proteins are tetrameric and have very similar enzymatic activities. In contrast, the subdomain migrates as a monomer or dimer on a gel filtration column and lacks enzymatic activity. Circular dichroism spectropolarimetry indicates that the core domain retains secondary structure very similar to full-length IMPDH, with 30% alpha-helix and 30% beta-sheet vs 33% alpha-helix and 29% beta-sheet for wild-type protein. Again, the subdomain protein is distinguished from both wild-type and core domain proteins by its content of secondary structure, with only 15% each of alpha-helix and beta-sheet. These studies demonstrate that the core domain of IMPDH expressed separately is both structurally intact and enzymatically active. The availability of the modules of IMPDH will aid in dissecting the architecture of this enzyme of the de novo purine nucleotide biosynthetic pathway, which is an important target for immunosuppressive and antiviral drugs.     

Prediction and Fourier transform infrared spectroscopy estimation of the secondary structure of a Bacillus licheniformis endo-beta-1,3-1,4-D-glucanase.

The secondary structure of a recombinant Bacillus licheniformis endo-beta-1,3-1,4-D-glucanase (EC.3.2.1.73) has been estimated by Fourier Transform Infrared Spectroscopy and also predicted by the algorithm of Chou and Fasman. From the curve fitting of the deconvolved IR spectrum, the most probable distribution of the secondary structural classes appears to be about 40% beta-sheet, 25% reverse turn, 24% non-ordered and 11% alpha-helix. From theoretical prediction of secondary structure the protein would present 37% beta-sheet, 31% reverse turn, 22% non-ordered and 10% alpha-helix.     

The conformation of T4 bacteriophage dihydrofolate reductase from circular dichroism

The secondary and tertiary structure of T4 bacteriophage dihydrofolate reductase is investigated by vacuum ultraviolet circular dichroism (CD) spectroscopy and probability analysis of the primary amino acid sequence. The far ultraviolet CD spectrum of the enzyme in the range of 260-178 nm is analyzed by the generalized inverse and variable selection methods developed by our laboratory. Variable selection yields an average content of 26% alpha-helix, 21% antiparallel beta-sheet, 10% parallel beta-sheet, 20% beta-turns, and 32% "other" structures within the T4 protein. The characteristic peaks of the CD spectrum indicate that the enzyme has a lot of antiparallel beta-sheet, which is typical of the alpha + beta tertiary class of globular proteins. The secondary structure of the protein is also analyzed by using four statistical methods on the amino acid sequence. Although the secondary structures predicted by each individual statistical method vary to a considerable extent, the fractions of each structure jointly predicted by a majority of the methods are in excellent agreement with our CD analysis. The alternating arrangement for some segments of alpha-helix and beta-sheet predicted from primary structure to be within the enzyme is characteristic of proteins containing parallel beta-sheet. This supports our conclusion that the protein contains both parallel and antiparallel beta-sheet structures, but finding both types of beta-sheet also means that the protein may have the variation on alpha/beta tertiary structure recently found in EcoRI endonuclease and thymidylate synthase. These observations, in conjunction with other physical properties of the T4 reductase, suggest that the enzyme perhaps shares an evolution in common with the dihydrofolate reductases derived from type I R-plasmids rather than with the host-cell protein.     

Determination of protein secondary structure using factor analysis of infrared spectra

A method is presented for determining the secondary structural composition of a protein in aqueous solution from its infrared spectrum. A factor analysis approach is used to analyze the infrared spectra of 18 proteins whose crystal structures are known from X-ray studies. Factor analysis followed by multiple linear regression identifies those eigenspectra that correlate with the variation in properties described by the calibration set. The properties of interest in this study are % alpha-helix, % beta-sheet, and % turns. In the analysis of an unknown, the factor loadings required to reproduce its spectrum are substituted in the regression equation for each property to predict its secondary structural composition. The accuracy of the method was determined by removing each standard, in turn, from the calibration set and using a calibration set generated from the remainder to predict its composition. By this method we obtain standard errors of prediction of 3.9% for alpha-helix, 8.3% for beta-sheet, and 6.6% for turns. The method may also be applied to the spectra of proteins in 2H2O. The method has important advantages over those currently in use for the quantitative analysis of the infrared spectra of proteins. Manipulation of the spectrum is kept to a minimum, no curve-fitting is necessary, and the several amide I band components need not be assigned.     

Estimation of protein secondary structure and error analysis from circular dichroism spectra

The estimation of protein secondary structure from circular dichroism spectra is described by a multivariate linear model with noise (Gauss-Markoff model). With this formalism the adequacy of the linear model is investigated, paying special attention to the estimation of the error in the secondary structure estimates. It is shown that the linear model is only adequate for the alpha-helix class. Since the failure of the linear model is most likely due to nonlinear effects, a locally linearized model is introduced. This model is combined with the selection of the estimate whose fractions of secondary structure summate to approximately one. Comparing the estimation from the CD spectra with the X-ray data (by using the data set of W.C. Johnson Jr., 1988, Annu. Rev. Biophys. Chem. 17, 145-166) the root mean square residuals are 0.09 (alpha-helix), 0.12 (anti-parallel beta-sheet), 0.08 (parallel beta-sheet), 0.07 (beta-turn), and 0.09 (other). These residuals are somewhat larger than the errors estimated from the locally linearized model. In addition to alpha-helix, in this model the beta-turn and "other" class are estimated adequately. But the estimation of the antiparallel and parallel beta-sheet class remains unsatisfactory. We compared the linear model and the locally linearized model with two other methods (S. W. Provencher and J. Gl?ckner, 1981, Biochemistry 20, 1085-1094; P. Manavalan and W. C. Johnson Jr., 1988, Anal. Biochem. 167, 76-85). The locally linearized model and the Provencher and Gl?ckner method provided the smallest residuals. However, an advantage of the locally linearized model is the estimation of the error in the secondary structure estimates.     

Minimal model for studying prion-like folding pathways

The Monte Carlo technique is used to simulate the energy landscape and the folding kinetics of a minimal prion-like protein model. We show that the competition between hydrogen-bonding and hydrophobic interactions yields two energetically favored secondary structures, an alpha-helix and a beta-hairpin. Folding simulations indicate that the probability of reaching the alpha-helix form from a denatured random conformation is much higher than the probability of reaching the beta-sheet form, even though the beta-sheet has a lower energy. The existence of a lower energy beta-sheet state gives the possibility for the normal alpha-helix structure to take a structural transformation into the beta-sheet structure under external influences.     

Genetic code redundancy and the evolutionary stability of protein secondary structure

The genetic code has an inherent bias towards some amino acids because of the variable number of synonymous codons per amino acid. The extent to which these biases are expressed in protein secondary structure is described through the analysis of the overall amino acid compositions of the alpha-helix, beta-sheet, beta-turn and random coil segments elucidated by X-ray crystallography. Given the concept of neutral mutation in proteins, the allocation of synonyms in the genetic code appears to protect secondary structures from amino acid changes and discourages the appearance of chemically complex residues. The level of protection is similar for each structural form, despite their clear preferences for certain amino acids. The organization of the code is therefore relevant to the preservation of conformation seen in the evolution of many protein families.     

Determination of alpha-helix and beta-sheet stability in the solid state: a solid-state NMR investigation of poly(L-alanine)

The relative stability of alpha-helix and beta-sheet secondary structure in the solid state was investigated using poly(L-alanine) (PLA) as a model system. Protein folding and stability has been well studied in solution, but little is known about solid-state environments, such as the core of a folded protein, where peptide packing interactions are the dominant factor in determining structural stability. (13)C cross-polarization with magic angle spinning (CPMAS) NMR spectroscopy was used to determine the backbone conformation of solid powder samples of 15-kDa and 21.4-kDa PLA before and after various sample treatments. Reprecipitation from helix-inducing solvents traps the alpha-helical conformation of PLA, although the method of reprecipitation also affects the conformational distribution. Grinding converts the secondary structure of PLA to a final steady-state mixture of 55% beta-sheet and 45% alpha-helix at room temperature regardless of the initial secondary structure. Grinding PLA at liquid nitrogen temperatures leads to a similar steady-state mixture with 60% beta-sheet and 40% alpha-helix, indicating that mechanical shear force is sufficient to induce secondary structure interconversion. Cooling the sample in liquid nitrogen or subjecting it to high pressure has no effect on secondary structure. Heating the sample without grinding results in equilibration of secondary structure to 50% alpha-helix/50% beta-sheet at 100 degrees C when starting from a mostly alpha-helical state. No change was observed upon heating a beta-sheet sample, perhaps due to kinetic effects and the different heating rate used in the experiments. These results are consistent with beta-sheet approximately 260 J/mol more stable than alpha-helix in solid-state PLA.     

Determination of the secondary structure content of proteins in aqueous solutions from their amide I and amide II infrared bands. Comparison between classical and partial least-squares methods

A method for estimating protein secondary structure from infrared spectra has been developed. The infrared spectra of H2O solutions of 13 proteins of known crystal structure have been recorded and corrected for the spectral contribution of water in the amide I and II region by using the algorithm of Dousseau et al. [Dousseau, F., Therrien, M., & Pézolet, M. (1989) Appl. Spectrosc. 43, 538-542]. This calibration set of proteins has been analyzed by using either a classical least-squares (CLS) method or the partial least-squares (PLS) method. The pure-structure spectra calculated by the classical least-squares method are in good agreement with spectra of poly(L-lysine) in the alpha-helix, beta-sheet, and undefined conformations. The results show that the best agreement between the secondary structure determined by X-ray crystallography and that predicted by infrared spectroscopy is obtained when both the amide I and II bands are used to generate the calibration set, when the PLS method is used, and when it is assumed that the secondary structure of proteins is composed of only four types of structure: ordered and disordered alpha-helices, beta-sheet, and undefined conformation. Attempts to include turns in the secondary structure estimation have led to a loss of accuracy. The standard deviation of the difference between X-ray and infrared secondary structure estimates with this method is 4.8% for the alpha-helix, 3.7% for the beta-sheet, and 5.1% for the undefined structure, whereas the regression coefficients are 0.95, 0.96, and 0.56, respectively. The spectra of the calibration proteins were also recorded in 2H2O solution.(ABSTRACT TRUNCATED AT 250 WORDS)     

Structural changes and aggregation process of Cu/containing amine oxidase in the presence of 2,2,2'-trifluoroethanol

Conformational and structural changes of lentil seedlings amine oxidase (LSAO) were studied in the presence of trifluoroethanol (TFE) by spectroscopic and analytical techniques. At TFE concentrations up to 5%, the induction of a structural transition from beta-sheet to alpha-helix and up to 10% TFE a structural transition from alpha-helix to beta-sheet as well as inactivation of the enzyme are observed. At TFE concentrations between 10-35%, LSAO proves to be prone to aggregation and beyond 35% TFE leads to a non-native protein structure with a high alpha-helix content. The obtained results revealed that the aggregation of LSAO is strongly linked to the nature of secondary structures.     

Effect of hydrophobic surfactant proteins SP-B and SP-C on phospholipid monolayers. Protein structure studied using 2D IR and beta correlation analysis

We have applied two-dimensional infrared (2D IR) and betanu correlation spectroscopy to in-situ IR spectroscopy of pulmonary surfactant proteins SP-B and SP-C in lipid-protein monolayers at the air-water interface. For both SP-B and SP-C, a statistical windowed autocorrelation method identified two separate surface pressure regions that contained maximum amide I intensity changes: 4-25 mN/m and 25-40 mN/m. For SP-C, 2D IR and betanu correlation analyses of these regions indicated that SP-C adopts a variety of secondary structure conformations, including alpha-helix, beta-sheet, and an intermolecular aggregation of extended beta-sheet structure. The main alpha-helix band split into two peaks at high surface pressures, indicative of two different helix conformations. At low surface pressures, all conformations of the SP-C molecule reacted identically to increasing surface pressure and reoriented in phase with each other. Above 25 mN/m, however, the increasing surface pressure selectively affected the coexisting protein conformations, leading to an independent reorientation of the protein conformations. The asynchronous 2D IR spectrum of SP-B showed the presence of two alpha-helix components, consistent with two separate populations of alpha-helix in SP-B-a hydrophobic fraction associated with the lipid chains and a hydrophilic fraction parallel to the membrane surface. The distribution of correlation intensity between the two alpha-helix cross peaks indicated that the more hydrophobic helix fraction predominates at low surface pressures whereas the more hydrophilic helix fraction predominates at high surface pressures. The different SP-B secondary structures reacted identically to increasing surface pressure, leading to a reorientation of all SP-B subunits in phase with one another.     

A comparative CD study of carbonic anhydrase isoenzymes with different number of tryptophans: impact on calculation of secondary structure content.

The CD spectra of human carbonic anhydrase I and II and bovine carbonic anhydrase III were recorded and analyzed. The 3D structures of these isoenzymes are known, showing very similar secondary structure and polypeptide-chain fold. The tryptophan content, however, differs between the isoenzymes, i.e., isoenzymes I, II, and III possess 6, 7, and 8 tryptophans, respectively. All of the tryptophans except the additional tryptophans in isoenzymes II and III, i.e., W245 and W47, are conserved. Despite the fact that X-ray structure determinations showed that the isoenzymes had highly similar secondary structure, the contents of alpha-helix and beta-sheet structure differed considerably when using different CD algorithms for estimation of the fractions of various secondary structural elements. This shows that aromatic amino acids also interfere in the wavelength region (far-UV) used to calculate the amount of secondary structure. Such interference is especially problematic when analyzing proteins like carbonic anhydrase, which consist mainly of beta-structure that gives rise to weak ellipticity bands, compared to the bands arising from alpha-helical structure.     

Deconvolution of the circular dichroism spectra of proteins: the circular dichroism spectra of the antiparallel beta-sheet in proteins.

A recently developed algorithm, called Convex Constraint Analysis (CCA), was successfully applied to determine the circular dichroism (CD) spectra of the pure beta-pleated sheet in globular proteins. On the basis of X-ray diffraction determined secondary structures, the original data set used (Perczel, A., Hollosi, M., Tusnady, G. Fasman, G.D. Convex constraint analysis: A natural deconvolution of circular dichroism curves of proteins, Prot. Eng., 4:669-679, 1991), was improved by the addition of proteins with high beta-pleated sheet content. The analysis yielded CD curves of the pure components of the main secondary structural elements (alpha-helix, antiparallel beta-pleated sheet, beta-turns, and unordered conformation), as well as a curve attributed to the "aromatic contribution" in the wavelength range of 195-240 nm. Upon deconvolution the curves obtained were assigned to various secondary structures. The calculated weights (percentages determining the contributions of each pure component curve in the measured CD spectra of a given protein) were correlated with the X-ray diffraction determined percentages in an assignment procedure and were evaluated. The Pearson product correlation coefficients (R) are significant for all five components. The new pure component curves, which were obtained through deconvolution of the protein CD spectra alone, are promising candidates for determining the percentages of the secondary structural components in globular proteins without the necessity of adopting an X-ray database. The CD spectrum of the CheY protein was interesting because it has the characteristic shape associated with the alpha-helical structure, but upon analysis yielded a considerable amount of beta-sheet in agreement with the X-ray structure.     

Secondary structure of the mammalian 70-kilodalton heat shock cognate protein analyzed by circular dichroism spectroscopy and secondary structure prediction

Heat shock proteins are rapidly synthesized when cells are exposed to stressful agents that cause protein damage. The 70-kDa heat shock induced proteins and their closely related constitutively expressed cognate proteins bind to unfolded and aberrant polypeptides and to hydrophilic peptides. The structural features of the 70-kDa heat shock proteins that confer the ability to associate with diverse polypeptides are unknown. In this study, we have used circular dichroism (CD) spectroscopy and secondary structure prediction to analyze the secondary structure of the mammalian 70-kDa heat shock cognate protein (hsc 70). The far-ultraviolet CD spectrum of hsc 70 indicates a large fraction of alpha-helix in the protein and resembles the spectra one obtains from proteins of the alpha/beta structural class. Analysis of the CD spectra with deconvolution methods yielded estimates of secondary structure content. The results indicate about 40% alpha-helix and 20% aperiodic structure within hsc 70 and between 16-41% beta-sheet and 21-0% beta-turn. The Garnier-Osguthorpe-Robson method of secondary structure prediction was applied to the rat hsc 70 amino acid sequence. The predicted estimates of alpha-helix and aperiodic structure closely matched the values derived from the CD analysis, whereas the predicted estimates of beta-sheet and beta-turn were midway between the CD-derived values. Present evidence suggests that the polypeptide ligand binding domain of the 70-kDa heat shock protein resides within the C-terminal 160 amino acids [Milarski, K. L., & Morimoto, R. I. (1989) J. Cell Biol. 109, 1947-1962].(ABSTRACT TRUNCATED AT 250 WORDS)