Protein Resonance Assignment Research Articles

Automatic 13C chemical shift reference correction for unassigned protein NMR spectra

Poor chemical shift referencing, especially for 13C in protein Nuclear Magnetic Resonance (NMR) experiments, fundamentally limits and even prevents effective study of biomacromolecules via NMR, including protein structure determination and analysis of protein dynamics. To solve this problem, we constructed a Bayesian probabilistic framework that circumvents the limitations of previous reference correction methods that required protein resonance assignment and/or three-dimensional protein structure. Our algorithm named Bayesian Model Optimized Reference Correction (BaMORC) can detect and correct 13C chemical shift referencing errors before the protein resonance assignment step of analysis and without three-dimensional structure. By combining the BaMORC methodology with a new intra-peaklist grouping algorithm, we created a combined method called Unassigned BaMORC that utilizes only unassigned experimental peak lists and the amino acid sequence. Unassigned BaMORC kept all experimental three-dimensional HN(CO)CACB-type peak lists tested within ± 0.4 ppm of the correct 13C reference value. On a much larger unassigned chemical shift test set, the base method kept 13C chemical shift referencing errors to within ± 0.45 ppm at a 90% confidence interval. With chemical shift assignments, Assigned BaMORC can detect and correct 13C chemical shift referencing errors to within ± 0.22 at a 90% confidence interval. Therefore, Unassigned BaMORC can correct 13C chemical shift referencing errors when it will have the most impact, right before protein resonance assignment and other downstream analyses are started. After assignment, chemical shift reference correction can be further refined with Assigned BaMORC. These new methods will allow non-NMR experts to detect and correct 13C referencing error at critical early data analysis steps, lowering the bar of NMR expertise required for effective protein NMR analysis.

A Slice-based 13C-detected NMR Spin System Forming and Resonance Assignment Method.

Nuclear magnetic resonance (NMR) spectroscopy is attracting more attention in the field of computational structural biology. Till recently, 1H-detected experiments are the dominant NMR technique used due to the high sensitivity of 1H nuclei. However, the current availability of high magnetic fields and cryogenically cooled probe heads allow researchers to overcome the low sensitivity of 13C nuclei. Consequently, 13C-detected experiments have become a popular technique in different NMR applications especially resonance assignment and structure determination of large proteins. In this paper, we propose the first spin system forming method for 13C-detected NMR spectra. Our method is able to accurately form spin systems based on as few as two 13C-detected spectra, CBCACON, and CBCANCO. Our method picks slices from the more trusted spectrum and uses them as feedback to direct the slice picking in the less trusted one. This feedback leads to picking the accurate slices that consequently helps to form better spin systems. We tested our method on a real dataset of 'Ubiquitin' and a benchmark simulated dataset consisting of 12 proteins. We fed our spin systems as inputs to a genetic algorithm to generate the chemical shift assignment, and obtained 92 percent correct chemical shift assignment for Ubiquitin. For the simulated dataset, we obtained an average recall of 86 percent and an average precision of 88 percent. Finally, our chemical shift assignment of Ubiquitin was given as an input to CS-ROSETTA server that generated structures close to the experimentally determined structure.

Mixed pyruvate labeling enables backbone resonance assignment of large proteins using a single experiment

Backbone resonance assignment is a critical first step in the investigation of proteins by NMR. This is traditionally achieved with a standard set of experiments, most of which are not optimal for large proteins. Of these, HNCA is the most sensitive experiment that provides sequential correlations. However, this experiment suffers from chemical shift degeneracy problems during the assignment procedure. We present a strategy that increases the effective resolution of HNCA and enables near-complete resonance assignment using this single HNCA experiment. We utilize a combination of 2-13C and 3-13C pyruvate as the carbon source for isotope labeling, which suppresses the one bond (1Jαβ) coupling providing enhanced resolution for the Cα resonance and amino acid-specific peak shapes that arise from the residual coupling. Using this approach, we can obtain near-complete (>85%) backbone resonance assignment of a 42 kDa protein using a single HNCA experiment.

Accelerating NMR-Based Structural Studies of Proteins by Combining Amino Acid Selective Unlabeling and Fast NMR Methods

In recent years, there has been a growing interest in fast acquisition and analysis of nuclear magnetic resonance (NMR) spectroscopy data for high throughput protein structure determination. Towards this end, rapid data collection techniques and methods to simplify the NMR spectrum such as amino acid selective unlabeling have been proposed recently. Combining these two approaches can speed up further the structure determination process. Based on this idea, we present three new two-dimensional (2D) NMR experiments, which together provide 15N, 1HN, 13Cα, 13Cβ, 13C′ chemical shifts for amino acid residues which are immediate C-terminal neighbors (i + 1) of residues that are selectively unlabeled. These experiments have high sensitivity and can be acquired rapidly using the methodology of G-matrix Fourier transform (GFT) NMR spectroscopy combined with non-uniform sampling (NUS). This is a first study involving the application of fast NMR methods to proteins samples prepared using a specific labeling scheme. Taken together, this opens up new avenues to using the method of selective unlabeling for rapid resonance assignment of proteins.

Solid‐State NMR H–N–(C)–H and H–N–C–C 3D/4D Correlation Experiments for Resonance Assignment of Large Proteins

Solid-state NMR spectroscopy can provide insight into protein structure and dynamics at the atomic level without inherent protein size limitations. However, a major hurdle to studying large proteins by solid-state NMR spectroscopy is related to spectral complexity and resonance overlap, which increase with molecular weight and severely hamper the assignment process. Here the use of two sets of experiments is shown to expand the tool kit of 1 H-detected assignment approaches, which correlate a given amide pair either to the two adjacent CO-CA pairs (4D hCOCANH/hCOCAcoNH), or to the amide 1 H of the neighboring residue (3D HcocaNH/HcacoNH, which can be extended to 5D). The experiments are based on efficient coherence transfers between backbone atoms using INEPT transfers between carbons and cross-polarization for heteronuclear transfers. The utility of these experiments is exemplified with application to assemblies of deuterated, fully amide-protonated proteins from approximately 20 to 60 kDa monomer, at magic-angle spinning (MAS) frequencies from approximately 40 to 55 kHz. These experiments will also be applicable to protonated proteins at higher MAS frequencies. The resonance assignment of a domain within the 50.4 kDa bacteriophage T5 tube protein pb6 is reported, and this is compared to NMR assignments of the isolated domain in solution. This comparison reveals contacts of this domain to the core of the polymeric tail tube assembly.

Detecting and accounting for multiple sources of positional variance in peak list registration analysis and spin system grouping

Peak lists derived from nuclear magnetic resonance (NMR) spectra are commonly used as input data for a variety of computer assisted and automated analyses. These include automated protein resonance assignment and protein structure calculation software tools. Prior to these analyses, peak lists must be aligned to each other and sets of related peaks must be grouped based on common chemical shift dimensions. Even when programs can perform peak grouping, they require the user to provide uniform match tolerances or use default values. However, peak grouping is further complicated by multiple sources of variance in peak position limiting the effectiveness of grouping methods that utilize uniform match tolerances. In addition, no method currently exists for deriving peak positional variances from single peak lists for grouping peaks into spin systems, i.e. spin system grouping within a single peak list. Therefore, we developed a complementary pair of peak list registration analysis and spin system grouping algorithms designed to overcome these limitations. We have implemented these algorithms into an approach that can identify multiple dimension-specific positional variances that exist in a single peak list and group peaks from a single peak list into spin systems. The resulting software tools generate a variety of useful statistics on both a single peak list and pairwise peak list alignment, especially for quality assessment of peak list datasets. We used a range of low and high quality experimental solution NMR and solid-state NMR peak lists to assess performance of our registration analysis and grouping algorithms. Analyses show that an algorithm using a single iteration and uniform match tolerances approach is only able to recover from 50 to 80% of the spin systems due to the presence of multiple sources of variance. Our algorithm recovers additional spin systems by reevaluating match tolerances in multiple iterations. To facilitate evaluation of the algorithms, we developed a peak list simulator within our nmrstarlib package that generates user-defined assigned peak lists from a given BMRB entry or database of entries. In addition, over 100,000 simulated peak lists with one or two sources of variance were generated to evaluate the performance and robustness of these new registration analysis and peak grouping algorithms.

NMR-Based Determination of the 3D Structure of the Ligand-Protein Interaction Site without Protein Resonance Assignment.

Molecular replacement in X-ray crystallography is the prime method for establishing structure-activity relationships of pharmaceutically relevant molecules. Such an approach is not available for NMR. Here, we establish a comparable method, called NMR molecular replacement (NMR(2)). The method requires experimentally measured ligand intramolecular NOEs and ligand-protein intermolecular NOEs as well as a previously known receptor structure or model. Our findings demonstrate that NMR(2) may open a new avenue for the fast and robust determination of the interaction site of ligand-protein complexes at atomic resolution.

Rapid NMR Assignments of Proteins by Using Optimized Combinatorial Selective Unlabeling

A new approach for rapid resonance assignments in proteins based on amino acid selective unlabeling is presented. The method involves choosing a set of multiple amino acid types for selective unlabeling and identifying specific tripeptides surrounding the labeled residues from specific 2D NMR spectra in a combinatorial manner. The methodology directly yields sequence specific assignments, without requiring a contiguously stretch of amino acid residues to be linked, and is applicable to deuterated proteins. We show that a 2D [(15) N,(1) H] HSQC spectrum with two 2D spectra can result in ∼50 % assignments. The methodology was applied to two proteins: an intrinsically disordered protein (12 kDa) and the 29 kDa (268 residue) α-subunit of Escherichia coli tryptophan synthase, which presents a challenging case with spectral overlaps and missing peaks. The method can augment existing approaches and will be useful for applications such as identifying active-site residues involved in ligand binding, phosphorylation, or protein-protein interactions, even prior to complete resonance assignments.

Protein resonance assignment at MAS frequencies approaching 100kHz: a quantitative comparison of J-coupling and dipolar-coupling-based transfer methods.

We discuss the optimum experimental conditions to obtain assignment spectra for solid proteins at magic-angle spinning (MAS) frequencies around 100 kHz. We present a systematic examination of the MAS dependence of the amide proton T 2' times and a site-specific comparison of T 2' at 93 kHz versus 60 kHz MAS frequency. A quantitative analysis of transfer efficiencies of building blocks, as they are used for typical 3D experiments, was performed. To do this, we compared dipolar-coupling and J-coupling based transfer steps. The building blocks were then combined into 3D experiments for sequential resonance assignment, where we evaluated signal-to-noise ratio and information content of the different 3D spectra in order to identify the best assignment strategy. Based on this comparison, six experiments were selected to optimally assign the model protein ubiquitin, solely using spectra acquired at 93 kHz MAS. Within 3 days of instrument time, the required spectra were recorded from which the backbone resonances have been assigned to over 96%.

Probing Fleeting Interactions in Large and Dynamic Nonribosomal Peptide Synthetases with Novel NMR Methods

Nonribosomal peptide synthetases (NRPSs) are modular, multi-domain, enzymatic assembly lines that combine simple substrates to synthesize complex natural products (e.g. penicillin, bacitracin, yersiniabactin). Substrates are each covalently attached to so-called carrier proteins (CP) via a post-translationally introduced phosphopantetheine arm. CPs present the substrates to catalytic domains in a series of transient, sequential domain/domain and domain/substrate interactions during biosynthesis. Thus, understanding the mechanism of NRPS assembly line synthesis necessitates both characterizing fleeting molecular interactions and unraveling the molecular determinants for the directionality of these interactions. Unfortunately, the large molecular weight of NRPS domains, the presence of dynamics between domains and within domains, and the lability of tethered substrates have precluded studying these processes at the molecular level. Here, we present methods that permitted us to overcome these limitations and we discuss our findings. First, by exploiting the non-invasive nature of NMR we found that chemical substrates interact with their carrier proteins and these interactions may modulate domain communication to provide directionality. Second, we present a method that facilitates resonance assignments in large proteins. By applying covariance to conventional spectra, simplified correlation maps permit to identify sequential residues or to assign side-chain resonances with a simple visual inspection. Hence, resonances of 37 kDa and 53 kDa NRPS domains could be assigned unambiguously, enabling us to describe molecular signatures of transient interactions involving these domains.

Accurate ab initio prediction of NMR chemical shifts of nucleic acids and nucleic acids/protein complexes.

NMR chemical shift predictions based on empirical methods are nowadays indispensable tools during resonance assignment and 3D structure calculation of proteins. However, owing to the very limited statistical data basis, such methods are still in their infancy in the field of nucleic acids, especially when non-canonical structures and nucleic acid complexes are considered. Here, we present an ab initio approach for predicting proton chemical shifts of arbitrary nucleic acid structures based on state-of-the-art fragment-based quantum chemical calculations. We tested our prediction method on a diverse set of nucleic acid structures including double-stranded DNA, hairpins, DNA/protein complexes and chemically-modified DNA. Overall, our quantum chemical calculations yield highly/very accurate predictions with mean absolute deviations of 0.3–0.6 ppm and correlation coefficients (r2) usually above 0.9. This will allow for identifying misassignments and validating 3D structures. Furthermore, our calculations reveal that chemical shifts of protons involved in hydrogen bonding are predicted significantly less accurately. This is in part caused by insufficient inclusion of solvation effects. However, it also points toward shortcomings of current force fields used for structure determination of nucleic acids. Our quantum chemical calculations could therefore provide input for force field optimization.

Automated resonance assignment of the 21 kDa stereo-array isotope labeled thioldisulfide oxidoreductase DsbA

The automated chemical shift assignment algorithm FLYA has been extended for use with stereo-array isotope labeled (SAIL) proteins to determine the sequence-specific resonance assignments of large proteins. Here we present the assignment of the backbone and sidechain chemical shifts of the 21kDa thioldisulfide oxidoreductase DsbA from Escherichia coli that were determined with the SAIL-FLYA algorithm in conjunction with automated peak picking. No manual corrections of peak lists or assignments were applied. The assignments agreed with manually determined reference assignments in 95.4% of the cases if 16 input spectra were used, 94.1% if only 3D 13C/15N-resolved NOESY, CBCA(CO)NH, and 2D [13C/15N,1H]-HSQC were used, and 86.8% if exclusively 3D 13C/15N-resolved NOESY spectra were used. Considering only the assignments that are classified as reliable by the SAIL-FLYA algorithm, the degrees of agreement increased to 97.5%, 96.5%, and 94.2%, respectively. With our approach it is thus possible to automatically obtain almost complete and correct assignments of proteins larger than 20kDa.

Pseudo 5D HN(C)N experiment to facilitate the assignment of backbone resonances in proteins exhibiting high backbone shift degeneracy

Assignment of protein backbone resonances is most routinely carried out using triple resonance three-dimensional NMR experiments involving amide 1H/15N resonances. However for intrinsically unstructured proteins, alpha-helical proteins or proteins containing several disordered fragments, the assignment becomes problematic because of high-degree of backbone shift degeneracy. In this backdrop, a novel reduced-dimensionality (RD) experiment –(5,3)D-hNCO-CANH- is presented to facilitate/validate the sequential backbone resonance assignment in such proteins. The proposed 3D NMR experiment makes use of the modulated amide 15N chemical shifts (resulting from the joint sampling along both its indirect dimensions) to resolve the ambiguity involved in connecting the neighboring amide resonances (i.e. HiNi and Hi−1Ni−1) for overlapping amide-NH peaks. The experiment -in combination with routine triple resonance 3D-NMR experiments involving backbone amide (1H/15N) and carbon (13Cα/13C′) chemical shifts- will serve as a powerful complementary tool to achieve the nearly complete assignment of protein backbone resonances in a time efficient manner.

BSH-CP based 3D solid-state NMR experiments for protein resonance assignment

We have recently presented band-selective homonuclear cross-polarization (BSH-CP) as an efficient method for CO-CA transfer in deuterated as well as protonated solid proteins. Here we show how the BSH-CP CO-CA transfer block can be incorporated in a set of three-dimensional (3D) solid-state NMR (ssNMR) pulse schemes tailored for resonance assignment of proteins at high static magnetic fields and moderate magic-angle spinning rates. Due to the achieved excellent transfer efficiency of 33 % for BSH-CP, a complete set of 3D spectra needed for unambiguous resonance assignment could be rapidly recorded within 1 week for the model protein ubiquitin. Thus we expect that BSH-CP could replace the typically used CO-CA transfer schemes in well-established 3D ssNMR approaches for resonance assignment of solid biomolecules.

Reduced dimensionality tailored HN(C)N experiments for facile backbone resonance assignment of proteins through unambiguous identification of sequential HSQC peaks

Two novel reduced dimensionality (RD) tailored HN(C)N [S.C. Panchal, N.S. Bhavesh, R.V. Hosur, Improved 3D triple resonance experiments, HNN and HN(C)N, for HN and 15N sequential correlations in (13C, 15N) labeled proteins: application to unfolded proteins, J. Biomol. NMR 20 (2001) 135–147] experiments are proposed to facilitate the backbone resonance assignment of proteins both in terms of its accuracy and speed. These experiments – referred here as (4,3)D-hNCOcaNH and (4,3)D-hNcoCANH – exploit the linear combination of backbone 15N and 13C′/13Cα chemical shifts simultaneously to achieve higher peak dispersion and randomness along their respective F1 dimensions. Simply, this has been achieved by modulating the backbone 15N(i) chemical shifts with that of 13C′ (i−1)/13Cα (i−1) spins following the established reduced dimensionality NMR approach [T. Szyperski, D.C. Yeh, D.K. Sukumaran, H.N. Moseley, G.T. Montelione, Reduced-dimensionality NMR spectroscopy for high-throughput protein resonance assignment, Proc. Natl. Acad. Sci. USA 99 (2002) 8009–8014]. Though the modification is simple it has resulted an ingenious improvement of HN(C)N both in terms of peak dispersion and easiness of establishing the sequential connectivities. The increased dispersion along F1 dimension solves two purposes here: (i) resolves the ambiguities arising because of degenerate 15N chemical shifts and (ii) reduces the signal overlap in F2(15N)–F3(1H) planes (an important requisite in HN(C)N based assignment protocol for facile and unambiguous identification of sequentially connected HSQC peaks). The performance of both these experiments and the assignment protocol has been demonstrated using bovine apo Calbindin-d9k (75 aa) and urea denatured UNC60B (a 152 amino acid ADF/cofilin family protein of Caenorhabditis elegans), as representatives of folded and unfolded protein systems, respectively.

Reliability of exclusively NOESY-based automated resonance assignment and structure determination of proteins

Protein structure determination by NMR can in principle be speeded up both by reducing the measurement time on the NMR spectrometer and by a more efficient analysis of the spectra. Here we study the reliability of protein structure determination based on a single type of spectra, namely nuclear Overhauser effect spectroscopy (NOESY), using a fully automated procedure for the sequence-specific resonance assignment with the recently introduced FLYA algorithm, followed by combined automated NOE distance restraint assignment and structure calculation with CYANA. This NOESY-FLYA method was applied to eight proteins with 63-160 residues for which resonance assignments and solution structures had previously been determined by the Northeast Structural Genomics Consortium (NESG), and unrefined and refined NOESY data sets have been made available for the Critical Assessment of Automated Structure Determination of Proteins by NMR project. Using only peak lists from three-dimensional (13)C- or (15)N-resolved NOESY spectra as input, the FLYA algorithm yielded for the eight proteins 91-98 % correct backbone and side-chain assignments if manually refined peak lists are used, and 64-96 % correct assignments based on raw peak lists. Subsequent structure calculations with CYANA then produced structures with root-mean-square deviation (RMSD) values to the manually determined reference structures of 0.8-2.0 Å if refined peak lists are used. With raw peak lists, calculations for 4 proteins converged resulting in RMSDs to the reference structure of 0.8-2.8 Å, whereas no convergence was obtained for the four other proteins (two of which did already not converge with the correct manual resonance assignments given as input). These results show that, given high-quality experimental NOESY peak lists, the chemical shift assignments can be uncovered, without any recourse to traditional through-bond type assignment experiments, to an extent that is sufficient for calculating accurate three-dimensional structures.

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 ((1)H, (15)N, (13)C(α) and (13)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 ((13)C(α) and (13)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 (15)N and backbone (13)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) (15)N and (13)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.

Out-and-back 13C–13C scalar transfers in protein resonance assignment by proton-detected solid-state NMR under ultra-fast MAS

We present here (1)H-detected triple-resonance H/N/C experiments that incorporate CO-CA and CA-CB out-and-back scalar-transfer blocks optimized for robust resonance assignment in biosolids under ultra-fast magic-angle spinning (MAS). The first experiment, (H)(CO)CA(CO)NH, yields (1)H-detected inter-residue correlations, in which we record the chemical shifts of the CA spins in the first indirect dimension while during the scalar-transfer delays the coherences are present only on the longer-lived CO spins. The second experiment, (H)(CA)CB(CA)NH, correlates the side-chain CB chemical shifts with the NH of the same residue. These high sensitivity experiments are demonstrated on both fully-protonated and 100%-H(N) back-protonated perdeuterated microcrystalline samples of Acinetobacter phage 205 (AP205) capsids at 60 kHz MAS.

Solution NMR resonance assignment strategies for β‐barrel membrane proteins

Membrane proteins in detergent micelles are large and dynamic complexes that present challenges for solution NMR investigations such as spectral overlap and line broadening. In this study, multiple methods are introduced to facilitate resonance assignment of β-barrel membrane proteins using Opa60 from Neisseria gonorrhoeae as a model system. Opa60 is an eight-stranded β-barrel with long extracellular loops (∼63% of the protein) that engage host receptors and induce engulfment of the bacterium. The NMR spectra of Opa60 in detergent micelles exhibits significant spectral overlap and resonances corresponding to the loop regions had variable line widths, which interfered with a complete assignment of the protein. To assign the β-barrel residues, trypsin cleavage was used to remove much of the extracellular loops while preserving the detergent solubilized β-barrel. The removal of the loop resonances significantly improved the assignment of the Opa60 β-barrel region (97% of the resonances corresponding to the β-barrel and periplasmic turns were assigned). For the loop resonance assignments, two strategies were implemented; modulating temperature and synthetic peptides. Lowering the temperature broadened many peaks beyond detection and simplified the spectra to only the most dynamic regions of the loops facilitating 27 loop resonances to be assigned. To further assign functionally important and unstructured regions of the extracellular loops, a synthetic 20 amino acid peptide was synthesized and had nearly complete spectral overlap with the full-length protein allowing 17 loop resonances to be assigned. Collectively, these strategies are effective tools that may accelerate solution NMR structure determination of β-barrel membrane proteins.

High‐Dimensional NMR Spectra for Structural Studies of Biomolecules

Recent developments in the acquisition and processing of NMR data sets facilitate the recording of ultra-high-resolution NMR spectra in a reasonable time. The new experiments allow easy resonance assignment for folded and unfolded proteins, as well as the precise determination of spectral parameters, for example, chemical shifts, NOE contacts, coupling constants or cross-correlated relaxation rates. Owing to exceptional resolution of 4D-6D spectroscopy, detailed studies of biomolecules of unprecedented complexity are now possible. Herein, the principles of acquisition and processing methods are presented. The main applications of high-dimensional NMR experiments, including backbone and side-chain resonance assignment in proteins, as well as heteronuclear edited NOE techniques are reviewed.