A new algorithm for reliable and general NMR resonance assignment

The new FLYA automated resonance assignment algorithm determines NMR chemical shift assignments on the basis of peak lists from any combination of multidimensional through-bond or through-space NMR experiments for proteins. Backbone and side-chain assignments can be determined. All experimental data are used simultaneously, thereby exploiting optimally the redundancy present in the input peak lists and circumventing potential pitfalls of assignment strategies in which results obtained in a given step remain fixed input data for subsequent steps. Instead of prescribing a specific assignment strategy, the FLYA resonance assignment algorithm requires only experimental peak lists and the primary structure of the protein, from which the peaks expected in a given spectrum can be generated by applying a set of rules, defined in a straightforward way by specifying through-bond or through-space magnetization transfer pathways. The algorithm determines the resonance assignment by finding an optimal mapping between the set of expected peaks that are assigned by definition but have unknown positions and the set of measured peaks in the input peak lists that are initially unassigned but have a known position in the spectrum. Using peak lists obtained by purely automated peak picking from the experimental spectra of three proteins, FLYA assigned correctly 96-99% of the backbone and 90-91% of all resonances that could be assigned manually. Systematic studies quantified the impact of various factors on the assignment accuracy, namely the extent of missing real peaks and the amount of additional artifact peaks in the input peak lists, as well as the accuracy of the peak positions. Comparing the resonance assignments from FLYA with those obtained from two other existing algorithms showed that using identical experimental input data these other algorithms yielded significantly (40-142%) more erroneous assignments than FLYA. The FLYA resonance assignment algorithm thus has the reliability and flexibility to replace most manual and semi-automatic assignment procedures for NMR studies of proteins.     

Influence of NMR Data Completeness on Structure Determinations of Homodimeric Proteins

The structure determination of homodimeric proteins by NMR using conventional NOESY experiments is still challenging due to the degeneracy of the chemical shifts in the identical monomers, which causes ambiguity in the NOE assignments. Residues involved in the interface between two monomers provide essential intermolecular NOEs for the structure determinations of homodimeric proteins. Hence NMR data, such as NOE peak lists and chemical shift assignments of these interface residues, play a crucial role for the successful structure determination of homodimeric proteins. This paper extends our previous report (Lin, Y.‐J.; Kirchner, D. K.; Güntert, P. J. Magn. Reson.­ 2012 , 222, 96) and investigates the influence of incomplete NOESY peak lists combined with incomplete ¹H chemical shift assignments of the interface residues on the structure determination of homodimeric proteins using the program CYANA. Data incompleteness was simulated by random omission of both NOESY cross peaks and interface ¹H chemical shifts. Our results for three proteins with different percentages of interface residues reveal that the algorithm can tolerate about 40–50% NOESY peak omission with complete interface chemical shift assignments, which indicates that partial NOESY peak omission does not cause severe problems when the interface chemical shifts are completely assigned. Combining NOESY peak omission with incomplete interface chemical shift assignments, the tolerance for interface chemical shift omission decreases with the extent of omitted NOESY peaks. The tolerance for unassigned interface side chain, methyl and aromatic chemical shifts is affected more strongly by NOESY peak omission than that for the omission of general interface ¹H chemical shifts including the backbone. In general about 10–30% peaks omission is tolerated in conjunction with missing chemical shift assignments. If more NOESY peaks are omitted calculations gradually become unstable and tend not to tolerate any missing interface chemical shifts. A large amount of omitted NOESY peaks, for instance 30% omission in our calculations, could decrease the tolerance for missing aromatic or methyl interface ¹H chemical shifts to as few as 2–4 missing chemical shifts, suggesting that complete aromatic and methyl ¹H chemical shift assignments are important when the NOESY peak data is significantly incomplete. Finally, for homodimeric proteins with a low percentage of interface residues, our results reveal that the omission of NOESY peaks, even at an extent of only 10%, can result in no tolerance against the omission of interface ¹H chemical shifts, suggesting that the completeness of both interface ¹H chemical shift assignments and NOESY peaks are important for the successful structure determination of proteins with a small homodimer interface.     

Sequence-specific resonance assignment of soluble nonglobular proteins by 7D APSY-NMR spectroscopy

Based on sequence-specific resonance assignments, NMR is the method of choice for obtaining atomic-resolution experimental data on soluble nonglobular proteins. So far, however, NMR assignment of unfolded polypeptides in solution has been a time-consuming task, mainly due to the small chemical shift dispersion, which has limited practical applications of the NMR approach. This paper presents an efficient, fully automated method for sequence-specific backbone and beta-carbon NMR assignment of soluble nonglobular proteins with sizes up to at least 150 residues. The procedure is based on new APSY (automated projection spectroscopy) experiments which benefit from the short effective rotational correlation times in soluble nonglobular polypeptides to record five- to seven-dimensional NMR data sets, which reliably resolves chemical shift degeneracies. Fully automated sequence-specific resonance assignments of the backbone nuclei and C(beta) are described for the uniformly (13)C,(15)N-labeled urea-denatured 148-residue outer membrane protein X (OmpX) from E. coli. The method is generally applicable to systems with similar spectroscopic properties as unfolded OmpX, and we anticipate that this paper may open the door for extensive atomic-resolution studies of chemical denaturant-unfolded proteins, as well as some classes of functional nonglobular polypeptides in solution.     

Protein structure determination by high-resolution solid-state NMR spectroscopy: application to microcrystalline ubiquitin

High-resolution solid-state NMR spectroscopy has become a promising method for the determination of three-dimensional protein structures for systems which are difficult to crystallize or exhibit low solubility. Here we describe the structure determination of microcrystalline ubiquitin using 2D (13)C-(13)C correlation spectroscopy under magic angle spinning conditions. High-resolution (13)C spectra have been acquired from hydrated microcrystals of site-directed (13)C-enriched ubiquitin. Inter-residue carbon-carbon distance constraints defining the global protein structure have been evaluated from 'dipolar-assisted rotational resonance' experiments recorded at various mixing times. Additional constraints on the backbone torsion angles have been derived from chemical shift analysis. Using both distance and dihedral angle constraints, the structure of microcrystalline ubiquitin has been refined to a root-mean-square deviation of about 1 A. The structure determination strategies for solid samples described herein are likely to be generally applicable to many proteins that cannot be studied by X-ray crystallography or solution NMR spectroscopy.     

Structure calculation from unambiguous long-range amide and methyl 1H-1H distance restraints for a microcrystalline protein with MAS solid-state NMR spectroscopy

Magic-angle spinning (MAS) solid-state NMR becomes an increasingly important tool for the determination of structures of membrane proteins and amyloid fibrils. Extensive deuteration of the protein allows multidimensional experiments with exceptionally high sensitivity and resolution to be obtained. Here we present an experimental strategy to measure highly unambiguous spatial correlations for distances up to 13 ?. Two complementary three-dimensional experiments, or alternatively a four-dimensional experiment, yield highly unambiguous cross-peak assignments, which rely on four encoded chemical shift dimensions. Correlations to residual aliphatic protons are accessible via synchronous evolution of the (15)N and (13)C chemical shifts, which encode valuable amide-methyl distance restraints. On average, we obtain six restraints per residue. Importantly, 50% of all restraints correspond to long-range distances between residues i and j with |i - j| > 5, which are of particular importance in structure calculations. Using ARIA, we calculate a high-resolution structure for the microcrystalline 7.2 kDa α-spectrin SH3 domain with a backbone precision of ～1.1 ?.     

Automated NMR resonance assignment of large proteins for protein-ligand interaction studies

The detection and structural characterization of protein-ligand interactions by solution NMR is central to functional biology research as well as to drug discovery. Here we present a robust and highly automated procedure for obtaining the resonance assignments necessary for studies of such interactions. The procedure relies on a combination of three automated projection spectroscopy (APSY) experiments, including the new 4D APSY-HNCACB, and the use of fractionally deuterated protein samples. This labeling pattern increases the experimental sensitivity on the one hand, but it leads to peak multiplets on the other hand. The latter complications are however overcome by the geometric APSY analysis of the projection spectra. The three APSY experiments thus provide high precision chemical shift correlations of the backbone and side chain methyl groups, allowing a reliable and robust assignment of the protein by suitable algorithms. The present approach doubles the molecular size limit of APSY-based assignments to 25 kDa, thus providing the basis for efficient characterization of protein-ligand interactions at atomic resolution by NMR, such as structure-based drug design. We show the application to two human proteins with molecular weights of 15 and 22 kDa, respectively, at concentrations of 0.4 mM and discuss the general applicability to studies of protein-protein and protein-nucleic acid complexes.     

Semiautomatic sequence-specific assignment of proteins based on the tertiary structure--the program st2nmr

The sequence-specific assignment of resonances is still the most time-consuming procedure that is necessary as the first step in high-resolution NMR studies of proteins. In many cases a reliable three-dimensional (3D) structure of the protein is available, for example, from X-ray spectroscopy or homology modeling. Here we introduce the st2nmr program that uses the 3D structure and Nuclear Overhauser Effect spectroscopy (NOESY) peak list(s) to evaluate and optimize trial sequence-specific assignments of spin systems derived from correlation spectra to residues of the protein. A distance-dependent target function that scores trial assignments based on the presence of expected NOESY crosspeaks is optimized in a Monte Carlo fashion. The performance of the program st2nmr is tested on real NMR data of an alpha-helical (cytochrome c) and beta-sheet (lipocalin) protein using homology models and/or X-ray structures; it succeeded in completely reproducing the correct sequence-specific assignments in most cases using 2D and/or 15N/13C Nuclear Overhauser Effect (NOE) data. Additionally to amino acid residues the program can also handle ligands that are bound to the protein, such as heme, and can be used as a complementary tool to fully automated assignment procedures.     

3D structure determination of the Crh protein from highly ambiguous solid-state NMR restraints

In a wide variety of proteins, insolubility presents a challenge to structural biology, as X-ray crystallography and liquid-state NMR are unsuitable. Indeed, no general approach is available as of today for studying the three-dimensional structures of membrane proteins and protein fibrils. We here demonstrate, at the example of the microcrystalline model protein Crh, how high-resolution 3D structures can be derived from magic-angle spinning solid-state NMR distance restraints for fully labeled protein samples. First, we show that proton-mediated rare-spin correlation spectra, as well as carbon-13 spin diffusion experiments, provide enough short, medium, and long-range structural restraints to obtain high-resolution structures of this 2 x 10.4 kDa dimeric protein. Nevertheless, the large number of 13C/15N spins present in this protein, combined with solid-state NMR line widths of about 0.5-1 ppm, induces substantial ambiguities in resonance assignments, preventing 3D structure determination by using distance restraints uniquely assigned on the basis of their chemical shifts. In the second part, we thus demonstrate that an automated iterative assignment algorithm implemented in a dedicated solid-state NMR version of the program ARIA permits to resolve the majority of ambiguities and to calculate a de novo 3D structure from highly ambiguous solid-state NMR data, using a unique fully labeled protein sample. We present, using distance restraints obtained through the iterative assignment process, as well as dihedral angle restraints predicted from chemical shifts, the 3D structure of the fully labeled Crh dimer refined at a root-mean-square deviation of 1.33 A.     

Structure determination of a membrane protein with two trans-membrane helices in aligned phospholipid bicelles by solid-state NMR spectroscopy

The structure of the membrane protein MerFt was determined in magnetically aligned phospholipid bicelles by solid-state NMR spectroscopy. With two trans-membrane helices and a 10-residue inter-helical loop, this truncated construct of the mercury transport membrane protein MerF has sufficient structural complexity to demonstrate the feasibility of determining the structures of polytopic membrane proteins in their native phospholipid bilayer environment under physiological conditions. PISEMA, SAMMY, and other double-resonance experiments were applied to uniformly and selectively (15)N-labeled samples to resolve and assign the backbone amide resonances and to measure the associated (15)N chemical shift and (1)H-(15)N heteronuclear dipolar coupling frequencies as orientation constraints for structure calculations. (1)H/(13)C/(15)N triple-resonance experiments were applied to selectively (13)C'- and (15)N-labeled samples to complete the resonance assignments, especially for residues in the nonhelical regions of the protein. A single resonance is observed for each labeled site in one- and two-dimensional spectra. Therefore, each residue has a unique conformation, and all protein molecules in the sample have the same three-dimensional structure and are oriented identically in planar phospholipid bilayers. Combined with the absence of significant intensity near the isotropic resonance frequency, this demonstrates that the entire protein, including the loop and terminal regions, has a well-defined, stable structure in phospholipid bilayers.     

Tryptophan chemical shift in peptides and proteins: a solid state carbon-13 nuclear magnetic resonance spectroscopic and quantum chemical investigation

We have obtained the carbon-13 nuclear magnetic resonance spectra of a series of tryptophan-containing peptides and model systems, together with their X-ray crystallographic structures, and used quantum chemical methods to predict the (13)C NMR shifts (or shieldings) of all nonprotonated aromatic carbons (C(gamma), C(delta 2) and C(epsilon 2). Overall, there is generally good accord between theory and experiment. The chemical shifts of Trp C(gamma) in several proteins, hen egg white lysozyme, horse myoglobin, horse heart cytochrome c, and four carbonmonoxyhemoglobins, are also well predicted. The overall Trp C(gamma) shift range seen in the peptides and proteins is 11.4 ppm, and individual shifts (or shieldings) are predicted with an rms error of approximately 1.4 ppm (R value = 0.86). Unlike C(alpha) and N(H) chemical shifts, which are primarily a function of the backbone phi,psi torsion angles, the Trp C(gamma) shifts are shown to be correlated with the side-chain torsion angles chi(1) and chi(2) and appear to arise, at least in part, from gamma-gauche interactions with the backbone C' and N(H) atoms. This work helps solve the problem of the chemical shift nonequivalences of nonprotonated aromatic carbons in proteins first identified over 30 years ago and opens up the possibility of using aromatic carbon chemical shift information in structure determination.     

Automated NMR assignment of protein side chain resonances using automated projection spectroscopy (APSY)

This paper describes an automated method for sequence-specific NMR assignment of the aliphatic resonances of protein side chains in small- and medium-sized globular proteins in aqueous solution. The method requires the recording of a five-dimensional (5D) automated projection spectroscopy (APSY-) NMR experiment and the subsequent analysis of the APSY peak list with the algorithm ALASCA (Algorithm for local and linear assignment of side chains from APSY data). The 5D APSY-HC(CC-TOCSY)CONH experiment yields 5D chemical shift correlations of aliphatic side chain C-H moieties with the backbone atoms H(N), N, and C'. A simultaneous variation of the TOCSY mixing times and the projection angles in this APSY-type TOCSY experiment gives access to all aliphatic C-H moieties in the 20 proteinogenic amino acids. The correlation peak list resulting from the 5D APSY-HC(CC-TOCSY)CONH experiment together with the backbone assignment of the protein under study is the sole input for the algorithm ALASCA that assigns carbon and proton resonances of protein side chains. The algorithm is described, and it is shown that the aliphatic parts of 17 of the 20 common amino acid side chains are assigned unambiguously, whereas the remaining three amino acids are assigned with a certainty of above 95%. The overall feasibility of the approach is demonstrated with the globular 116-residue protein TM1290, for which reference assignments are known. For this protein, 97% of the expected side chain carbon atoms and 87% of the expected side chain protons were detected with the 5D APSY-HC(CC-TOCSY)CONH experiment in 24 h of spectrometer time, and all these resonances were correctly assigned by ALASCA. Based on the experience with TM1290, we expect that the approach presented in this work is routinely applicable to globular proteins with sizes up to at least 120 amino acids.     

Protein NMR Resonance Assignment without Spectral Analysis: 5D SOlid‐State Automated Projection SpectroscopY (SO‐APSY)

Narrow proton signals, high sensitivity, and efficient coherence transfers provided by fast magic‐angle spinning at high magnetic fields make automated projection spectroscopy feasible for the solid‐state NMR analysis of proteins. We present the first ultrahigh dimensional implementation of this approach, where 5D peak lists are reconstructed from a number of 2D projections for protein samples of different molecular sizes and aggregation states, which show limited dispersion of chemical shifts or inhomogeneous broadenings. The resulting datasets are particularly suitable to automated analysis and yield rapid and unbiased assignments of backbone resonances.     

Isotropic chemical shifts in magic-angle spinning NMR spectra of proteins

Here we examine the effect of magic-angle spinning (MAS) rate upon lineshape and observed peak position for backbone carbonyl (C') peaks in NMR spectra of uniformly-(13)C,15N-labeled (U-(13)C,15N) solid proteins. 2D N-C' spectra of U-(13)C,15N microcrystalline protein GB1 were acquired at six MAS rates, and the site-resolved C' lineshapes were analyzed by numerical simulations and comparison to spectra from a sparsely labeled sample (derived from 1,3-(13)C-glycerol). Spectra of the U-(13)C,15N sample demonstrate large variations in the signal-to-noise ratio and peak positions, which are absent in spectra of the sparsely labeled sample, in which most 13C' sites do not possess a directly bonded 13CA. These effects therefore are a consequence of rotational resonance, which is a well-known phenomenon. Yet the magnitude of this effect pertaining to chemical shift assignment has not previously been examined. To quantify these effects in high-resolution protein spectra, we performed exact numerical two- and four-spin simulations of the C' lineshapes, which reproduced the experimentally observed features. Observed peak positions differ from the isotropic shift by up to 1.0 ppm, even for MAS rates relatively far (a few ppm) from rotational resonance. Although under these circumstances the correct isotropic chemical shift values may be determined through simulation, systematic errors are minimized when the MAS rate is equivalent to approximately 85 ppm for 13C. This moderate MAS condition simplifies spectral assignment and enables data sets from different labeling patterns and spinning rates to be used most efficiently for structure determination.     

Dipolar waves as NMR maps of protein structure

The anisotropy of nuclear spin interactions results in a unique mapping of structure to the resonance frequencies and split tings observed in NMR spectra, however, the determination of molecular structure from experimentally measured spectral parameters is complicated by angular ambiguities resulting from the symmetry properties of dipole-dipole and chemical shift interactions. This issue can be addressed through the periodicity inherent in secondary structure elements, which can be used as an index of topology. Distinctive wheel-like patterns are observed in two-dimensional 1H-15N heteronuclear dipolar/15N chemical shift PISEMA (polarization inversion spin-exchange at the magic angle) spectra of helical membrane proteins in highly aligned lipid bilayer samples. One-dimensional dipolar waves are an extension of two-dimensional PISA (polarity index slant angle) wheels to map protein structure in NMR spectra of both highly and weakly aligned samples. Dipolar waves describe the periodic wavelike variations of the magnitudes of the static heteronuclear dipolar couplings as a function of residue number in the absence of chemical shift effects. Weakly aligned samples of proteins display these same effects, primarily as residual dipolar couplings (RDCs), in solution NMR spectra. The corresponding properties of the RDCs in solution NMR spectra of weakly aligned helices represent a convergence of solid-state and solution NMR approaches to structure determination.     

酞侧基聚芳醚酮的核磁共振研究(Ⅰ)──链结构及NMR谱的归属

用1D NMR方法研究酞侧基聚芳醚酮(PEK-C)链结构,用二维同核化学位移相关与二维异核化学位移相关实验方法对1D NMR谱峰进行归属,探讨了二维异核远程相关实验在缩聚高分子研究中的应用,为PEK-C修饰机理以及共混相容机理的研究提供重要信息。溶液NMR谱图数据表明,PEK-C具有较规整的链结构。     

Direct use of unassigned resonances in NMR structure calculations with proxy residues

We present a method that significantly enhances the robustness of (automated) NMR structure determination by allowing the NOE data corresponding to unassigned NMR resonances to be used directly in the calculations. The unassigned resonances are represented by additional atoms or groups of atoms that have no interaction with the regular protein atoms except through distance restraints. These so-called "proxy" residues can be used to generate NOE-based distance restraints in a similar fashion as for the assigned part of the protein. If sufficient NOE information is available, the restraints are expected to place the proxies at positions close to the correct atoms for the unassigned resonance, which can facilitate subsequent assignment. Convergence can be further improved by supplying additional information about the possible identities of the unassigned resonances. We have implemented this approach in the widely used automated assignment and structure calculation protocols ARIA and CANDID. We find that it significantly increases the robustness of structure calculations with regard to missing assignments and yields structures of higher quality. Our approach is still able to find correctly folded structures with up to 30% randomly missing resonance assignments, and even when only backbone and beta resonances are present! This should be of significant value to NMR-based structural proteomics initiatives.     

Assignment of the backbone resonances for microcrystalline ubiquitin

Site-specific assignments for the solid-state NMR spectra of uniformly 13C,15N-enriched ubiquitin are described. The assignments are derived from three three-dimensional 15N-13C-13C correlation spectra collected at 400 MHz on microcrystalline material. A few residues (the loop near Threonine 9 and the C-terminal fragment) were missing and correspond to regions previously reported to be mobile on the basis of X-ray crystallography and solution NMR studies. A few additional sites exhibit shifts that differ from previously reported solution NMR assignments. Nonetheless, these de novo assignments indicate close agreement between the chemical shifts observed in solution and those in microcrystalline or precipitated solids. The methods utilized are likely to be generally applicable for other noncrystalline, nonsoluble proteins.     

Solid-State NMR of a Large Membrane Protein by Paramagnetic Relaxation Enhancement

Membrane proteins play an important role in many biological functions. Solid-state NMR spectroscopy is uniquely suited for studying structure and dynamics of membrane proteins in a membranous environment. The major challenge to obtain high quality solid-state NMR spectra of membrane proteins is sensitivity, due to limited quantities of labeled high-molecular-weight proteins. Here we demonstrate the incorporation of paramagnetic metal (Cu(2+)) ions, through either EDTA or a chelator lipid, into membrane protein samples for rapid data collection under fast magic-angle spinning (MAS) and low power (1)H decoupling. Spectral sensitivity of DsbB (20 kDa), an integral membrane protein, more than doubles in the same experimental time due to (1)H T(1) relaxation enhancement by Cu(2+) ions, with DsbB native fold and active site intact. This technique can be implemented to acquire multidimensional solid-state NMR spectra for chemical shift assignments and structure elucidation of large membrane proteins with small sample quantities.     

The development of an NMR chemical shift prediction application with the accuracy necessary to grade proton NMR spectra for identity

We have developed an NMR chemical shift prediction system that enables high throughput automatic grading of NMR spectra. In support of high throughput synthetic efforts for our drug discovery program, a rapid and accurate analysis for identity was needed. The system was designed and implemented to take advantage of the NMR assignments that had been tabulated on internally generated research compounds. The system has been operational for four years and has been used in conjunction with an internally written grading program to successfully analyze several hundred thousand samples based only on their 1D ¹H spectrum. A focused test of the system's accuracy on 1006 molecules demonstrated the ability to estimate the proton chemical shift with an average error of +/?0.16 ppm. This level of chemical shift accuracy allows for reliable structure confirmation by automated analysis using only proton NMR. Copyright © 2009 John Wiley & Sons, Ltd.