ExSer: A standalone tool to mine protein data bank (PDB) for secondary structural elements

Detailed structural analysis of protein necessitates investigation at primary, secondary and tertiary levels, respectively. Insight into protein secondary structures pave way for understanding the type of secondary structural elements involved (α-helices, β-strands etc.), the amino acid sequence that encode the secondary structural elements, number of residues, length and, percentage composition of the respective elements in the protein. Here we present a standalone tool entitled "ExSer" which facilitate an automated extraction of the amino acid sequence that encode for the secondary structural regions of a protein from the protein data bank (PDB) file. AVAILABILITY: ExSer is freely downloadable from http://code.google.com/p/tool-exser/     

Determination of protein fold class from Raman or Raman optical activity spectra using random forests

Knowledge of the fold class of a protein is valuable because fold class gives an indication of protein function and evolution. Fold class can be accurately determined from a crystal structure or NMR structure, though these methods are expensive, time-consuming, and inapplicable to all proteins. In contrast, vibrational spectra [infra-red, Raman, or Raman optical activity (ROA)] are rapidly obtained for proteins under wide range of biological molecules under diverse experimental and physiological conditions. Here, we show that the fold class of a protein can be determined from Raman or ROA spectra by converting a spectrum into data of 10 cm⁻¹ bin widths and applying the random forest machine learning algorithm. Spectral data from 605 and 1785 cm⁻¹ were analyzed, as well as the amide I, II, and III regions in isolation and in combination. ROA amide II and III data gave the best performance, with 33 of 44 proteins assigned to one of the correct four top-level structural classification of proteins (SCOP) fold class (all α, all β, α and β, and disordered). The method also shows which spectral regions are most valuable in assigning fold class.     

A new catalog of protein beta-sheets

Systematic protein folding studies depend on protein three-dimensional structure annotation, the assignment of amino acid structural types from atomic coordinates. Significant stabilizing factors between adjacent beta-sheet peptide chains have recently been characterized and were not considered during the development of previously published annotation methods. To produce an accurate beta-sheet domain catalog and to encompass the full beta-sheet spectacle, we developed a method, beta-Spider, which evaluates a packing energy between adjacent peptide chains in accordance with the newly discovered stabilizing factors. While considering important energetic factors, our approach also minimizes the use of subjective criteria, such as (phi,psi) boundaries and sets of H-bonding motifs that are used in other existing methods. As a result of the application of beta-Spider to a set of available high-resolution X-ray crystal structures, we present here a new beta-sheet catalog that differs considerably from the one produced by the most acclaimed DSSP method. The catalog includes new H-bonding motifs that were never reported.     

Conformational change from antiparallel beta-sheet to alpha-helix in a series of depsipeptide, -(Leu-Leu-Lac)(n)-: syntheses, spectroscopic studies, and crystal structures of Boc-Leu-Lac-OEt and Boc-(Leu-Leu-Lac)(n)-OEt (n = 1, 2)

The depsipeptides Boc-Leu-Lac-OEt (1) and Boc-(Leu-Leu-Lac)(n)-OEt (n = 1, 2) (2 and 3, respectively) (Boc = tert-butyloxycarbonyl, Lac = L-lactic acid residue) has been synthesized and studied by crystallographic, CD spectroscopic, and ESI-MS analyses. In the packing cells, those three compounds adopt beta-strand conformations. Each molecule is linked into a dimer (1) or an infinite assembly (2 and 3) by tight hydrogen bonds of the type NH...O==C. Interestingly, the hexamer, 3 shows the first example of antiparallel pleated beta-sheet crystal structure for a depsipeptide molecule. In the packing cells, especially for 3, the ester groups O--C==O are perpendicularly oriented to the amide groups NH--C==O and beta-sheet planes to avoid the interaction between --O--(ester) and O==C. Therefore, when the chain length become longer, the O...O==C repulsion interaction works as a beta-sheet breaker and hence promotes an alpha-helical structure as observed for Boc-(Leu-Leu-Lac)(3)-Leu-Leu-OEt (4) (Oku et al. Biopolymers 2004, 75, 242-254) and Boc-(Leu-Leu-Lac)(n)-OEt (n = 4-6) (5-7) (Katakai et al., Biopolymers 1996, 38, 285-290), in which the O...O==C repulsion does not cause significant structural changes in alpha-helical main chains. Therefore from the structural and spectroscopic analyses, we have found governing factors for the specificity in the beta-sheet and alpha-helix decision in this series of depsipeptides, -(Leu-Leu-Lac)(n)-.     

Addition of a peptide fragment on an alpha-helical depsipeptide induces alpha/3(10)-conjugated helix: synthesis, crystal structure, and CD spectra of Boc-Leu-Leu-Ala-(Leu-Leu-Lac)3-Leu-Leu-OEt

The depsipeptide Boc(1)-Leu(2)-Leu(3)-Ala(4)-Leu(5)-Leu(6)-Lac(7)-Leu(8)-Leu(9)-Lac(10)-Leu(11)-Leu(12)-Lac(13)-Leu(14)-Leu(15)-OEt(16) (1) (Boc = tert-butyloxycarbonyl, Lac = L-lactic acid residue) has been synthesized from the peptide Boc-Leu-Leu-Ala-OEt (2) and a depsipeptide, Boc-(Leu-Leu-Lac)(3)-Leu-Leu-OEt (3). Single crystals of 1 were successfully obtained and the structure has been solved by direct methods (such as Sir2002 and Shake-and-Bake). Interestingly, 1 adopts an alpha/3(10)-conjugated helix containing a kink at the junction of peptide and depsipeptide segments, Leu3-Lac7. This is significantly different from the conformation of 3, which has a straight alpha-helical structure with standard phi and psi angles. Microcrystalline CD spectra were also studied to compare structural properties of 1 and 3. The differences between alpha/3(10)- and alpha-helices appear in these CD spectra.     

gamma-Zein secondary structure in solution by circular dichroism

The proline-rich N-terminal domain of gamma-zein has been reported in relevant processes, which include its ability to cross the cell membranes. Evidences indicate that synthetic hexapeptide (PPPVHL), naturally found in N-terminal portion of gamma-zein, can adopt the polyproline II (PPII) conformation in aqueous solution. The secondary structure of gamma-zein in maize protein bodies had been analyzed by solid state Fourier transform infrared and nuclear magnetic resonance spectroscopies. However, it was not possible to measure PPII content in physiological environment since the beta-sheet and PPII signals overlap in both solid state techniques. Here, the secondary structure of gamma-zein has been analyzed by circular dichroism in SDS aqueous solution with and without ditiothreitol (DTT), and in 60% of 2-propanol and water with DTT. The results show that gamma-zein has high helical content in all solutions. The PPII conformation was present at about 7% only in water/DTT solution.     

The human prion protein alpha2 helix: a thermodynamic study of its conformational preferences

We have synthesized both free and terminally-blocked peptide corresponding to the second helical region of the globular domain of normal human prion protein, which has recently gained the attention of structural biologists because of a possible role in the nucleation process and fibrillization of prion protein. The profile of the circular dichroism spectrum of the free peptide was that typical of alpha-helix, but was converted to that of beta-structure in about 16 h. Instead, below 2.1 x 10(-5) M, the spectrum of the blocked peptide exhibited a single band centered at 200 nm, unequivocally associated to random conformations, which did not evolve even after 24 h. Conformational preferences of this last peptide have been investigated as a function of temperature, using trifluoroethanol or low-concentration sodium dodecyl sulfate as alpha- or beta-structure inducers, respectively. Extrapolation of free energy data to zero concentration of structuring agent highlighted that the peptide prefers alpha-helical to beta-type organization, in spite of results from prediction algorithms. However, the free energy difference between the two forms, as obtained by a thermodynamic cycle, is subtle (roughly 5-8 kJ mol(-1) at any temperature from 280 K to 350 K), suggesting conformational ambivalence. This result supports the view that, in the prion protein, the structural behavior of the peptide is governed by the cellular microenvironment.     

HRMAS NMR observation of beta-sheet secondary structure on a water swollen solid support.

In this paper HRMAS NMR was used to investigate whether peptides on a peptidyl resin swollen in aqueous solution can adopt an intramolecular beta-sheet structure. A model peptide YQNPDGSQA, that was previously shown to adopt such a secondary structure in solution, (Blanco et al, J. Am. Chem. Soc., 1993) was grafted onto three different solid supports that swell in aqueous solution to examine the influence of the resin on the structure. Both parameters of resin loading and pH inside the swollen peptidyl resin proved to be important for the physicochemical behaviour of the peptide on the support.     

Contribution of secondary structure to the heat capacity and enthalpy distribution of the unfolded state in proteins.

We have recently shown that one can construct the enthalpy distribution for protein molecules from experimental knowledge of the temperature dependence of the heat capacity. For many proteins the enthalpy distribution evaluated at the midpoint of the denaturation transition (corresponding to the maximum in the heat capacity vs temperature curve) is broad and biphasic, indicating two different populations of molecules (native and unfolded) with distinctly different enthalpies. At temperatures above the denaturation point, the heat capacity for the unfolded state in many proteins is quite large and using the analysis just mentioned, we obtain a gaussian-like enthalpy distribution that is very broad. A large value of the heat capacity indicates that there are structural changes going on in the unfolded state above the transition temperature. In the present paper we investigate the origin of this large heat capacity by considering the presence of changing amounts of secondary structure (specifically, alpha-helix) in the unfolded state. For this purpose we use the empirical estimates of the Zimm-Bragg sigma and s factors for all of the native amino acids in water as determined by Scheraga and co-workers. Using myoglobin as an example, we calculate probability profiles and distribution functions for the total number of helix states in the specific-sequence molecule. Given the partition function for the specific-sequence molecule, we can then calculate a set of enthalpy moments for the molecule from which we obtain a good estimate of the enthalpy distribution in the unfolded state. This distribution turns out to be quite narrow when compared with the distribution obtained from the raw heat capacity data. We conclude that there must be other major structural changes (backbone and solvent) that are not accounted for by the inclusion of alpha-helix in the unfolded state.     

Disruption of the beta-sheet structure of a protected pentapeptide, related to the beta-amyloid sequence 17-21, induced by a single, helicogenic C(alpha)-tetrasubstituted alpha-amino acid.

Fifteen years ago it was shown that an alpha-aminoisobutyric acid (Aib) residue is significantly more effective than an L-Pro or a D-amino acid residue in inducing beta-sheet disruption in short model peptides. As this secondary structure element is known to play a crucial role in the neuropathology of Alzheimer's disease, it was decided to check the effect of Aib (and other selected, helix inducer, C(alpha)-tetrasubstituted alpha-amino acids) on the beta-sheet conformation adopted by a protected pentapeptide related to the sequence 17-21 of the beta-amyloid peptide. By use of FT-IR absorption and 1H NMR techniques it was found that the strong self-association characterizing the pentapeptide molecules in weakly polar organic solvents is completely abolished by replacing a single residue with Aib or one of its congeners.     

Knowledge-based protein secondary structure assignment

We have developed an automatic algorithm STRIDE for protein secondary structure assignment from atomic coordinates based on the combined use of hydrogen bond energy and statistically derived backbone torsional angle information. Parameters of the pattern recognition procedure were optimized using designations provided by the crystallographers as a standard-of-truth. Comparison to the currently most widely used technique DSSP by Kabsch and Sander (Biopolymers 22:2577-2637, 1983) shows that STRIDE and DSSP assign secondary structural states in 58 and 31% of 226 protein chains in our data sample, respectively, in greater agreement with the specific residue-by-residue definitions provided by the discoverers of the structures while in 11% of the chains, the assignments are the same. STRIDE delineates every 11th helix and every 32nd strand more in accord with published assignments. © 1995 Wiley-Liss, Inc.     

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

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

Application of ATR‐FTIR spectroscopy to the study of thermally induced changes in secondary structure of protein molecules in solid state

FTIR spectroscopy in combination with ATR sampling technique is the most accessible analytical technique to study secondary structure of proteins both in solid and aqueous solution. Although several studies have demonstrated the applications of ATR‐FTIR to study conformational changes of solid dried proteins due to dehydration, there are no reports that demonstrate the application of ATR‐FTIR in the study of thermally induced changes of secondary structure of biomolecules directly on the solid state. In this study, four biomolecules of pharmaceutical interest, lysozyme, myoglobine, chymotripsin and human growth hormone (hGH), were studied on the solid state before and after different thermal treatments in order to relate changes of secondary structure to partial or total thermal denaturation processes. The results obtained provide experimental evidence that protein thermal denaturation in the solid state can be detected by displacement of carbonyl bands which correspond to conformational transformations between α–helix to β‐sheet or intermolecular β‐sheet; the molecules studied undergo this transformation when exposed to a temperature close to their denaturation temperature which may become irreversible depending on the extent of the heating treatment. These findings demonstrate that ATR‐FTIR is an effective and time efficient technique that allows the monitoring of the protein thermal denaturation process of solid samples without further reconstitution or prior sample preparation. © 2015 Wiley Periodicals, Inc. Biopolymers 103: 574–584, 2015.     

CD studies on films of amyloid proteins and polypeptides: quantitative g-factor analysis indicates a common folding motif

Irrespective of the constituent protein, all amyloid fibrils show similar morphology in the electron microscope and x-ray diffraction patterns characteristic of a "cross-beta" structure, with extended beta-strands perpendicular to the fibril's long axis. Little is known about the amount or type of this structure. I have measured CD spectra of films formed from a number of amyloid proteins and polypeptides, and estimated their contents of extended secondary structure, by analysis of their g-factor spectra, the ratio of the CD and absorbance signals (P. McPhie, Analytical Biochemistry, 2001, Vol. 293, pp. 109-119). Amyloid films of Abeta-(1-40) peptide, beta-2-microglobulin, insulin, and three homopolypeptides show very intense CD spectra, compatible with the presence of a beta-helix-like structure, arranged in a common framework in the fibrils. The extent of this structure was estimated as 45-80% in the protein fibrils and 30-80% in the polypeptide fibrils.     

Analysis of interactive packing of secondary structural elements in alpha/beta units in proteins

An alpha-helix and a beta-strand are said to be interactively packed if at least one residue in each of the secondary structural elements loses 10% of its solvent accessible contact area on association with the other secondary structural element. An analysis of all such 5,975 nonidentical alpha/beta units in protein structures, defined at < or = 2.5 A resolution, shows that the interaxial distance between the alpha-helix and the beta-strand is linearly correlated with the residue-dependent function, log[(V/nda)/n-int], where V is the volume of amino acid residues in the packing interface, nda is the normalized difference in solvent accessible contact area of the residues in packed and unpacked secondary structural elements, and n-int is the number of residues in the packing interface. The beta-sheet unit (beta u), defined as a pair of adjacent parallel or antiparallel hydrogen-bonded beta-strands, packing with an alpha-helix shows a better correlation between the interaxial distance and log(V/nda) for the residues in the packing interface. This packing relationship is shown to be useful in the prediction of interaxial distances in alpha/beta units using the interacting residue information of equivalent alpha/beta units of homologous proteins. It is, therefore, of value in comparative modeling of protein structures.     

FTIR microspectroscopic analysis of the effects of certain drugs on oxidative stress and brain protein structure

This study reports the changes in lipids and proteins of different brain areas of nicotine, D+-amphetamine, and nicotine and D+-amphetamine treated rats by monitoring lipid peroxidation and protein beta-sheet formation using infrared microspectroscopy. Compared with the untreated brain samples, the peroxide level is relatively higher in the amphetamine-treated brain sections, both in the cortex and hippocampus area. However, this peroxide increase is attenuated when administering amphetamine plus nicotine. Analogous drug-dependent trends for protein beta-sheet content are observed, which suggests a connection between lipid oxidation involved in oxidative stress and beta-sheet protein structure generally present in neurodegenerative diseases. The above property of nicotine is of interest, in the sense that it might reduce the production of beta-amyloid proteins in Alzheimer's disease.     

Structures of intact glycoproteins from vibrational circular dichroism

Vibrational circular dichroism (VCD) spectra for the glycoproteins alpha1-acid glycoprotein (AGP) and bovine submaxillary mucin (BSM), have been measured in D2O solutions and for the films prepared from aqueous (H2O) buffer solutions in the 1800 to 900 cm(-1) region. The solution VCD results revealed that AGP has beta-sheet structure, along with a significant amount of alpha-helix as evidenced from a W pattern in the amide I region. The VCD of BSM solution suggested a polyproline II type structure, characterized by the appearance of strong negative couplet in the amide I region. The film VCD results on AGP and BSM suggested that the secondary structures of polypeptide fold in the film state are similar to those in the solution. The absence of any significant film VCD in the low frequency region (1200-900 cm(-1)), suggested that the dominant linkage for carbohydrate residues is likely to be a beta linkage. VCD spectroscopy gains importance in the secondary structural analysis of polypeptide fold in glycoproteins due to the absence of interfering VCD from the carbohydrate residues in the conformationally sensitive amide I region. Also, film VCD studies permit measurements in the low wavenumber region (1200-900 cm(-1)) that reveal the dominant type of linkage for carbohydrate residues. Such clear structural information is unlike that from ECD, where ECD bands of acylated amino sugar residues interfere with those of polypeptide backbone in the conformationally sensitive far-UV region.     

Crystal-state 3D-structural characterization of novel, Aib-based, turn and helical peptides.

The crystal-state conformations of the hexapeptide amide Pht-(Aib)(6)-NH-C(CH(3))(2)-O-OtBu (7), the hexapeptide Ac-L-aIle-(Aib)(5)-OtBu (6), the pentapeptide Z-(Aib)(3)-L-Glu(OtBu)-Aib-O-(CH(2))(2)-(1)Nap (5), the tetrapeptides Z-(Aib)(2)-L-His(N(tau)-Trt)-Aib-OMe (4 I) and Z-(Aib)(2)-L-Nva-Aib-OtBu (4 II), the tripeptide Pyr-(Aib)(3)-OtBu (3 I), the dipeptide amides Pyr-(Aib)(2)-(4)NH-TEMPO (3 II) and Piv-(Aib)(2)-NH-C(CH(3))(2)-O-OtBu (3 III), and the dipeptides Pht-Aib-betaAc(6)c-OtBu (2 I), Pht-Aib-NH-C(CH(3))(2)-O-OtBu (2 II) and Boc-gGly-mAib-OH (2 III) have been determined by X-ray diffraction analyses. All peptides investigated are characterized by one or more turn/helix forming Aib residues. Except the three short dipeptides, all are folded into C==O...H--N intramolecularly H-bonded 3(10)-helices, or into various types of beta-turns. In the structure of 6, two independent molecules of opposite screw sense were observed in the asymmetric unit, generating diastereomeric 3(10)-helices.     

Thermal folding and mechanical unfolding pathways of protein secondary structures

Mechanical stretching of secondary structures is studied through molecular dynamics simulations of a Go-like model. Force versus displacement curves are studied as a function of the stiffness and velocity of the pulling device. The succession of stretching events, as measured by the order in which contacts are ruptured, is compared to the sequencing of events during thermal folding and unfolding. Opposite cross-correlations are found for an alpha-helix and a beta-hairpin structure. In a tandem of two alpha-helices, the two constituent helices unravel nearly simultaneously. A simple condition for simultaneous versus sequential unraveling of repeat units is presented.