Analyzing protein circular dichroism spectra for accurate secondary structures.

We have developed an algorithm to analyze the circular dichroism of proteins for secondary structure. Its hallmark is tremendous flexibility in creating the basis set, and it also combines the ideas of many previous workers. We also present a new basis set containing the CD spectra of 22 proteins with secondary structures from high quality X-ray diffraction data. High flexibility is obtained by doing the analysis with a variable selection basis set of only eight proteins. Many variable selection basis sets fail to give a good analysis, but good analyses can be selected without any a priori knowledge by using the following criteria: (1) the sum of secondary structures should be close to 1.0, (2) no fraction of secondary structure should be less than -0.03, (3) the reconstructed CD spectrum should fit the original CD spectrum with only a small error, and (4) the fraction of alpha-helix should be similar to that obtained using all the proteins in the basis set. This algorithm gives a root mean square error for the predicted secondary structure for the proteins in the basis set of 3.3% for alpha-helix, 2.6% for 3(10)-helix, 4.2% for beta-strand, 4.2% for beta-turn, 2.7% for poly(L-proline) II type 3(1)-helix, and 5.1% for other structures when compared with the X-ray structure.     

An infrared and circular dichroism combined approach to the analysis of protein secondary structure.

Selected regions of infarred (ir) and circular dichroism (CD) spectral data from 10 proteins were combined and analyzed by a factor analysis method. The regions consisted of the area normalized amide I region from 1700 to 1600 cm-1 for the ir spectra and from 178 to 240 nm for the CD spectra. Each CD spectrum was scaled by a factor of 0.5 before appending the data to the ir spectral data. The scaling factor was deemed necessary to account for relative intensity differences between the ir and CD data and provided nearly optimum agreement between secondary structure estimated by the combined approach to secondary structure determined by X-ray crystallography. The ir/CD combined approach to estimation of helix, beta-sheet, beta-turn, and other or undefined secondary structure agreed with X-ray crystallographic determined structure better than estimation using data from either method alone. Correlation coefficients between X-ray and ir/CD combined secondary structure determinations were 0.99 for helix, 0.90 for beta-sheet, 0.70 for beta-turn, and 0.78 for other structure. The four most significant eigenvectors or basis spectra from eigenanalysis of the ir/CD data are presented as well as generalized inverse spectra for four secondary structures.     

The circular dichroism of tumor necrosis factor-α: Measurement into the vacuum UV and analysis for secondary structure

The circular dichroism (CD) spectrum of tumor necrosis factor-α has been measured into the vacuum UV to 168 nm. Analysis of the CD for secondary structure is in good agreement with X-ray diffraction results, but the analysis is somewhat unstable. Adding the CD of this protein together with its X-ray determined secondary structure to the basis set should improve subsequent analyses of CD spectra for other all-β proteins.     

Analysis of protein circular dichroism spectra for secondary structure using a simple matrix multiplication

Inverse circular dichroism (CD) spectra are presented for each of the five major secondary structures of proteins: alpha-helix, antiparallel and parallel beta-sheet, beta-turn, and other (random) structures. The fraction of the each secondary structure in a protein is predicted by forming the dot product of the corresponding inverse CD spectrum, expressed as a vector, with the CD spectrum of the protein digitized in the same way. We show how this method is based on the construction of the generalized inverse from the singular value decomposition of a set of CD spectra corresponding to proteins whose secondary structures are known from X-ray crystallography. These inverse spectra compute secondary structure directly from protein CD spectra without resorting to least-squares fitting and standard matrix inversion techniques. In addition, spectra corresponding to the individual secondary structures, analogous to the CD spectra of synthetic polypeptides, are generated from the five most significant CD eigenvectors.     

Determination of secondary structure of globular proteins using circular dichroism spectra

A ridge regression method is presented for prediction of the secondary structure of proteins by the circular dichroism spectra (CD) from 190 to 236 nm. Eight types of the secondary structure were calculated on a microcalculator. The method is based on the X-ray data of Kabsh and Sander. The teaching rule is constructed on CD spectra of 30 proteins of all structural classes of the globular proteins (alpha, alpha/beta, alpha + beta and beta-proteins). The errors of the methods are analysed by removing each protein from the reference set and analyzing its structure in terms of the remaining proteins. Correlation coefficients and root-mean square deviations between CD and X-ray data were: 0.99 and 0.03 for alpha-helix, 0.86 and 0.02 for 3(10)-helix, 0.92 and 0.06 for antiparallel beta-sheet, 0.86 and 0.03 for parallel beta-sheet, 0.94 and 0.01 for T3 beta-turn, 0.85 and 0.02 for other beta-turn, 0.84 and 0.03 for S-bends, 0.83 and 0.04 for "random" structure.     

Estimation of protein secondary structure from circular dichroism spectra: inclusion of denatured proteins with native proteins in the analysis

We have expanded our reference set of proteins used in the estimation of protein secondary structure by CD spectroscopy from 29 to 37 proteins by including 3 additional globular proteins with known X-ray structure and 5 denatured proteins. We have also modified the self-consistent method for analyzing protein CD spectra, SELCON3, by including a new selection criterion developed by W. C. Johnson, Jr. (Proteins Struct. Funct. Genet. 35, 307-312, 1999). The secondary structure corresponding to the denatured proteins was approximated to be 90% unordered, owing to the spectral similarity of the denatured proteins and unordered structures. We examined the thermal denaturation of ribonuclease T1 by CD using both the original and expanded sets of reference proteins and obtained more consistent results with the expanded set. The expanded set of reference proteins will be helpful for the determination of protein secondary structure from protein CD spectra with higher reliability, especially of proteins with significant unordered structure content and/or in the course of denaturation.     

Molecular Dynamics Ensemble Refinement of Intrinsically Disordered Peptides According to Deconvoluted Spectra from Circular Dichroism

We have developed a computational method of atomistically refining the structural ensemble of intrinsically disordered peptides (IDPs) facilitated by experimental measurements using circular dichroism spectroscopy (CD). A major challenge surrounding this approach stems from the deconvolution of experimental CD spectra into secondary structure features of the IDP ensemble. Currently available algorithms for CD deconvolution were designed to analyze the spectra of proteins with stable secondary structures. Herein, our work aims to minimize any bias from the peptide deconvolution analysis by implementing a non-negative linear least-squares fitting algorithm in conjunction with a CD reference data set that contains soluble and denatured proteins (SDP48). The non-negative linear least-squares method yields the best results for deconvolution of proteins with higher disordered content than currently available methods, according to a validation analysis of a set of protein spectra with Protein Data Bank entries. We subsequently used this analysis to deconvolute our experimental CD data to refine our computational model of the peptide secondary structure ensemble produced by all-atom molecular dynamics simulations with implicit solvent. We applied this approach to determine the ensemble structures of a set of short IDPs, that mimic the calmodulin binding domain of calcium/calmodulin-dependent protein kinase II and its 1-amino-acid and 3-amino-acid mutants. Our study offers a, to our knowledge, novel way to solve the ensemble secondary structures of IDPs in solution, which is important to advance the understanding of their roles in regulating signaling pathways through the formation of complexes with multiple partners.     

Prediction of protein secondary structure from circular dichroism using theoretically derived spectra

Circular dichroism (CD) is a spectroscopic technique commonly used to investigate the structure of proteins. Major secondary structure types, alpha‐helices and beta‐strands, produce distinctive CD spectra. Thus, by comparing the CD spectrum of a protein of interest to a reference set consisting of CD spectra of proteins of known structure, predictive methods can estimate the secondary structure of the protein. Currently available methods, including K2D2, use such experimental CD reference sets, which are very small in size when compared to the number of tertiary structures available in the Protein Data Bank (PDB). Conversely, given a PDB structure, it is possible to predict a theoretical CD spectrum from it. The methodological framework for this calculation was established long ago but only recently a convenient implementation called DichroCalc has been developed. In this study, we set to determine whether theoretically derived spectra could be used as reference set for accurate CD based predictions of secondary structure. We used DichroCalc to calculate the theoretical CD spectra of a nonredundant set of structures representing most proteins in the PDB, and applied a straightforward approach for predicting protein secondary structure content using these theoretical CD spectra as reference set. We show that this method improves the predictions, particularly for the wavelength interval between 200 and 240 nm and for beta‐strand content. We have implemented this method, called K2D3, in a publicly accessible web server at http://www. ogic.ca/projects/k2d3 . Proteins 2012. © 2011 Wiley Periodicals, Inc.     

Automatic amide I frequency selection for rapid quantification of protein secondary structure from Fourier transform infrared spectra of proteins

Here we report the development of a new neural network based approach for rapid quantification of protein secondary structure from Fourier transform infrared (FTIR) spectra of proteins. A technique for efficiently reducing the amount of spectral data by almost 90% is suggested to facilitate faster neural network analysis. Additionally, an automatic procedure is introduced for selecting only those regions within the amide I band of protein FTIR spectra, which can be best related to secondary structure contents by subsequent neural network analysis. Based on a given reference set of FTIR spectra from proteins with known secondary structure, a subset of merely 29 out of 101 amide I absorbance values could be identified, which lead to an improved prediction accuracy. The average prediction accuracy achieved for helix, sheet, turn, bend, and other is 4.96% which is better than that achieved by alternative methods that have been previously reported indicating the significant potential of this approach. Our suggested automatic amide I frequency selection procedure may be easily extended to identify promising regions from spectral data recorded by other spectroscopic techniques, like for example circular dichroism spectroscopy.     

Accuracy of protein secondary structure determination from circular dichroism spectra based on immunoglobulin examples

Strong contribution of the aromatic amino acid side chain chromophores to the far-UV circular dichroism (CD) spectra substantially distorts a relatively weak CD signal originating from beta sheet, the main type of immunoglobulin secondary structure. In this study we compared the secondary structure calculated from the far-UV CD spectra with the X-ray data for three antibody Fab fragments. Calculations were performed with three different algorithms, using two sets of reference proteins. Low standard deviations between all six estimates indicate stable mathematical solutions. Despite pronounced differences in the shape and amplitude of the CD spectra, we found a strong correlation between CD and X-ray data in the secondary structure for every protein studied. The number and average length of the secondary structure elements estimated from the CD spectra closely resemble those of the X-ray data. Agreement between spectroscopic and crystallographic results demonstrates that modern methods of secondary structure calculation are resilient to distortions of the far-UV CD spectra of immunoglobulins caused by aromatic side chain chromophores.     

Protein secondary structure from Fourier transform infrared spectroscopy: a data base analysis

An infrared (ir) method to determine the secondary structure of proteins in solution using the amide I region of the spectrum has been devised. The method is based on the circular dichroism (CD) matrix method for secondary structure analysis given by Compton and Johnson (L. A. Compton and W. C. Johnson, 1986, Anal. Biochem. 155, 155-167). The infrared data matrix was constructed from the normalized Fourier transform infrared spectra from 1700 to 1600 cm-1 of 17 commercially available proteins. The secondary structure matrix was constructed from the X-ray data of the seventeen proteins with secondary structure elements of helix, beta-sheet, beta-turn, and other (random). The CD and ir methods were compared by analyzing the proteins of the CD and ir databases as unknowns. Both methods produce similar results compared to structures obtained by X-ray crystallographic means with the CD slightly better for helix conformation, and the ir slightly better for beta-sheet. The relatively good ir analysis for concanavalin A and alpha-chymotrypsin indicate that the ir method is less affected by the presence of aromatic groups. The concentration of the protein and the cell path length need not be known for the ir analysis since the spectra can be normalized to the total ir intensity in the amide I region. The ir spectra for helix, beta-sheet, beta-turn, and other, as extracted from the data-base, agree with the literature band assignments. The ir data matrix and the inverse matrix necessary to analyze unknown proteins are presented.     

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.     

Comparison of the conformation and orientation of alamethicin and melittin in lipid membranes

The secondary structure of alamethicin in lipid membranes below and above the lipid phase transition temperature Tt is determined by Raman spectroscopy and circular dichroism (CD) measurements. In both cases structural data are obtained by fitting the experimental spectra by a superposition of the spectra of 15 reference proteins of known three-dimensional structure. According to the Raman experiments, in a lipid bilayer above Tt alamethicin is helical from residue 1 to 12, whereas below Tt the helix extends from residue 1 to 16. The remaining C-terminal part is nonhelical up to the end residue 20 both above and below Tt. A considerable lower helix content is derived from CD, namely, 38% and 46% above and below Tt, respectively, in agreement with several reported values for CD in the literature. It is shown that the commonly used set of CD spectra of water-soluble reference proteins is unsuitable to describe the CD spectra of alamethicin correctly. Therefore the secondary structure of alamethicin as derived from CD measurements is at the present state of analysis unreliable. In contrast to the case of alamethicin, the CD spectra of melittin in lipid membranes are correctly described by the reference protein spectra. The helix content of melittin is determined thereby to be 72% in lipid membranes above Tt and 75% below Tt. The data are in accord with a structure where the hydrophobic part of melittin adopts a bent helix as determined recently by Raman spectroscopy [Vogel, H., & J?hnig, F. (1986) Biophys. J. 50, 573]. The orientational order parameters of the helical parts of alamethicin and of melittin in a lipid membrane are deduced from the difference between a corresponding CD spectrum of a polypeptide in planar multibilayers and that in lipid vesicles. The presented method for determining helix order parameters is new and may be generally applicable to other membrane proteins. The orientation of the helical part of both polypeptides depends on the physical state of the lipid bilayer at maximal membrane hydration and in the ordered lipid state furthermore on the degree of membrane hydration. Under conditions where alamethicin and melittin are incorporated in an aggregated form in a fluid lipid membrane at maximal water content the helical segments are oriented preferentially parallel to the membrane normal. Cooling such lipid membranes to a temperature below Tt changes the orientation of the helical part of alamethicin as well as melittin toward the membrane plane.(ABSTRACT TRUNCATED AT 400 WORDS)     

Analysis of the circular dichroism spectrum of proteins using the convex constraint algorithm: a practical guide.

Due to the time scale of circular dichroism (CD) measurements, it is theoretically possible to deconvolute such a spectrum if the pure CD spectra differ significantly from one another. In the last decade several methods have been published aiming at obtaining the conformational weights, or percentages (which are the coefficients for a linear combination) of the so-called typical secondary structural elements making up the three-dimensional structure of proteins. Two methods that can be used to determine the secondary structures of proteins are described here. The first method, called LINCOMB, is a simple algorithm based on a least-squares fit with a set of reference spectra representing the known secondary structures and yielding an estimation of weights attributed to alpha-helix, beta-pleated sheet (mainly antiparallel), beta-turns, unordered form, and aromatic/disulfide (or nonpeptide) contributions of the protein being analyzed. This method requires a "template" or reference curve set, which was obtained from the second method. The second method, "convex constraint analysis," is a general deconvolution method for a CD spectra set of any variety of conformational type. The algorithm, based on a set of three constraints, is able to deconvolute a set of CD curves to its common "pure"-component curves and conformational weights. To analyze a single CD spectrum with this method, the spectrum is appended to the data set used as a reference data set. As a way to determine the reliability of the algorithm and provide a guideline to its usage, some applications are presented.     

Native-like secondary structure of molten globules

The most common evidence for the existence of secondary structure in a globular protein is the presence of a strongly pronounced far-UV circular dichroism (CD) spectrum. Although CD spectra of native proteins are well described and their quantitative analysis is widely used, similar studies for denatured proteins have still to be done. Far-UV CD spectra of nine proteins in the native and the pH-induced molten globule states were acquired and analyzed. Singular value decomposition showed that the spectra of molten globules could be described as a superposition of at least three independent components (most likely alpha-, beta- and irregular structure). A self-consistent procedure of CD spectra analysis revealed the existence of a clear correlation between the shape of the molten globule spectra and the content of secondary structure elements in the corresponding native proteins, as determined from X-ray data. A mathematical expression of this correlation in terms of the Pierson coefficient amounts to the value of 0.9 for both the alpha-helix and the beta-structure. Thus, the secondary structure of proteins in the molten globule state is close to that in the native state.     

Secondary structural analysis of retrovirus integrase: characterization by circular dichroism and empirical prediction methods

The retrovirus integrase (IN) protein is essential for integration of viral DNA into host DNA. The secondary structure of the purified IN protein from avian myeloblastosis virus was investigated by both circular dichroism (CD) spectroscopy and five empirical prediction methods. The secondary structures determined from the resolving of CD spectra through a least-squares curve fitting procedure were compared with those predicted from four statistical methods, e.g., the Chou-Fasman, Garnier-Osguthorpe-Robson, Nishikawa-Ooi, and a JOINT scheme which combined all three of these methods, plus a pure a priori one, the Ptitsyn-Finkelstein method. Among all of the methods used, the Nishikawa-Ooi prediction gave the closest match in the composition of secondary structure to the CD result, although the other methods each correctly predicted one or more secondary structural group. Most of the alpha-helix and beta-sheet states predicted by the Ptitsyn-Finkelstein method were in accord with the Nishikawa-Ooi method. Secondary structural predictions by the Nishikawa-Ooi method were extended further to include IN proteins from four phylogenetic distinct retroviruses. The structural relationships between the four most conserved amino acid blocks of these IN proteins were compared using sequence homology and secondary structure predictions.     

On the analysis of membrane protein circular dichroism spectra

Analysis of circular dichroism spectra of proteins provides information about protein secondary structure. Analytical methods developed for such an analysis use structures and spectra of a set of reference proteins. The reference protein sets currently in use include soluble proteins with a wide range of secondary structures, and perform quite well in analyzing CD spectra of soluble proteins. The utility of soluble protein reference sets in analyzing membrane protein CD spectra, however, has been questioned in a recent study that found current reference protein sets to be inadequate for analyzing membrane proteins. We have examined the performance of reference protein sets available in the CDPro software package for analyzing CD spectra of 13 membrane proteins with available crystal structures. Our results indicate that the reference protein sets currently available for CD analysis perform reasonably well in analyzing membrane protein CD spectra, with performance indices comparable to those for soluble proteins. Soluble + membrane protein reference sets, which were constructed by combining membrane proteins with soluble protein reference sets, gave improved performance in both soluble and membrane protein CD analysis.     

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.     

Determination of the secondary structure of proteins from circular dichroism spectra. IV. Contribution of aromatic amino acid residues into circular dichroism spectra of proteins in the peptide region

A new method for determination of the secondary protein structure from the CD spectra taking into account the contribution of aromatic amino acid residues is proposed. New proteins reference CD spectra for five secondary structures (alpha-helices, antiparallel and parallel beta-structures, beta-bends and irregular form) without contribution of aromatic residues are obtained. By means of this new method the secondary structure of sixteen different proteins was analysed. There is a good correlation of these results with the X-ray data.     

Distinct circular dichroism spectroscopic signatures of polyproline II and unordered secondary structures: Applications in secondary structure analyses

Circular dichroism (CD) spectroscopy is a valuable method for defining canonical secondary structure contents of proteins based on empirically‐defined spectroscopic signatures derived from proteins with known three‐dimensional structures. Many proteins identified as being “Intrinsically Disordered Proteins” have a significant amount of their structure that is neither sheet, helix, nor turn; this type of structure is often classified by CD as “other”, “random coil”, “unordered”, or “disordered”. However the “other” category can also include polyproline II (PPII)‐type structures, whose spectral properties have not been well‐distinguished from those of unordered structures. In this study, synchrotron radiation circular dichroism spectroscopy was used to investigate the spectral properties of collagen and polyproline, which both contain PPII‐type structures. Their native spectra were compared as representatives of PPII structures. In addition, their spectra before and after treatment with various conditions to produce unfolded or denatured structures were also compared, with the aim of defining the differences between CD spectra of PPII and disordered structures. We conclude that the spectral features of collagen are more appropriate than those of polyproline for use as the representative spectrum for PPII structures present in typical amino acid‐containing proteins, and that the single most characteristic spectroscopic feature distinguishing a PPII structure from a disordered structure is the presence of a positive peak around 220nm in the former but not in the latter. These spectra are now available for inclusion in new reference data sets used for CD analyses of the secondary structures of soluble proteins.