Empirical correlation between protein backbone 15N and 13C secondary chemical shifts and its application to nitrogen chemical shift re-referencing

Predicted and Experimental NMR Chemical Shifts at Variable Temperatures: The Effect of Protein Conformational Dynamics

NMR chemical shifts provide a sensitive probe of protein structure and dynamics but remain challenging to predict and interpret. We examine the effect of protein conformational distributions on 15N chemical shifts for dihydrofolate reductase (DHFR), comparing QM/MM predicted shifts with experimental shifts in solution as well as frozen distributions. Representative snapshots from MD trajectories exhibit variation in predicted 15N chemical shifts of up to 25 ppm. The average over the fluctuations is in significantly better agreement with room temperature solution experimental values than the prediction for any single optimal conformations. Meanwhile, solid-state NMR (SSNMR) measurements of frozen solutions at 105 K exhibit broad lines whose widths agree well with the widths of distributions of predicted shifts for samples from the trajectory. The backbone torsion angle ψi-1 varies over 60° on the picosecond time scale, compensated by φi. These fluctuations can explain much of the shift variation.

NMR-Based Characterization of the Interaction between Yeast Oxa1-CTD and Ribosomes

In mitochondria, the major subunits of oxidative phosphorylation complexes are translated by the mitochondrial ribosome (mito-ribosome). The correct insertion and assembly of these subunits into the inner mitochondrial membrane (IMM) are facilitated by mitochondrial oxidase assembly protein 1 (Oxa1) during the translation process. This co-translational insertion process involves an association between the mito-ribosome and the C-terminus of Oxa1 (Oxa1-CTD) Nuclear magnetic resonance (NMR) methods were mainly used to investigate the structural characterization of yeast Oxa1-CTD and its mode of interaction with the E. coli 70S ribosome. Oxa1-CTD forms a transient α-helical structure within the residues P342-Q385, which were reported to form an α-helix when combining with the ribosome. Two conserved contact sites that could interact with the ribosome were further identified. The first site was located on the very end of the N-terminus (V321-I327), and the second one encompassed a stretch of amino acid residues I348-Q370. Based on our discoveries and previous reports, a model has been proposed in which Oxa1-CTD interacts with ribosomes, accompanied by transient-to-stable transitions at the second contact site. These observations may enhance our understanding of the potential role of Oxa1-CTD in facilitating the assembly of oxidative phosphorylation complexes and provide insight into the structural characteristics of Oxa1-CTD.

The influence of random-coil chemical shifts on the assessment of structural propensities in folded proteins and IDPs

In studying secondary structural propensities of proteins by nuclear magnetic resonance (NMR) spectroscopy, secondary chemical shifts (SCSs) serve as the primary atomic scale observables. For SCS calculation, the selection of an appropriate random coil chemical shift (RCCS) dataset is a crucial step, especially when investigating intrinsically disordered proteins (IDPs). The scientific literature is abundant in such datasets, however, the effect of choosing one over all the others in a concrete application has not yet been studied thoroughly and systematically. Hereby, we review the available RCCS prediction methods and to compare them, we conduct statistical inference by means of the nonparametric sum of ranking differences and comparison of ranks to random numbers (SRD-CRRN) method. We try to find the RCCS predictors best representing the general consensus regarding secondary structural propensities. The existence and the magnitude of resulting differences on secondary structure determination under varying sample conditions (temperature, pH) are demonstrated and discussed for globular proteins and especially IDPs.

Predicted and Experimental NMR Chemical Shifts at Variable Temperatures: The Effect of Protein Conformational Dynamics

NMR chemical shifts provide a sensitive probe of protein structure and dynamics. Prediction of shifts, and therefore interpretation of shifts, particularly for the frequently measured amidic 15 N sites, remains a tall challenge. We demonstrate that protein 15 N chemical shift prediction from QM/MM predictions can be improved if conformational variation is included via MD sampling, focusing on the antibiotic target, E. coli Dihydrofolate reductase (DHFR). Variations of up to 25 ppm in predicted 15 N chemical shifts are observed over the trajectory. For solution shifts the average of fluctuations on the low picosecond timescale results in a superior prediction to a single optimal conformation. For low temperature solid state measurements, the histogram of predicted shifts for locally minimized snapshots with specific solvent arrangements sampled from the trajectory explains the heterogeneous linewidths; in other words, the conformations and associated solvent are 'frozen out' at low temperatures and result in inhomogeneously broadened NMR peaks. We identified conformational degrees of freedom that contribute to chemical shift variation. Backbone torsion angles show high amplitude fluctuations during the trajectory on the low picosecond timescale. For a number of residues, including I60, ψ varies by up to 60º within a conformational basin during the MD simulations, despite the fact that I60 (and other sites studied) are in a secondary structure element and remain well folded during the trajectory. Fluctuations in ψ appear to be compensated by other degrees of freedom in the protein, including φ of the succeeding residue, resulting in "rocking" of the amide plane with changes in hydrogen bonding interactions. Good agreement for both room temperature and low temperature NMR spectra provides strong support for the specific approach to conformational averaging of computed chemical shifts.

The residual structure of acid-denatured β2 -microglobulin is relevant to an ordered fibril morphology

β₂ -Microglobulin (β2m) forms amyloid fibrils in vitro under acidic conditions. Under these conditions, the residual structure of acid-denatured β2m is relevant to seeding and fibril extension processes. Disulfide (SS) bond-oxidized β2m has been shown to form rigid, ordered fibrils, whereas SS bond-reduced β2m forms curvy, less-ordered fibrils. These findings suggest that the presence of an SS bond affects the residual structure of the monomer, which subsequently influences the fibril morphology. To clarify this process, we herein performed NMR experiments. The results obtained revealed that oxidized β2m contained a residual structure throughout the molecule, including the N- and C-termini, whereas the residual structure of the reduced form was localized and other regions had a random coil structure. The range of the residual structure in the oxidized form was wider than that of the fibril core. These results indicate that acid-denatured β2m has variable conformations. Most conformations in the ensemble cannot participate in fibril formation because their core residues are hidden by residual structures. However, when hydrophobic residues are exposed, polypeptides competently form an ordered fibril. This conformational selection phase may be needed for the ordered assembly of amyloid fibrils.

NMR Spectroscopy in the Conformational Analysis of Peptides: An Overview

BACKGROUND

NMR spectroscopy is one of the most powerful tools to study the structure and interaction properties of peptides and proteins from a dynamic perspective. Knowing the bioactive conformations of peptides is crucial in the drug discovery field to design more efficient analogue ligands and inhibitors of protein-protein interactions targeting therapeutically relevant systems.

OBJECTIVE

This review provides a toolkit to investigate peptide conformational properties by NMR.

METHODS

Articles cited herein, related to NMR studies of peptides and proteins were mainly searched through PubMed and the web. More recent and old books on NMR spectroscopy written by eminent scientists in the field were consulted as well.

RESULTS

The review is mainly focused on NMR tools to gain the 3D structure of small unlabeled peptides. It is more application-oriented as it is beyond its goal to deliver a profound theoretical background. However, the basic principles of 2D homonuclear and heteronuclear experiments are briefly described. Protocols to obtain isotopically labeled peptides and principal triple resonance experiments needed to study them, are discussed as well.

CONCLUSION

NMR is a leading technique in the study of conformational preferences of small flexible peptides whose structure can be often only described by an ensemble of conformations. Although NMR studies of peptides can be easily and fast performed by canonical protocols established a few decades ago, more recently we have assisted to tremendous improvements of NMR spectroscopy to investigate instead large systems and overcome its molecular weight limit.

A transient α-helix in the N-terminal RNA recognition motif of polypyrimidine tract binding protein senses RNA secondary structure

The polypyrimidine tract binding protein (PTB) is a multi-domain protein involved in alternative splicing, mRNA localization, stabilization, polyadenylation and translation initiation from internal ribosome entry sites (IRES). In this latter process, PTB promotes viral translation by interacting extensively with complex structured regions in the 5′-untranslated regions of viral RNAs at pyrimidine-rich targets located in single strand and hairpin regions. To better understand how PTB recognizes structured elements in RNA targets, we solved the solution structure of the N-terminal RNA recognition motif (RRM) in complex with an RNA hairpin embedding the loop sequence UCUUU, which is frequently found in IRESs of the picornovirus family. Surprisingly, a new three-turn α3 helix C-terminal to the RRM, folds upon binding the RNA hairpin. Although α3 does not mediate any contacts to the RNA, it acts as a sensor of RNA secondary structure, suggesting a role for RRM1 in detecting pyrimidine tracts in the context of structured RNA. Moreover, the degree of helix formation depends on the RNA loop sequence. Finally, we show that the α3 helix region, which is highly conserved in vertebrates, is crucial for PTB function in enhancing Encephalomyocarditis virus IRES activity.

Temperature dependence of NMR chemical shifts: Tracking and statistical analysis

Isotropic chemical shifts measured by solution nuclear magnetic resonance (NMR) spectroscopy offer extensive insights into protein structure and dynamics. Temperature dependences add a valuable dimension; notably, the temperature dependences of amide proton chemical shifts are valuable probes of hydrogen bonding, temperature‐dependent loss of structure, and exchange between distinct protein conformations. Accordingly, their uses include structural analysis of both folded and disordered proteins, and determination of the effects of mutations, binding, or solution conditions on protein energetics. Fundamentally, these temperature dependences result from changes in the local magnetic environments of nuclei, but correlations with global thermodynamic parameters measured via calorimetric methods have been observed. Although the temperature dependences of amide proton and nitrogen chemical shifts are often well approximated by a linear model, deviations from linearity are also observed and may be interpreted as evidence of fast exchange between distinct conformational states. Here, we describe computational methods, accessible via the Shift‐T web server, including an automated tracking algorithm that propagates initial (single temperature) ¹H—¹⁵N cross peak assignments to spectra collected over a range of temperatures. Amide proton and nitrogen temperature coefficients (slopes determined by fitting chemical shift vs. temperature data to a linear model) are subsequently calculated. Also included are methods for the detection of systematic, statistically significant deviation from linearity (curvature) in the temperature dependences of amide proton chemical shifts. The use and utility of these methods are illustrated by example, and the Shift‐T web server is freely available at http://meieringlab.uwaterloo.ca/shiftt.

NMRFx Processor: a cross-platform NMR data processing program

NMRFx Processor is a new program for the processing of NMR data. Written in the Java programming language, NMRFx Processor is a cross-platform application and runs on Linux, Mac OS X and Windows operating systems. The application can be run in both a graphical user interface (GUI) mode and from the command line. Processing scripts are written in the Python programming language and executed so that the low-level Java commands are automatically run in parallel on computers with multiple cores or CPUs. Processing scripts can be generated automatically from the parameters of NMR experiments or interactively constructed in the GUI. A wide variety of processing operations are provided, including methods for processing of non-uniformly sampled datasets using iterative soft thresholding. The interactive GUI also enables the use of the program as an educational tool for teaching basic and advanced techniques in NMR data analysis.

Conformationally selective multidimensional chemical shift ranges in proteins from a PACSY database purged using intrinsic quality criteria

We have determined refined multidimensional chemical shift ranges for intra-residue correlations ((13)C-(13)C, (15)N-(13)C, etc.) in proteins, which can be used to gain type-assignment and/or secondary-structure information from experimental NMR spectra. The chemical-shift ranges are the result of a statistical analysis of the PACSY database of >3000 proteins with 3D structures (1,200,207 (13)C chemical shifts and >3 million chemical shifts in total); these data were originally derived from the Biological Magnetic Resonance Data Bank. Using relatively simple non-parametric statistics to find peak maxima in the distributions of helix, sheet, coil and turn chemical shifts, and without the use of limited "hand-picked" data sets, we show that ~94% of the (13)C NMR data and almost all (15)N data are quite accurately referenced and assigned, with smaller standard deviations (0.2 and 0.8 ppm, respectively) than recognized previously. On the other hand, approximately 6% of the (13)C chemical shift data in the PACSY database are shown to be clearly misreferenced, mostly by ca. -2.4 ppm. The removal of the misreferenced data and other outliers by this purging by intrinsic quality criteria (PIQC) allows for reliable identification of secondary maxima in the two-dimensional chemical-shift distributions already pre-separated by secondary structure. We demonstrate that some of these correspond to specific regions in the Ramachandran plot, including left-handed helix dihedral angles, reflect unusual hydrogen bonding, or are due to the influence of a following proline residue. With appropriate smoothing, significantly more tightly defined chemical shift ranges are obtained for each amino acid type in the different secondary structures. These chemical shift ranges, which may be defined at any statistical threshold, can be used for amino-acid type assignment and secondary-structure analysis of chemical shifts from intra-residue cross peaks by inspection or by using a provided command-line Python script (PLUQin), which should be useful in protein structure determination. The refined chemical shift distributions are utilized in a simple quality test (SQAT) that should be applied to new protein NMR data before deposition in a databank, and they could benefit many other chemical-shift based tools.

NMR-based modeling and refinement of protein 3D structures

NMR is a well-established method to characterize the structure and dynamics of biomolecules in solution. High-quality structures can now be produced thanks to both experimental advances and computational developments that incorporate new NMR parameters and improved protocols and force fields in the structure calculation and refinement process. In this chapter, we give a short overview of the various types of NMR data that can provide structural information, and then focus on the structure calculation methodology itself. We discuss and illustrate with tutorial examples "classical" structure calculation, refinement, and structure validation approaches.

An overview of tools for the validation of protein NMR structures

Biomolecular structures at atomic resolution present a valuable resource for the understanding of biology. NMR spectroscopy accounts for 11% of all structures in the PDB repository. In response to serious problems with the accuracy of some of the NMR-derived structures and in order to facilitate proper analysis of the experimental models, a number of program suites are available. We discuss nine of these tools in this review: PROCHECK-NMR, PSVS, GLM-RMSD, CING, Molprobity, Vivaldi, ResProx, NMR constraints analyzer and QMEAN. We evaluate these programs for their ability to assess the structural quality, restraints and their violations, chemical shifts, peaks and the handling of multi-model NMR ensembles. We document both the input required by the programs and output they generate. To discuss their relative merits we have applied the tools to two representative examples from the PDB: a small, globular monomeric protein (Staphylococcal nuclease from S. aureus, PDB entry 2kq3) and a small, symmetric homodimeric protein (a region of human myosin-X, PDB entry 2lw9).

Conformational transitions in the Ramachandran space of amino acids using the dynamic rotational isomeric state (DRIS) model

The dynamic rotational isomeric state model is applied to predict the internal dynamics of the 20 amino acids. Transition rates between rotational isomeric states are calculated from molecular dynamics simulations of Gly-Gly-X-Gly-Gly peptides where X represents one of the 20 amino acids. Predicted relaxation times are in good agreement with fluorescence quenching rate measurements.

Recommendations of the wwPDB NMR Validation Task Force

As methods for analysis of biomolecular structure and dynamics using nuclear magnetic resonance spectroscopy (NMR) continue to advance, the resulting 3D structures, chemical shifts, and other NMR data are broadly impacting biology, chemistry, and medicine. Structure model assessment is a critical area of NMR methods development, and is an essential component of the process of making these structures accessible and useful to the wider scientific community. For these reasons, the Worldwide Protein Data Bank (wwPDB) has convened an NMR Validation Task Force (NMR-VTF) to work with wwPDB partners in developing metrics and policies for biomolecular NMR data harvesting, structure representation, and structure quality assessment. This paper summarizes the recommendations of the NMR-VTF, and lays the groundwork for future work in developing standards and metrics for biomolecular NMR structure quality assessment.

Vivaldi: visualization and validation of biomacromolecular NMR structures from the PDB

We describe Vivaldi (VIsualization and VALidation DIsplay; http://pdbe.org/vivaldi), a web-based service for the analysis, visualization, and validation of NMR structures in the Protein Data Bank (PDB). Vivaldi provides access to model coordinates and several types of experimental NMR data using interactive visualization tools, augmented with structural annotations and model-validation information. The service presents information about the modeled NMR ensemble, validation of experimental chemical shifts, residual dipolar couplings, distance and dihedral angle constraints, as well as validation scores based on empirical knowledge and databases. Vivaldi was designed for both expert NMR spectroscopists and casual non-expert users who wish to obtain a better grasp of the information content and quality of NMR structures in the public archive. © Proteins 2013. © 2012 Wiley Periodicals, Inc.

HASH: a program to accurately predict protein Hα shifts from neighboring backbone shifts

Chemical shifts provide not only peak identities for analyzing nuclear magnetic resonance (NMR) data, but also an important source of conformational information for studying protein structures. Current structural studies requiring H(α) chemical shifts suffer from the following limitations. (1) For large proteins, the H(α) chemical shifts can be difficult to assign using conventional NMR triple-resonance experiments, mainly due to the fast transverse relaxation rate of C(α) that restricts the signal sensitivity. (2) Previous chemical shift prediction approaches either require homologous models with high sequence similarity or rely heavily on accurate backbone and side-chain structural coordinates. When neither sequence homologues nor structural coordinates are available, we must resort to other information to predict H(α) chemical shifts. Predicting accurate H(α) chemical shifts using other obtainable information, such as the chemical shifts of nearby backbone atoms (i.e., adjacent atoms in the sequence), can remedy the above dilemmas, and hence advance NMR-based structural studies of proteins. By specifically exploiting the dependencies on chemical shifts of nearby backbone atoms, we propose a novel machine learning algorithm, called HASH, to predict H(α) chemical shifts. HASH combines a new fragment-based chemical shift search approach with a non-parametric regression model, called the generalized additive model, to effectively solve the prediction problem. We demonstrate that the chemical shifts of nearby backbone atoms provide a reliable source of information for predicting accurate H(α) chemical shifts. Our testing results on different possible combinations of input data indicate that HASH has a wide rage of potential NMR applications in structural and biological studies of proteins.

An automated system designed for large scale NMR data deposition and annotation: application to over 600 assigned chemical shift data entries to the BioMagResBank from the Riken Structural Genomics/Proteomics Initiative internal database

Biomolecular NMR chemical shift data are key information for the functional analysis of biomolecules and the development of new techniques for NMR studies utilizing chemical shift statistical information. Structural genomics projects are major contributors to the accumulation of protein chemical shift information. The management of the large quantities of NMR data generated by each project in a local database and the transfer of the data to the public databases are still formidable tasks because of the complicated nature of NMR data. Here we report an automated and efficient system developed for the deposition and annotation of a large number of data sets including (1)H, (13)C and (15)N resonance assignments used for the structure determination of proteins. We have demonstrated the feasibility of our system by applying it to over 600 entries from the internal database generated by the RIKEN Structural Genomics/Proteomics Initiative (RSGI) to the public database, BioMagResBank (BMRB). We have assessed the quality of the deposited chemical shifts by comparing them with those predicted from the PDB coordinate entry for the corresponding protein. The same comparison for other matched BMRB/PDB entries deposited from 2001-2011 has been carried out and the results suggest that the RSGI entries greatly improved the quality of the BMRB database. Since the entries include chemical shifts acquired under strikingly similar experimental conditions, these NMR data can be expected to be a promising resource to improve current technologies as well as to develop new NMR methods for protein studies.

Combining NMR ensembles and molecular dynamics simulations provides more realistic models of protein structures in solution and leads to better chemical shift prediction

While chemical shifts are invaluable for obtaining structural information from proteins, they also offer one of the rare ways to obtain information about protein dynamics. A necessary tool in transforming chemical shifts into structural and dynamic information is chemical shift prediction. In our previous work we developed a method for 4D prediction of protein (1)H chemical shifts in which molecular motions, the 4th dimension, were modeled using molecular dynamics (MD) simulations. Although the approach clearly improved the prediction, the X-ray structures and single NMR conformers used in the model cannot be considered fully realistic models of protein in solution. In this work, NMR ensembles (NMRE) were used to expand the conformational space of proteins (e.g. side chains, flexible loops, termini), followed by MD simulations for each conformer to map the local fluctuations. Compared with the non-dynamic model, the NMRE+MD model gave 6-17% lower root-mean-square (RMS) errors for different backbone nuclei. The improved prediction indicates that NMR ensembles with MD simulations can be used to obtain a more realistic picture of protein structures in solutions and moreover underlines the importance of short and long time-scale dynamics for the prediction. The RMS errors of the NMRE+MD model were 0.24, 0.43, 0.98, 1.03, 1.16 and 2.39 ppm for (1)Hα, (1)HN, (13)Cα, (13)Cβ, (13)CO and backbone (15)N chemical shifts, respectively. The model is implemented in the prediction program 4DSPOT, available at http://www.uef.fi/4dspot.

Orphan macrodomain protein (human C6orf130) is an O-acyl-ADP-ribose deacylase: solution structure and catalytic properties

Post-translational modification of proteins/histones by lysine acylation has profound effects on the physiological function of modified proteins. Deacylation by NAD(+)-dependent sirtuin reactions yields as a product O-acyl-ADP-ribose, which has been implicated as a signaling molecule in modulating cellular processes. Macrodomain-containing proteins are reported to bind NAD(+)-derived metabolites. Here, we describe the structure and function of an orphan macrodomain protein, human C6orf130. This unique 17-kDa protein is a stand-alone macrodomain protein that occupies a distinct branch in the phylogenic tree. We demonstrate that C6orf130 catalyzes the efficient deacylation of O-acetyl-ADP-ribose, O-propionyl-ADP-ribose, and O-butyryl-ADP-ribose to produce ADP-ribose (ADPr) and acetate, propionate, and butyrate, respectively. Using NMR spectroscopy, we solved the structure of C6orf130 in the presence and absence of ADPr. The structures showed a canonical fold with a deep ligand (ADPr)-binding cleft. Structural comparisons of apo-C6orf130 and the ADPr-C6orf130 complex revealed fluctuations of the β(5)-α(4) loop that covers the bound ADPr, suggesting that the β(5)-α(4) loop functions as a gate to sequester substrate and offer flexibility to accommodate alternative substrates. The ADPr-C6orf130 complex identified amino acid residues involved in substrate binding and suggested residues that function in catalysis. Site-specific mutagenesis and steady-state kinetic analyses revealed two critical catalytic residues, Ser-35 and Asp-125. We propose a catalytic mechanism for deacylation of O-acyl-ADP-ribose by C6orf130 and discuss the biological implications in the context of reversible protein acylation at lysine residues.

Quantitative use of chemical shifts for the modeling of protein complexes

Despite recent advances in the modeling of protein-protein complexes by docking, additional information is often required to identify the best solutions. For this purpose, NMR data deliver valuable restraints that can be used in the sampling and/or the scoring stage, like in the data-driven docking approach HADDOCK that can make use of NMR chemical shift perturbation (CSP) data to define the binding site on each protein and drive the docking. We show here that a quantitative use of chemical shifts (CS) in the scoring stage can help to resolve ambiguities. A quantitative CS-RMSD score based on (1) H(α) ,(13) C(α) , and (15) N chemical shifts ranks the best solutions always at the top, as demonstrated on a small benchmark of four complexes. It is implemented in a new docking protocol, CS-HADDOCK, which combines CSP data as ambiguous interaction restraints in the sampling stage with the CS-RMSD score in the final scoring stage. This combination of qualitative and quantitative use of chemical shifts increases the reliability of data-driven docking for the structure determination of complexes from limited NMR data.

Overview on the use of NMR to examine protein structure

Any protein structure determination process contains several steps, starting from obtaining a suitable sample, then moving on to acquiring data and spectral assignment, and lastly to the final steps of structure determination and validation. This unit describes all of these steps, starting with the basic physical principles behind NMR and some of the most commonly measured and observed phenomena such as chemical shift, scalar and residual coupling, and the nuclear Overhauser effect. Then, in somewhat more detail, the process of spectral assignment and structure elucidation is explained. Furthermore, the use of NMR to study protein-ligand interaction, protein dynamics, or protein folding is described.

Time-resolved NMR: extracting the topology of complex enzyme networks

The use of nondestructive NMR spectroscopy for enzymatic studies offers unique opportunities to identify nearly all enzymatic byproducts and detect unstable short-lived products or intermediates at the molecular level; however, numerous challenges must be overcome before it can become a widely used tool. The biosynthesis of acetyl-coenzyme A (acetyl-CoA) by acetyl-CoA synthetase is used here as a case study for the development of an analytical NMR-based time-course assay platform. We describe an algorithm to deconvolve superimposed spectra into spectra for individual molecules, and further develop a model to simulate the acetyl-CoA synthetase enzyme reaction network using the data derived from time-course NMR. Simulation shows indirectly that synthesis of acetyl-CoA is mediated via an enzyme-bound intermediate (possibly acetyl-AMP) and is accompanied by a nonproductive loss from an intermediate. The ability to predict enzyme function based on partial knowledge of the enzymatic pathway topology is also discussed.

A probabilistic approach for validating protein NMR chemical shift assignments

It has been estimated that more than 20% of the proteins in the BMRB are improperly referenced and that about 1% of all chemical shift assignments are mis-assigned. These statistics also reflect the likelihood that any newly assigned protein will have shift assignment or shift referencing errors. The relatively high frequency of these errors continues to be a concern for the biomolecular NMR community. While several programs do exist to detect and/or correct chemical shift mis-referencing or chemical shift mis-assignments, most can only do one, or the other. The one program (SHIFTCOR) that is capable of handling both chemical shift mis-referencing and mis-assignments, requires the 3D structure coordinates of the target protein. Given that chemical shift mis-assignments and chemical shift re-referencing issues should ideally be addressed prior to 3D structure determination, there is a clear need to develop a structure-independent approach. Here, we present a new structure-independent protocol, which is based on using residue-specific and secondary structure-specific chemical shift distributions calculated over small (3-6 residue) fragments to identify mis-assigned resonances. The method is also able to identify and re-reference mis-referenced chemical shift assignments. Comparisons against existing re-referencing or mis-assignment detection programs show that the method is as good or superior to existing approaches. The protocol described here has been implemented into a freely available Java program called "Probabilistic Approach for protein Nmr Assignment Validation (PANAV)" and as a web server ( http://redpoll.pharmacy.ualberta.ca/PANAV ) which can be used to validate and/or correct as well as re-reference assigned protein chemical shifts.

4D prediction of protein (1)H chemical shifts

A 4D approach for protein (1)H chemical shift prediction was explored. The 4th dimension is the molecular flexibility, mapped using molecular dynamics simulations. The chemical shifts were predicted with a principal component model based on atom coordinates from a database of 40 protein structures. When compared to the corresponding non-dynamic (3D) model, the 4th dimension improved prediction by 6-7%. The prediction method achieved RMS errors of 0.29 and 0.50 ppm for Halpha and HN shifts, respectively. However, for individual proteins the RMS errors were 0.17-0.34 and 0.34-0.65 ppm for the Halpha and HN shifts, respectively. X-ray structures gave better predictions than the corresponding NMR structures, indicating that chemical shifts contain invaluable information about local structures. The (1)H chemical shift prediction tool 4DSPOT is available from http://www.uku.fi/kemia/4dspot .

CheckShift improved: fast chemical shift reference correction with high accuracy

The construction of a consistent protein chemical shift database is an important step toward making more extensive use of this data in structural studies. Unfortunately, progress in this direction has been hampered by the quality of the available data, particularly with respect to chemical shift referencing, which is often either inaccurate or inconsistently annotated. Preprocessing of the data is therefore required to detect and correct referencing errors. In an earlier study we developed CheckShift, a program for performing this task automatically. Now we spent substantial effort in improving the running time of the CheckShift algorithm, which resulted in an running time decrease of 90%, thereby achieving equivalent quality to the former version of CheckShift. The reason for the running time decrease is twofold. Firstly we improved the search for the optimal re-referencing offset considerably. Secondly, as CheckShift is based on a secondary structure prediction from the amino acid sequence (formally PsiPred was used), we evaluated a wide range of available secondary structure prediction programs focusing on the special needs of the CheckShift algorithm. The results of this evaluation prove empirically that we can use faster secondary structure prediction programs than PsiPred without sacrificing CheckShift's accuracy. Very recently Wang and Markley (2009) gave a small list of extreme outliers of the former version of the CheckShift web-server. Those were due to the empirical reduction of the search space implemented in the old version. The new version of CheckShift now gives very similar results to RefDB and LACS for all outliers mentioned in Table 1 of Wang and Markley (2009).