Dual transformation of homonuclear solid-state NMR spectra—an option to decrease measuring time

Long measurement times due to low sensitivity are a prime concern in solid-state NMR and limit the application of multidimensional experiments severely. One possibility to address this problem could be post-experimental suppression of noise and a reduction of the number of increments needed for higher dimensional data sets. This can be achieved by a hybrid approach based on the combination of separately Fourier transformed and covariance processed datasets. The method is applied to synthetic sets as well as to experimental two-dimensional homonuclear solid-state NMR spectra of peptide samples. It is demonstrated that a reduction in experiment time by a factor of 4 can be achieved for the case of a ¹³C-¹³C correlation spectrum on the nonapeptide bradykinin.     

Site-specific ¹⁹F NMR chemical shift and side chain relaxation analysis of a membrane protein labeled with an unnatural amino acid

Site‐specific ¹⁹F chemical shift and side chain relaxation analysis can be applied on large size proteins. Here, one‐dimensional ¹⁹F spectra and T₁, T₂ relaxation data were acquired on a SH3 domain in aqueous buffer containing 60% glycerol, and a nine‐transmembrane helices membrane protein diacyl‐glycerol kinase (DAGK) in dodecyl phosphochoine (DPC) micelles. The high quality of the data indicates that this method can be applied to site‐specifically analyze side chain internal mobility of membrane proteins or large size proteins.     

Classification of Brain Tumors by Ex Vivo 1H NMR Spectroscopy

Abstract: Ex vivo biopsy samples (n = 42) from human brain tumors and normal brain have been examined by high-resolution proton magnetic resonance spectroscopy. Parameters from one-dimensional ¹H spectra, two-dimensional COSY spectra, and transverse relaxation time (T₂) data were used to classify the tumors according to the histopathological diagnoses. The ratio of the area between 3.4 and 3.1 ppm to that between 1.5 and 1.1 ppm distinguished glioblastomas from astrocytomas and normal brain, and appeared to be indicative of malignant potential. In support of the one-dimensional data, cross-peaks in the COSY spectra of brain specimens classified glioblastomas and metastases into one group and the more benign tumors, meningiomas, astrocytomas, and normal brain into a second group. The transverse relaxation of the resonance at 1.3 ppm was fitted by a model with two T₂ values. The longer T₂ value could be used to distinguish glioblastomas from normal brain, the latter having a much longer long T₂ value. Astrocytomas showed a continuum of T₂ values between glioblastomas and normal brain, with the grade of the astrocytoma correlating roughly with the value of the long T₂ component.     

Automated solvent artifact removal and base plane correction of multidimensional NMR protein spectra by AUREMOL-SSA

Strong solvent signals lead to a disappearance of weak protein signals close to the solvent resonance frequency and to base plane variations all over the spectrum. AUREMOL-SSA provides an automated approach for solvent artifact removal from multidimensional NMR protein spectra. Its core algorithm is based on singular spectrum analysis (SSA) in the time domain and is combined with an automated base plane correction in the frequency domain. The performance of the method has been tested on synthetic and experimental spectra including two-dimensional NOESY and TOCSY spectra and a three-dimensional ¹H,¹³C-HCCH-TOCSY spectrum. It can also be applied to frequency domain spectra since an optional inverse Fourier transformation is included in the algorithm.     

Spin-state selection for increased confidence in cross-correlation rates measurements

A new approach is described for measuring chemical shift anisotropy (CSA)/dipolar cross-correlated relaxation (CCR) rates based on the selection of the individual ¹⁵N doublet components prior to the relaxation period. The method uses the spin-state-selective element (S³E) of Sørensen and co-authors [Meissner et al. (1997) J. Mag. Reson., 128, 92–97]. The main advantage of the new method compared to other J-resolved experiments is that it does not create problems of additional signal overlap encountered in coupled spectra. At the same time, this approach allows a simpler control of magnetization pathways than the indirect methods. The method is demonstrated for the B3 domain of protein G.Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1007/s10858-004-7562-8     

Fully automated sequence-specific resonance assignments of hetero- nuclear protein spectra

Full automation of the analysis of spectra is a prerequisite for high-throughput NMR studies in structural or functional genomics. Sequence-specific assignments often form the major bottleneck. Here, we present a procedure that yields nearly complete backbone and side chain resonance assignments starting from a set of heteronuclear three-dimensional spectra. Neither manual intervention, e.g., to correct lists obtained from peak picking before feeding these to an assignment program, nor protein-specific information, e.g., structures of homologous proteins, were required. By combining two earlier published procedures, AUTOPSY [Koradi et al. (1998) J. Magn. Reson., 135, 288–297] and GARANT [Bartels et al. (1996) J. Biomol. NMR, 7, 207–213], with a new program, PICS, all necessary steps from spectra analyses to sequence-specific assignments were performed fully automatically. Characteristic features of the present approach are a flexible design allowing as input almost any combination of NMR spectra, applicability to side chains, robustness with respect to parameter choices (such as noise levels) and reproducibility. In this study, automated resonance assignments were obtained for the 14 kD blue copper protein azurin from P. aeruginosa using five spectra: HNCACB, HNHA, HCCH-TOCSY, ¹⁵N-NOESY-HSQC and ¹³C-NOESY-HSQC. Peaks from these three-dimensional spectra were filtered and calibrated with the help of two two-dimensional spectra: ¹⁵N-HSQC and ¹³C-HSQC. The rate of incorrect assignments is less than 1.5% for backbone nuclei and about 3.5% when side chain protons are also considered.     

AUREMOL-RFAC-3D,combination of R-factors and their use for automated quality assessment of protein solution structures

We present here the computer program AUREMOL-RFAC-3D that is a generalization of the previously published program RFAC for the fully automated estimation of residual indices (R-factors) from 2D NOESY spectra. It is part of the larger AUREMOL software package (www.auremol.de). RFAC-3D calculates R-factors directly from two-dimensional homonuclear NOESY spectra as well as from three-dimensional ¹⁵N or ¹³C edited NOESY-HSQC spectra and thus extends the application range to larger proteins. The fully automated method includes automated peak picking and integration, a Bayesian noise and artifact recognition and the use of the complete relaxation matrix formalism. To enhance the reliability of the calculated R-factors the method is also generalized to calculate combined R-factors from a set of 2D and 3D-spectra. For an optimal combination of the information derived from different sources a plausible formalism had to be derived. In addition, we present a novel direct R-factors based measure that correlates an R-factors as defined in this paper to the root mean square deviation of the actual structure from the optimal structure. The new program has been successfully tested on the histidine containing phosphocarrier protein (HPr) from Staphylococcus carnosus and on the Ras-binding domain (RBD) of the Ral guanine-nucleotide dissociation stimulation factor (RalGDS).     

Direct measurement of the 15N CSA/dipolar relaxation interference from coupled HSQC spectra

Here we propose a method for the measurement of the ¹⁵N CSA/dipolar relaxation interference based on direct comparison of the ¹⁵N doublet components observed in a ¹H-coupled ¹H-¹⁵N HSQC-type spectrum. This allows the determination of the cross-correlation rates with no need for correction factors associated with other methods. The signal overlap problem of coupled HSQC spectra is addressed here by using the IPAP scheme (Ottiger et al., 1998). The approach is applied to the B3 domain of protein G to show that the method provides accurate measurements of the ¹⁵N CSA/dipolar cross-correlation rates.     

Computer assignment of the backbone resonances of labelled proteins using two-dimensional correlation experiments

Summary We present ALPS (Assignment for Labelled Protein Spectra), a flexible computer program for the automatic assignment of backbone NMR resonances of ¹⁵N/¹³C-labelled proteins. The program constructs pseudoresidues from peak-picking lists of a set of two-dimensional triple resonance experiments and uses either a systematic search or a simulated annealing-based optimization to perform the assignment. This method has been successfully tested on two-dimensional triple resonance spectra of Rhodobacter capsulatus ferrocytochrome c ₂ (116 amino acids).     

The relative orientation of the fibronectin 6F11F2 module pair: A 15N NMR relaxation study

The structure of a pair of modules (⁶F1¹F2), that forms part of the collagen-binding region of fibronectin, is refined using heteronuclear relaxation data. A structure of the pair was previously derived from ¹H-¹H NOE and ³ J _HHN data [Bocquier et al. (1999) Structure, 7, 1451–1460] and a weak module–module interface, comprising Leu19 and Leu28, in ⁶F1, and Tyr68 in ²F1, was identified. In this study, the definition of the average relative orientation of the two modules is improved using the dependence of ¹⁵N relaxation on rotational diffusion anisotropy. This structure refinement is based on the selection of a subset of structures from sets calculated with NOE and ³ J _HHN data alone, using the quality of the fits to the relaxation data as the selection criterion. This simple approach is compared to a refinement strategy where ¹⁵N relaxation data are included in the force field as additional restraints [Tjandra et al. (1997) Nat. Struct. Biol., 4, 443–449].     

Addressing the overlap problem in the quantitative analysis of two dimensional NMR spectra: Application to <Superscript>15</Superscript>N relaxation measurements

A quantitative analysis of 2D ¹H-¹⁵N spectra is often complicated by resonance overlap. Here a simple method is presented for resolving overlapped correlations by recording 2D projection planes from HNCO data sets. Applications are presented involving the measurement of ¹⁵N T₁ relaxation rates in a high molecular weight protein, malate synthase G, and in a system that exchanges between folded and unfolded states, the drkN SH3 domain. By supplementing relaxation data recorded in the conventional way as a series of 2D ¹H-¹⁵N data sets with a series of a pair of projection planes the number of dynamics probes is increased significantly for both systems studied.     

Computer-assisted assignment of 2D1H NMR spectra of proteins: Basic algorithms and application to phoratoxin B

Summary A suite of computer programs (CLAIRE) is described which can be of assistance in the process of assigning 2D¹H NMR spectra of proteins. The programs embody a software implementation of the sequential assignment approach first developed by Wüthrich and co-workers (K. Wüthrich. G. Wider, G. Wagner and W. Braun (1982)J. Mol. Biol. 155, 311). After data-abstraction (peakpicking), the software can be used to detect patterns (spin systems), to find cross peaks between patterns in 2D NOE data sets and to generate assignments that are consistent with all available data and which satisfy a number of constraints imposed by the user. An interactive graphics program calledCONPAT is used to control the entire assignment process as well as to provide the essential feedback from the experimental NMR spectra. The algorithms are described in detail and the approach is demonstrated on a set of spectra from the mistletoe protein phoratoxin B, a homolog of crambin. The results obtained compare well with those reported earlier based entirely on a manual assignment process.     

MUNIN: A new approach to multi-dimensional NMR spectra interpretation

A new method, MUNIN (Multi-dimensional NMR spectra interpretation), is introduced for the automated interpretation of three-dimensional NMR spectra. It is based on a mathematical concept referred to as three-way decomposition. An NMR spectrum is decomposed into a sum of components, with each component corresponding to one or a group of peaks. Each component is defined as the direct product of three one-dimensional shapes. A consequence is reduction in dimensionality of the spectral data used in further analysis. The decomposition may be applied to frequency-domain or time-domain data, or to a mixture of these. Features of MUNIN include good resolution in crowded regions and the absence of assumptions about line shapes. Uniform sampling of time-domain data, a prerequisite for discrete Fourier transform, is not required. This opens an avenue for the processing of NMR data that do not follow oscillating behaviour, e.g. from relaxation measurements. The application of MUNIN is illustrated for a ¹H-¹⁵N-NOESY-HSQC, where each component is defined as the set of all NOE peaks formed by a given amide group. As a result, the extraction of structural information simply consists of one-dimensional peak picking of the shape along the NOE-axis obtained for each amide group.     

A general Monte Carlo/simulated annealing algorithm for resonance assignment in NMR of uniformly labeled biopolymers

We describe a general computational approach to site-specific resonance assignments in multidimensional NMR studies of uniformly ¹⁵N,¹³C-labeled biopolymers, based on a simple Monte Carlo/simulated annealing (MCSA) algorithm contained in the program MCASSIGN2. Input to MCASSIGN2 includes lists of multidimensional signals in the NMR spectra with their possible residue-type assignments (which need not be unique), the biopolymer sequence, and a table that describes the connections that relate one signal list to another. As output, MCASSIGN2 produces a high-scoring sequential assignment of the multidimensional signals, using a score function that rewards good connections (i.e., agreement between relevant sets of chemical shifts in different signal lists) and penalizes bad connections, unassigned signals, and assignment gaps. Examination of a set of high-scoring assignments from a large number of independent runs allows one to determine whether a unique assignment exists for the entire sequence or parts thereof. We demonstrate the MCSA algorithm using two-dimensional (2D) and three-dimensional (3D) solid state NMR spectra of several model protein samples (α-spectrin SH3 domain and protein G/B1 microcrystals, HET-s_218–289 fibrils), obtained with magic-angle spinning and standard polarization transfer techniques. The MCSA algorithm and MCASSIGN2 program can accommodate arbitrary combinations of NMR spectra with arbitrary dimensionality, and can therefore be applied in many areas of solid state and solution NMR.     

Quantitative internuclear distancesvia two-dimensional nuclear magnetic resonance spectra: A test case and a DNA octamer duplex

Two-dimensional proton nuclear magnetic resonance nuclear Overhauser effect experiments have been performed at a series of mixing times on proflavine and on a DNA octamer duplex [d-(GGAATTCC)]₂ in solution. Using the complete matrix approach recently explored theoretically (Keepers and James, 1984), proton-proton internuclear distances were determined quantitatively for proflavine from the two-dimensional nuclear Overhauser effect results. Since proflavine is a rigid molecule with X-ray crystal structure determined, interproton distances obtained from the two-dimensional nuclear Overhauser effect experiments in solution can be compared with those for the crystalline compound agreement is better than 10 %. Experimental two-dimensional nuclear Overhauser effect spectral data for [d-(GGAATTCC)]₂ were analyzed by comparison with theoretical two-dimensional nuclear Overhauser effect spectra at each mixing time calculated using the complete 70 × 70 relaxation matrix. The theoretical spectra were calculated using two structures: a standard B-form DNA structure and an energy-minimized structure based on similarity of the octamer's six internal residues with those of [d-(CGCGAATTCGCG)]₂, for which the crystal structure has been determined. Neither the standard B-DNA nor the energy-minimized structure yield theoretical two-dimensional nuclear Overhauser effect spectra which accurately reproduce all experimental peak intensities. But many aspects of the experimental spectra can be represented by both the B-DNA and the energy-minimized structure. In general, the energy-minimized structure yields theoretical two-dimensional nuclear Overhauser effect spectra which mimic many, if not all, features of the experimental, spectra including structural characteristics at the purine-pyrimidine junction.     

Model-free analysis for large proteins at high magnetic field strengths

Protein backbone dynamics is often characterized using model-free analysis of three sets of ¹⁵N relaxation data: longitudinal relaxation rate (R ₁), transverse relaxation rate (R ₂), and ¹⁵N–{H} NOE values. Since the experimental data is limited, a simplified model-free spectral density function is often used that contains one Lorentzian describing overall rotational correlation but not one describing internal motion. The simplified spectral density function may be also used in estimating the overall rotational correlation time, by making the R ₂/R ₁ largely insensitive to internal motions, as well as used as one of the choices in the model selection protocol. However, such approximation may not be valid for analysis of relaxation data of large proteins recorded at high magnetic field strengths since the contribution to longitudinal relaxation from the Lorentzian describing the overall rotational diffusion of the molecule is comparably small relative to that describing internal motion. Here, we quantitatively estimate the errors introduced by the use of the simplified spectral density in model-free analysis for large proteins at high magnetic field strength.     

Improved TROSY-HNCA experiment with suppression of conformational exchange induced relaxation

A general method for improving of the sensitivity of the TROSY-type triple resonance experiments in the presence of conformational exchange-induced (CSX) relaxation is proposed based on the use of CPMG-INEPT (Müller et al., J. Am. Chem. Soc., 1995, 117, 11043–11048) during the N–C polarization transfer periods. Significantly improved sensitivity is demonstrated for the majority of cross-peaks in the new [¹⁵N,¹H]-TROSY-XY-HNCA experiment, measured with partially folded RNase AS-Protein, with negligible loss of sensitivity for resonances unaffected by CSX relaxation. In addition, a comparison of cross-peak amplitudes in [¹⁵N,¹N]-TROSY-XY-HNCA and conventional [¹⁵N,¹H]-TROSY-HNCA spectra provides a quick and sensitive estimation of the CSX relaxation contribution.     

Side-chain to backbone correlations from solid-state NMR of perdeuterated proteins through combined excitation and long-range magnetization transfers

Proteins with excessive deuteration give access to proton detected solid-state NMR spectra of extraordinary resolution and sensitivity. The high spectral quality achieved after partial proton back-exchange has been shown to start a new era for backbone assignment, protein structure elucidation, characterization of protein dynamics, and access to protein parts undergoing motion. The large absence of protons at non-exchangeable sites, however, poses a serious hurdle for characterization of side chains, which play an important role especially for structural understanding of the protein core and the investigation of protein–protein and protein–ligand interactions, e.g. This has caused the perdeuteration approach to almost exclusively be amenable to backbone characterization only. In this work it is shown that a combination of isotropic ¹³C mixing with long-range ¹H/¹³C magnetization transfers can be used effectively to also access complete sets of side-chain chemical shifts in perdeuterated proteins and correlate these with the protein backbone with high unambiguity and resolution. COmbined POlarization from long-Range transfers And Direct Excitation (COPORADE) allows this strategy to yield complete sets of aliphatic amino acid resonances with reasonable sensitivity.     

Biomass production of site selective 13C/15N nucleotides using wild type and a transketolase E. coli mutant for labeling RNA for high resolution NMR

Characterization of the structure and dynamics of nucleic acids by NMR benefits significantly from position specifically labeled nucleotides. Here an E. coli strain deficient in the transketolase gene (tktA) and grown on glucose that is labeled at different carbon sites is shown to facilitate cost-effective and large scale production of useful nucleotides. These nucleotides are site specifically labeled in C1′ and C5′ with minimal scrambling within the ribose ring. To demonstrate the utility of this labeling approach, the new site-specific labeled and the uniformly labeled nucleotides were used to synthesize a 36-nt RNA containing the catalytically essential domain 5 (D5) of the brown algae group II intron self-splicing ribozyme. The D5 RNA was used in binding and relaxation studies probed by NMR spectroscopy. Key nucleotides in the D5 RNA that are implicated in binding Mg²⁺ ions are well resolved. As a result, spectra obtained using selectively labeled nucleotides have higher signal-to-noise ratio compared to those obtained using uniformly labeled nucleotides. Thus, compared to the uniformly ¹³C/¹⁵N-labeled nucleotides, these specifically labeled nucleotides eliminate the extensive ¹³C–¹³C coupling within the nitrogenous base and ribose ring, give rise to less crowded and more resolved NMR spectra, and accurate relaxation rates without the need for constant-time or band-selective decoupled NMR experiments. These position selective labeled nucleotides should, therefore, find wide use in NMR analysis of biologically interesting RNA molecules.     

High resolution <Superscript>13</Superscript>C-detected solid-state NMR spectroscopy of a deuterated protein

High resolution ¹³C-detected solid-state NMR spectra of the deuterated beta-1 immunoglobulin binding domain of the protein G (GB1) have been collected to show that all ¹⁵N, ¹³C′, ¹³Cα and ¹³Cβ sites are resolved in ¹³C–¹³C and ¹⁵N–¹³C spectra, with significant improvement in T ₂ relaxation times and resolution at high magnetic field (750 MHz). The comparison of echo T ₂ values between deuterated and protonated GB1 at various spinning rates and under different decoupling schemes indicates that ¹³Cα T ₂′ times increase by almost a factor of two upon deuteration at all spinning rates and under moderate decoupling strength, and thus the deuteration enables application of scalar-based correlation experiments that are challenging from the standpoint of transverse relaxation, with moderate proton decoupling. Additionally, deuteration in large proteins is a useful strategy to selectively detect polar residues that are often important for protein function and protein–protein interactions.