A polynomial-time nuclear vector replacement algorithm for automated NMR resonance assignments.

High-throughput NMR structural biology can play an important role in structural genomics. We report an automated procedure for high-throughput NMR resonance assignment for a protein of known structure, or of a homologous structure. These assignments are a prerequisite for probing protein-protein interactions, protein-ligand binding, and dynamics by NMR. Assignments are also the starting point for structure determination and refinement. A new algorithm, called Nuclear Vector Replacement (NVR) is introduced to compute assignments that optimally correlate experimentally measured NH residual dipolar couplings (RDCs) to a given a priori whole-protein 3D structural model. The algorithm requires only uniform( 15)N-labeling of the protein and processes unassigned H(N)-(15)N HSQC spectra, H(N)-(15)N RDCs, and sparse H(N)-H(N) NOE's (d(NN)s), all of which can be acquired in a fraction of the time needed to record the traditional suite of experiments used to perform resonance assignments. NVR runs in minutes and efficiently assigns the (H(N),(15)N) backbone resonances as well as the d(NN)s of the 3D (15)N-NOESY spectrum, in O(n(3)) time. The algorithm is demonstrated on NMR data from a 76-residue protein, human ubiquitin, matched to four structures, including one mutant (homolog), determined either by x-ray crystallography or by different NMR experiments (without RDCs). NVR achieves an assignment accuracy of 92-100%. We further demonstrate the feasibility of our algorithm for different and larger proteins, using NMR data for hen lysozyme (129 residues, 97-100% accuracy) and streptococcal protein G (56 residues, 100% accuracy), matched to a variety of 3D structural models. Finally, we extend NVR to a second application, 3D structural homology detection, and demonstrate that NVR is able to identify structural homologies between proteins with remote amino acid sequences using a database of structural models.     

13C NMR spectroscopy of labeled pyridoxal 5'-phosphate. Model studies, D-serine dehydratase, and L-glutamate decarboxylase.

Pyridoxal 5'-phosphate labeled to the extent of 90% with 13C in the 4' (aldehyde) and 5' (methylene) positions has been synthesized. 13C NMR spectra of this material and of natural abundance pyridoxal 5'-phosphate are reported, as well as 13C NMR spectra of the Schiff base formed by reaction of pyridoxal 5'-phosphate with n-butylamine, the secondary amine formed by reduction of this Schiff base, the thiazolidine formed by reaction of pyridoxal 5'-phosphate with cysteine, the hexahydropyrimidine formed by reaction of pyridoxal 5'-phosphate with 1,3-diaminobutane, and pyridoxamine 5'-phosphate. The range of chemical shifts for carbon 4' in these compounds is more than 100 ppm, and thus this chemical shift is expected to be a sensitive indicator of structure in enzyme-bound pyridoxal 5'-phosphate. The chemical shift of carbon 5', on the other hand, is insensitive to these structure changes. 13C NMR spectra have been obtained at pH 7.8 and 9.4 for D-serine dehydratase (Mr = 46,000) containing natural abundance pyridoxal 5'-phosphate and containing 13C-enriched pyridoxal 5'-phosphate. The enriched material contains two new resonances not present in the natural abundance material, one at 167.7 ppm with a linewidth of approximately 24 Hz, attributed to carbon 4' of the Schiff base in the bound coenzyme, and one at 62.7 Hz with a linewidth of approximately 48 Hz attributed to carbon 5' of the bound Schiff base. A large number of resonances due to individual amino acids are assigned. The NMR spectrum changes only slightly when the pH is raised to 9.4. The widths of the two enriched coenzyme resonances indicate that the coenzyme is rather rigidly bound to the enzyme but probably has limited motional freedom relative to the protein. 13C NMR spectra have been obtained for L-glutamate decarboxylase containing natural abundance pyridoxal 5'-phosphate and 13C-enriched pyridoxal 5'-phosphate. Under conditions where the two enriched 13C resonances are clearly visible in D-serine dehydratase, no resonances are visible in enriched L-glutamate decarboxylase, presumably because the coenzyme is rigidly bound to the protein and the 300,000 molecular weight of this enzyme produces very short relaxation times for the bound coenzyme and thus very broad lines.     

Structural and orientational constraints of bacteriorhodopsin in purple membranes determined by oriented-sample solid-state NMR spectroscopy

We report for the first time, oriented-sample solid-state NMR experiments, specifically polarization inversion spin exchange at the magic angle (PISEMA) and 1H-15N heteronuclear chemical shift correlation (HETCOR), applied to an integral seven-transmembrane protein, bacteriorhodopsin (bR), in natural membranes. The spectra of [15N]Met-bR revealed clearly distinguishable signals from the helical and loop regions. By deconvolution of the helix resonances, it was possible to establish constraints for some helix tilt angles. It was estimated that the extracellular section of helix B has a tilt of less than 5 degrees from the membrane normal, while the tilt of helix A was estimated to be 18-22 degrees , both of which are in agreement with most crystal structures. Comparison of the experimental PISEMA spectrum with simulated spectra based on crystal structures showed that PISEMA and HETCOR experiments are extremely sensitive to the polytopic protein structure, and the solid-state NMR spectra for membrane-embedded bR matched most favorably with the recent 1FBB electron crystallography structure. These results suggest that this approach has the potential to yield structural and orientational constraints for large integral polytopic proteins whilst intercalated and functionally competent in a natural membrane.     

Synthesis of Methyl (5Z,9Z,17R)-17-Methylnonadeca-5,9-dienoate,the (R)-Enantiomer of the Structure Proposed for a Metabolite of the Philippine Sponge Plakinastrella sp.

The (R)-enantiomer (1) of methyl (5Z,9Z)-17-methyl-nonadeca-5,9-dienoate, the structure proposed for a metabolite of the Philippine sponge, Plakinastrella sp., was synthesized. The ¹H- and ¹³C-NMR spectra of the synthetic material were different from those reported for the natural product. The proposed structure 1 is therefore incorrect.     

Development of a method for reconstruction of crowded NMR spectra from undersampled time-domain data

NMR is a unique methodology for obtaining information about the conformational dynamics of proteins in heterogeneous biomolecular systems. In various NMR methods, such as transferred cross-saturation, relaxation dispersion, and paramagnetic relaxation enhancement experiments, fast determination of the signal intensity ratios in the NMR spectra with high accuracy is required for analyses of targets with low yields and stabilities. However, conventional methods for the reconstruction of spectra from undersampled time-domain data, such as linear prediction, spectroscopy with integration of frequency and time domain, and analysis of Fourier, and compressed sensing were not effective for the accurate determination of the signal intensity ratios of the crowded two-dimensional spectra of proteins. Here, we developed an NMR spectra reconstruction method, “conservation of experimental data in analysis of Fourier” (Co-ANAFOR), to reconstruct the crowded spectra from the undersampled time-domain data. The number of sampling points required for the transferred cross-saturation experiments between membrane proteins, photosystem I and cytochrome b ₆ f, and their ligand, plastocyanin, with Co-ANAFOR was half of that needed for linear prediction, and the peak height reduction ratios of the spectra reconstructed from truncated time-domain data by Co-ANAFOR were more accurate than those reconstructed from non-uniformly sampled data by compressed sensing.     

Identification of metabolites from 2D 1H-13C HSQC NMR using peak correlation plots

Background

Identification of individual components in complex mixtures is an important and sometimes daunting task in several research areas like metabolomics and natural product studies. NMR spectroscopy is an excellent technique for analysis of mixtures of organic compounds and gives a detailed chemical fingerprint of most individual components above the detection limit. For the identification of individual metabolites in metabolomics, correlation or covariance between peaks in ¹H NMR spectra has previously been successfully employed. Similar correlation of 2D ¹H-¹³C Heteronuclear Single Quantum Correlation spectra was recently applied to investigate the structure of heparine. In this paper, we demonstrate how a similar approach can be used to identify metabolites in human biofluids (post-prostatic palpation urine).

Results

From 50 ¹H-¹³C Heteronuclear Single Quantum Correlation spectra, 23 correlation plots resembling pure metabolites were constructed. The identities of these metabolites were confirmed by comparing the correlation plots with reported NMR data, mostly from the Human Metabolome Database.

Conclusions

Correlation plots prepared by statistically correlating ¹H-¹³C Heteronuclear Single Quantum Correlation spectra from human biofluids provide unambiguous identification of metabolites. The correlation plots highlight cross-peaks belonging to each individual compound, not limited by long-range magnetization transfer as conventional NMR experiments.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-014-0413-z) contains supplementary material, which is available to authorized users.     

3-(4-Methyl-3-pentenyl)-2(5H)-furanone, alpha,alpha-acariolide and 4-(4-methyl-3-pentenyl)-2(5H)-furanone, alpha,beta-acariolide: new monoterpene lactones from the astigmatid mites, Schwiebea araujoae and Rhizoglyphus sp. (Astigmata: Acaridae)

A new monoterpene lactone from the acarid mite, Schwiebea araujoae, was elucidated without its isolation by GC/FT-IR and GC/MS analyses to be 3-(4-methyl-3-pentenyl)-2(5H)-furanone (1) and tentatively named as alpha,alpha-acariolide. The structure of 1 was identified by its synthesis from alpha-bromo-gamma-butyrolactone via 4 reaction steps. The synthesized compound gave the same GC/MS and GC/FT-IR spectra as those of the natural product. The other monoterpene lactone was likewise elucidated from the unidentified Rhizoglyphus mite to be 4-(4-methyl-3-pentenyl)-2(5H)-furanone (2) and named as alpha,beta-acariolide; it was also identified by its synthesis in 5 reaction steps from the same butyrolactone as the starting material. GC/MS and GC/FT-IR spectra of the preparation were identical to those of the natural product.     

Recent progress in heteronuclear long-range NMR of complex carbohydrates: 3D H2BC and clean HMBC

The new NMR experiments 3D H2BC and clean HMBC are explored for challenging applications to a complex carbohydrate at natural abundance of ¹³C. The 3D H2BC experiment is crucial for sequential assignment as it yields heteronuclear one- and two-bond together with COSY correlations for the ¹H spins, all in a single spectrum with good resolution and non-informative diagonal-type peaks suppressed. Clean HMBC is a remedy for the ubiquitous problem of strong coupling induced one-bond correlation artifacts in HMBC spectra of carbohydrates. Both experiments work well for one of the largest carbohydrates whose structure has been determined by NMR, not least due to the enhanced resolution offered by the third dimension in 3D H2BC and the improved spectral quality due to artifact suppression in clean HMBC. Hence these new experiments set the scene to take advantage of the sensitivity boost achieved by the latest generation of cold probes for NMR structure determination of even larger and more complex carbohydrates in solution.     

Synthesis of methyl (5Z,9Z,17R)-17-methylnonadeca-5,9-dienoate, the (R)-enantiomer of the structure proposed for a metabolite of the Philippine sponge Plakinastrella sp.

The (R)-enantiomer (1) of methyl (5Z,9Z)-17-methylnonadeca-5,9-dienoate, the structure proposed for a metabolite of the Philippine sponge, Plakinastrella sp., was synthesized. The 1H- and 13C-NMR spectra of the synthetic material were different from those reported for the natural product. The proposed structure 1 is therefore incorrect.     

RNAPack: an integrated NMR approach to RNA structure determination

Over the last decade, a vast number of useful nuclear magnetic resonance (NMR) experiments have been developed and successfully employed to determine the structure and dynamics of RNA oligonucleotides. Despite this progress, high-resolution RNA structure determination by NMR spectroscopy still remains a lengthy process and requires programming and extensive calibrations to perform NMR experiments successfully. To accelerate RNA structure determination by NMR spectroscopy, we have designed and programmed a package of RNA NMR experiments, called RNAPack. The user-friendly package contains a set of semiautomated single, double, and triple resonance NMR experiments, which are fully optimized for high-resolution RNA solution structure determination on Varian NMR spectrometers. RNAPack provides an autocalibration feature that allows rapid calibration of all NMR experiments in a single step and thereby speeds up the NMR data collection and eliminates user errors. In our laboratory, we have successfully employed this technology to solve RNA solution structures of domains of the internal ribosome entry site of the genomic hepatitis C viral RNA in less than 3 months. RNAPack therefore makes NMR spectroscopy an attractive and rapid structural tool and allows integration of atomic resolution structural information into biochemical studies of large RNA systems.     

An NMR method towards the routine chiral determination of natural products

State-of-the-art structure elucidation and dereplication of natural products is incomplete without the determination of enantiomeric purity, especially when compounds are to be biologically evaluated. An NMR procedure is presented in order to distinguish and determine enantiomers in natural product samples. The method is also of value in the structure elucidation process by providing information, which is otherwise of a non-routine nature. Using enantiomeric 1-acetoxychavicol acetates and carvones as model compounds, this study presents a chiral NMR procedure that allows distinction and determination of chiral antipodes of natural products in a routine set-up.     

The J-UNIO protocol for automated protein structure determination by NMR in solution

The J-UNIO (JCSG protocol using the software UNIO) procedure for automated protein structure determination by NMR in solution is introduced. In the present implementation, J-UNIO makes use of APSY-NMR spectroscopy, 3D heteronuclear-resolved [(1)H,(1)H]-NOESY experiments, and the software UNIO. Applications with proteins from the JCSG target list with sizes up to 150 residues showed that the procedure is highly robust and efficient. In all instances the correct polypeptide fold was obtained in the first round of automated data analysis and structure calculation. After interactive validation of the data obtained from the automated routine, the quality of the final structures was comparable to results from interactive structure determination. Special advantages are that the NMR data have been recorded with 6-10 days of instrument time per protein, that there is only a single step of chemical shift adjustments to relate the backbone signals in the APSY-NMR spectra with the corresponding backbone signals in the NOESY spectra, and that the NOE-based amino acid side chain chemical shift assignments are automatically focused on those residues that are heavily weighted in the structure calculation. The individual working steps of J-UNIO are illustrated with the structure determination of the protein YP_926445.1 from Shewanella amazonensis, and the results obtained with 17 JCSG targets are critically evaluated.     

The solution structure of a [3Fe-4S] ferredoxin: oxidised ferredoxin II from Desulfovibrio gigas

The use of standard 2D NMR experiments in combination with 1D NOE experiments allowed the assignment of 51 of the 58 spin systems of oxidised [3Fe-4S] ferredoxin isolated from Desulfovibrio gigas. The NMR solution structure was determined using data from 1D NOE and 2D NOESY spectra, as distance constraints, and information from the X-ray structure for the spin systems not detected by NMR in torsion angle dynamics calculations to produce a family of 15 low target function structures. The quality of the NMR family, as judged by the backbone r.m.s.d. values, was good (0.80?Å), with the majority of φ/ψ angles falling within the allowed region of the Ramachandran plot. A comparison with the X-ray structure indicated that the overall global fold is very similar in solution and in the solid state. The determination of the solution structure of ferredoxin II (FdII) in the oxidised state (FdII_ox) opens the way for the determination of the solution structure of the redox intermediate state of FdII (FdII_int), for which no X-ray structure is available.     

Attenuated total reflection IR spectroscopy as a tool to investigate the structure, orientation and tertiary structure changes in peptides and membrane proteins

During the last few years, attenuated total reflection Fourier transform infrared spectroscopy (ATR-FTIR) has become one of the most powerful methods to determine the structure of biological materials and in particular of components of biological membranes, like proteins that cannot be studied by x-ray crystallography and NMR. ATR-FTIR requires a little amount of material (1-100 microg) and spectra are recorded in a matter of minutes. The environment of the molecules can be modulated so that their conformation can be studied as a function of temperature, pressure, pH, as well as in the presence of specific ligands. For instance, replacement of amide hydrogen by deuterium is extremely sensitive to environmental changes and the kinetics of exchange can be used to detect tertiary conformational changes in the protein structure. Moreover, in addition to the conformational parameters that can be deduced from the shape of the infrared spectra, the orientation of various parts of the molecule can be estimated with polarized IR. This allows more precise analysis of the general architecture of the membrane molecules within the biological membranes. The present review focuses on ATR-IR as an experimental approach of special interest for the study of the structure, orientation, and tertiary structure changes in peptides and membrane proteins.     

VirtualSpectrum,a tool for simulating peak list for multi-dimensional NMR spectra

NMR spectroscopy is a widely used technique for characterizing the structure and dynamics of macromolecules. Often large amounts of NMR data are required to characterize the structure of proteins. To save valuable time and resources on data acquisition, simulated data is useful in the developmental phase, for data analysis, and for comparison with experimental data. However, existing tools for this purpose can be difficult to use, are sometimes specialized for certain types of molecules or spectra, or produce too idealized data. Here we present a fast, flexible and robust tool, VirtualSpectrum, for generating peak lists for most multi-dimensional NMR experiments for both liquid and solid state NMR. It is possible to tune the quality of the generated peak lists to include sources of artifacts from peak overlap, noise and missing signals. VirtualSpectrum uses an analytic expression to represent the spectrum and derive the peak positions, seamlessly handling overlap between signals. We demonstrate our tool by comparing simulated and experimental spectra for different multi-dimensional NMR spectra and analyzing systematically three cases where overlap between peaks is particularly relevant; solid state NMR data, liquid state NMR homonuclear ¹H and ¹⁵N-edited spectra, and 2D/3D heteronuclear correlation spectra of unstructured proteins. We analyze the impact of protein size and secondary structure on peak overlap and on the accuracy of structure determination based on data of different qualities simulated by VirtualSpectrum.     

Strategies for dereplication of natural compounds using high-resolution tandem mass spectrometry

Complete structural elucidation of natural products is commonly performed by nuclear magnetic resonance spectroscopy (NMR), but annotating compounds to most likely structures using high-resolution tandem mass spectrometry is a faster and feasible first step. The CASMI contest 2016 (Critical Assessment of Small Molecule Identification) provided spectra of eighteen compounds for the best manual structure identification in the natural products category. High resolution precursor and tandem mass spectra (MS/MS) were available to characterize the compounds. We used the Seven Golden Rules, Sirius2 and MS-FINDER software for determination of molecular formulas, and then we queried the formulas in different natural product databases including DNP, UNPD, ChemSpider and REAXYS to obtain molecular structures. We used different in-silico fragmentation tools including CFM-ID, CSI:FingerID and MS-FINDER to rank these compounds. Additional neutral losses and product ion peaks were manually investigated. This manual and time consuming approach allowed for the correct dereplication of thirteen of the eighteen natural products.     

Two-dimensional NMR studies of staphylococcal nuclease. 1. Sequence-specific assignments of hydrogen-1 signals and solution structure of the nuclease H124L-thymidine 3',5'-bisphosphate-Ca2+ ternary complex

Staphylococcal nuclease H124L is a recombinant protein produced in Escherichia coli whose sequence is identical with that of the nuclease produced by the V8 variant of Staphylococcus aureus. The enzyme-metal ion activator-nucleotide inhibitor ternary complex, nuclease H124L-thymidine 3',5'-bisphosphate-Ca2+, was investigated by two-dimensional (2D) NMR techniques. Efficient overproduction of the enzyme facilitated the production of random fractionally deuterated protein, which proved essential for detailed NMR analysis. 1H NMR spin systems were analyzed by conventional 2D 1H[1H] methods: COSY, relayed COSY, HOHAHA, and NOESY. Assignments obtained by 1H NMR experiments were confirmed and extended by 1H-13C and 1H-15N heteronuclear NMR experiments [Wang, J., Hinck, A. P., Loh, S. N., & Markley, J. L. (1990) Biochemistry (following paper in this issue)]. Spectra of the ternary complexes prepared with protein at natural abundance and at 50% random fractional deuteration provided the information needed for sequence-specific assignments of 121 of the 149 amino acid residues. Short- and intermediate-range NOE connectivities allowed the determination of secondary structural features of the ternary complex: three alpha-helical domains and three antiparallel beta-pleated sheets with several reverse turns. A number of nonsequential long-range HN-HN and H alpha-HN connectivities revealed additional information about the spatial arrangement of these secondary structural elements. The solution structure of this ternary complex shows a close correspondence to the crystal structure of the nuclease wt-thymidine 3',5'-bisphosphate-Ca2+ ternary complex [Cotton, F. A., Hazen, E. E., & Legg, M. J. (1979) Proc. Natl. Acad. Sci. U.S.A. 76, 2551-2555].     

Structure determination of membrane proteins in five easy pieces

Rotational Alignment (RA) solid-state NMR provides the basis for a general method for determining the structures of membrane proteins in phospholipid bilayers under physiological conditions. Membrane proteins are high priority targets for structure determination, and are challenging for existing experimental methods. Because membrane proteins reside in liquid crystalline phospholipid bilayer membranes it is important to study them in this type of environment. The RA solid-state NMR approach we have developed can be summarized in five steps, and incorporates methods of molecular biology, biochemistry, sample preparation, the implementation of NMR experiments, and structure calculations. It relies on solid-state NMR spectroscopy to obtain high-resolution spectra and residue-specific structural restraints for membrane proteins that undergo rotational diffusion around the membrane normal, but whose mobility is otherwise restricted by interactions with the membrane phospholipids. High resolution spectra of membrane proteins alone and in complex with other proteins and ligands set the stage for structure determination and functional studies of these proteins in their native, functional environment.     

What are the structures of phakellistatin 2 and phakellistatin 4?

Summary Phakellistatins 2 and 4 are cyclic peptides with cancer cell growth inhibitory properties isolated from the ocean marine spongesPhakellia carteri andPhakellia costata. To verify the proposed structures of natural phakellistatins 2 and 4, the given sequential structures were synthesized and their NMR spectra compared with the natural product. A completely different spectral pattern indicates that the structure of the natural compound should be different in both cases. The synthetic compounds according to the formula of phakellistatin 2 and 4 turned out to be inactive against several human cancer cell lines.     

Biochemical and NMR spectroscopy evidence for a new tertiary A-U base pair in lupin ribosomal 5 S RNA structure.

The new model for the tertiary structure of ribosomal 5 S rRNA from plants recently proposed by some of us has been already supported by RNase H digestions in the presence of complementary oligodeoxynucleotides. These results are confirmed now by the new biochemical and NMR spectroscopy data. Diethylpyrocarbonate (DEP) and monoperphthalic acid (MPA) are the reagents with the high specificity toward single-stranded adenosine residues. Our experiments clearly show that under native conditions adenosine 100 (A100) of lupin 5 S rRNA is not available for reaction toward these reagents. However under denaturing conditions this residue reacts with DEP and MPA. The detailed analysis of the lupin 5 S rRNA by NMR spectra provide the data on the specific interaction of A100-U53. Thus, we have seen for the first time the NMR signal due to the A100-U53 tertiary base pair, which as we believe, stabilizes interactions between loops B and E.