Method of analyzing 3D NMR spectrum

Abstract

There is disclosed a method of symmetrizing a three-dimensional NMR spectrum which is obtained by applying a preparation pulse or pulse sequence, sampling the resulting free induction decay signals with a sequence having two evolution periods and a detection period to obtain data S (t.sub.1, t.sub.2, t.sub.3), and taking the Fourier transform of the data to create the three-dimensional NMR spectrum .alpha. (.omega..sub.1, .omega..sub.2, .omega..sub.3). The 3D NMR spectrum is sliced perpendicularly to the frequency axis .omega..sub.3. Data about the 3D NMR spectrum is arithmetically processed in such a way that peaks having no symmetrical partners at symmertical positions determined also taking account of the direction of the frequency axis .omega..sub.3 are found and erased to obtain a spectrum .alpha..sub.1 (.omega..sub.1, .omega..sub.2, .omega..sub.3). The data about this spectrum is arithmetically processed in such a way that peaks complementary to the peaks having no symmetrical partners in the sliced planes of the spectrum .alpha..sub.1 (.omega..sub.1, .omega..sub.2, .omega..sub.3 ) are added to obtain a symmetrized spectrum .alpha..sub.2 (.omega..sub.1, .omega..sub.2, .omega..sub.3).

Description

FIELD OF THE INVENTION

The present invention relates to three-dimensional NMR spectroscopy and, more particularly, to a method of analyzing the three-dimensional spectrum for symmetrical peaks by first symmetrizing a three-dimensional NMR spectrum.

BACKGROUND OF THE INVENTION

In recent years, two-dimensional NMR has enjoyed wide acceptance. Two-dimensional NMR spectra are easy to analyze. Especially, the interaction between nuclear spins can be elucidated by two-dimensional NMR spectroscopy.

The prior art 2D NMR experiment is now described in detail by referring to FIG. 1, where a pulse sequence with two 90.degree. pulses is used. A general process of measurement by 2D NMR spectroscopy consists of a preparation period preceding the first 90.degree. pulse, an evolution period of t.sub.1, and a detection period of t.sub.2. The preparation period is necessary to maintain the nuclear spin magnetization in its appropriate initial condition. The preparation pulse, or the first 90.degree. pulse, brings the magnetization into nonequilibrium state. This state is caused to evolve in the evolution period t.sub.1. In the detection period t.sub.2 subsequent to the application of the detection period, or the second 90.degree. pulse, the resulting free induction decay signal is detected. The phase and the amplitude of this decay signal reflect the behavior of the magnetization taken in the evolution period t.sub.1. The period t.sub.1 is varied so as to assume n discrete values. The free induction decay signal is detected at each of the n values. As a result, a set of data S (t.sub.1, t.sub.2) is obtained from n free induction decay signals. The obtained data includes information about the behavior of the magnetization in the evolution period t.sub.1, as well as information about the behavior of the magnetization in the detection period t.sub.2. Data about a 2D NMR spectrum .alpha.(.omega..sub.1, .omega..sub.2) is derived by taking the 2D Fourier transform of the data S (t.sub.1, t.sub.2) with respect to t.sub.2 and t.sub.1. FIG. 2 is one example of this 2D NMR spectrum, and in which the spectrum is plotted on the .omega..sub.1 -.omega..sub.2 plane. In 2D NMR spectra, diagonal peaks and cross-peaks appear symmetrically. The relations between the peaks can be known from the symmetrical peaks.

If numerous superimposed peaks appear in a 2D NMR spectrum, it is difficult to understand which of the peaks are coupled to each other. Accordingly, 3D NMR has been proposed which separates the superimposed peaks and permits one to understand them.

FIG. 3 shows a simplest pulse sequence used in 3D NMR spectroscopy. This pulse sequence consists of three 90.degree. pulses. A preparation period, a first evolution period of t.sub.1, a second evolution period of t.sub.2, and a detection period of t.sub.3 are separated by the three pulses.

The first evolution period t.sub.1 is varied so as to take n discrete values. The second evolution period t.sub.2 is changed so as to assume m discrete values. As a result, n.times.m free induction decay signals are produced, and all of them are detected. Thus, a set of data S (t.sub.1, t.sub.2, t.sub.3) is derived from the n.times.m free induction decay signals. The data includes information about the behavior of the magnetization in the period t.sub.3. In addition, the data includes information about the behavior of the magnetization in the periods t.sub.1 and t.sub.3. A three-dimensional NMR spectrum .alpha. (.omega..sub.1, .omega..sub.2, .omega..sub.3) having three frequency parameters .omega..sub.1, .omega..sub.2, .omega..sub.3 is obtained by taking the three-dimensional Fourier transform of the data S (t.sub.1, t.sub.2, t.sub.3) with respect to t.sub.3, t.sub.2, t.sub.1. In this case, .omega..sub.1, .omega..sub.2 and .omega..sub.3 vary in the same domain. However, in one of the most useful 3D experiments as shown in FIG. 9, .omega..sub.1 and .omega..sub.2 vary in the same domain and .omega..sub.3 varies in another domain. 3D experiments developed from homonuclear 2D experiments by adding a dimension of another kind of nucleus belong to this class. Homonuclear 2D spectra are symmetrical and 3D spectra also have this property.

If this 3D NMR spectrum is sliced (.alpha.(.omega..sub.1, .omega..sub.2, .omega..sub.3)) at right angles to the frequency axis .omega..sub.3 at a certain value of .omega..sub.3, for example .omega..sub.3 =.omega..sub.30, then one 2D NMR spectrum is obtained. This 2D NMR spectrum differs greatly in symmetry from 2D NMR spectra obtained by ordinary 2D NMR spectroscopy. In particular, a 2D NMR spectrum obtained by 2D NMR spectroscopy is essentially symmetrical, i.e., .alpha.(.omega..sub.1, .omega..sub.2)=.alpha.(.omega..sub.2, .omega..sub.1). However, this equation does not hold for a 2D NMR spectrum derived by slicing a 3D NMR spectrum.

More specifically, two nuclear spins interact with each other if they are chemically directly coupled to each other or if they are close to each other. As an example, as shown in FIG. 4(a), we now discuss the case in which protons .sup.1 H are spatially close to each other and chemically coupled to .sup.15 N. Let .delta..sub.1, .delta..sub.2, .delta..sub.3, and .delta..sub.4 be the chemical shifts of these four nuclear spins .sup.1 H and .sup.15 N. Since couplings A and B are symmetrical, the corresponding peaks in the 2D NMR spectrum appear symmetrically. However, the chemical shift caused by .sup.1 H differs from .delta..sub.3 and .delta..sub.4. In a 3D NMR spectrum in which these chemical shifts can be distinguished, a pair of symmetrical peaks appear in different sliced planes that are shifted along the frequency axis .omega..sub.3. In each sliced plane, or in each 2D NMR spectrum, the spectrum is asymmetrical. Therefore, in order to search sliced planes for symmetrical peaks, sliced planes or 2D NMR spectra which differ in .omega..sub.3 must be scrutinized one after another. It is very cumbersome and time-consuming to perform this work.

As shown in FIG. 4(b), where carbon nucleus C is coupled to proton .sup.1 H causing chemical shift .delta..sub.1, coupling B has no symmetrical partner. For this reason, in a 3D NMR spectrum, only a peak corresponding to coupling B appears; any peak forming a symmetrical partner does not appear. In this way, information that is useful for making analysis, based on symmetrical peaks in a 3D NMR spectrum, is not contained in those peaks which have no symmetrical partners. Rather, such peaks complicate spectra and can be said to be harmful.

SUMMARY OF THE INVENTION

It is an object of the present invention to provide a method of easily and quickly searching a three-dimensional NMR spectrum for symmetrical peaks by symmetrizing the spectrum in such a way that symmetrical peaks appear in the same sliced plane at all times.

It is another object of the invention to provide a method of easily and quickly searching the spectrum for the symmetrical peaks by symmetrizing the spectrum in such a way that peaks having no symmetrical partners are removed.

The above objects are achieved by method of symmetrizing a three-dimensional NMR spectrum .alpha.(.omega..sub.1, .omega..sub.2, .omega..sub.3) which is obtained by applying a preparation pulse or pulse sequence, sampling the resulting free induction decay signals with a sequence having two evolution periods and a detection period to obtain data S (t.sub.1, t.sub.2, t.sub.3), and taking the Fourier transform of the data to create the three-dimensional NMR spectrum .alpha.(.omega..sub.1, .omega..sub.2, .omega..sub.3), said method comprising the steps of: slicing the three-dimensional NMR spectrum perpendicularly to the frequency axis .omega..sub.3 ; arithmetically processing data about the three-dimensional NMR spectrum in such a way that peaks having no symmetrical partners at symmetrical positions determined, also taking account of the direction of the frequency axis .omega..sub.3, are found and erased to obtain a spectrum .alpha..sub.1 (.omega..sub.1, .omega..sub.2, .omega..sub.3); adding peaks complementary to the peaks each existing alone in each sliced plane of the spectrum .alpha..sub.1 (.omega..sub.1, .omega..sub.2, .omega..sub.3) to obtain a symmetrized spectrum .alpha..sub.2 (.omega..sub.1, .omega..sub.2, .omega..sub.3).

Other objects and features of the invention will appear in the course of the description thereof which follows.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a diagram showing a fundamental pulse sequence used for two-dimensional NMR spectroscopy;

FIG. 2 is a two-dimensional NMR spectrum;

FIG. 3 is a diagram showing a fundamental pulse sequence used for three-dimensional NMR spectroscopy;

FIGS. 4(a) and 4(b) are diagrams showing couplings of protons .sup.1 H to .sup.15 N;

FIG. 5 is a flowchart illustrating a sequence in which the novel method is carried out;

FIG. 6 is a visual representation of data about a three-dimensional NMR spectrum;

FIG. 7 is a diagram showing the relation of .alpha.(.omega..sub.11, .omega..sub.21, .omega..sub.31) to .alpha.(.omega..sub.21, .omega..sub.11, X).vertline.X.epsilon..OMEGA..sub.3 when arbitrary values .omega..sub.11, .omega..sub.21, .omega..sub.31 are substituted into .omega..sub.1, .omega..sub.2, .omega..sub.3, respectively;

FIG. 8 is a diagram shown in the chemical structure of a sample that is a protein labeled with .sup.15 N;

FIG. 9 is a diagram of a pulse sequence used for three-dimensional NMR experiments; and

FIG. 10 shows changes in a three-dimensional NMR spectrum processed in accordance with the invention, and analyzing the symmetrized spectrum .alpha..sub.2 (.omega..sub.1, .omega..sub.2, .omega..sub.3) for symmetrical peaks.

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT

FIG. 5 is a flowchart illustrating the fundamental procedure for analyzing a three-dimensional NMR spectrum for symmetrical peaks by symmetrizing a three-dimensional NMR spectrum in accordance with this invention. We now express data about the 3D NMR spectrum by means of .alpha.(.omega..sub.1, .omega..sub.2, .omega..sub.3); where .alpha. is a real function; and .omega..sub.1, .omega..sub.2, .omega..sub.3 are real numbers. It is assumed that .omega..sub.1 and .omega..sub.2 have the same domain. It is also assumed that some peaks are scattered in the 3D NMR spectrum and that these peaks are given by delta function .delta..

At step 1, 3D NMR experiments are made, and data S (t.sub.1, t.sub.2, t.sub.3) is obtained.

At step 2, data about a three-dimensional NMR spectrum .alpha.(.omega..sub.1, .omega..sub.2, .omega..sub.3) consisting of three frequency parameters .omega..sub.1, .omega..sub.2, .omega..sub.3 is obtained by taking the three-dimensional Fourier transform of the data S (t.sub.1, t.sub.2, t.sub.3) with respect to t.sub.3, t.sub.2, t.sub.1. A visual representation of the data .alpha.(.omega..sub.1, .omega..sub.2, .omega..sub.3) is shown in FIG. 6.

At step 3, the obtained data .alpha.(.omega..sub.1, .omega..sub.2, .omega..sub.3) is arithmetically processed, using max function and min function given by equation (1), to produce data about a spectrum .alpha..sub.1 (.omega..sub.1, .omega..sub.2, .omega..sub.3) from which all the peaks having no symmetrical partners have been removed.

.alpha..sub.1 (.omega..sub.1, .omega..sub.2, .omega..sub.3)=min [.alpha.(.omega..sub.1, .omega..sub.2, .omega..sub.3), max {.alpha.(.omega..sub.2, .omega..sub.1, X).vertline.X.epsilon..OMEGA..sub.3 }] (1) where .OMEGA..sub.3 is the domain of the frequency .omega..sub.3.

FIG. 7 shows the relation of .alpha.(.omega..sub.11, .omega..sub.21, .omega..sub.31) to .alpha.(.omega..sub.21, .omega..sub.11, X) .vertline.X .epsilon..OMEGA..sub.3 when arbitrary values .omega..sub.11, .omega..sub.21, .omega..sub.31 are substituted into .omega..sub.1, .omega..sub.2, .omega..sub.3, respectively. In FIG. 7, point P is given by .alpha.(.omega..sub.11, .omega..sub.21, .omega..sub.31). Point P and P' are symmetrical with respect to the central plane .omega..sub.1 =.omega..sub.2. A straight line L which passes through the point P' and runs parallel to the axis .omega..sub.3 is given by .alpha.(.omega..sub.21, .omega..sub.11, X) .vertline.X .epsilon..OMEGA..sub.3.

Function max {.alpha.(.omega..sub.21, .omega..sub.11, X) .vertline.X .epsilon..OMEGA..sub.3 } is used to find the maximum value of the straight line L. If any peak exists on the line L, the value of the function max {.alpha.(.omega..sub.21, .omega..sub.11, X).vertline.X.epsilon..OMEGA..sub.3 } is the intensity of the peak. If no peaks exist, the value of the function max {.alpha.(.omega..sub.21, .omega..sub.11, X) .vertline.X.epsilon..OMEGA..sub.3 } is zero.

Equation (1) means that the smaller one of the values of the function max {.alpha.(.omega..sub.21, .omega..sub.11, X) .vertline.X.epsilon..OMEGA..sub.3 } and the value of .alpha.(.omega..sub.11, .omega..sub.21, .omega..sub.31) is left and taken as .alpha..sub.1 (.omega..sub.11, .omega..sub.21, .omega..sub.31). Therefore, if no peaks are present on the line L, then the relation .alpha..sub.1 (.omega..sub.11, .omega..sub.21, .omega..sub.31)=0 holds. In summary, by performing the processing given by equation (1), data about a new spectrum .alpha..sub.1 (.omega..sub.1, .omega..sub.2, .omega..sub.3) is obtained. In this new spectrum, those peaks which have no symmetrical partners have been erased.

At step 4, The obtained data .alpha..sub.1 (.omega..sub.1, .omega..sub.2, .omega..sub.3) is arithmetically processed, using max function given by equation (2), to give rise to data .alpha..sub.2 (.omega..sub.1, .omega..sub.2, .omega..sub.3).

.alpha..sub.2 (.omega..sub.1, .omega..sub.2, .omega..sub.3)=max {.alpha..sub.1 (.omega..sub.1, .omega..sub.2, .omega..sub.3), .alpha..sub.1 (.omega..sub.2, .omega..sub.1, .omega..sub.3)}(2) In this processing, in each 2D NMR spectrum given by a sliced plane vertical to the axis .omega..sub.3, data .alpha..sub.1 (.omega..sub.1, .omega..sub.2, .omega..sub.3) is compared with the data .alpha..sub.1 (.omega..sub.2, .omega..sub.1, .omega..sub.3) at the symmetrical position. The strength of the larger data is taken as .alpha..sub.2 (.omega..sub.1, .omega..sub.2, .omega..sub.3). The relation .alpha..sub.2 (.omega..sub.1, .omega..sub.2, .omega..sub.3 =c holds, as long as .alpha..sub.1 (.omega..sub.2, .omega..sub.1, .omega..sub.3)=c (c.noteq.0) holds, even if .alpha..sub.1 (.omega..sub.1, .omega..sub.2, .omega..sub.3)=0. If any peak has no symmetrical partner on the opposite side of the axis of symmetry (.omega..sub.1 =.omega..sub.2) in each 2D spectrum, the symmetrical peak of the same strength is supplemented by this processing.

At step 5, the obtained data .alpha..sub.2 (.omega..sub.1, .omega..sub.2, .omega..sub.3) represent a 3D NMR spectrum in which peaks having no symmetrical partners have been removed by the processing given by equation (1). In addition, in each individual 2D NMR spectrum constituting the 3D NMR spectrum, symmetrical partners are complemented. These two sets of processing given by equations (1) and (2) are next described in further detail.

We now discuss a sample consisting of a protein labeled with .sup.15 N. As shown in FIG. 8, this sample has long peptide chains. Adjacent protons .sup.1 H are spatially close to each other. Protons .sup.1 H are chemically coupled to nitrogen atoms .sup.15 N.

Three-dimensional NMR experiments are made on this sample with a pulse sequence as shown in FIG. 9. These 3D NMR experiments are a combination of 2D NOESY experiments for collecting information regarding the correlation of adjacent protons .sup.1 H and 2D HMQC experiments for collecting information regarding chemical coupling between protons .sup.1 H and nitrogen atoms .sup.15 N.

Where the sample consists of three peptide chains as shown in FIG. 8, a 3D NMR spectrum having three sliced planes as shown in FIG. 10 (a), (b), (c) is obtained. Two-dimensional NMR spectra appear on these sliced planes which are perpendicular to the frequency axis .omega..sub.3. In the past, it has been known that three peaks arising from .sup.15 N-.sup.1 H coupling are associated, by tracing symmetrical peaks along the broken lines from FIG. 10 (a). As already described, observation of numerous sliced planes along the axis .omega..sub.3 is needed for this purpose. Furthermore, in order to know that a peak indicated by .DELTA. in FIG. 10 (a) has no symmetrical partner, all other sliced planes must be checked.

On the other hand, in accordance with the present invention, spectra shown in FIG. 10 (a), (b), (c) are changed into spectra shown in FIG. 10 (a'), (b'), (c'), respectively, by the processing given by equation (1). It can be seen that the peak existing alone in FIG. 10 (a) has disappeared in FIG. 10 (a'). Spectra shown in FIG. 10 (a'), (b'), (c') are changed into spectra shown in FIG. 10 (a"), (b"), (c"), respectively, by the processing given by equation (2). It can be seen that symmetrical partners are complemented on the opposite side of the axis of symmetry in the sliced planes shown in FIG. 10 (a"), (b"), (c").

Where no complementary peaks exist as shown in FIG. 10 (a), (b), (c), one must search for the symmetrical peaks while estimating symmetrical locations. In accordance with the invention, it is only necessary to search other sliced planes for peaks which coincide in position with the complemented peaks. Hence, it is easy to search the other sliced planes for the symmetrical peaks precisely. In this case, vertical alignment of the proper peaks with the complemented peaks can be known by drawing the spectra carried by the sliced planes on transparent sheets, superimposing all the sheets, and watching the spectra from just or slantly above.

Having thus described our invention with the detail and particularity required by the Patent Laws, what is desired to be protected by Letters Patent is set forth in the following claims.

Claims

1. A three-dimensional NMR pulsed spectroscopy method comprising the steps for:

applying a preparation pulse or pulse sequence to a specimen,

sampling the resulting free induction decay signals with a sequence having two evolution periods and a detection period to obtain data S(t.sub.1, t.sub.2, t.sub.3) wherein t.sub.1 is the duration of of the first evolution period, t.sub.2 is the duration of the second evolution period and t.sub.3 is the time from the start of the detection period;

taking the Fourier transform of the data to create the three-dimensional NMR spectrum .alpha.(.omega..sub.1, .omega..sub.2, .omega..sub.3);

slicing the three-dimensional NMR spectrum perpendicularly to the frequency axis .omega..sub.3 to obtain a plurality of two-dimensional slice spectrums;

arithmetically processing data in each two-dimensional slice spectrum in such a way that peaks having no symmetrical partners at symmetrical positions in any two-dimensional slice spectrum are found and erased to obtain a new three-dimensional spectrum .alpha..sub.1 (.omega..sub.1, .omega..sub.2, .omega..sub.3);

adding peaks complementary to the peaks existing alone in each tow-dimensional slice spectrum to obtain a symmetrized three-dimensional spectrum .alpha..sub.2 (.omega..sub.1, .omega..sub.2, .omega..sub.3); and

reviewing the symmetrized spectrum for tracing symmetrical peaks indicative of specimen structure.

2. A three-dimensional NMR pulsed spectroscopy method comprising the steps for:

applying a preparation pulse or pulse sequence to a specimen;

sampling the resulting free induction decay signals with a sequence having two evolution periods and a detection period to obtain data S(t.sub.1, t.sub.2, t.sub.3) wherein t.sub.1 is the duration of the first evolution period, t.sub.2 is the duration of the second evolution period and t.sub.3 is the time from the start of the detection period;

taking the Fourier transform of the data to create the three-dimensional NMR spectrum .alpha.(.omega..sub.1, .omega..sub.2, .omega..sub.3);

slicing the three-dimensional NMR spectrum perpendicularly to the frequency axis .omega..sub.3 to obtain a plurality of two-dimensional slice spectrums;

arithmetically processing data in each two-dimensional slice spectrum according to the equation given by

.alpha..sub. (.omega..sub.1, .omega..sub.2, .omega..sub.3)=min [.alpha.(.omega..sub.1, .omega..sub.2, .omega..sub.3), max {.alpha.(.omega..sub.2, .omega..sub.1, X).vertline.X.epsilon..OMEGA..sub.3 }]

to obtain a spectrum .alpha..sub.1 (.omega..sub.1, .omega..sub.2, .omega..sub.3), where .OMEGA..sub.3 is the domain of frequency .omega..sub.3 ;

arithmetically processing data in each two-dimensional slice spectrum according to the equation given by

.alpha..sub.2 (.omega..sub.1, .omega..sub.2, .omega..sub.3)=max {.alpha..sub.1 (.omega..sub.1, .omega..sub.2, .omega..sub.3), .alpha..sub.1 (.omega..sub.2, .omega..sub.1, .omega..sub.3)}

to obtain a symmetrized three-dimensional spectrum .alpha..sub.2 (.omega..sub.1, .omega..sub.2, .omega..sub.3); and

reviewing the symmetrized spectrum for tracing symmetrical peaks indicative of specimen structure.