Secondary structural analysis of proteins based on 13C chemical shift assignments in unresolved solid-state NMR spectra enhanced by fragmented structure database

Magic-angle-spinning solid-state ¹³C NMR spectroscopy is useful for structural analysis of non-crystalline proteins. However, the signal assignments and structural analysis are often hampered by the signal overlaps primarily due to minor structural heterogeneities, especially for uniformly-¹³C,¹⁵N labeled samples. To overcome this problem, we present a method for assigning ¹³C chemical shifts and secondary structures from unresolved two-dimensional ¹³C–¹³C MAS NMR spectra by spectral fitting, named reconstruction of spectra using protein local structures (RESPLS). The spectral fitting was conducted using databases of protein fragmented structures related to ¹³C^α, ¹³C^β, and ¹³C′ chemical shifts and cross-peak intensities. The experimental ¹³C–¹³C inter- and intra-residue correlation spectra of uniformly isotope-labeled ubiquitin in the lyophilized state had a few broad peaks. The fitting analysis for these spectra provided sequence-specific C^α, C^β, and C′ chemical shifts with an accuracy of about 1.5 ppm, which enabled the assignment of the secondary structures with an accuracy of 79 %. The structural heterogeneity of the lyophilized ubiquitin is revealed from the results. Test of RESPLS analysis for simulated spectra of five different types of proteins indicated that the method allowed the secondary structure determination with accuracy of about 80 % for the 50–200 residue proteins. These results demonstrate that the RESPLS approach expands the applicability of the NMR to non-crystalline proteins exhibiting unresolved ¹³C NMR spectra, such as lyophilized proteins, amyloids, membrane proteins and proteins in living cells.     

New 13C-detected experiments for the assignment of intrinsically disordered proteins

NMR assignment of intrinsically disordered proteins (IDPs) by conventional HN-detected methods is hampered by the small dispersion of the amide protons chemical shifts and exchange broadening of amide proton signals. Therefore several alternative assignment strategies have been proposed in the last years. Attempting to seize that dispersion of ¹³C′ and ¹⁵N chemical shifts holds even in IDPs, we recently proposed two ¹³C-detected experiments to directly correlate the chemical shifts of two consecutive ¹³C′–¹⁵N groups in proteins, i.e. without mediation of other nuclei. Main drawback of these experiments is the interruption of the connection at prolines. Here we present new ¹³C-detected experiments to correlate consecutive ¹³C′–¹⁵N groups in IDPs, hacacoNcaNCO and hacaCOncaNCO, that overcome this limitation. Moreover, the experiments provide recognition of glycine residues, thereby facilitating the assignment process.     

Partial site-specific assignment of a uniformly C, N enriched membrane protein, light-harvesting complex 1 (LH1), by solid state NMR

Partial site-specific assignments are reported for the solid state NMR spectra of light-harvesting complex 1, a 160 kDa integral membrane protein. The assignments were derived from 600 MHz ¹⁵N-¹³CO-¹³Cα and ¹⁵N-¹³Cα-¹³CX correlation spectra, using uniformly ¹³C, ¹⁵N enriched hydrated material, in an intact and precipitated form. Sequential assignments were verified using characteristic ¹⁵N-¹³Cα-¹³Cβ side chain chemical shifts observed in 3D experiments. Tertiary contacts found in 2D DARR spectra of the selectively ¹³C enriched sample provided further confirmatory evidence for the assignments. The assignments include the region of the Histidine ligands binding the Bacteriochlorophyll chromophore. The chemical shifts of Cα and Cβ resonances indicated the presence of typical α-helical secondary structure, consistent with previous studies.     

Reduced dimensionality (4,3)D-hnCOCANH experiment: an efficient backbone assignment tool for NMR studies of proteins

Sequence specific resonance assignment of proteins forms the basis for variety of structural and functional proteomics studies by NMR. In this context, an efficient standalone method for rapid assignment of backbone (¹H, ¹⁵N, ¹³C^α and ¹³C′) resonances of proteins has been presented here. Compared to currently available strategies used for the purpose, the method employs only a single reduced dimensionality experiment—(4,3)D-hnCOCANH and exploits the linear combinations of backbone (¹³C^α and ¹³C′) chemical shifts to achieve a dispersion relatively better compared to those of individual chemical shifts (see the text). The resulted increased dispersion of peaks—which is different in sum (CA + CO) and difference (CA ? CO) frequency regions—greatly facilitates the analysis of the spectrum by resolving the problems (associated with routine assignment strategies) arising because of degenerate amide ¹⁵N and backbone ¹³C chemical shifts. Further, the spectrum provides direct distinction between intra- and inter-residue correlations because of their opposite peak signs. The other beneficial feature of the spectrum is that it provides: (a) multiple unidirectional sequential (i→i + 1) ¹⁵N and ¹³C correlations and (b) facile identification of certain specific triplet sequences which serve as check points for mapping the stretches of sequentially connected HSQC cross peaks on to the primary sequence for assigning the resonances sequence specifically. On top of all this, the F ₂–F ₃ planes of the spectrum corresponding to sum (CA + CO) and difference (CA ? CO) chemical shifts enable rapid and unambiguous identification of sequential HSQC peaks through matching their coordinates in these two planes (see the text). Overall, the experiment presented here will serve as an important backbone assignment tool for variety of structural and functional proteomics and drug discovery research programs by NMR involving well behaved small folded proteins (MW < 15 kDa) or a range of intrinsically disordered proteins.      

Accurate measurements of the effects of deuteration at backbone amide positions on the chemical shifts of 15N, 13Cα, 13Cβ, 13CO and 1Hα nuclei in proteins

An approach towards accurate NMR measurements of deuterium isotope effects on the chemical shifts of all backbone nuclei in proteins (¹⁵N, ¹³C_α, ¹³CO, ¹H_α) and ¹³C_β nuclei arising from ¹H-to-D substitutions at amide nitrogen positions is described. Isolation of molecular species with a defined protonation/deuteration pattern at successive backbone nitrogen positions in the polypeptide chain allows quantifying all deuterium isotope shifts of these nuclei from the first to the fourth order. Some of the deuterium isotope shifts measured in the proteins ubiquitin and GB1 can be interpreted in terms of backbone geometry via empirical relationships describing their dependence on (φ; ψ) backbone dihedral angles. Because of their relatively large variability and notable dependence on the protein secondary structure, the two- and three-bond ¹³C_α isotope shifts, ²ΔC_α(N_iD) and ³ΔC_α(N_i+1D), and three-bond ¹³C_β isotope shifts, ³ΔC_β(N_iD), are useful reporters of the local geometry of the protein backbone.     

Practical use of chemical shift databases for protein solid-state NMR: 2D chemical shift maps and amino-acid assignment with secondary-structure information

We introduce a Python-based program that utilizes the large database of ¹³C and ¹⁵N chemical shifts in the Biological Magnetic Resonance Bank to rapidly predict the amino acid type and secondary structure from correlated chemical shifts. The program, called PACSYlite Unified Query (PLUQ), is designed to help assign peaks obtained from 2D ¹³C–¹³C, ¹⁵N–¹³C, or 3D ¹⁵N–¹³C–¹³C magic-angle-spinning correlation spectra. We show secondary-structure specific 2D ¹³C–¹³C correlation maps of all twenty amino acids, constructed from a chemical shift database of 262,209 residues. The maps reveal interesting conformation-dependent chemical shift distributions and facilitate searching of correlation peaks during amino-acid type assignment. Based on these correlations, PLUQ outputs the most likely amino acid types and the associated secondary structures from inputs of experimental chemical shifts. We test the assignment accuracy using four high-quality protein structures. Based on only the Cα and Cβ chemical shifts, the highest-ranked PLUQ assignments were 40–60 % correct in both the amino-acid type and the secondary structure. For three input chemical shifts (CO–Cα–Cβ or N–Cα–Cβ), the first-ranked assignments were correct for 60 % of the residues, while within the top three predictions, the correct assignments were found for 80 % of the residues. PLUQ and the chemical shift maps are expected to be useful at the first stage of sequential assignment, for combination with automated sequential assignment programs, and for highly disordered proteins for which secondary structure analysis is the main goal of structure determination.     

An NMR method for the determination of protein binding interfaces using TEMPOL-induced chemical shift perturbations

Background

The determination of protein–protein interfaces is of crucial importance to understand protein function and to guide the design of compounds. To identify protein–protein interface by NMR spectroscopy, ¹³C NMR paramagnetic shifts induced by freely diffusing 4-hydroxy-2, 2, 6, 6-tetramethyl-piperidine-1-oxyl (TEMPOL) are promising, because TEMPOL affects distinct ¹³C NMR chemical shifts of the solvent accessible nuclei belonging to proteins of interest, while ¹³C nuclei within the interior of the proteins may be distinguished by a lack of such shifts.

Method

We measured the ¹³C NMR paramagnetic shifts induced by TEMPOL by recording ¹³C–¹³C TOCSY spectra for ubiquitin in the free state and the complex state with yeast ubiquitin hydrolase1 (YUH1).

Results

Upon complexation of ubiquitin with YUH1, ¹³C NMR paramagnetic shifts associated with the protein binding interface were reduced by 0.05 ppm or more. The identified interfacial atoms agreed with the prior X-ray crystallographic data.

Conclusions

The TEMPOL-induced ¹³C chemical shift perturbation is useful to determine precise protein–protein interfaces.

General significance

The present method is a useful method to determine protein–protein interface by NMR, because it has advantages in easy sample preparations, simple data analyses, and wide applicabilities.     

Direct correlation of consecutive C′–N groups in proteins: a method for the assignment of intrinsically disordered proteins

Two novel 3D ¹³C-detected experiments, hNcocaNCO and hnCOcaNCO, are proposed to facilitate the resonance assignment of intrinsically disordered proteins. The experiments correlate the ¹⁵N and ¹³C′ chemical shifts of two consecutive amide moieties without involving other nuclei, thus taking advantage of the good dispersion shown by the ¹⁵N–¹³C′ correlations, even for proteins that lack a well defined tertiary structure. The new pulse sequences were successfully tested using Nupr1, an intrinsically disordered protein of 93 residues.     

RefDB: a database of uniformly referenced protein chemical shifts

RefDB is a secondary database of reference-corrected protein chemical shifts derived from the BioMagResBank (BMRB). The database was assembled by using a recently developed program (SHIFTX) to predict protein ¹H, ¹³C and ¹⁵N chemical shifts from X-ray or NMR coordinate data of previously assigned proteins. The predicted shifts were then compared with the corresponding observed shifts and a variety of statistical evaluations performed. In this way, potential mis-assignments, typographical errors and chemical referencing errors could be identified and, in many cases, corrected. This approach allows for an unbiased, instrument-independent solution to the problem of retrospectively re-referencing published protein chemical shifts. Results from this study indicate that nearly 25% of BMRB entries with ¹³C protein assignments and 27% of BMRB entries with ¹⁵N protein assignments required significant chemical shift reference readjustments. Additionally, nearly 40% of protein entries deposited in the BioMagResBank appear to have at least one assignment error. From this study it evident that protein NMR spectroscopists are increasingly adhering to recommended IUPAC ¹³C and ¹⁵N chemical shift referencing conventions, however, approximately 20% of newly deposited protein entries in the BMRB are still being incorrectly referenced. This is cause for some concern. However, the utilization of RefDB and its companion programs may help mitigate this ongoing problem. RefDB is updated weekly and the database, along with its associated software, is freely available at http://redpoll.pharmacy.ualberta.ca and the BMRB website.     

Deuterium isotope shifts for backbone 1H, 15N and 13C nuclei in intrinsically disordered protein α-synuclein

Intrinsically disordered proteins (IDPs) are abundant in nature and characterization of their potential structural propensities remains a widely pursued but challenging task. Analysis of NMR secondary chemical shifts plays an important role in such studies, but the output of such analyses depends on the accuracy of reference random coil chemical shifts. Although uniform perdeuteration of IDPs can dramatically increase spectral resolution, a feature particularly important for the poorly dispersed IDP spectra, the impact of deuterium isotope shifts on random coil values has not yet been fully characterized. Very precise ²H isotope shift measurements for ¹³C^??, ¹³C^??, ¹³C??, ¹⁵N, and ¹H^N have been obtained by using a mixed sample of protonated and uniformly perdeuterated ??-synuclein, a protein with chemical shifts exceptionally close to random coil values. Decomposition of these isotope shifts into one-bond, two-bond and three-bond effects as well as intra- and sequential residue contributions shows that such an analysis, which ignores conformational dependence, is meaningful but does not fully describe the total isotope shift to within the precision of the measurements. Random coil ²H isotope shifts provide an important starting point for analysis of such shifts in structural terms in folded proteins, where they are known to depend strongly on local geometry.     

Automated 1H and 13C chemical shift prediction using the BioMagResBank

A computer program has been developed to accurately and automatically predict the ¹H and ¹³C chemical shifts of unassigned proteins on the basis of sequence homology. The program (called SHIFTY) uses standard sequence alignment techniques to compare the sequence of an unassigned protein against the BioMagResBank – a public database containing sequences and NMR chemical shifts of nearly 200 assigned proteins [Seavey et al. (1991) J. Biomol. NMR, 1, 217–236]. From this initial sequence alignment, the program uses a simple set of rules to directly assign or transfer a complete set of ¹H or ¹³C chemical shifts (from the previously assigned homologues) to the unassigned protein. This homologous assignment protocol takes advantage of the simple fact that homologous proteins tend to share both structural similarity and chemical shift similarity. SHIFTY has been extensively tested on more than 25 medium-sized proteins. Under favorable circumstances, this program can predict the ¹H or ¹³C chemical shifts of proteins with an accuracy far exceeding any other method published to date. With the expo- nential growth in the number of assigned proteins appearing in the literature (now at a rate of more than 150 per year), we believe that SHIFTY may have widespread utility in assigning individual members in families of related proteins, an endeavor that accounts for a growing portion of the protein NMR work being done today.     

NMR study of silk I structure of Bombyx mori silk fibroin with 15N- and 13C-NMR chemical shift contour plots

The metastable state silk I structures of Bombyx mori silk fibroin in the solid state were studied on the basis of ¹⁵N- and ¹³C-nmr chemical shifts of Ala, Ser, and Gly residues. The ¹⁵N cross-polarization magic angle spinning (CP/MAS) nmr spectra of the precipitated fraction after chymotrypsin hydrolysis of B. mori silk fibroin with the silk I and silk II forms were measured to determine the ¹⁵N chemical shifts of Gly, Ala, and Ser residues. For comparison, ¹⁵N CP/MAS nmr chemical shifts of Ala were measured for [¹⁵N] Ala Philosamia cynthia ricini silk fibroin with antiparallel β-sheet and α-helix forms. The ¹³C CP/MAS nmr chemical shifts of Ala, Ser, and Gly residues of B. mori silk fibroin with the silk I and silk II forms, as well as ¹³C CP/MAS nmr chemical shifts of Ala residue of P. c. ricini silk fibroin with β-sheet and α-helix forms, are used for the examination of the silk I structure. Both silk I and α-helix peaks are shifted to a lower field than silk II (β-sheet) for the Cα carbons of the Ala residues, while both Cβ carbon peaks are shifted to higher field. However, the silk I peak of the ¹⁵N nucleus of the Ala residue is shifted to lower field than the silk II peak, but the α-helix peak is shifted to high field. Thus, the difference in the structure between the silk I and α-helix is reflected in a different manner between the ¹³C and ¹⁵N chemical shifts. The Cα and Cβ chemical shift contour plots for Ala and Ser residues, and the Cα plot for the Gly residue, were prepared from the Protein Data Bank data obtained for 12 proteins and used for discussing the silk I structure quantitatively from the conformation-dependent chemical shifts. The plots reported by Le and Oldfield for ¹⁵N chemical shifts were also used for the purpose. All these chemical shift data support Fossey's model (Ala: ϕ = −80°, φ = 150°, Gly: ϕ = −150°, φ = 80°) and do not support Lotz and Keith's model (Ala: ϕ = −104.6°, φ = 112.2°, Gly: ϕ = 79.8°, φ = 49.7° or Ala: ϕ = −124.5°, φ = 88.2°, Gly: ϕ = −49.8°, φ = −76.1°) as the silk I structure. © 1997 John Wiley & Sons, 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.     

CPMG relaxation dispersion NMR experiments measuring glycine 1Hα and 13Cα chemical shifts in the ‘invisible’ excited states of proteins

Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion NMR experiments are extremely powerful for characterizing millisecond time-scale conformational exchange processes in biomolecules. A large number of such CPMG experiments have now emerged for measuring protein backbone chemical shifts of sparsely populated (>0.5%), excited state conformers that cannot be directly detected in NMR spectra and that are invisible to most other biophysical methods as well. A notable deficiency is, however, the absence of CPMG experiments for measurement of ¹H^α and ¹³C^α chemical shifts of glycine residues in the excited state that reflects the fact that in this case the ¹H^α, ¹³C^α spins form a three-spin system that is more complex than the AX ¹H^α–¹³C^α spin systems in the other amino acids. Here pulse sequences for recording ¹H^α and ¹³C^α CPMG relaxation dispersion profiles derived from glycine residues are presented that provide information from which ¹H^α, ¹³C^α chemical shifts can be obtained. The utility of these experiments is demonstrated by an application to a mutant of T4 lysozyme that undergoes a millisecond time-scale exchange process facilitating the binding of hydrophobic ligands to an internal cavity in the protein.     

Comparative analysis of 13C chemical shifts of β-sheet amyloid proteins and outer membrane proteins

Cross-β amyloid fibrils and membrane-bound β-barrels are two important classes of β-sheet proteins. To investigate whether there are systematic differences in the backbone and sidechain conformations of these two families of proteins, here we analyze the ¹³C chemical shifts of 17 amyloid proteins and 7 β-barrel membrane proteins whose high-resolution structures have been determined by NMR. These 24 proteins contain 373 β-sheet residues in amyloid fibrils and 521 β-sheet residues in β-barrel membrane proteins. The ¹³C chemical shifts are shown in 2D ¹³C–¹³C correlation maps, and the amino acid residues are categorized by two criteria: (1) whether they occur in β-strand segments or in loops and turns; (2) whether they are water-exposed or dry, facing other residues or lipids. We also examine the abundance of each amino acid in amyloid proteins and β-barrels and compare the sidechain rotameric populations. The ¹³C chemical shifts indicate that hydrophobic methyl-rich residues and aromatic residues exhibit larger static sidechain conformational disorder in amyloid fibrils than in β-barrels. In comparison, hydroxyl- and amide-containing polar residues have more ordered sidechains and more ordered backbones in amyloid fibrils than in β-barrels. These trends can be explained by steric zipper interactions between β-sheet planes in cross-β fibrils, and by the interactions of β-barrel residues with lipid and water in the membrane. These conformational trends should be useful for structural analysis of amyloid fibrils and β-barrels based principally on NMR chemical shifts.

     

CSI 2.0: a significantly improved version of the Chemical Shift Index

Protein chemical shifts have long been used by NMR spectroscopists to assist with secondary structure assignment and to provide useful distance and torsion angle constraint data for structure determination. One of the most widely used methods for secondary structure identification is called the Chemical Shift Index (CSI). The CSI method uses a simple digital chemical shift filter to locate secondary structures along the protein chain using backbone ¹³C and ¹H chemical shifts. While the CSI method is simple to use and easy to implement, it is only about 75–80 % accurate. Here we describe a significantly improved version of the CSI (2.0) that uses machine-learning techniques to combine all six backbone chemical shifts (¹³C_α, ¹³C_β, ¹³C, ¹⁵N, ¹HN, ¹H_α) with sequence-derived features to perform far more accurate secondary structure identification. Our tests indicate that CSI 2.0 achieved an average identification accuracy (Q3) of 90.56 % for a training set of 181 proteins in a repeated tenfold cross-validation and 89.35 % for a test set of 59 proteins. This represents a significant improvement over other state-of-the-art chemical shift-based methods. In particular, the level of performance of CSI 2.0 is equal to that of standard methods, such as DSSP and STRIDE, used to identify secondary structures via 3D coordinate data. This suggests that CSI 2.0 could be used both in providing accurate NMR constraint data in the early stages of protein structure determination as well as in defining secondary structure locations in the final protein model(s). A CSI 2.0 web server (http://csi.wishartlab.com) is available for submitting the input queries for secondary structure identification.     

Relationship between1H and13C NMR chemical shifts and the secondary and tertiary structure of a zinc finger peptide

Summary Essentially complete assignments have been obtained for the¹H and protonated¹³C NMR spectra of the zinc finger peptide Xfin-31 in the presence and absence of zinc. The patterns observed for the¹H and¹³C chemical shifts of the peptide in the presence of zinc, relative to the shifts in the absence of zinc, reflect the zinc-mediated folding of the unstructured peptide into a compact globular structure with distinct elements of secondary structure. Chemical shifts calculated from the 3D solution structure of the peptide in the presence of zinc and the observed shifts agree to within ca. 0.2 and 0.6 ppm for the backbone C^aH and NH protons, respectively. In addition, homologous zinc finger proteins exhibit similar correlations between their¹H chemical shifts and secondary structure.     

Chemical shift assignments of a minimal Rna14p/Rna15p heterodimer from the yeast cleavage factor IA complex

The two yeast proteins Rna14p and Rna15p form part of the cleavage/polyadenylation factor IA (CF IA) complex that is involved in the 3′ processing of pre-mRNA. Association of the two proteins is mediated by a small C-terminal peptide from Rna14p and a region in Rna15p that corresponds to the hinge domain first identified within the human orthologue. Here I report the ¹H, ¹³C and ¹⁵N spectral assignments for a bacterially co-expressed heterodimer of Rna14p/Rna15p. Further analysis of secondary chemical shifts reveals that both peptides are predominantly α-helical within the complex.     

Site-selective <Superscript>13</Superscript>C labeling of histidine and tryptophan using ribose

Nuclear magnetic resonance spectroscopy studies of ever larger systems have benefited from many different forms of isotope labeling, in particular, site specific isotopic labeling. Site specific ¹³C labeling of methyl groups has become an established means of probing systems not amenable to traditional methodology. However useful, methyl reporter sites can be limited in number and/or location. Therefore, new complementary site specific isotope labeling strategies are valuable. Aromatic amino acids make excellent probes since they are often found at important interaction interfaces and play significant structural roles. Aromatic side chains have many of the same advantages as methyl containing amino acids including distinct ¹³C chemical shifts and multiple magnetically equivalent ¹H positions. Herein we report economical bacterial production and one-step purification of phenylalanine with ¹³C incorporation at the Cα, Cγ and Cε positions, resulting in two isolated ¹H-¹³C spin systems. We also present methodology to maximize incorporation of phenylalanine into recombinantly overexpressed proteins in bacteria and demonstrate compatibility with ILV-methyl labeling. Inexpensive, site specific isotope labeled phenylalanine adds another dimension to biomolecular NMR, opening new avenues of study.