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.     

Stereospecific assignments of protein NMR resonances based on the tertiary structure and 2D/3D NOE data

In many cases of protein structure determination by NMR a high-quality structure is required. An important contribution to structural precision is stereospecific assignment of magnetically nonequivalent prochiral methylene and methyl groups, eliminating the need for introducing pseudoatoms and pseudoatom corrections in distance restraint lists. Here, we introduce the stereospecific assignment program that uses the resonance assignment, a preliminary 3D structure and 2D and/or 3D nuclear Overhauser effect spectroscopy peak lists for stereospecific assignment. For each prochiral group the algorithm automatically calculates a score for the two different stereospecific assignment possibilities, taking into account the presence and intensity of the nuclear Overhauser effect (NOE) peaks that are expected from the local environment of each prochiral group (i.e., the close neighbors). The performance of the algorithm has been tested and used on NMR data of alpha-helical and beta-sheet proteins using homology models and/or X-ray structures. The program produced no erroneous stereospecific assignments provided the NOEs were carefully picked and the 3D model was sufficiently accurate. The set of NOE distance restraints produced by nmr2st using the results of the SSA module was superior in generating good-quality ensembles of NMR structures (low deviations from upper limits in conjunction with low root-mean-square-deviation values) in the first round of structure calculations. The program uses a novel approach that employs the entire 3D structure of the protein to obtain stereospecific assignment; it can be used to speed up the NMR structure refinement and to increase the quality of the final NMR ensemble even when no scalar or residual dipolar coupling information is available.     

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.     

Deuterated protein folds obtained directly from unassigned nuclear overhauser effect data

We demonstrate the feasibility of determining the global fold of a highly deuterated protein from unassigned experimental NMR nuclear Overhauser effect (NOE) data only. The method relies on the calculation of a spatial configuration of covalently unconnected protons-a "cloud"-directly from unassigned distance restraints derived from 13C- and 15N-edited NOESY spectra. Each proton in the cloud, labeled by its chemical shift and that of the directly bound 13C or 15N, is subsequently mapped to specific atoms in the protein. This is achieved via graph-theoretical protocols that search for connectivities in graphs that encode the structural information within the cloud. The peptidyl HN chain is traced by seeking for all possible routes and selecting the one that yields the minimal sum of sequential distances. Complete proton identification in the cloud is achieved by linking the side-chain protons to proximal main-chain HNs via bipartite graph matching. The identified protons automatically yield the NOE assignments, which in turn are used for structure calculation with RosettaNMR, a protocol that incorporates structural bias derived from protein databases. The method, named Sparse-Constraint CLOUDS, was applied to experimental NOESY data on the 58-residue Z domain of staphylococcal protein A. The generated structures are of similar accuracy to those previously reported, which were derived via a conventional approach involving a larger NMR data set. Additional tests were performed on seven reported protein structures of various folds, using restraint lists simulated from the known atomic coordinates.     

Completely automated, highly error-tolerant macromolecular structure determination from multidimensional nuclear overhauser enhancement spectra and chemical shift assignments

The major rate-limiting step in high-throughput NMR protein structure determination involves the calculation of a reliable initial fold, the elimination of incorrect nuclear Overhauser enhancement (NOE) assignments, and the resolution of NOE assignment ambiguities. We present a robust approach to automatically calculate structures with a backbone coordinate accuracy of 1.0-1.5 A from datasets in which as much as 80% of the long-range NOE information (i.e., between residues separated by more than five positions in the sequence) is incorrect. The current algorithm differs from previously published methods in that it has been expressly designed to ensure that the results from successive cycles are not biased by the global fold of structures generated in preceding cycles. Consequently, the method is highly error tolerant and is not easily funnelled down an incorrect path in either three-dimensional structure or NOE assignment space. The algorithm incorporates three main features: a linear energy function representation of the NOE restraints to allow maximization of the number of simultaneously satisfied restraints during the course of simulated annealing; a method for handling the presence of multiple possible assignments for each NOE cross-peak which avoids local minima by treating each possible assignment as if it were an independent restraint; and a probabilistic method to permit both inactivation and reactivation of all NOE restraints on the fly during the course of simulated annealing. NOE restraints are never removed permanently, thereby significantly reducing the likelihood of becoming trapped in a false minimum of NOE assignment space. The effectiveness of the algorithm is demonstrated using completely automatically peak-picked experimental NOE data from two proteins: interleukin-4 (136 residues) and cyanovirin-N (101 residues). The limits of the method are explored using simulated data on the 56-residue B1 domain of Streptococcal protein G.     

Conformation and absolute configuration of nocathiacin I determined by NMR spectroscopy and chiral capillary electrophoresis

Nocathiacin I (BMS-249524) is a highly cross-linked thiazolyl peptide that displays potent activity against Gram-positive bacteria, including a number of antibiotic-resistant strains. This natural product contains 10 chiral centers. NMR studies have been performed to characterize the solution structure of nocathiacin I. A uniformly 13C,15N-labeled sample was used to obtain NMR assignments. Restrained simulated annealing calculations were performed by using accurately determined NOE distance restraints. All of the chiral centers were allowed to float during the simulated annealing protocol. Two clusters of structures were obtained that satisfy the NOE restraints very well and that are reasonably consistent with vicinal J-coupling constants. Within each cluster, all 10 chiral centers are uniquely defined. The two clusters are effectively mirror images of each other: all chiral centers that have the R(S) configuration in one cluster have the S(R) configuration in the other. The single threonine residue in nocathiacin I was subsequently determined to be l-threonine by chiral capillary electrophoresis, allowing the absolute configurations of all 10 chiral centers to be defined.     

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.     

Automated protein structure determination from NMR spectra

Fully automated structure determination of proteins in solution (FLYA) yields, without human intervention, three-dimensional protein structures starting from a set of multidimensional NMR spectra. Integrating existing and new software, automated peak picking over all spectra is followed by peak list filtering, the generation of an ensemble of initial chemical shift assignments, the determination of consensus chemical shift assignments for all (1)H, (13)C, and (15)N nuclei, the assignment of NOESY cross-peaks, the generation of distance restraints, and the calculation of the three-dimensional structure by torsion angle dynamics. The resulting, preliminary structure serves as additional input to the second stage of the procedure, in which a new ensemble of chemical shift assignments and a refined structure are calculated. The three-dimensional structures of three 12-16 kDa proteins computed with the FLYA algorithm coincided closely with the conventionally determined structures. Deviations were below 0.95 A for the backbone atom positions, excluding the flexible chain termini. 96-97% of all backbone and side-chain chemical shifts in the structured regions were assigned to the correct residues. The purely computational FLYA method is suitable for substituting all manual spectra analysis and thus overcomes a main efficiency limitation of the NMR method for protein structure determination.     

Protein-ligand NOE matching: a high-throughput method for binding pose evaluation that does not require protein NMR resonance assignments

Given the three-dimensional (3D) structure of a protein, the binding pose of a ligand can be determined using distance restraints derived from assigned intra-ligand and protein-ligand nuclear Overhauser effects (NOEs). A primary limitation of this approach is the need for resonance assignments of the ligand-bound protein. We have developed an approach that utilizes data from 3D 13C-edited, 13C/15N-filtered HSQC-NOESY spectra for evaluating ligand binding poses without requiring protein NMR resonance assignments. Only the 1H NMR assignments of the bound ligand are essential. Trial ligand binding poses are generated by any suitable method (e.g., computational docking). For each trial binding pose, the 3D 13C-edited, 13C/15N-filtered HSQC-NOESY spectrum is predicted, and the predicted and observed patterns of protein-ligand NOEs are matched and scored using a fast, deterministic bipartite graph matching algorithm. The best scoring (lowest "cost") poses are identified. Our method can incorporate any explicit restraints or protein assignment data that are available, and many extensions of the basic procedure are feasible. Only a single sample is required, and the method can be applied to both slowly and rapidly exchanging ligands. The method was applied to three test cases: one complex involving muscle fatty acid-binding protein (mFABP) and two complexes involving the leukocyte function-associated antigen 1 (LFA-1) I-domain. Without using experimental protein NMR assignments, the method identified the known binding poses with good accuracy. The addition of experimental protein NMR assignments improves the results. Our "NOE matching" approach is expected to be widely applicable; i.e., it does not appear to depend on a fortuitous distribution of binding pocket residues.     

Impact of 15N R2/R1 relaxation restraints on molecular size, shape, and bond vector orientation for NMR protein structure determination with sparse distance restraints

(15)N R(2)/R(1) relaxation data contain information on molecular shape and size as well as on bond vector orientations relative to the diffusion tensor. Since the diffusion tensor can be directly calculated from the molecular coordinates, direct inclusion of (15)N R(2)/R(1) restraints in NMR structure calculations without any a priori assumptions is possible. Here we show that (15)N R(2)/R(1) restraints are particularly valuable when only sparse distance restraints are available. Using three examples of proteins of varying size, namely, GB3 (56 residues), ubiquitin (76 residues), and the N-terminal domain of enzyme I (EIN, 249 residues), we show that incorporation of (15)N R(2)/R(1) restraints results in large and significant increases in coordinate accuracy that can make the difference between being able or unable to determine an approximate global fold. For GB3 and ubiquitin, good coordinate accuracy was obtained using only backbone hydrogen-bond restraints supplemented by (15)N R(2)/R(1) relaxation restraints. For EIN, the global fold could be determined using sparse nuclear Overhauser enhancement (NOE) distance restraints involving only NH and methyl groups in conjunction with (15)N R(2)/R(1) restraints. These results are of practical significance in the study of larger and more complex systems, where the increasing spectral complexity and number of chemical shift degeneracies reduce the number of unambiguous NOE assignments that can be readily obtained, resulting in progressively reduced NOE coverage as the size of the protein increases.     

Structures of protein-protein complexes are docked using only NMR restraints from residual dipolar coupling and chemical shift perturbations

NMR structures of protein-protein and protein-ligand complexes rely heavily on intermolecular NOEs. Recent work has shown that if no significant conformational changes occur upon complex formation residual dipolar coupling can replace most of the NOE restraints in protein-protein complexes, while restraints derived from chemical shift perturbations can largely replace intermolecular NOEs in protein-ligand structures. By combining restraints from chemical shift perturbations with orientation restraints derived from measurements of residual dipolar couplings, we show that the structure of the EIN-HPr complex can be calculated without NOE restraints. The final structure, built from the crystal structures of EIN and HPr in their uncomplexed form and docked only with NMR restraints, places HPr within 2.5 A of the position determined from the mean NMR structure of the complex.     

Protein NMR recall, precision, and F-measure scores (RPF scores): structure quality assessment measures based on information retrieval statistics

One of the most important challenges in modern protein NMR is the development of fast and sensitive structure quality assessment measures that can be used to evaluate the "goodness-of-fit" of the 3D structure with NOESY data, to indicate the correctness of the fold and accuracy of the resulting structure. Quality assessment is especially critical for automated NOESY interpretation and structure determination approaches. This paper describes new NMR quality assessment scores, including Recall, Precision, and F-measure scores (referred to here are "NMR RPF" scores), which quickly provide global measures of the goodness-of-fit of the 3D structures with NOESY peak lists using methods from information retrieval statistics. The sensitivity of the F-measure is improved using a scaled Fold Discriminating Power (DP) score. These statistical RPF scores are quite rapid to compute since NOE assignments and complete relaxation matrix calculations are not required. A graphical method for site-specific assessment of structure quality based on the Precision statistic is also described. These statistical measures are demonstrated to be valuable for assessing protein NMR structure accuracy. Their relationships to other proposed NMR "R-factors" and structure quality assessment scores are also discussed.     

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.     

Solution and solid state13C-NMR of barley straw lignins

Grass lignins are formed by the polymerization of phenoxy radicals and contain a variety of carbon-carbon and carbon-oxygen bonds. They are similar to the hardwood lignins, but differ by containing a substantial proportion of esterified cinnamic acids. Detailed nuclear magnetic resonance studies in conjunction with chemical analysis have given new information on the structure of grass lignins. Milled straw lignins (MSL) from barley were examined by both solution and solid-state (CP/MAS) NMR before and after acetylation. The assignment of the carbon-13 (100 MHz) solution spectra was achieved using model compound data, nuclear Overhauser enhancement (NOE) suppression, and insensitive nuclei enhanced by polarization transfer (INEPT) techniques. The NOE suppression permitted quantitative analysis of lignin, giving information on the ratio of specific carbon atoms. Use of the relaxation agent, chromium acetylacetonate, enabled accumulation of sufficient spectral data to give a spectrum suitable for integration after 90 h. The INEPT technique, which had not previously been used for lignin analysis, was successfully applied to acetylated MSL. This technique increased signal intensities 3–4-fold and simplified the spectrum by inverting methylene carbons and eliminating or inverting quaternary carbons. Comparison of this spectrum with the normal spectrum permitted accurate assignment of quaternary and methine carbons. The solid-state carbon-13 CP/MAS NMR was used to examinein situ lignin and the isolated MSL. The¹³C-CP/MAS spectrum ofin situ lignin shows that cellulose and hemicellulose resonances dominate with little evidence ot the aromatic structure of lignin. the¹³C-CP/MAS of MSL shows reduced carbohydrate resonances and increased aromatic resonances. The extent of modification to the barley straw was estimated and results indicate the presence oflignin-carbohydrate complexes. Detailed information on the nature of the linkages between lignin components and between lignin and carbohydrate components has been obtained from these spectra.     

Structure determination by restrained molecular dynamics using NMR pseudocontact shifts as experimentally determined constraints

The structure of a DNA octamer d(TTGGCCAA)(2) complexed to chromomycin-A(3) and a single divalent cobalt ion has been solved by using the pseudocontact shifts due to the unpaired electrons on the cobalt. A protocol was developed and critically evaluated for using the pseudocontact shifts in structure determination. The pseudocontact shifts were input as experimental restraints in molecular dynamics simulations with or without NOE constraints. Both the magnitude and orientation of the susceptibility anisotropy tensor required for the shift calculations were determined during the simulations by iterative refinement. The pseudocontact shifts could be used to define the structure to a very high precision and accuracy compared with a corresponding NOE-determined structure. Convergence was obtained from different starting structures and tensors. A structure determination using both NOE's and pseudocontact shifts revealed a general agreement between the two data sets. However, some evidence for a discrepancy between NOE's and pseudocontact shifts was observed in the backbone and terminal base pairs of the DNA. Violations in shift or NOE restraints remaining in the final structures were examined and may be a reflection of motional averaging of the constraints and evidence for flexibility. This work demonstrates that pseudocontact shifts are a powerful tool for NMR structure determination.     

Site-directed parallel spin-labeling and paramagnetic relaxation enhancement in structure determination of membrane proteins by solution NMR spectroscopy

A major challenge for the structure determination of integral membrane proteins by solution NMR spectroscopy is the limited number of NOE restraints in these systems stemming from extensive deuteration. Paramagnetic relaxation enhancement (PRE) by means of nitroxide spin-labels can provide valuable long-range distance information but, in practice, has limits in its application to membrane proteins because spin-labels are often incompletely reduced in highly apolar environments. Using the integral membrane protein OmpA as a model system, we introduce a method of parallel spin-labeling with paramagnetic and diamagnetic labels and show that distances in the range 15-24 Angstroms can be readily determined. The protein was labeled at 11 water-exposed and lipid-covered sites, and 320 PRE distance restraints were measured. The addition of these restraints resulted in significant improvement of the calculated backbone structure of OmpA. Structures of reasonable quality can even be calculated with PRE distance restraints only, i.e., in the absence of NOE distance restraints.     

Three‐dimensional structure of cyclic antibiotic teicoplanin aglycone using NMR distance and dihedral angle restraints in a DMSO solvation model

The three‐dimensional solution conformation of teicoplanin aglycone was determined using NMR spectroscopy. A combination of NOE and dihedral angle restraints in a DMSO solvation model was used to calculate an ensemble of structures having a root mean square deviation of 0.17 Å. The structures were generated using systematic searches of conformational space for optimal satisfaction of distance and dihedral angle restraints. Comparison of the NMR‐derived structure of teicoplanin aglycone with the X‐ray structure of a teicoplanin aglycone analog revealed a common backbone conformation with deviation of two aromatic side chain substituents. Experimentally determined backbone ¹³C chemical shifts showed good agreement with those computed at the density functional level of theory, providing a cross validation of the backbone conformation. The flexible portion of the molecule was consistent with the region that changes conformation to accommodate protein binding. The results showed that a hydrogen‐bonded DMSO molecule in combination with NMR‐derived restraints together enabled calculation of structures that satisfied experimental data. Copyright © 2015 John Wiley & Sons, Ltd.     

The occurrence of head-to-head butene units in ethylene-butene copolymers and corresponding ethyl branches in low-density (branched) polyethylenes

A high resolution carbon-13 NMR study of an ethylene-butene copolymer has yielded a spectrum whose resonances could be assigned to isolated ethyl branches, 1,3-diethyl branching as well as 1,2-diethyl branches resulting from head-to-head butene polymerization. All these structures are present in the same sample. An examination of all the published spectra of low-density (branched) polyethylenes has revealed many examples of heretofore unassigned resonances which can be associated with this latter type of diethyl branching in these polymers. Spin-lattice relaxation times have been determined, when feasible, for the carbon atoms in each of the three branched structures.     

Fully automatic assignment of small molecules' NMR spectra without relying on chemical shift predictions

We present a method for the automatic assignment of small molecules' NMR spectra. The method includes an automatic and novel self‐consistent peak‐picking routine that validates NMR peaks in each spectrum against peaks in the same or other spectra that are due to the same resonances. The auto‐assignment routine used is based on branch‐and‐bound optimization and relies predominantly on integration and correlation data; chemical shift information may be included when available to fasten the search and shorten the list of viable assignments, but in most cases tested, it is not required in order to find the correct assignment. This automatic assignment method is implemented as a web‐based tool that runs without any user input other than the acquired spectra. Copyright © 2015 John Wiley & Sons, Ltd.