Comparison of and limits of accuracy for statistical analyses of vibrational and electronic circular dichroism spectra in terms of correlations to and predictions of protein secondary structure

Inhibition of Human Colorectal Cancer by a Natural Product 7-Acetylhorminone and Interactions with BSA/HSA: Multispectral Analysis and In Silico and In Vitro Studies

We have semi-synthesized a natural product 7-acetylhorminone from crude extract of Premna obtusifolia (Indian headache tree), which is active against colorectal cancer after probation through computational screening methods as it passed through the set parameters of pharmacokinetics (most important nonblood-brain barrier permeant) and drug likeliness (e.g., Lipinski's, Ghose's, Veber's rule) which most other phytoconstituents failed to pass combined with docking with EGFR protein which is highly upregulated in the colorectal carcinoma cell. The structure of 7-acetylhorminone was confirmed by single crystal X-ray diffraction studies and 1H NMR, 13C NMR, and COSY studies. To validate the theoretical studies, first, in vitro experiments were carried out against human colorectal carcinoma cell lines (HCT116) which revealed the potent cytotoxic efficacy of 7-acetylhorminone and verified preliminary investigation. Second, the drugability of 7-acetylhorminone interaction with serum albumin proteins (HSA and BSA) is evaluated both theoretically and experimentally via steady-state fluorescence spectroscopic studies, circular dichroism, isothermal titration calorimetry, and molecular docking. In summary, this study reveals the applicability of 7-acetylhorminone as a potent drug candidate or as a combinatorial drug against colorectal cancer.

Farnesylation and lipid unsaturation are critical for the membrane binding of the C-terminal segment of G-Protein Receptor Kinase 1

Many proteins are modified by the covalent addition of different types of lipids, such as myristoylation, palmitoylation and prenylation. Lipidation is expected to promote membrane association of proteins. Visual phototransduction involves many lipid-modified proteins. The G-Protein-coupled receptor of rod photoreceptors, rhodopsin, is inactivated by G-Protein-coupled Receptor Kinase 1 (GRK1). The C-terminus of GRK1 is farnesylated and its truncation has been shown to result in a very high decrease of its enzymatic activity, most likely because of the loss of its membrane localization. Little information is available on the membrane binding of GRK1 as well as of most prenylated proteins. Measurements of the membrane binding of the non-farnesylated and farnesylated C-terminal segment of GRK1 were thus performed using lipids typical of those found in rod outer segment disk membranes. Their random coil secondary structure was determined using circular dichroism and infrared spectroscopy. The non-farnesylated C-terminal segment of GRK1 has no surface activity. In contrast, the farnesylated C-terminal segment of GRK1 shows a particularly strong binding to lipid monolayers bearing at least one unsaturated fatty acyl chain. No binding is observed in the presence of monolayers of saturated phospholipids, in agreement with the low affinity of farnesylated Ras proteins for lipids in the liquid-ordered state. Altogether, these data demonstrate that the farnesyl group of the C-terminal segment of GRK1 is mandatory for its membrane binding, which is favored by particular lipids or lipid mixtures. This information will also be useful for the understanding of the membrane binding of other prenylated proteins.

Mini review: Instrumentation for vibrational circular dichroism spectroscopy, still a role for dispersive instruments

Vibrational circular dichroism (VCD) has become a standard method for determination of absolute stereochemistry, particularly now that reliable commercial instrumentation has become available. These instruments use a now well-documented Fourier transform infrared-based approach to measure VCD that has virtually displaced initial dispersive infrared-based designs. Nonetheless, many papers have appeared reporting dispersive VCD data, especially for biopolymers. Instrumentation designed with these original methods, particularly after more recent updates optimizing performance in selected spectral regions, has been shown still to have advantages for specific applications. This article presents a mini-review of dispersive VCD instrument designs and includes sample spectra obtained for various biopolymer (particularly peptide) samples. Complementary reviews of Fourier transform-VCD designs are broadly available.

Comparison between the behavior of different hydrophobic peptides allowing membrane anchoring of proteins

Membrane binding of proteins such as short chain dehydrogenase reductases or tail-anchored proteins relies on their N- and/or C-terminal hydrophobic transmembrane segment. In this review, we propose guidelines to characterize such hydrophobic peptide segments using spectroscopic and biophysical measurements. The secondary structure content of the C-terminal peptides of retinol dehydrogenase 8, RGS9-1 anchor protein, lecithin retinol acyl transferase, and of the N-terminal peptide of retinol dehydrogenase 11 has been deduced by prediction tools from their primary sequence as well as by using infrared or circular dichroism analyses. Depending on the solvent and the solubilization method, significant structural differences were observed, often involving α-helices. The helical structure of these peptides was found to be consistent with their presumed membrane binding. Langmuir monolayers have been used as membrane models to study lipid-peptide interactions. The values of maximum insertion pressure obtained for all peptides using a monolayer of 1,2-dioleoyl-sn-glycero-3-phospho-ethanolamine (DOPE) are larger than the estimated lateral pressure of membranes, thus suggesting that they bind membranes. Polarization modulation infrared reflection absorption spectroscopy has been used to determine the structure and orientation of these peptides in the absence and in the presence of a DOPE monolayer. This lipid induced an increase or a decrease in the organization of the peptide secondary structure. Further measurements are necessary using other lipids to better understand the membrane interactions of these peptides.

Alanine scan of core positions in ubiquitin reveals links between dynamics, stability, and function

Mutations at solvent-inaccessible core positions in proteins can impact function through many biophysical mechanisms including alterations to thermodynamic stability and protein dynamics. As these properties of proteins are difficult to investigate, the impacts of core mutations on protein function are poorly understood for most systems. Here, we determined the effects of alanine mutations at all 15 core positions in ubiquitin on function in yeast. The majority (13 of 15) of alanine substitutions supported yeast growth as the sole ubiquitin. Both the two null mutants (I30A and L43A) were less stable to temperature-induced unfolding in vitro than wild type (WT) but were well folded at physiological temperatures. Heteronuclear NMR studies indicated that the L43A mutation reduces temperature stability while retaining a ground-state structure similar to WT. This structure enables L43A to bind to common ubiquitin receptors in vitro. Many of the core alanine ubiquitin mutants, including one of the null variants (I30A), exhibited an increased accumulation of high-molecular-weight species, suggesting that these mutants caused a defect in the processing of ubiquitin-substrate conjugates. In contrast, L43A exhibited a unique accumulation pattern with reduced levels of high-molecular-weight species and undetectable levels of free ubiquitin. When conjugation to other proteins was blocked, L43A ubiquitin accumulated as free ubiquitin in yeast. Based on these findings, we speculate that ubiquitin's stability to unfolding may be required for efficient recycling during proteasome-mediated substrate degradation.

Near-infrared and mid-infrared Fourier transform vibrational circular dichroism of proteins in aqueous solution

Vibrational circular dichroism (VCD) of a series of proteins in H(2)O solution with differing secondary structure are reported for the first time in the near-infrared (NIR) region as well as the NH-stretching region. The Fourier transform (FT) near-infrared (NIR) measurements were carried out between 6000 to 4000 cm(-1). FT-VCD measurements were simultaneously carried out for the mid-infrared (mid-IR) region from 2000 to 800 cm(-1) for direct comparison to VCD in the NIR region. The NIR VCD spectra of proteins show distinct spectral features for different protein structural motifs, indicating a valuable new method to study protein conformations. The principal VCD transitions in the NIR region are two combination bands, the amide A-II and B-II bands, of the amide A and B fundamentals with the amide II fundamental, and the second overtone of the amide II, referred to as the amide 3 x II band. VCD in the amide A and B band region consisting primarily of NH stretching motions were successfully obtained in H(2)O for the first time for an insulin fibril sample. Similar to the enhanced VCD signal observed in amide I and II regions, the amide A and B VCD of insulin fibril shows strong intensity enhancements, providing an additional valuable probe of protein fibril growth and development in solution. The relative sensitivities of the mid-IR, N-H stretching, and NIR regions are discussed.

Bioinspired interactions studied by vibrational circular dichroism

Vibrational circular dichroism (VCD) spectra are reliable indicators of the spatial structure of chiral molecules. The specific and characteristic feature of vibrational spectroscopy, and therefore also of VCD, where the energy of some vibrational modes is predominantly focused to a specific part of the molecule, enables monitoring both the structure of the molecule dissolved in different solvents and under different physicochemical conditions and molecular interactions. This minireview deals with recent contributions covering structural information on the bioinspired interactions obtained by means of VCD, especially in the following areas: interaction of DNA with biomolecules and biogenic metals, guanine tetramers and quadruplexes, biointeractions of bile pigments, and polypeptide and protein interactions with other biomolecules.

Unique physical properties and interactions of the domains of methylated DNA binding protein 2

Methylated DNA binding protein 2 (MeCP2) is a methyl CpG binding protein whose key role is the recognition of epigenetic information encoded in DNA methylation patterns. Mutation or misregulation of MeCP2 function leads to Rett syndrome as well as a variety of other autism spectrum disorders. Here, we have analyzed in detail the properties of six individually expressed human MeCP2 domains spanning the entire protein with emphasis on their interactions with each other, with DNA, and with nucleosomal arrays. Each domain contributes uniquely to the structure and function of the full-length protein. MeCP2 is approximately 60% unstructured, with nine interspersed alpha-molecular recognition features (alpha-MoRFs), which are polypeptide segments predicted to acquire secondary structure upon forming complexes with binding partners. Large increases in secondary structure content are induced in some of the isolated MeCP2 domains and in the full-length protein by binding to DNA. Interactions between some MeCP2 domains in cis and trans seen in our assays likely contribute to the structure and function of the intact protein. We also show that MeCP2 has two functional halves. The N-terminal portion contains the methylated DNA binding domain (MBD) and two highly disordered flanking domains that modulate MBD-mediated DNA binding. One of these flanking domains is also capable of autonomous DNA binding. In contrast, the C-terminal portion of the protein that harbors at least two independent DNA binding domains and a chromatin-specific binding domain is largely responsible for mediating nucleosomal array compaction and oligomerization. These findings led to new mechanistic and biochemical insights regarding the conformational modulations of this intrinsically disordered protein, and its context-dependent in vivo roles.

Rett syndrome-causing mutations in human MeCP2 result in diverse structural changes that impact folding and DNA interactions

Most cases of Rett syndrome (RTT) are caused by mutations in the methylated DNA-binding protein, MeCP2. Here, we have shown that frequent RTT-causing missense mutations (R106W, R133C, F155S, T158M) located in the methylated DNA-binding domain (MBD) of MeCP2 have profound and diverse effects on its structure, stability, and DNA-binding properties. Fluorescence spectroscopy, which reports on the single tryptophan in the MBD, indicated that this residue is strongly protected from the aqueous environment in the wild type but is more exposed in the R133C and F155S mutations. In the mutant proteins R133C, F155S, and T158M, the thermal stability of the domain was strongly reduced. Thermal stability of the wild-type protein was increased in the presence of unmethylated DNA and was further enhanced by DNA methylation. DNA-induced thermal stability was also seen, but to a lesser extent, in each of the mutant proteins. Circular dichroism (CD) of the MBD revealed differences in the secondary structure of the four mutants. Upon binding to methylated DNA, the wild type showed a subtle but reproducible increase in alpha-helical structure, whereas the F155S and R106W did not acquire secondary structure with DNA. Each of the mutant proteins studied is unique in terms of the properties of the MBD and the structural changes induced by DNA binding. For each mutation, we examined the extent to which the magnitude of these differences correlated with the severity of RTT patient symptoms.

K2D2: estimation of protein secondary structure from circular dichroism spectra

Background

Circular dichroism spectroscopy is a widely used technique to analyze the secondary structure of proteins in solution. Predictive methods use the circular dichroism spectra from proteins of known tertiary structure to assess the secondary structure contents of a protein with unknown structure given its circular dichroism spectrum.

Results

We developed K2D2, a method with an associated web server to estimate protein secondary structure from circular dichroism spectra. The method uses a self-organized map of spectra from proteins with known structure to deduce a map of protein secondary structure that is used to do the predictions.

Conclusion

The K2D2 server is publicly accessible at . It accepts as input a circular dichroism spectrum and outputs the estimated secondary structure content (alpha-helix and beta-strand) of the corresponding protein, as well as an estimated measure of error.

Secondary structure of a truncated form of lecithin retinol acyltransferase in solution and evidence for its binding and hydrolytic action in monolayers

Lecithin retinol acyltransferase (LRAT) is a 230 amino acids membrane-associated protein which catalyzes the esterification of all-trans-retinol into all-trans-retinyl ester. The enzymatic activity of a truncated form of LRAT (tLRAT) which contains the residues required for catalysis but which is lacking N- and C-terminal hydrophobic segments has been shown to depend on the detergent used for its solubilization. Moreover, it is unknown whether tLRAT can bind membranes in the absence of these hydrophobic segments. The present study has allowed to measure the membrane binding and hydrolytic action of tLRAT in lipid monolayers by use of polarization modulation infrared reflection absorption spectroscopy and Brewster angle microscopy. Moreover, the proportion of the secondary structure components of tLRAT was determined in three different detergents by infrared absorption spectroscopy, vibrational circular dichroism and electronic circular dichroism which allowed to explain its detergent dependent activity. In addition, the secondary structure of tLRAT in the absence of detergent was very similar to that in Triton X-100 thus suggesting that, compared to the other detergents assayed, the secondary structure of this protein is very little perturbed by this detergent.

Concentration-independent estimation of protein secondary structure by circular dichroism: a comparison of methods

Estimation of a protein's secondary structure from its circular dichroism spectrum usually requires accurate knowledge of the concentration and pathlength of the sample. Two recently described methods avoid this problem by analysis of g-factor spectra (McPhie, Anal. Biochem. 293, 109-119) or scaling of relative intensities (Raussens et al., Anal. Biochem. 319, 114-121). Application of the two methods to the same samples shows that they can have similar efficacies. Calculation with the latter method is more rapid, but the performance of the former is maintained over reduced wavelength ranges.

Flexibility of the cytoplasmic domain of the phototaxis transducer II from Natronomonas pharaonis

Chemo- and phototaxis systems in bacteria and archaea serve as models for more complex signal transduction mechanisms in higher eukaryotes. Previous studies of the cytoplasmic fragment of the phototaxis transducer (pHtrII-cyt) from the halophilic archaeon Natronomonas pharaonis showed that it takes the shape of a monomeric or dimeric rod under low or high salt conditions, respectively. CD spectra revealed only approximately 24% helical structure, even in 4 M KCl, leaving it an open question how the rod-like shape is achieved. Here, we conducted CD, FTIR, and NMR spectroscopic studies under different conditions to address this question. We provide evidence that pHtrII-cyt is highly dynamic with strong helical propensity, which allows it to change from monomeric to dimeric helical coiled-coil states without undergoing dramatic shape changes. A statistical analysis of predicted disorder for homologous sequences suggests that structural flexibility is evolutionarily conserved within the methyl-accepting chemotaxis protein family.

An infrared spectroscopic study of the conformational transition of elastin-like polypeptides

The infrared spectroscopy of elastin-like polypeptides and the relation to the inverse thermal transition are discussed. To correlate the spectroscopic observations with structure a density function theory model was created that captures the essential hydrogen bonding and packing of the beta-spiral structure proposed for elastin and elastin-like polypeptides. The infrared spectrum was calculated using periodic boundary conditions and a method for estimating the difference dipole moment permits both frequencies and intensities to be obtained for the modeling of spectra. The two observed amide I bands at 1615 cm(-1) and 1656 cm(-1) are shown to arise from the beta-spiral structure. The increase in intensity of these bands with increasing salt concentration and temperature is assigned to the closer association of strands of the beta-spiral. The sharp inverse temperature transition is observed within 1 degrees C and involves a change in secondary structure that involves formation of interstrand beta-sheets for approximately 25% of the amino acids. This conclusion is consistent with available data and simulations that have been reported to date.

Prediction of protein secondary structure content using amino acid composition and evolutionary information

Knowing protein structure and inferring its function from the structure are one of the main issues of computational structural biology, and often the first step is studying protein secondary structure. There have been many attempts to predict protein secondary structure contents. Previous attempts assumed that the content of protein secondary structure can be predicted successfully using the information on the amino acid composition of a protein. Recent methods achieved remarkable prediction accuracy by using the expanded composition information. The overall average error of the most successful method is 3.4%. Here, we demonstrate that even if we only use the simple amino acid composition information alone, it is possible to improve the prediction accuracy significantly if the evolutionary information is included. The idea is motivated by the observation that evolutionarily related proteins share the similar structure. After calculating the homolog-averaged amino acid composition of a protein, which can be easily obtained from the multiple sequence alignment by running PSI-BLAST, those 20 numbers are learned by a multiple linear regression, an artificial neural network and a support vector regression. The overall average error of method by a support vector regression is 3.3%. It is remarkable that we obtain the comparable accuracy without utilizing the expanded composition information such as pair-coupled amino acid composition. This work again demonstrates that the amino acid composition is a fundamental characteristic of a protein. It is anticipated that our novel idea can be applied to many areas of protein bioinformatics where the amino acid composition information is utilized, such as subcellular localization prediction, enzyme subclass prediction, domain boundary prediction, signal sequence prediction, and prediction of unfolded segment in a protein sequence, to name a few.

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.

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.

The optimization of protein secondary structure determination with infrared and circular dichroism spectra

We have used the circular dichroism and infrared spectra of a specially designed 50 protein database [Oberg, K.A., Ruysschaert, J.M. & Goormaghtigh, E. (2003) Protein Sci. 12, 2015-2031] in order to optimize the accuracy of spectroscopic protein secondary structure determination using multivariate statistical analysis methods. The results demonstrate that when the proteins are carefully selected for the diversity in their structure, no smaller subset of the database contains the necessary information to describe the entire set. One conclusion of the paper is therefore that large protein databases, observing stringent selection criteria, are necessary for the prediction of unknown proteins. A second important conclusion is that only the comparison of analyses run on circular dichroism and infrared spectra independently is able to identify failed solutions in the absence of known structure. Interestingly, it was also found in the course of this study that the amide II band has high information content and could be used alone for secondary structure prediction in place of amide I.

Optical spectroscopic differentiation of various equilibrium denatured states of horse cytochrome c

Thermal unfolding of cytochrome c (cyt c) from several states has been studied using equilibrium spectroscopic techniques. CD in the uv, vibrational circular dichroism, infrared, and uv-vis absorption spectra measured at various temperatures, pHs, salt concentrations, and GuHCl concentrations are used to show the conformational as well as heme structural differences between native and various denatured states. The difference in thermal denaturation behaviors of cyt c starting from acid denatured, molten globule (MG), and the A and native states are explored. Different final high temperature states were observed for cytochrome c unfolding from four different initial states (native, MG, A, and acid denatured state) by electronic CD, Fourier transform infrared (FTIR), and vibrational CD (VCD). Consistent with this, different thermal unfolding pathways for the MG and A states are suggested by the FTIR and VCD data for this process.

Individual tyrosine side-chain contributions to circular dichroism of ribonuclease

Experimental and theoretical studies using site-directed mutants of ribonuclease A (RNase A) offer more extensive information on the tyrosine side-chain contributions to the circular dichroism (CD) of the enzyme. Bovine pancreatic RNase A has three exposed tyrosine residues (Tyr73, Tyr76, and Tyr115) and three buried tyrosine residues (Tyr25, Tyr92 and Tyr97). The difference CD spectra between the wild type and the mutants at pH 7.0 (Deltaepsilon(277,wt) - Deltaepsilon(277,mut)) show bands with more negative DeltaDeltaepsilon(277) values for Y73F and Y115F than those for Y25F and Y92F and bands with positive DeltaDeltaepsilon(277) values for Y76F and Y97F. The theoretical calculations are in good semiquantitative agreement for all the mutants. The pH difference spectrum (pH 11.3-7.0) for the wild type shows a negative band at 295 nm and an enhanced positive band at 245 nm. The three mutants at buried tyrosine sites and one mutant at an exposed tyrosine site (Y76F) exhibit pH-difference spectra that are similar to that of the wild type. In contrast, two mutants at exposed tyrosine sites (Y73F and Y115F) exhibit diminished 295-nm negative bands and, instead of positive bands at 245 nm, negative bands are observed. Our results indicate that Tyr73 and Tyr115, two of the exposed tyrosine residues, are the largest contributors to the 277- and 245-nm CD bands of RNaseA, but the buried tyrosine residues and the one remaining exposed residue also contribute to these bands. Disulfide contributions to the 277- and 240-nm bands and the peptide contribution to the 240-nm band are confirmed theoretically.

Calculations of protein circular dichroism from first principles

Understanding the relationship between the amino acid sequence of a protein and its unique, compact 3D structure is one of the grand challenges in molecular biophysics. One particularly exciting approach is time-resolved electronic circular dichroism (CD) spectroscopy, which offers resolution on a nanosecond (or faster) time scale, although it does not provide the spatial resolution of techniques like X-ray crystallography or NMR. The thrust of our work is to underpin fast time scale spectroscopic studies of protein folding with a stronger theoretical foundation. Ultimately, we seek to use molecular dynamics simulations to study the influence of conformational dynamics and conformational transitions on the electronic CD spectra of proteins. We discuss how improved quantum chemical models of individual chromophores, including aromatic sidechains, can be incorporated into calculations of the electronic structure of proteins and their CD.

Rationally selected basis proteins: a new approach to selecting proteins for spectroscopic secondary structure analysis

Protein basis sets have been extensively used as reference data for the determination of protein structure with optical methods such as circular dichroism and infrared spectroscopies. We have taken a new approach to basis protein selection by utilizing three crystal structure classification databases: CATH, SCOP, and PDB_SELECT. Through the use of the information available in these and other online resources, we identified 115 commercially available proteins as potential basis set candidates. By carefully screening the quality of the crystal structures and commercial protein preparations, we obtained a final set of 50 rationally selected proteins (RaSP50) that has been optimized for use in spectroscopic protein structure determination studies. These proteins span the full range of known protein folds as well as alpha-helix and beta-sheet contents, and they represent a more comprehensive variety of fold types than any previous reference set. This report includes a detailed presentation of the reasoning behind the rational protein selection process, a description of the properties of the RaSP50 set, and a discussion of the types of structural and spectral variations that are represented in the set.

Optical spectroscopic investigations of model beta-sheet hairpins in aqueous solution

In this contribution we report optical spectroscopic data on a series of designed beta hairpins previously shown by NMR to contain a substantial population of beta-sheet structure. These models contain a designed hydrophobic cluster and a (D)Pro-Gly sequence to promote formation of a turn geometry. FTIR, electronic and vibrational CD (ECD and VCD) spectra for these small peptides are comparable to expected bandshapes for peptides of high beta-sheet content. The (D)Pro-Gly sequence provides a better turn motif than Asn-Gly as measured by its beta-sheet spectral characteristics. IR and VCD spectra are in qualitative agreement with theoretical simulations based on transfer of parameters from ab initio quantum mechanical force field and intensity computations for the turn and strands. These calculations provide assignments for some distinguishing modes in both IR and VCD spectra. Increased sheet structure can be induced in these hairpins by use of mixed solvent conditions. Thermal denaturation studies reveal that these hairpins undergo very broad unfolding transitions. Guanidine hydrochloride unfolding transitions for the selected hairpin models are similarly broad. However, the "end-states" of temperature and chaotropic denaturation are spectroscopically differentiable.

Structural composition of betaI- and betaII-proteins

Circular dichroism spectra of proteins are sensitive to protein secondary structure. The CD spectra of alpha-rich proteins are similar to those of model alpha-helices, but beta-rich proteins exhibit CD spectra that are reminiscent of CD spectra of either model beta-sheets or unordered polypeptides. The existence of these two types of CD spectra for beta-rich proteins form the basis for their classification as betaI- and betaII-proteins. Although the conformation of beta-sheets is largely responsible for the CD spectra of betaI-proteins, the source of betaII-protein CD, which resembles that of unordered polypeptides, is not completely understood. The CD spectra of unordered polypeptides are similar to that of the poly(Pro)II helix, and the poly(Pro)II-type (P2) structure forms a significant fraction of the unordered conformation in globular proteins. We have compared the beta-sheet and P2 structure contents in beta-rich proteins to understand the origin of betaII-protein CD. We find that betaII-proteins have a ratio of P2 to beta-sheet content greater than 0.4, whereas for betaI-proteins this ratio is less than 0.4. The beta-sheet content in betaI-proteins is generally higher than that in betaII-proteins. The origin of two classes of CD spectra for beta-rich proteins appears to lie in their relative beta-sheet and P2 structure contents.

Discriminating 3(10)- from alpha-helices: vibrational and electronic CD and IR absorption study of related Aib-containing oligopeptides

Model peptides based on -(Aib-Ala)(n)-, and (Aib)(n)-Leu-(Aib)(2) sequences, which have varying amounts of 3(10)-helical character, were studied by use of vibrational and electronic circular dichroism (VCD and ECD) and Fourier transform infrared (FTIR) absorption spectroscopies to test the correlation of spectral response and conformation. The data indicate that these peptides, starting from a length of about four to six residues, predominantly adopt a 3(10)-helical conformation at room temperature. The longest model peptides, depending on the series, may evidence some alpha-helical contribution to the spectra, while the shorter ones, with less than six residues, have much less order. The IR absorption spectra (as supported by theory) showed only small frequency changes between 3(10)- and alpha-helices. By contrast, solvent effects are a source of much bigger perturbations. The ECD results show that the intensity ratio for the approximately 222-nm to approximately 208-nm bands, while useful for distinguishing between these two helical types in some sequences, may have a narrower range of application than VCD. However, the VCD data presented here continue to support the proposed discrimination between alpha- and 3(10)-helices based on qualitative amide I and II bandshape differences. The present study shows the intensities of the 3(10)-helical amide I (peak-to-peak) to its amide II VCD to be of the same order and useful for discriminating them from alpha-helices, whose amide I dominates the amide II in intensity. This qualitative result is experimentally independent of the amount of alphaMe-substituted residues in the sequence. These experimental VCD results are consistent in detail with theoretical spectral simulations for Ac-(Ala)(8)-NH(2), Ac-(Aib-Ala)(4)-NH(2), and Ac-(Aib)(8)-NH(2) in 3(10)- and alpha-helical conformations.

Purification and structural characterization of the central hydrophobic domain of oleosin

The oil bodies of rapeseeds contain a triacylglycerol matrix surrounded by a monolayer of phospholipids embedded with abundant structural alkaline proteins termed oleosins and some other minor proteins. Oleosins are unusual proteins because they contain a 70-80-residue uninterrupted nonpolar domain flanked by relatively polar C- and N-terminal domains. Although the hydrophilic N-terminal domain had been studied, the structural feature of the central hydrophobic domain remains unclear due to its high hydrophobicity. In the present study, we reported the generation, purification, and characterization of a 9-kDa central hydrophobic domain from rapeseed oleosin (19 kDa). The 9-kDa central hydrophobic domain was produced by selectively degrading the N and C termini with enzymes and then purifying the digest by SDS-PAGE and electroelution. We have also reconstituted the central domain into liposomes and synthetic oil bodies to determine the secondary structure of the domain using CD and Fourier transform infrared (FTIR) spectroscopy. The spectra obtained from CD and FTIR were analyzed with reference to structural information of the N-terminal domain and the full-length rapeseed oleosin. Both CD and FTIR analysis revealed that 50-63% of the domain was composed of beta-sheet structure. Detailed analysis of the FTIR spectra indicated that 80% of the beta-sheet structure, present in the central domain, was arranged in parallel to the intermolecular beta-sheet structure. Therefore, interactions between adjacent oleosin proteins would give rise to a stable beta-sheet structure that would extend around the surface of the seed oil bodies stabilizing them in emulsion systems. The strategies used in our present study are significant in that it could be generally used to study difficult proteins with different independent structural domains, especially with long hydrophobic domains.

Pretransitional structural changes in the thermal denaturation of ribonuclease S and S protein

Two mechanisms have been proposed for the thermal unfolding of ribonuclease S (RNase S). The first is a sequential partial unfolding of the S peptide/S protein complex followed by dissociation, whereas the second is a concerted denaturation/dissociation. The thermal denaturation of ribonuclease S and its fragment, the S protein, were followed with circular dichroism and infrared spectra. These spectra were analyzed by the principal component method of factor analysis. The use of multiple spectral techniques and of factor analysis monitored different aspects of the denaturation simultaneously. The unfolding pathway was compared with that of the parent enzyme ribonuclease A (RNase A), and a model was devised to assess the importance of the dissociation in the unfolding. The unfolding patterns obtained from the melting curves of each protein imply the existence of multiple intermediate states and/or processes. Our data provide evidence that the pretransition in the unfolding of ribonuclease S is due to partial unfolding of the S protein/S peptide complex and that the dissociation occurs at higher temperature. Our observations are consistent with a sequential denaturation mechanism in which at least one partial unfolding step comes before the main conformational transition, which is instead a concerted, final unfolding/dissociation step.

Protein and peptide secondary structure and conformational determination with vibrational circular dichroism

Vibrational circular dichroism (VCD) provides alternative views of protein and peptide conformation with advantages over electronic (UV) CD (ECD) or IR spectroscopy. VCD is sensitive to short-range order, allowing it to discriminate beta-sheet and various helices as well as disordered structure. Quantitative secondary structure analyses use protein VCD bandshapes, but are best combined with ECD and IR for balance. Much recent work has focused on empirical and theoretical VCD analyses of peptides, with detailed prediction of helix, sheet and hairpin spectra and site-specific application of isotopic substitution for structure and folding.

Thermal unfolding of ribonuclease A in phosphate at neutral pH: deviations from the two-state model

The thermal denaturation of ribonuclease A (RNase A) in the presence of phosphate at neutral pH was studied by differential scanning calorimetry (DSC) and a combination of optical spectroscopic techniques to probe the existence of intermediate states. Fourier transform infrared (FTIR) spectra of the amide I' band and far-uv circular dichroism (CD) spectra were used to monitor changes in the secondary structure. Changes in the tertiary structure were monitored by near-uv CD. Spectral bandshape changes with change in temperature were analyzed using factor analysis. The global unfolding curves obtained from DSC confirmed that structural changes occur in the molecule before the main thermal denaturation transition. The analysis of the far-uv CD and FTIR spectra showed that these lower temperature-induced modifications occur in the secondary structure. No pretransition changes in the tertiary structure (near-uv CD) were observed. The initial changes observed in far-uv CD were attributed to the fraying of the helical segments, which would explain the loss of spectral intensity with almost no modification of spectral bandshape. Separate analyses of different regions of the FTIR amide I' band indicate that, in addition to alpha-helix, part of the pretransitional change also occurs in the beta-strands.

SOMCD: method for evaluating protein secondary structure from UV circular dichroism spectra

This article presents SOMCD, an improved method for the evaluation of protein secondary structure from circular dichroism spectra, based on Kohonen's self-organizing maps (SOM). Protein circular dichroism (CD) spectra are used to train a SOM, which arranges the spectra on a two-dimensional map. Location in the map reflects the secondary structure composition of a protein. With SOMCD, the prediction of beta-turn has been included. The number of spectra in the training set has been increased, and it now includes 39 protein spectra and 6 reference spectra. Finally, SOM parameters have been chosen to minimize distortion and make the network produce clusters with known properties. Estimation results show improvements compared with the previous version, K2D, which, in addition, estimated only three secondary structure components; the accuracy of the method is more uniform over the different secondary structures.

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.

Simulations of oligopeptide vibrational CD: effects of isotopic labeling

Simulated ir absorption and vibrational CD (VCD) spectra of four alanine-based octapeptides, each having its main chain constrained to a different secondary structure conformation, were analyzed and compared with experimental results for several different peptides. The octapeptide simulations were based on transfer of property tensors from a series of ab initio calculations for a short L-alanine based segment containing 3 peptide bonds with relative straight phi, psi angles fixed to those appropriate for alpha-helix, 3(10)-helix, ProII-like helix, and beta-sheet-like strand. The tripeptide force field (FF) and atomic polar tensors were obtained with density functional theory techniques at the BPW91/6-31G** level and the atomic axial tensor at the mixed BPW91/6-31G**/HF/6-31G level. Allowing for frequency correction due to the FF limitations, the octapeptide results obtained are qualitatively consistent with experimental observations for ir and VCD spectra of polypeptides and oligopeptides in established conformations. In all cases, the correct VCD sign patterns for the amide I and II bands were predicted, but the intensities did have some variation from the experimental patterns. Predicted VCD changes upon deuteration of either the peptide or side-chains as well as for (13)C isotopic labeling of the amide C=O at specific sites in the peptide chain were computed for analysis of experimental observations. A combination of theoretical modeling with experimental data for labeled compounds leads both to enhanced resolution of component transitions and added conformational applicability of the VCD spectra.

Enhanced prediction accuracy of protein secondary structure using hydrogen exchange Fourier transform infrared spectroscopy

A novel equilibrium hydrogen exchange Fourier transform IR (HX-FTIR) spectroscopy method for predicting secondary structure content was employed using spectra obtained for a training set of 23 globular proteins. The IR bandshape and frequency changes resulting from controlled levels of H-D exchange were observed to be protein-dependent. Their analysis revealed these variations to be partly correlated to secondary structure. For each protein, a set of 6 spectra was measured with a systematic variation of the solvent H-D ratio and was subjected to factor analysis. The most significant component spectra for each protein, representing independent aspects of the spectral response to deuteration, were each subjected to a second factor analysis over the entire training set. Restricted multiple regression (RMR) analysis using the loadings of the principal components from 19 of these H-D analyses revealed an improvement in prediction accuracy compared with conventional bandshape-based analyses of FTIR data. Nearly a factor of 2 reduction in error for prediction of helix fractions was found using s1, the average spectral response for the H-D set. In some cases, significant error reduction for prediction of minor components was found using higher factors. Using the same analytical methods, prediction errors with this new deuteration-response-FTIR method were shown to be even better than those obtained by use of electronic circular dichroism (ECD) data for helix predictions and to be significantly lower for ECD-based sheet prediction, making these the best secondary structure predictions obtained with the RMR method. Tests of a limited variable selection scheme showed further improvements, consistent with previous results of this approach using ECD data.

Novel matrix descriptor for secondary structure segments in proteins: demonstration of predictability from circular dichroism spectra

An extension to standard protein secondary structure predictions using optical spectra that encompasses the number and average lengths of segments of uniform secondary structure in the sequence is demonstrated. The connectivity and numbers of segments can be described by a matrix descriptor [sij] (i, j representing segment types such as helix and beta-sheet strands). Independent knowledge of the fractional concentration of each secondary structure type and of the total number of residues in the protein then with [sij] yields the average segment length of each type. The physical background for prediction of this extended structural descriptor from spectral data is summarized, rules for its generation from reference X-ray structures are defined, and formal variants of its form are discussed. Using a novel neural network approach to analyze a training set of electronic circular dichroism (ECD) and vibrational circular dichroism (VCD) spectra for 23 proteins, matrix descriptors encompassing helix, sheet, and other forms are predicted. The results show that the matrix descriptor can be predicted to an accuracy comparable to that of conventionally predicted average fractional secondary structures. In this respect the ECD predictions of [sij] were significantly more accurate than the VCD ones, which may result from the longer range length dependence of the ECD bandshape and intensity. Summary results for a parallel analysis using Fourier transform infrared spectra indicate somewhat lower reliability than those for VCD.

Predictions of protein secondary structures using factor analysis on Fourier transform infrared spectra: effect of Fourier self-deconvolution of the amide I and amide II bands

Fourier self-deconvolution (FSD) was performed on protein amide I and II Fourier transform infrared (FTIR) spectra to test if the resultant increased band shape variation would lead to improvements in protein secondary structure prediction with our factor analysis based restricted multiple regression (RMR) methods. FTIR spectra of 23 proteins dissolved in H2O were measured and normalized to a constant amide I peak absorbance. The deconvolved spectra were renormalized by area so that the deconvolved spectra sets had the same area as before. Principal component analysis of the deconvolved spectra sets was carried out, which was followed by a selective multiple linear regression (RMR) analysis of the principal component loadings with regard to the fractional components (FC) of secondary structure. As compared to analyses based on the original spectra set, helix and sheet predictions were not noticeably improved by FSD; but, if a very large number of component spectra (16) were retained in the pool to select which loadings to be used in the RMR optimization, better predictions of turn and "other" resulted. The prediction quality varied depending on the deconvolution parameters used.

A new attenuated total reflectance Fourier transform infrared spectroscopy method for the study of proteins in solution

An attenuated total reflectance Fourier transform infrared method has been developed that allows collection of spectra from proteins in solution. This method eliminates any structural perturbations induced by the internal reflection element (IRE), and thus the spectra reflect the solution conformation of the protein. A key feature of the method is subtraction of the signal from any protein adsorbed to the IRE. The advantages of this method include the small amount of sample required and the high sampling rate. Attenuated total reflectance (ATR)-Fourier transform infrared spectroscopy (FTIR) is more versatile than transmission FTIR because it is possible to collect spectra of nontransparent samples, to use samples of very low protein concentration (< or = 0.3 mg/ml), and to study proteins in the presence of strongly absorbing solutes (such as denaturants). The experimental procedures and data processing routines developed were evaluated by collecting spectra from a set of 13 proteins and evaluating their accuracy with a partial least-squares analysis. The relative mean and standard deviation errors for the basis set analysis were 6.3% for alpha-helix, 5.9% for beta-sheet/extended structure, and 4.4% for turn, which are similar to values from comparable analyses of transmission FTIR spectra. In addition, a detailed comparison between this solution ATR method and the hydrated thin-film ATR technique is presented.