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.     

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.     

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.     

Analyses of circular dichroism spectra of membrane proteins

Circular dichroism (CD) spectroscopy is a valuable technique for the determination of protein secondary structures. Many linear and nonlinear algorithms have been developed for the empirical analysis of CD data, using reference databases derived from proteins of known structures. To date, the reference databases used by the various algorithms have all been derived from the spectra of soluble proteins. When applied to the analysis of soluble protein spectra, these methods generally produce calculated secondary structures that correspond well with crystallographic structures. In this study, however, it was shown that when applied to membrane protein spectra, the resulting calculations produce considerably poorer results. One source of this discrepancy may be the altered spectral peak positions (wavelength shifts) of membrane proteins due to the different dielectric of the membrane environment relative to that of water. These results have important consequences for studies that seek to use the existing soluble protein reference databases for the analyses of membrane proteins.     

Improved estimation of the secondary structures of proteins by vacuum-ultraviolet circular dichroism spectroscopy

The vacuum-ultraviolet circular dichroism (VUVCD) spectra of 16 globular proteins (insulin, lactate dehydrogenase, glucose isomerase, lipase, conalbumin, transferrin, catalase, subtilisin A, alpha-amylase, staphylococcal nuclease, papain, thioredoxin, carbonic anhydrase, elastase, avidin, and xylanase) were successfully measured in aqueous solutions at 25 degrees C from 260 to 160 nm under a high vacuum using a synchrotron-radiation VUVCD spectrophotometer. These proteins exhibited characteristic CD spectra below 190 nm that were related to their different secondary structures, which could not be detected with a conventional CD spectrophotometer. The component spectra of alpha-helices, beta-strands, turns, and unordered structures were obtained by deconvolution analysis of the VUVCD spectra of 31 reference proteins including the 15 proteins reported in our previous paper [Matsuo, K. et al. (2004) J. Biochem. 135, 405-411]. Prediction of the secondary-structure contents using the SELCON3 program was greatly improved, especially for alpha-helices, by extending the short-wavelength limit of CD spectra to 160 nm and by increasing the number of reference proteins. The numbers of alpha-helix and beta-strand segments, which were calculated from the distorted alpha-helix and beta-strand contents, were close to those obtained on X-ray crystallography. These results demonstrate the usefulness of synchrotron-radiation VUVCD spectroscopy for the secondary structure analysis of proteins.     

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.     

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.     

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.     

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.     

Using circular dichroism spectra to estimate protein secondary structure

Circular dichroism (CD) is an excellent tool for rapid determination of the secondary structure and folding properties of proteins that have been obtained using recombinant techniques or purified from tissues. The most widely used applications of protein CD are to determine whether an expressed, purified protein is folded, or if a mutation affects its conformation or stability. In addition, it can be used to study protein interactions. This protocol details the basic steps of obtaining and interpreting CD data, and methods for analyzing spectra to estimate the secondary structural composition of proteins. CD has the advantage that measurements may be made on multiple samples containing < or =20 microg of proteins in physiological buffers in a few hours. However, it does not give the residue-specific information that can be obtained by x-ray crystallography or NMR.     

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.     

Bioinformatics analyses of circular dichroism protein reference databases

MOTIVATION: Circular dichroism (CD) spectroscopy has become established as a key method for determining the secondary structure contents of proteins which has had a significant impact on molecular biology. Many excellent mathematical protocols have been developed for this purpose and their quality is above question. However, reference database sets of proteins, with CD spectra matched to secondary structure components derived from X-ray structures, provide the key resource for this task. These databases were created many years ago, before most CD spectrophotometers became standardized and before it was commonplace to validate X-ray structures prior to publication. The analyses presented here were undertaken to investigate the overall quality of these reference databases in light of their extensive usage in determining protein secondary structure content from CD spectra. RESULTS: The analyses show that there are a number of significant problems associated with the CD reference database sets in current use. There are disparities between CD spectra for the same protein collected by different groups. These include differences in magnitudes, peak positions or both. However, many current reference sets are now amalgamations of spectra from these groups, introducing inconsistencies that can lead to inaccuracies in the determination of secondary structure components from the CD spectra. A number of the X-ray structures used fall short on the validation criteria now employed as standard for structure determination. Many have substantial percentages of residues in the disallowed regions of the Ramachandran plot. Hence their calculated secondary structure components, used as a foundation for the reference databases, are likely to be in error. Additionally, the coverage of secondary structure space in the reference datasets is poorly correlated to the secondary structure components found in the Protein Data Bank. A conclusion is that a new reference CD database with cross-correlated, machine-independent CD spectra and validated X-ray structures that cover more secondary structure components, including diverse protein folds, is now needed. However, that reasonably accurate values for the secondary structure content of proteins can be determined from spectra is a testament to CD spectroscopy being a very powerful technique.     

Absorption flattening in the circular dichroism spectra of small membrane fragments

The inhomogeneous distribution of chromophore occurring in a particulate suspension can result in a reduction in the apparent molar ellipticity recorded in circular dichroism (CD) spectra. The possibility of such a systematic error has often been a matter of concern when CD spectra of cell membrane proteins are recorded. The recent publication of CD spectra for bacteriorhodopsin in native and sonicated membranes, in detergent-solubilized form, and reconstituted into small unilamellar vesicles [Mao, D., & Wallace, B. A. (1984) Biochemistry 23, 2667-2673] gives a unique opportunity to apply the theoretical analysis of Gordon and Holzwarth [Gordon, D. J., & Holzwarth, G. (1971) Arch. Biochem. Biophys. 142, 481-488] so as to provide a definitive answer to the question of whether absorption flattening is significant for membrane particles. We show here that the data of Mao and Wallace can be combined with the theoretical analysis of Gordon and Holzwarth to rule out significant absorption flattening effects over the range 200-240 nm for submicrometer-sized membranes. In addition, the results show that absorption flattening can be disregarded even at 190 nm for membranous material in the size range below 100 nm. The demonstration that there are no major flattening effects in the CD spectra of bacteriorhodopsin, particularly in the region of 200-240 nm, means that the experimental spectra are incompatible with the proposal that this transmembrane protein contains seven transmembrane helices.     

Convex constraint analysis: a natural deconvolution of circular dichroism curves of proteins

A new algorithm, called convex constraint analysis, has been developed to deduce the chiral contribution of the common secondary structures directly from experimental CD curves of a large number of proteins. The analysis is based on CD data reported by Yang, J.T., Wu, C.-S.C. and Martinez, H.M. [Methods Enzymol., 130, 208-269 (1986)]. Application of the decomposition algorithm for simulated protein data sets resulted in component spectra [B (lambda, i)] identical to the originals and weights [C (i, k)] with excellent Pearson correlation coefficients (R) [Chang, C.T., Wu, C.-S.C. and Yang, J.T. (1978) Anal. Biochem., 91, 12-31]. Test runs were performed on sets of simulated protein spectra created by the Monte Carlo technique using poly-L-lysine-based pure component spectra. The significant correlational coefficients (R greater than 0.9) demonstrated the high power of the algorithm. The algorithm, applied to globular protein data, independent of X-ray data, revealed that the CD spectrum of a given protein is composed of at least four independent sources of chirality. Three of the computed component curves show remarkable resemblance to the CD spectra of known protein secondary structures. This approach yields a significant improvement in secondary structural evaluations when compared with previous methods, as compared with X-ray data, and yields a realistic set of pure component spectra. The new method is a useful tool not only in analyzing CD spectra of globular proteins but also has the potential for the analysis of integral membrane proteins.     

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.     

Differentiation between transmembrane helices and peripheral helices by the deconvolution of circular dichroism spectra of membrane proteins.

The interpretation of the circular dichroism (CD) spectra of proteins to date requires additional secondary structural information of the proteins to be analyzed, such as X-ray or NMR data. Therefore, these methods are inappropriate for a CD database whose secondary structures are unknown, as in the case of the membrane proteins. The convex constraint analysis algorithm (Perczel, A., Hollósi, M., Tusnády, G., & Fasman, G. D., 1991, Protein Eng. 4, 669-679), on the other hand, operates only on a collection of spectral data to extract the common spectral components with their spectral weights. The linear combinations of these derived "pure" CD curves can reconstruct the original data set with great accuracy. For a membrane protein data set, the five-component spectra so obtained from the deconvolution consisted of two different types of alpha helices (the alpha helix in the soluble domain and the alpha T helix, for the transmembrane alpha helix), a beta-pleated sheet, a class C-like spectrum related to beta turns, and a spectrum correlated with the unordered conformation. The deconvoluted CD spectrum for the alpha T helix was characterized by a positive red-shifted band in the range 195-200 nm (+95,000 deg cm2 dmol-1), with the intensity of the negative band at 208 nm being slightly less negative than that of the 222-nm band (-50,000 and -60,000 deg cm2 dmol-1, respectively) in comparison with the regular alpha helix, with a positive band at 190 nm and two negative bands at 208 and 222 nm with magnitudes of +70,000, -30,000, and -30,000 deg cm2 dmol-1, respectively.     

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.     

Wavelength shifts in solid-state circular dichroism spectra: a possible explanation

We have devised an artificial sampling approach generating "absorption flattening" (AF) on the UV-Vis spectrum of a solution of a chiral compound: Tris (ethylendiamine)cobalt(III) chloride (Lambda-Co(en)(3)Cl(3)). We have observed a concomitant red shift of the CD maxima. Related CD and absorption spectra have been calculated from spectra recorded on diluted homogeneous solutions, thus monitoring the effect of AF on both types of data. Experimental data are in good agreement with calculated spectra. Simulations with suitable bandshapes show that the red shift of the CD spectrum is due to AF. On the basis of these results, we conclude that AF is an important cause of distortions in CD spectra for inhomogeneous samples. Plans to compensate or at least to take into account this effect are presented.     

Experimental and calculated circular dichroism spectra of monoaza[5]helicenes

Circular dichroism (CD) spectra have been measured in the range of 400-200 nm on CH₃OH solutions of both enantiomers for the almost complete series of monoaza[5]helicenes, namely the molecules where the hetero N atom occupies positions 1, 3, 4, 5, 6, and 7, respectively (the 2 isomer is missing due to difficulties in the synthesis). CD spectra recorded at controlled room temperature allow one to define precise racemization rates, that are nicely interpreted on the basis of DFT molecular orbital calculations. Time-dependent DFT methods provide us with calculated CD and UV spectra, that are compared with the corresponding experimental data. We discuss the role of the N atom in determining the height of the racemization barrier and in shaping the appearance of the CD spectra.     

Estimation of protein secondary structure from circular dichroism spectra: comparison of CONTIN, SELCON, and CDSSTR methods with an expanded reference set

We have expanded the reference set of proteins used in SELCON3 by including 11 additional proteins (selected from the reference sets of Yang and co-workers and Keiderling and co-workers). Depending on the wavelength range and whether or not denatured proteins are included in the reference set, five reference sets were constructed with the number of reference proteins varying from 29 to 48. The performance of three popular methods for estimating protein secondary structure fractions from CD spectra (implemented in software packages CONTIN, SELCON3, and CDSSTR) and a variant of CONTIN, CONTIN/LL, that incorporates the variable selection method in the locally linearized model in CONTIN, were examined using the five reference sets described here, and a 22-protein reference set. Secondary structure assignments from DSSP were used in the analysis. The performances of all three methods were comparable, in spite of the differences in the algorithms used in the three software packages. While CDSSTR performed the best with a smaller reference set and larger wavelength range, and CONTIN/LL performed the best with a larger reference set and smaller wavelength range, the performances for individual secondary structures were mixed. Analyzing protein CD spectra using all three methods should improve the reliability of predicted secondary structural fractions. The three programs are provided in CDPro software package and have been modified for easier use with the different reference sets described in this paper. CDPro software is available at the website: http://lamar.colostate.edu/ approximately sreeram/CDPro.