Six- and seven-dimensional experiments by combination of sparse random sampling and projection spectroscopy dedicated for backbone resonance assignment of intrinsically disordered proteins.

Two novel six- and seven-dimensional NMR experiments are proposed. The new experiments employ non-uniform sampling that enables achieving high resolution in four indirectly detected dimensions and synchronous sampling in the additional dimensions using projection spectroscopy principle. The resulted data sets could be processed as five-dimensional data using existing software. The experiments facilitate resonance assignment of intrinsically disordered proteins. The novel experiments were successfully tested using 1 mM sample of α-synuclein on 600 and 800 MHz NMR spectrometers equipped with standard room temperature probes. The experiments allowed backbone assignment from a 1-day acquisition.

Introduction

In recent years, intrinsically disordered proteins (IDPs) have attracted widespread interest in molecular biology research (Wright and Dyson 1999; Habchi et al. 2014). Their properties like structural propensity, dynamics and interactions could be effectively studied by solution NMR spectroscopy methods. However, the intrinsic disorder results in a fast conformational dynamics which causes an effective averaging of chemical shifts. Thus, very poor peak separation makes resonance assignment difficult, even for a relatively small disordered protein fragments. As a consequence, the signal overlap in the spectra strongly limits the possibility of standard NMR experiments to provide sufficient resolution for IDPs characterization. On the other hand, the fast local dynamics significantly decreases transverse relaxation rates, enabling the application of long multidimensional pulse sequences, and making possible to achieve long evolution times in order to increase the resolution and facilitate assignment of resonances.

The projection NMR spectroscopy methods are founded on the idea of Accordion Spectroscopy (Bodenhausen and Ernst 1981), and employ synchronous incrementing of two or more evolution periods in the pulse sequence. The particularly important approach of this type is simultaneous sampling of more than one chemical shift evolution, which is referred to as radial sampling. Such an option is utilized in projection spectroscopy (Kupče and Freeman 2005; Coggins et al. 2004, 2010), and requires the algebraic decoding of peak frequencies (Kim and Szyperski 2003; Koźmiński and Zhukov 2003; Hiller et al. 2005), or the reconstruction of multidimensional spectrum (Kupče and Freeman 2003).

To date, several high dimensional strategies have been proposed for backbone resonance assignment of IDPs (Reddy and Hosur 2014; Goradia et al. 2015; Yoshimura et al. 2015), including direct 13C detection (Bermel et al. 2006, 2009; Nováček et al. 2011; Pantoja-Uceda and Santoro 2013, 2014) and HA detection (Mäntylahti et al. 2010, 2011; Yao et al. 2014) which allow for revealing of proline residues. The number of dimensions, i.e. number of simultaneously correlated frequencies, could be increased by employing automated projection spectroscopy (APSY) (Hiller et al. 2007; Narayanan et al. 2010; Yao et al. 2014), or the non-uniformly sampled (NUS) high-dimensional experiments (Kazimierczuk et al. 2010; Zawadzka-Kazimierczuk et al. 2012).

APSY approach enabled the acquisition of experiments with dimensionality up to 7, by recording a set of 2D projections of high dimensionality spectrum at different angles. The high-dimensional peak list is achieved by algebraic reconstruction from automatically collected peak coordinates in the projections. The possible disadvantage of this technique is severe peak overlap on 2D spectral projections, which may influence the precision of frequency determination.

The advantage of non-uniformly sampled high-dimensionality experiments, followed by SMFT (sparse multidimensional Fourier transform) processing (Kazimierczuk et al. 2009), is the ability for visual inspection of the spectrum of the full dimensionality. Such examination of 2D cross-sections helps to assign peaks suffering from partial overlap.

SMFT processing procedure with the output in the form of a series of 2D cross-sections, requires, however, a so-called base spectrum of dimensionality of N-2, which becomes impractical for N > 5. Therefore, we propose here new 6 and 7D experiments that feature high resolution and high dimensionality resulting from the use of non-uniform sampling in the indirectly detected dimensions combined with projection spectroscopy principle. The experiments utilize NUS in four indirectly sampled dimensions, while additional frequencies are added by synchronous incrementing with t1. Thus, the previously described 5D SMFT processing with a 3D base spectrum could be used for the proposed experiments, resulting in a 5D spectrum with additional dimensions projected into F1.

Since the N and CO chemical shifts were found to be best dispersed in IDPs (Piai et al. 2014; Nowakowski et al. 2015), the proposed experiments are derived from 5D (H)NCO(NCA)CONH (Zawadzka-Kazimierczuk et al. 2012), by adding HN (6D) or HN and CA (7D) chemical shifts. The experimental verification was performed on a 1 mM sample of α-synuclein, employing 600 and 800 MHz NMR spectrometers with standard RT triple resonance probes.

Materials and methods

Pulse sequence

The proposed pulse sequence is based on 5D (H)NCO(NCA)CONH presented in (Zawadzka-Kazimierczuk et al. 2012) and is depicted in the Fig. 1.

Fig. 1

HNCO(N)CACONH experiment. Top scheme shows the coherence pathway involved in the experiment. Bottom is the pulse sequence scheme. Rectangles represent hard pulses. Filled and empty symbols represent 90° and 180° pulses, respectively. 1H and 15N composite pulse decouplings are performed with WALTZ-16 (Shaka et al. 1983), at γB1/2π of 5.4 and 1.14 kHz at the 800 MHz spectrometer, and 4.0 and 0.97 kHz at 600 MHz, respectively. Simultaneous inversion of CA and CO spins was archived using 6-element composite pulse (Shaka 1985). Selective CA and CO 13C 90° (180°) pulses were applied with rf field strength adjusted to |ΔΩCA-CO|/√15 (√3). At the 800 MHz spectrometer 90° and 180°, rectangular and sinc-shaped pulses (bell-shaped at the scheme) of the duration of 40.1 (35.9) μs, and 65.8 (58.8) μs, respectively, were used. Whereas, at the 600 MHz spectrometer of, 53.5 (47.9) μs, and 87.8 (78.4) μs, respectively. Off-resonance pulses were applied using phase modulation of the carrier. The amplitude (Tm−1) of G1–G9 PFG pulses were set to: 0.212, 0.154, 0.137, 0.0926, 0.0820, 0.347, 0.154, 0.0579, 0.352 at the 800 MHz spectrometer and at 600 MHz: 0.222, 0.162, 0.143, 0.0970, 0.0858, 0.364, 0.162, 0.0606, 0.368, respectively. The PFG duration of G1–G5, G8 of 0.5 ms, G6 and G7 of 2.0 ms and G9 of 0.2 ms, were used. Evolutions for H, N, CA were in semi-constant-time mode: a  = (t  + Δ)/2; b  = t (1 − Δ/t imax)/2; c  = Δ(1 − t /t imax)/2 or in constant-time mode: a  = (t  + Δ)/2; b  = 0; c  = (Δ − t )/2 where Δ stands for ΔN–H, ΔN–CO, ΔN–CA. Evolution for CO in t 2 is in real-time mode. Delays were set as follows: ΔN–H = 5.4 ms, ΔN–CO = 28 ms, ΔCO–N–CA = 28 ms, ΔN–CA = 54 ms, ΔCO–CA = 9.1 ms. The four step phase cycle was used: ϕ1 = x, -x; ϕ2 = 2x, 2(-x) and ϕrec = x, 2(-x), x = ϕ1 + ϕ2. In t 1, t 2, t 3, t 4 dimensions quadrature was accomplished using States-TPPI method, by incrementing ϕ1, ϕ2, ϕ3, ϕ4 phases, respectively. In t 5, t 6 dimensions quadrature was accomplished using States method by incrementing ϕ5, ϕ6 phases and adding additional increments to phase ϕ1 = ϕ1 + ϕ5 + ϕ6, thus, the ϕ1 phase was incremented by additional 90° for each sine modulation in t 1, t 5 and t 6. The phase ψ = x was inverted simultaneously with the last gradient (G9) pulse to achieve echo-antiecho coherence transfer selection in the indirect dimension. The coherence selection gradients (G6 and G9) were applied at magic angle (600 MHz) or along z-axis (800 MHz). 90° and 180° water 1.2 ms sinc-shaped flipback pulses were used for ϕ5 phase equal to x and y, respectively. For 5D HNCO(N)CACONH experiment t 5, t 6 and ϕ5, ϕ6 were set to 0 s and x, respectively. Additional dimensions are achieved by setting evolution times (t 5 and/or t 6) proportional to t 1. The ratio of maximum evolution times determined the respective projection angle. Note, that two different 6D and one 7D experiments could be acquired according to presented scheme

The original pulse sequence was designed to obtain both sequential and auto-correlation peaks by setting ΔN–CA delay to 28.6 ms. It surely allows easier visual inspection of the spectra. On the other hand, the doubled number of peaks increases sampling artefacts level. Moreover, lack of auto-correlation peaks present in the spectrum reduces spectral overlap, which may affect accuracy of measured peaks positions. This is especially important while studying IDPs. The coherence transfer amplitudes are described by Eqs. (1) and (2):

Assuming , values of 11, 7 and 35 Hz (Sattler et al. 1999), respectively, the optimal choice of ΔN–CA delay length for both high intensity of sequential signal (Ni−1(ω1), COi−2(ω2), COi−1(ω3), Ni(ω4), HiN(ω5)) and the suppression of auto-correlation signal (Ni(ω1), COi−1(ω2), COi−1(ω3), Ni(ω4), HiN(ω5)) can be obtained with the value of 54 ms (see Fig. 2), similar approach was also presented by Nováček et al. 2011. In the present work the auto-correlation peaks suppressed version was used.

Fig. 2

Transfer efficiencies for auto-correlation and sequential peaks as a function of ΔN–CA delay length

Quadrature for t1 dimensions is achieved using a standard States-TPPI procedure. Quadratures in coevolved dimensions require recording and storing of two data sets, cosine and sine modulated, respectively. It is obtained by 90 degree shift of phases ϕ5, and ϕ5, ϕ6 (for 6D and 7D versions respectively) (Koźmiński and Zhukov 2003; Kim and Szyperski 2003). To allow standard SMFT data processing the π/2 phase shifts for sine modulated components are compensated by simultaneous π/2 shifting of ϕ1 phase, for t1 dimension. Therefore, such procedure of processing of spectra after SMFT approach is limited to a simple co-addition to obtain spectra differing in frequency signs in coevolved dimensions.

Data processing and inspection

The presented experiments are, in fact, projections of six- and seven-dimensional experiments to the five- dimensional frequency space. Therefore, both of the presented experiments were processed as a usual 5D spectrum using SMFT approach. Last three dimensions (CO, N, HN) were fixed with chemical shifts obtained from 3D HNCO experiment. This lead to a set of 2D cross-sections showing correlations of coevolved dimensions (t1, t5, t6) with CO (t2). A simple co-addition of differently modulated data sets result with two data sets (in a case of 6D), which differ with signs of involved frequencies (7D version results in four different combinations), this procedure is schematically illustrated in the Figs. 3 and 4.

Fig. 3

Pictorial representation of resulting peak pattern when two dimensions are coevolved. Two modulations lead to spectra I and II containing 2 peaks each. Sum and difference of obtained spectra gives two different spectra containing only one peak each, for which frequencies are a linear combinations of frequencies from coevolved dimensions. a coefficient is equal to t 2/t 1. Empty and filled circles represent positive and negative signal intensities, respectively

Fig. 4

Pictorial representation of resulting peak pattern when three dimensions are coevolved. Four modulations lead to spectra I, II, III and IV containing 4 peaks each. Co-addition of obtained spectra gives four different spectra containing only one peak each, for which frequencies are a linear combinations of frequencies from coevolved dimensions. a, and b coefficients are equal to t 2/t 1, and t 3/t 1, respectively. Empty and filled circles represent positive and negative signal intensities, respectively

All 5D spectra were processed using ToASTD (Kazimierczuk et al. 2006) and reduced (Kazimierczuk et al. 2009) programs which take peak list from 3D HNCO as the input for SMFT routine. For 3D HNCO sampling artefacts were removed from the spectra using Signal Separation Algorithm (Stanek and Koźmiński 2010), no such procedure was performed for 5D spectra. All used programs are available at http://nmr.cent3.uw.edu.pl/software. After processing, all spectra were inspected using Sparky (Goddard and Kneller 2000).

Four distinct peaks frequencies in the projected dimension: , …, on the transformations for 6D experiment are given by the Eq. (3):In the case of 7D data the analogous expressions for peak frequencies: , …, are given by the Eq. (4),, are the values of resonance frequencies of N, H and CA nuclei of the residue with respect to the corresponding carrier offsets. , , are maximum evolution times in N, H and CA dimensions, respectively. As a result, for 6D spectrum, a system of four equations (two of them are linearly independent) with two unknown is obtained. In the case of 7D spectrum system of 8 equations (four linearly independent) with three unknown is created. Solving these systems of equations results in , , frequencies. The first system of equations (Eq. 3) is unambiguously defined, however, the second one (Eq. 4) is overdetermined. Three equations are sufficient to solve it but in order to increase the precision of a result, all four equations with equal weights were always used.

NMR spectroscopy

NMR sample contained 1 mM of 13C, 15N-labeled α-synuclein in 20 mM sodium phosphate buffer, pH 6.5, 200 mM NaCl. NMR spectra were recorded at 288 K on Agilent 800 MHz and 600 MHz spectrometers both equipped with room temperature probes. 7D HNCO(N)CACONH was acquired on the 800 MHz spectrometer in 31 h with evolution times set to 20, 30, 25, 52, 25, 30 ms for HiN, Ni, COi−1, CAi−1, COi, Ni+1 dimensions, respectively. 6D HNCO(NCA)CONH was acquired on both, 800 MHz and 600 MHz, spectrometers in 13 and 23 h, respectively. Evolution times for the measurement on the 800 MHz spectrometer were set to: 20, 30, 25, 25, 30 ms, for HiN, Ni, COi−1, COi, Ni+1 dimensions, respectively, and to sustain high enough resolution on the lower field evolution times were set to: 20, 40, 40, 40, 40 ms, respectively on the 600 MHz spectrometer.

Results

In order to verify results obtained from 6D and 7D projection experiments full backbone assignment of the α-synuclein was performed using 5D HN(CA)CONH (Kazimierczuk et al. 2010), 5D (H)NCO(NCA)CONH (Zawadzka-Kazimierczuk et al. 2012) and 5D HabCabCONH (Kazimierczuk et al. 2010) experiments. To avoid any discrepancies, exactly the same sample as for 6D and 7D measurements was used.

Analysis of both recorded 6D spectra led to identification of all but two expected sequential peaks. Two missing residues were the first two residues in the protein sequence which suggest that those signals were broadened beyond detection limit due to amide proton exchange. Similar result was obtained from 7D spectrum, the same set of sequential signals was identified.

Performing sequential assignment of α-synuclein due to its disordered nature can be regarded as a challenging task. Despite α-synuclein’s moderate size (140 a.a.), full potential of the “HNCO-HNCO” strategy can be seen, for example, in the case of G36CO-V37N and G73CO-V74N signals which are hardly distinguishable using COi−1Ni connectivity (see Fig. 5). Application of 6D experiment quickly resolves such ambiguity without using any additional experiment as the HN chemical shifts of aforementioned residues vary by almost 0.2 ppm (see Fig. 6). Moreover, peak positions in the resulting spectra are further differentiated when seventh dimension (CAi−1) is introduced (see Fig. 7). In addition to a better signal dispersion CA chemical shift provide partial information about residue type.

Fig. 5

Two 2D cross-sections from 5D (H)NCO(NCA)CONH experiment (on the left) corresponding to HNCO peak positions of G36CO-V37N and G73CO-V74N signals showing COi−1Ni sequential peaks. 2D NCO projection from 3D HNCO (on the right) showing highly crowded glycine region. Unambiguous assignment is prevented due to signal overlap

Fig. 6

Four 2D cross-sections from 6D HNCO(NCA)CONH experiment (at the top) (++ and +− corresponds to the signs in the co-addition of obtained raw data sets) corresponding to the HNCO peak positions of G36CO-V37N and G73CO-V74N signals showing COi−1(N H ) sequential peaks. Small difference in frequencies in the measured peak positions encodes additional information of HN chemical shifts. Two slices from 3D HNCO spectra obtained at HN the positions calculated from 6D experiment data (at the bottom). Unambiguous assignment is possible owing to the difference in HN chemical shifts of V37 and V74. Note that y-axis at the cross-sections from the 6D HNCO(NCA)CONH spectrum is labelled with 15N chemical shift scale, however, the peak frequencies in this dimension are given by the Eq. (3)

Fig. 7

Eight 2D cross-sections from 7D HNCO(N)CACONH experiment (++++ , ++−, +−+− and +− + correspond to the signs in the co-addition of obtained raw data sets) corresponding to the HNCO peak positions of G36CO-V37N and G73CO-V74N signals showing CA i−1COi−1(N H N) sequential peaks. While coevolution of HN leads to a more robust connectivity adding of CA frequency further differentiate signal positions due to the difference in CA chemical shifts of V37 and V74. Note that y-axis at the cross-sections from the 7D HNCO(N)CACONH spectrum is labelled with 15N chemical shift scale, however, the peak frequencies in this dimension are given by the Eq. (4)

Discussion

The most important advantage of presented approach is the establishment of highly robust connectivities due to the expansion of sequential peaks frequencies by addition of Hi−1N frequency. Similar approach was already proposed in the literature which involved the combination of two different 5D experiments (Kazimierczuk et al. 2013; Piai et al. 2014) having shared sequential nitrogen dimension. Here it was presented that such “HNCO-HNCO” strategy can be successfully implemented within one experiment. Additional expansion to 7D experiment giving CA chemical shift can further limit number of overlapping signals as well as provide partial information needed for identifying amino acid type.

True benefit of proposed approach is the possibility of a visual inspection of resulting spectra which allows to use the expertise and the experience of a spectroscopist in resolving most difficult, severely overlapped cases. As a result, all of non-proline residues of α-synuclein were successfully identified and assigned from a single experiment (except aforementioned first two residues).

Surprisingly, even example of a moderate size α-synuclein already justifies expanding “CON-CON” strategy. Presented approach can be especially beneficial in performing the resonance assignment of proteins of a larger size than α-synuclein, since the spectra in such a case will be more crowded. What is more, sensitivity of presented experiments is high enough to record good quality spectra in a relatively short time compared to other six- and seven-dimensional approaches (Hiller et al. 2007; Yao et al. 2014), even without the use of cryogenically cooled probe.

We have demonstrated the principles of the new experiment and data inspection protocol. Although, it could be modified in several ways. The proposed experiment is fully compatible with BEST approach (Lescop et al. 2007; Solyom et al. 2013) and co-solute paramagnetic relaxation enhancement (Oktaviani et al. 2015), which could be used for acceleration of the signal repetition rate, and thus increase the number of data points acquired in a given time. Amide hydrogen exchange process, particularly effective in the IDP molecules, may cause a severe signal loss at conditions close to physiological (temperature and pH). Therefore, the pulse sequence could be modified using the aliphatic proton excitation combined with the 13C or HA detection (Mäntylahti et al. 2010). Moreover, aliphatic protons or 13C detection would allow observation of additional cross-peaks from proline residues (Bermel et al. 2006; Hellman et al. 2014). Additionally, the amide protons exchange causes further decay of 15N–1H antiphase coherences. Hence, the refocusing of 15N–1H couplings and 1H composite pulse decoupling applied during the main part of pulse sequence is used to maximize the sensitivity.

Conclusions

We presented six- and seven-dimensional experiments for backbone assignment of intrinsically disordered proteins. Described techniques combine projection spectroscopy and SMFT data processing that allow access to full dimensional spectra information. 6D experiment correlates HiN, Ni, COi−1, Hi−1N, Ni−1, COi−2. 7D adds additional correlation with CAi−1 dimension. We proved that these experiments can be successfully run on a standard room-temperature probe, even on a medium field NMR spectrometer (600 MHz).