Gifa V. 4: A complete package for NMR data set processing

Summary The Gifa program is designed for processing, displaying and analysing 1D, 2D and 3D NMR data sets. It has been constructed in a modular fashion, based on three independent modules: a set of commands that perform all the basic processing operations such as apodisation functions, a complete set of Fourier Transforms, phasing and baseline correction, peak-picking and line fitting, linear prediction and maximum entropy processing; a set of command language primitives that permit the execution of complex macro commands; and a set of graphic commands that permit to build a complete graphic user interface, allowing the user to interact easily with the program. We have tried to create a versatile program that can be easily extended according to the user's requirements and that is adapted to a novice as well as an experienced user. The program runs on any UNIX computer, with or without graphic display, in interactive or batch mode.     

Anatomic Twist to a Straightforward Ablation

Atrioventricular (AV) junction ablation for treatment of refractory atrial fibrillation is a well defined, standardized procedure and the simplest of commonly performed radiofrequency ablations in the field of cardiac electrophysiology. We report successful AV junction ablation using an inferior approach in a case of inferior vena cava interruption. Inability during the procedure to initially pass the ablation catheter into the right ventricle, combined with low amplitude electrograms, led to suspicion of an anatomic abnormality. This was determined to be a heterotaxy syndrome with inferior vena cava interruption and azygos continuation, draining in turn into the superior vena cava. Advancing Schwartz right 0 (SRO) sheath through the venous abnormality into the right atrium allowed adequate catheter stability to successfully induce complete AV block with radiofrequency energy.     

Analytical solution to the coupled evolution of multidimensional NMR data

A substantial time savings in the collection of multidimensional NMR data can be achieved by coupling the evolution of nuclei in the indirect dimensions. In order to save time, the sampling of the indirect dimensions is inherently incomplete. Therefore, many algorithms and samplings schemes have been developed aimed at separating the coevolved frequencies into analyzable data with limited artifacts. This paper extends the use of circulant matrices to describe coupled evolution with convolutions. By understanding the data in terms of convolutions, there is an exact solution to the inversion problem of extracting the orthogonal vectors from the coupled dimensions. Previously, this inversion problem has been solved using peak coordinates extracted from spectra. In contrast, the method described here uses spectra directly. This solution suggests a simple sampling scheme of collecting N orthogonal spectra, and N + 1 projections at specific projection angles, however, the theory developed can be extended generally to arbitrary projection angles. The circulant matrix methodology is demonstrated for simulated and real data. Further, an algorithm for separating overlapped signals in the detected dimension is presented. The algorithm involves the forward calculation of the coupled spectra from the orthogonal spectra, followed by back calculation of the orthogonal spectra from the coupled spectra, thus permitting rigorous cross-validation. This algorithm is shown to be robust in that erroneous solutions give rise to large artifacts. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

A nomenclature and data model to describe NMR experiments

Despite ongoing efforts in organising NMR information, there is no consistent and well-described generic standard for naming NMR experiments. The main reason for the absence of a universal naming system is that the information content of the coherence pathways is difficult to describe in full detail. We propose a system that describes the common and generic elements of the coherence pathways produced by pulse sequences. The system itself is formalised by an ‘NMR experiment protocol’ model, which is described in the Universal Modelling Language (UML) as part of the CCPN data model. Furthermore, normalized experiment names can be derived from this proposed model. We hope this article will stimulate discussion to organise the wealth of NMR experiments, and that by bringing this discussion into the public domain we can improve and expand our proposed system to include as much information and as many NMR experiments as possible. Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users.Joint first authors     

An automated system designed for large scale NMR data deposition and annotation: application to over 600 assigned chemical shift data entries to the BioMagResBank from the Riken Structural Genomics/Proteomics Initiative internal database

Biomolecular NMR chemical shift data are key information for the functional analysis of biomolecules and the development of new techniques for NMR studies utilizing chemical shift statistical information. Structural genomics projects are major contributors to the accumulation of protein chemical shift information. The management of the large quantities of NMR data generated by each project in a local database and the transfer of the data to the public databases are still formidable tasks because of the complicated nature of NMR data. Here we report an automated and efficient system developed for the deposition and annotation of a large number of data sets including (1)H, (13)C and (15)N resonance assignments used for the structure determination of proteins. We have demonstrated the feasibility of our system by applying it to over 600 entries from the internal database generated by the RIKEN Structural Genomics/Proteomics Initiative (RSGI) to the public database, BioMagResBank (BMRB). We have assessed the quality of the deposited chemical shifts by comparing them with those predicted from the PDB coordinate entry for the corresponding protein. The same comparison for other matched BMRB/PDB entries deposited from 2001-2011 has been carried out and the results suggest that the RSGI entries greatly improved the quality of the BMRB database. Since the entries include chemical shifts acquired under strikingly similar experimental conditions, these NMR data can be expected to be a promising resource to improve current technologies as well as to develop new NMR methods for protein studies.     

NMR View: A computer program for the visualization and analysis of NMR data

Summary NMR View is a computer program designed for the visualization and analysis of NMR data. It allows the user to interact with a practically unlimited number of 2D, 3D and 4D NMR data files. Any number of spectral windows can be displayed on the screen in any size and location. Automatic peak picking and facilitated peak analysis features are included to aid in the assignment of complex NMR spectra. NMR View provides structure analysis features and data transfer to and from structure generation programs, allowing for a tight coupling between spectral analysis and structure generation. Visual correlation between structures and spectra can be done with the Molecular Data Viewer, a molecular graphics program with bidirectional communication to NMR View. The user interface can be customized and a command language is provided to allow for the automation of various tasks.Inquiries concerning the availability of NMR View and the Molecular Data Viewer should be sent via email to johnsonb@merck.com or to Bruce A. Johnson, Merck Research Laboratories, RY80Y-103, P.O. Box 2000, Rahway, NJ 07065, U.S.A.     

EZ-ASSIGN, a program for exhaustive NMR chemical shift assignments of large proteins from complete or incomplete triple-resonance data

For several of the proteins in the BioMagResBank larger than 200 residues, 60 % or fewer of the backbone resonances were assigned. But how reliable are those assignments? In contrast to complete assignments, where it is possible to check whether every triple-resonance Generalized Spin System (GSS) is assigned once and only once, with incomplete data one should compare all possible assignments and pick the best one. But that is not feasible: For example, for 200 residues and an incomplete set of 100 GSS, there are 1.6 × 10²⁶⁰ possible assignments. In “EZ-ASSIGN”, the protein sequence is divided in smaller unique fragments. Combined with intelligent search approaches, an exhaustive comparison of all possible assignments is now feasible using a laptop computer. The program was tested with experimental data of a 388-residue domain of the Hsp70 chaperone protein DnaK and for a 351-residue domain of a type III secretion ATPase. EZ-ASSIGN reproduced the hand assignments. It did slightly better than the computer program PINE (Bahrami et al. in PLoS Comput Biol 5(3):e1000307, 2009) and significantly outperformed SAGA (Crippen et al. in J Biomol NMR 46:281–298, 2010), AUTOASSIGN (Zimmerman et al. in J Mol Biol 269:592–610, 1997), and IBIS (Hyberts and Wagner in J Biomol NMR 26:335–344, 2003). Next, EZ-ASSIGN was used to investigate how well NMR data of decreasing completeness can be assigned. We found that the program could confidently assign fragments in very incomplete data. Here, EZ-ASSIGN dramatically outperformed all the other assignment programs tested.     

Improved methods for the acquisition and interpretation of NMR metabolomic data

The high spectral congestion typically observed in one-dimensional (1D) 1H nuclear magnetic resonance (NMR) spectra of tissue extracts and biofluids limits the metabolic information that can be extracted. This study evaluates the application of two-dimensional J-resolved (JRES) spectroscopy for metabolomics, which can provide proton-decoupled projected 1D spectra (p-JRES). This approach is illustrated by an investigation of embryogenesis in Japanese medaka (Oryzias latipes), an established fish model for developmental toxicology. When combined with optimized spectral pre-processing,(2) including a 0.005-ppm bin width for data segmentation and a logarithmic transformation, the reduced congestion in the p-JRES spectra increases the likelihood that a specific metabolite can be accurately integrated and thus increases the extractable information content of the spectra. Principal components analysis of the p-JRES spectra reveals the concept of a developmental trajectory that summarizes the changes in the NMR-visible metabolome throughout medaka embryogenesis. Advantages and potential disadvantages of the p-JRES approach are discussed.     

Estimating the number of protein folds and families from complete genome data

Using the data on proteins encoded in complete genomes, combined with a rigorous theory of the sampling process, we estimate the total number of protein folds and families, as well as the number of folds and families in each genome. The total number of folds in globular, water- soluble proteins is estimated at about 1000, with structural information currently available for about one-third of the number. The sequenced genomes of unicellular organisms encode from approximately 25%, for the minimal genomes of the Mycoplasmas, to 70-80% for larger genomes, such as Escherichia coli and yeast, of the total number of folds. The number of protein families with significant sequence conservation was estimated to be between 4000 and 7000, with structures available for about 20% of these.     

Rapid prediction of multi-dimensional NMR data sets

We present a computational environment for Fast Analysis of multidimensional NMR DAta Sets (FANDAS) that allows assembling multidimensional data sets from a variety of input parameters and facilitates comparing and modifying such ??in silico?? data sets during the various stages of the NMR data analysis. The input parameters can vary from (partial) NMR assignments directly obtained from experiments to values retrieved from in silico prediction programs. The resulting predicted data sets enable a rapid evaluation of sample labeling in light of spectral resolution and structural content, using standard NMR software such as Sparky. In addition, direct comparison to experimental data sets can be used to validate NMR assignments, distinguish different molecular components, refine structural models or other parameters derived from NMR data. The method is demonstrated in the context of solid-state NMR data obtained for the cyclic nucleotide binding domain of a bacterial cyclic nucleotide-gated channel and on membrane-embedded sensory rhodopsin II. FANDAS is freely available as web portal under WeNMR (http://www.wenmr.eu/services/FANDAS).     

Surface fractality of proteins from theory and NMR data

Different approaches to study protein surface fractality are considered. An approach based on analysis of surface versus molecular weight dependence is shown to be an informative tool for investigation of protein surface behaviour. An evidence for protein surface fractality, obtained with the use of this analysis from the data of both NMR measurements in protein solutions and computer analysis of protein structures, is presented. Obtained value of fractal dimension of protein surface (ds congruent to 2.2) is in a good agreement with the results of conventional approach (with variation of yardstick length) to protein surface fractality. A conclusion is made that surface enlargement due to the rise of protein molecular weight is accompanied by the increase of maximum scale of irregularities on protein surface. Possible effect of surface fractality on hydrodynamic characteristics of protein molecules in solution is discussed.     

Molecular mechanics and dynamics of an abasic frameshift in DNA and comparison to NMR data

In a previous publication (Ph. Cuniasse, L.C. Sowers, R. Eritja, B. Kaplan, M.F. Goodman, J.A.H. Cognet, M. Le Bret, W. Guschlbauer and G.V. Fazakerley, Biochemistry 28, 2018 (1989), we determined by two dimensional NMR studies and molecular mechanics calculations the three-dimensional structure of a non-selfcomplementary oligonucleotide: [sequence; see text] where dr, at the center of the first strand, is a model abasic site. In order to explain all the results arising from NMR measurements, we found that an equilibrium between two conformations was necessary. These conformations differ mainly by the sugar pucker of G5 which is C2' endo or C3' endo. The latter is stabilized by addition of counterions between phosphate residues P3 and P4. In this paper, we have constructed systematically, all possible structures as a function of torsion angles delta of dr4 and of G5 by molecular mechanics in the presence or absence of counterions. Since these conformations were not forced with NMR distance measurements, this method allows detailed comparisons between all possible conformations and NMR data. Maps of contour lines of the potential energy, of fits to NMR distance measurements, and of helical twist as a function of torsion angles delta of dr4 and of G5 unravel the difficulties associated with the study of the G5 sugar pucker conformation equilibrium. Sugar puckers and proton distances are very sensitive criteria to monitor molecular dynamics. Relying on these experimental criteria, we have tested many molecular dynamics preparation phases and we propose a new warm-up and equilibration procedure for molecular dynamics. Thus we show with a 290 ps molecular dynamic run that G5 is in conformational equilibrium and that all NMR data are well reproduced.     

MUNIN: Application of three-way decomposition to the analysis of heteronuclear NMR relaxation data**

MUNIN (Multidimensional NMR Spectra Interpretation), a recently introduced approach exploiting the mathematical concept of three-way decomposition, is proposed for separation and quantitative relaxation measurements of strongly overlapped resonances in sets of heteronuclear two-dimensional spectra that result from typical relaxation experiments. The approach is general and may also be applied to sets of two-dimensional spectra with arbitrary modulation along the third dimension (e.g., J-coupling, diffusion). Here, the method is applied for the analysis of ¹⁵N rotating frame relaxation data.     

KUJIRA,a package of integrated modules for systematic and interactive analysis of NMR data directed to high-throughput NMR structure studies

The recent expansion of structural genomics has increased the demands for quick and accurate protein structure determination by NMR spectroscopy. The conventional strategy without an automated protocol can no longer satisfy the needs of high-throughput application to a large number of proteins, with each data set including many NMR spectra, chemical shifts, NOE assignments, and calculated structures. We have developed the new software KUJIRA, a package of integrated modules for the systematic and interactive analysis of NMR data, which is designed to reduce the tediousness of organizing and manipulating a large number of NMR data sets. In combination with CYANA, the program for automated NOE assignment and structure determination, we have established a robust and highly optimized strategy for comprehensive protein structure analysis. An application of KUJIRA in accordance with our new strategy was carried out by a non-expert in NMR structure analysis, demonstrating that the accurate assignment of the chemical shifts and a high-quality structure of a small protein can be completed in a few weeks. The high completeness of the chemical shift assignment and the NOE assignment achieved by the systematic analysis using KUJIRA and CYANA led, in practice, to increased reliability of the determined structure.     

SESAME: A least-squares approach to the evaluation of protein structures computed from NMR data

Summary A method is proposed for defining a probability distribution on an ensemble of protein conformations from a 2D NOE spectrum, while at the same time back-calculating the experimental spectrum from the ensemble. This enables one to assess the relative quality and significance of the conformations, and to test the consistency of the ensemble as a whole with the experimental spectrum. The method eliminates the need to integrate the cross-peak intensities and is surprisingly insensitive to random noise in the spectrum. In this communication, these advantages are demonstrated by applying the method to simulated data, for which the correct result is already known.Abbreviations NOE nuclear Overhauser effect - BPTI basic pancreatic trypsin inhibitor - rmsd root-mean-square deviation