A probabilistic model for secondary structure prediction from protein chemical shifts

Protein chemical shifts encode detailed structural information that is difficult and computationally costly to describe at a fundamental level. Statistical and machine learning approaches have been used to infer correlations between chemical shifts and secondary structure from experimental chemical shifts. These methods range from simple statistics such as the chemical shift index to complex methods using neural networks. Notwithstanding their higher accuracy, more complex approaches tend to obscure the relationship between secondary structure and chemical shift and often involve many parameters that need to be trained. We present hidden Markov models (HMMs) with Gaussian emission probabilities to model the dependence between protein chemical shifts and secondary structure. The continuous emission probabilities are modeled as conditional probabilities for a given amino acid and secondary structure type. Using these distributions as outputs of first‐ and second‐order HMMs, we achieve a prediction accuracy of 82.3%, which is competitive with existing methods for predicting secondary structure from protein chemical shifts. Incorporation of sequence‐based secondary structure prediction into our HMM improves the prediction accuracy to 84.0%. Our findings suggest that an HMM with correlated Gaussian distributions conditioned on the secondary structure provides an adequate generative model of chemical shifts. Proteins 2013; © 2012 Wiley Periodicals, Inc.     

The effect of ring currents on carbon chemical shifts in cytochromes

Calculations suggest that some carbon chemical shifts in proteinsshould have large ring current shifts (>1 ppm). We present¹³C, ¹⁵N and ¹H assignments forcytochrome c₂ from Rhodospirillum rubrum, compare these withshifts for other cytochromes c, and show that the calculated ring currentshifts are similar to experimentally observed shifts, but that there remainsubstantial conformation-dependent shifts of side-chain carbons. Ringcurrent shifts as large as 6 ppm are observed. We show that the ring currenteffects do not seriously affect the Chemical Shift Index method fordelineating secondary structure, but may have an impact on more precisemethods for generating structural constraints.     

Predicting the redox state and secondary structure of cysteine residues in proteins using NMR chemical shifts

Different programs and methods were employed to superimpose protein structures, using members of four very different protein families as test subjects, and the results of these efforts were compared. Algorithms based on human identification of key amino acid residues on which to base the superpositions were nearly always more successful than programs that used automated techniques to identify key residues. Among those programs automatically identifying key residues, MASS could not superimpose all members of some families, but was very efficient with other families. MODELLER, MultiProt, and STAMP had varying levels of success. A genetic algorithm program written for this project did not improve superpositions when results from neighbor-joining and pseudostar algorithms were used as its starting cases, but it always improved superpositions obained by MODELLER and STAMP. A program entitled PyMSS is presented that includes three superposition algorithms featuring human interaction.     

Diversity of functions of proteins with internal symmetry in spatial arrangement of secondary structural elements.

We carry out a systematic analysis of the correlation between similarity of protein three-dimensional structures and their evolutionary relationships. The structural similarity is quantitatively identified by an all-against-all comparison of the spatial arrangement of secondary structural elements in nonredundant 967 representative proteins, and the evolutionary relationship is judged according to the definition of superfamily in the SCOP database. We find the following symmetry rule: a protein pair that has similar folds but belong to different superfamilies has (with a very rare exception) certain internal symmetry in its common similar folds. Possible reasons behind the symmetry rule are discussed.     

Prediction algorithm for amino acid types with their secondary structure in proteins (PLATON) using chemical shifts

The algorithm PLATON is able to assign sets of chemical shifts derived from a single residue to amino acid types with its secondary structure (amino acid species). A subsequent ranking procedure using optionally two different penalty functions yields predictions for possible amino acid species for the given set of chemical shifts. This was demonstrated in the case of the -spectrin SH3 domain and applied to 9 further protein data sets taken from the BioMagRes database. A database consisting of reference chemical shift patterns (reference CSPs) was generated from assigned chemical shifts of proteins with known 3D-structure. This reference CSP database is used in our approach for extracting distributions of amino acid types with their most likely secondary structure elements (namely -helix, -sheet, and coil) for single amino acids by comparison with query CSPs. Results obtained for the 10 investigated proteins indicates that the percentage of correct amino acid species in the first three positions in the ranking list, ranges from 71.4% to 93.2% for the more favorable penalty function. Where only the top result of the ranking list for these 10 proteins is considered, 36.5% to 83.1% of the amino acid species are correctly predicted. The main advantage of our approach, over other methods that rely on average chemical shift values is the ability to increase database content by incorporating newly derived CSPs, and therefore to improve PLATON's performance over time.     

The 13C chemical shifts of amino acids in aqueous solution containing organic solvents: Application to the secondary structure characterization of peptides in aqueous trifluoroethanol solution

Summary The ¹³C chemical shifts for all of the protonated carbons of the 20 common amino acid residues in the protected linear pentapeptide Gly-Gly-X-Gly-Gly have been obtained in water at low pH as well as in aqueous solution containing 10, 20 and 30% acetonitrile or trifluoroethanol. Dioxane was used as an internal reference and its carbon chemical shift value was found to be 66.6 ppm relative to external TMS in water. Comparison of the different referencing methods for ¹³C chemical shifts in organic cosolvent mixtures showed that an external standard (either TMS or TSP capillary) was the most appropriate. In the present study, external TSP was chosen to define the 0 ppm of the ¹³C chemical shift scale. When the difference in referencing the dioxane carbon resonance is taken into account, the carbon chemical shift values of the amino acids in aqueous solution are similar to those previously reported (Richarz and Wüthrich (1978) Biopolymers, 17, 2133–2141; Howarth and Lilley (1979) Prog. NMR Spectrosc., 12, 1–40). The pentapeptides studied were assumed to be in a random coil conformation and the measured ¹³C chemical shifts were used as reference values to correlate carbon chemical shifts with the secondary structure of two well-characterized peptides, bombesin and the 1–29 amino acid fragment of Nle²⁷ human growth hormone-releasing factor. In both cases, the C chemical shifts exhibited a characteristic positive deviation from the random coil values, which indicates the presence of -helices.     

Cooperative effects in hydrogen-bonding of protein secondary structure elements: a systematic analysis of crystal data using Secbase

A systematic analysis of the hydrogen-bonding geometry in helices and beta sheets has been performed. The distances and angles between the backbone carbonyl O and amide N atoms were correlated considering more than 1500 protein chains in crystal structures determined to a resolution better than 1.5 A. They reveal statistically significant trends in the H-bond geometry across the different secondary structural elements. The analysis has been performed using Secbase, a modular extension of Relibase (Receptor Ligand Database) which integrates information about secondary structural elements assigned to individual protein structures with the various search facilities implemented into Relibase. A comparison of the mean hydrogen-bond distances in alpha helices and 3(10) helices of increasing length shows opposing trends. Whereas in alpha helices the mean H-bond distance shrinks with increasing helix length and turn number, the corresponding mean dimension in 3(10) helices expands in a comparable series. Comparing similarly the hydrogen-bond lengths in beta sheets there is no difference to be found between the mean H-bond length in antiparallel and parallel beta sheets along the strand direction. In contrast, an interesting systematic trend appears to be given for the hydrogen bonds perpendicular to the strands bridging across an extended sheet. With increasing number of accumulated strands, which results in a growing number of back-to-back piling hydrogen bonds across the strands, a slight decrease of the mean H-bond distance is apparent in parallel beta sheets whereas such trends are obviously not given in antiparallel beta sheets. This observation suggests that cooperative effects mutually polarizing spatially well-aligned hydrogen bonds are present either in alpha helices and parallel beta sheets whereas such influences seem to be lacking in 3(10) helices and antiparallel beta sheets.