METHOD FOR SELECTIVELY DETECTING DOPAMINE BASED ON MAGNETIC RESONANCE NUCLEAR SPIN SINGLET STATE

Abstract

The present invention discloses a method for selectively detecting dopamine based on magnetic resonance nuclear spin singlet state. The method uses the nuclear spin singlet of three hydrogen atoms on the dopamine benzene ring to achieve selective detection of dopamine signals in a complicated system. The present invention is based on magnetic resonance technology to detect dopamine, has good accuracy, sensitivity and selectivity, can accurately detect the signal of dopamine from the complicated system, and the interference of signals of other substances are well eliminated. Meanwhile, the present invention further has the advantages of simple operation and non-intervention, can be used for monitoring the content and distribution of the dopamine in a living body, and has important application value in the fields of biology and medicine.

Claims

1. A method for selectively detecting dopamine molecules based on magnetic resonance nuclear (NMR) spin singlet comprising: step A: obtain the chemical shifts and J coupling values of the three .sup.1H spins on the benzene ring of dopamine in the test sample; step B: according to the spin characteristics of dopamine, design a pulse sequence that utilizes the nuclear spin singlet states to realize the selection of dopamine signals; step C: by the designed pulse sequence, prepare the spin singlet states of the three .sup.1H spin coupling system on the benzene ring of dopamine; during the existence of the spin singlets, the gradient field pulse is applied to remove the other signals except the spin singlet in the sample, resulting in the targeting selection of dopamine signal.

2. The method of claim 1, wherein the purpose of step A is to obtain the nuclear spin coupling characteristics of dopamine molecules by analyzing the NMR spectrum of the test sample.

3. The method of claim 1, wherein the designed pulse sequence of step B contains the preparation module of the spin singlets of the nuclear spin coupling system of dopamine, and the module to transform the spin singlets to the observable signals.

4. The method of claim 3, wherein, in the module, firstly a 90-degree hard pulse is applied along the x direction; after a delay time of τ.sub.1, a 180-degree hard pulse is applied along the x direction; then, a delay time of (τ.sub.1+τ.sub.2) is applied, followed by a 90-degree hard pulse along the y direction, and then a delay time of τ.sub.3 is applied; from the start of the pulse sequence to the end of τ.sub.3, this pulse sequence module is used to prepare the spin singlets of the three .sup.1H spin system on the benzene ring of dopamine; after τ.sub.3, a delay time of τ.sub.4 is given, followed by a 90-degree hard pulse along the y direction and then a delay time of τ.sub.5 is applied; this pulse sequence module is used to transform the spin singlets into the observable signals; in the designed pulse sequence, the values of τ.sub.1, τ.sub.2, τ.sub.3 can affect the efficiency of the singlet preparation of the three .sup.1H spins on the benzene ring of dopamine; in order to maximize the singlet preparation efficiency, optimization of τ.sub.1, τ.sub.2, τ.sub.3 is required; τ.sub.4 and τ.sub.5 are related to the dopamine signals for the final detection; in order to maximize the final detection signals, τ.sub.4 and τ.sub.5 need to be optimized too.

5. The method of claim 1, wherein the designed pulse sequence of step B comprises the pulse sequence module in which the gradient field pulses are applied during the evolution of the nuclear spin singlets to remove all other NMR signals except the nuclear spin singlets; from the start of the pulse sequence to the end of τ.sub.3 is the pulse sequence module for the singlet preparation; after this pulse sequence module, the first gradient field pulse, g.sub.1, is applied, then there is the singlet evolution stage; in the singlet evolution stage, a decoupling pulse can be applied; the purpose of the decoupling pulse is to minimize the influence of the environmental factors on the spin singlets evolution; to better remove all of the other NMR signals except the spin singlets, a second gradient field pulse g.sub.2 can be applied after the decoupling period.

6. The method of claim 1, wherein the method specifically comprises the following steps: steps 1: put the D.sub.2O aqueous solution with a mass fraction of 2%-5% dopamine in the magnetic resonance instrument, and apply a 90-degree hard pulse to the D.sub.2O aqueous solution of dopamine, which makes the .sup.1H signals excited, then obtain the .sup.1H spectrum of dopamine, and in turn the chemical shifts of three .sup.1H on the benzene ring of dopamine and the J coupling values among the protons; steps 2: according to the pulse form of preparing the singlets from the two-spin system of the weakly coupled system, for the three-spin system consisting of three .sup.1H atoms on the benzene ring of dopamine, the pulse parameters for the preparation and detection of dopamine singlet in the three-spin system are calculated based on chemical shifts and J coupling values obtained from step 1 by using MATLAB; as a result, the pulse sequence to prepare the nuclear spin singlet states of three .sup.1H atoms on the dopamine benzene ring with the maximum efficiency can be obtained; steps 3: the complete pulse sequence is obtained by combining the pulse form to prepare the singlet of the two-spin system of the weakly coupled system and the pulse parameters required for the preparation and detection of dopamine singlet in the three-spin system calculated in step 2, which is applied to the dopamine D.sub.2O aqueous solution to prepare and detect the singlet of the three-spin system consisting of three .sup.1H on the benzene ring; steps 4: on the basis of preparation and detection of dopamine singlet in step 3, two gradient field pulses with different amplitudes and a continuous wave (CW) decoupling pulse are applied between the pulses for the singlet preparation and the pulses for the signal detection to form a new pulse sequence; the function of the new pulse sequence can be divided into three parts: the first part is to obtain the singlets of the three .sup.1H on the benzene ring of dopamine; the second part is to keep the singlet states of three .sup.1H on the benzene ring of the dopamine and filter the other non-single state signals because the nuclear spin singlet state is not affected by the gradient field pulses and the CW decoupling pulse; the third part is to detect the singlet states of three .sup.1H on the benzene ring; in the end, only the three .sup.1H signals on the benzene ring are kept, achieving the purpose of selective signal filtering; in this process, it is necessary to continuously optimize the time of the CW pulse to achieve the best filtering efficiency.

7. The method of claim 6, as for the .sup.1H spectrum of dopamine described in step 1, the three signals on the left side represent the three .sup.1H signals on the benzene ring, the single peak in the middle is the water signal, and the ethyl signal of dopamine is on the right side; the J coupling values and chemical shifts between the three .sup.1H atoms on the benzene ring of dopamine are obtained from the .sup.1H spectrum of dopamine.

8. The method of claim 6, wherein step 2 is as follows: firstly, a 90-degree hard pulse is applied along x direction, after a delay time of τ.sub.1, a 180-degree hard pulse is applied along the x direction; then, a delay time of (τ.sub.1+τ.sub.2) is applied, followed by a 90-degree hard pulse along the y direction, and then a delay time of τ.sub.2/2 namely τ.sub.3 is applied; the function of this pulse is to prepare the singlet states of three .sup.1H on the benzene ring of dopamine, which is called preparation pulse; because the singlet states cannot be directly detected, another pulse is needed to detect the singlets of the three .sup.1H of three-spin system on the benzene ring of dopamine, which is called detection pulse; the form of detection pulse is as follows, after a delay time of τ.sub.4, followed by a 90-degree hard pulse along the y direction and then a delay time of τ.sub.5 is applied; next, take ADC sample immediately until the sampling signal decay is completed; in this process, the values of τ.sub.1 and τ.sub.2 impact the efficiency of the singlet states consisting of three .sup.1H atoms on the benzene ring of dopamine; in order to maximize the efficiency of the singlet states, MATLAB software is used to calculate the values of τ.sub.1 and τ.sub.2; first of all, 64 basic operators of the three-spin system are constructed in the MATLAB script, and then the Hamiltonian of the three-spin system consisting of three .sup.1H atoms on the benzene ring of dopamine is written; finally, the operating operators corresponding to the 90-degree hard pulse and the 180-degree hard pulse are obtained; then, the system is continuously evolving from the thermal equilibrium signals under the operating operator and Hamiltonian of the hard pulse, and the evolution time τ.sub.1 and τ.sub.2 are continuously optimized to maximize the singlet states preparation efficiency; similarly, on the basis of generating singlet states, evolution time τ.sub.4 and τ.sub.5 are optimized to maximize the singlet states detection efficiency; finally, the complete pulse to prepare and detect three-spin system singlet states of dopamine is obtained by combining the pulse form of the two-spin system singlet of the weakly coupled system and the calculated pulse parameters.

9. The method of claim 6, as described in step 3, the singlet state of the three-spin system consisting of three .sup.1H atoms on the benzene ring of dopamine is prepared and detected; specifically, firstly, the complete pulse obtained in step 2 is written into the computer by the NMR instrument language; secondly, a D.sub.2O aqueous solution of dopamine is put into the magnetic resonance instrument, and then the field-locking, field-shimming, matching, and tuning are performed; finally, the radio frequency center of the transmitter is set to the three .sup.1H on the benzene ring of dopamine, and the complete pulse written into the computer is applied to prepare and detect the singlet states of the dopamine.

10. The method of claim 6, as described in step 4, the CW decoupling pulse and two gradient field pulses with different amplitudes are applied between the pulses for the singlet preparation and the pulses for the signal detection; this forms a new pulse module; specifically, the duration of the CW pulse is between 50 ms and 1 s; the amplitude varies from 1 watt to 15 watts; the duration of each of the two gradient field pulses with different amplitudes along the z direction is between 1 ms and 5 ms, with the amplitude varying from 5 to 10 Gauss/cm; this new pulse module including the CW decoupling pulse and the two gradient field pulses with different amplitudes is written into the computer by the NMR instrument language; then, the routine procedures such as field-locking, field-shimming, matching, and tuning are performed; finally, the radio frequency center of the transmitter is set to the three .sup.1H on the dopamine benzene ring, and the new pulse sequence which has been written into the computer including the CW decoupling pulse and the two gradient field pulses with different amplitudes is applied to prepare and detect the singlet states of dopamine.

11. The method of claim 1, wherein the three .sup.1H spin coupling system on the benzene ring of dopamine is prepared into the nuclear spin singlets in step C; there are several multiple spin systems in a dopamine molecule; the spin system formed by the three .sup.1H spins on the benzene ring of dopamine can be effectively prepared for the nuclear spin singlets.

12. The method of claim 5, wherein the method specifically comprises the following steps: steps 1: put the D.sub.2O aqueous solution with a mass fraction of 2%-5% dopamine in the magnetic resonance instrument, and apply a 90-degree hard pulse to the D.sub.2O aqueous solution of dopamine, which makes the .sup.1H signals excited, then obtain the .sup.1H spectrum of dopamine, and in turn the chemical shifts of three .sup.1H on the benzene ring of dopamine and the J coupling values among the protons; steps 2: according to the pulse form of preparing the singlets from the two-spin system of the weakly coupled system, for the three-spin system consisting of three .sup.1H atoms on the benzene ring of dopamine, the pulse parameters for the preparation and detection of dopamine singlet in the three-spin system are calculated based on chemical shifts and J coupling values obtained from step 1 by using MATLAB; as a result, the pulse sequence to prepare the nuclear spin singlet states of three .sup.1H atoms on the dopamine benzene ring with the maximum efficiency can be obtained; steps 3: the complete pulse sequence is obtained by combining the pulse form to prepare the singlet of the two-spin system of the weakly coupled system and the pulse parameters required for the preparation and detection of dopamine singlet in the three-spin system calculated in step 2, which is applied to the dopamine D.sub.2O aqueous solution to prepare and detect the singlet of the three-spin system consisting of three .sup.1H on the benzene ring; steps 4: on the basis of preparation and detection of dopamine singlet in step 3, two gradient field pulses with different amplitudes and a continuous wave (CW) decoupling pulse are applied between the pulses for the singlet preparation and the pulses for the signal detection to form a new pulse sequence; the function of the new pulse sequence can be divided into three parts: the first part is to obtain the singlets of the three .sup.1H on the benzene ring of dopamine; the second part is to keep the singlet states of three .sup.1H on the benzene ring of the dopamine and filter the other non-single state signals because the nuclear spin singlet state is not affected by the gradient field pulses and the CW decoupling pulse; the third part is to detect the singlet states of three .sup.1H on the benzene ring; in the end, only the three .sup.1H signals on the benzene ring are kept, achieving the purpose of selective signal filtering; in this process, it is necessary to continuously optimize the time of the CW pulse to achieve the best filtering efficiency.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0031] FIG. 1 is the molecular structure of dopamine and the .sup.1H NMR spectrum of the dopamine aqueous solution sample.

[0032] FIG. 2 is the .sup.1H NMR spectrum of the three .sup.1H spins on the benzene ring of dopamine aqueous solution sample to prepare the dopamine singlet with the spectrum range from 6.61 ppm to 6.89 ppm.

[0033] FIG. 3 is the scheme of the pulse sequence for the singlet preparation of dopamine. The black rectangle represents a 90-degree hard pulse, and the white square represents a 180-degree hard pulse. x and y represent the phases of the radio frequency (RF) pulses. τ.sub.1-τ.sub.5 is the time for the spin evolution. In FIG. 3, the pulse sequence module from the start of the pulse sequence to the end of τ.sub.3 is the preparation module of the spin singlets of the nuclear spin coupling system of dopamine. The pulse sequence module from the beginning of τ.sub.4 to the end of τ.sub.5 is the module to transform the spin singlets to the observable signals.

[0034] FIG. 4 is the .sup.1H NMR spectrum obtained through preparing and detecting the singlet states of the dopamine in Example 1. The pulse sequence in FIG. 3 is used to acquire the spectrum.

[0035] FIG. 5 is the scheme of the pulse after adding the gradient field pulse and the high power CW decoupling pulse on the basis of preparing and detecting singlet. The selected pulse sequence of dopamine signals is realized by using the nuclear spin singlet. The black rectangle represents a 90-degree pulse and the white square represents a 180-degree pulse. x and y represent the phases of the radio frequency (RF) pulses. τ.sub.5 is the time for the spin evolution. τ.sub.m is the decoupling pulse time, and g.sub.1 and g.sub.2 represent the gradient field pulses.

[0036] FIG. 6 is the .sup.1H NMR spectrum of the dopamine aqueous solution sample. In Example 1 of the invention, after applying the gradient field pulses and high power CW decoupling pulse on the basis of the original pulse sequence, the .sup.1H NMR spectrum is obtained by preparing and detecting the spin singlets of dopamine; the pulse sequence in FIG. 5 is used to acquire the spectrum. The water signal in the spectrum was significantly suppressed.

[0037] FIG. 7 is the .sup.1H NMR spectrum of a mixture aqueous solution sample. The solute in this sample is a mixture consisting of dopamine, creatine, inositol, and glutamine.

[0038] FIG. 8 is the .sup.1H NMR spectrum of a mixture aqueous solution sample. The solute in this sample is a mixture consisting of dopamine, creatine, inositol, and glutamine. After using the pulse sequence which has adding the gradient field pulses and the high power CW decoupling pulse on the basis of the original pulse sequence, the .sup.1H NMR spectrum can be obtained by preparing and detecting the spin singlets of dopamine in the mixture aqueous solution sample; The pulse sequence in FIG. 5 is used to acquire the spectrum. In this spectrum, except the signals from the benzene ring of dopamine, all other signals are greatly suppressed.

[0039] FIG. 9 is the .sup.1H NMR spectrum is obtained by preparing and detecting the spin singlets of dopamine in the complex system in Example 2 of this invention.

DESCRIPTION OF THE PREFERRED EMBODIMENTS

[0040] The present invention will be further described in detail with the reference to the following specific examples and figures.

Example 1

[0041] NMR Instrument: Bruker AVANCE III 500 NMR spectroscopy.
The example was carried out as follows:

[0042] Step A: Analyze the .sup.1H NMR spectrum of dopamine and make the signal assignment. In the .sup.1H NMR spectrum of dopamine (see FIG. 1 and FIG. 2), if setting the center frequency between peaks 1 and 3 to zero, the resonance frequencies and J coupling values of peaks 1, 2 and 3 (i.e., the characteristics of the three .sup.1H spin system of dopamine) are:

ω.sub.1=36.5 Hz,ω.sub.2=−7.8 Hz,ω.sub.3=−36.5 Hz;

J.sub.12=0 Hz,J.sub.13=8.14 Hz,J.sub.23=2.18 Hz.

[0043] Step B: Design a pulse sequence that can be used to selectively detect the signals of dopamine molecules based on the nuclear spin singlets. According to the characteristics of the spin system of dopamine, the pulse sequence in FIG. 5 was designed. The core steps and ideas of the design of the pulse sequence are as follows: Firstly, the Pauli operator is used to construct the base operator of the three-spin system. Then according to the characteristics of the three .sup.1H spin system of dopamine obtained in Step A (i.e., the three .sup.1H resonance frequencies and the J coupling values), the spin Hamiltonians and the evolution operators of the three-spin system are as follows:

Ĥ=ω.sub.1Î.sub.1z+ω.sub.2Î.sub.2z+ω.sub.3Î.sub.3z+2πJ.sub.12Î.sub.1zÎ.sub.2z+2πJ.sub.13Î.sub.1zÎ.sub.3z+2πÎ.sub.2zÎ.sub.3z

û(τ)=e.sup.−iĤτ

û*(τ)=e.sup.iĤτ

[0044] wherein, ω.sub.1, ω.sub.2, ω.sub.3 represent the resonance frequencies of the three .sup.1H spins on the benzene ring of dopamine. Î.sub.1z,Î.sub.2z,Î.sub.3z are the base operators along the z axis of the three-spin system. τ is the evolution time, û(τ) is the evolution operator after τ time, û*(τ) is the transposed complex conjugate of ŷ(τ), and Ĥ is the spin system Hamiltonian.

[0045] The operator of the φ degree hard pulse along the x direction is as follows:

û.sub.x(φ)=e.sup.−i*(Î.sup.1x.sup.+Î.sup.2x.sup.+Î.sup.3x.sup.)*φ

û*.sub.x(φ)=e.sup.i*(Î.sup.1x.sup.+Î.sup.2x.sup.+Î.sup.3x.sup.)*φ

[0046] wherein, Î.sub.1x,Î.sub.2x,Î.sub.3x are the base operators along the x axis of the three-spin system, and û*.sub.x(φ) is the transposed complex conjugate of û.sub.x(φ).

[0047] The three-spin system is in the thermal equilibrium at room temperature, and thus the density operator is as follows:

{circumflex over (ρ)}.sub.0=Î.sub.1z+Î.sub.2z+Î.sub.3z

The evolution process of the density operator corresponding to the pulse sequence in FIG. 5 is as follows:
1. After applying a 90-degree hard pulse along the x direction, the density operator can be written as:

[00001]${\hat{\rho }}_1={\hat{u}}_x\left(\frac{\pi }{2}\right){\hat{\rho }}_0{\hat{u}}_x^{*}\left(\frac{\pi }{2}\right)$

2. After the evolution through time of τ.sub.1, the density operator becomes:

{circumflex over (ρ)}.sub.2=û(τ.sub.1){circumflex over (ρ)}.sub.1û*(τ.sub.1)

3. After applying a 180-degree hard pulse along the x direction, the density operator becomes:

{circumflex over (ρ)}.sub.3=û.sub.x(π){circumflex over (ρ)}.sub.2û*.sub.x(π)

4. After the evolution through time of τ.sub.1+τ.sub.2, the density operator becomes:

{circumflex over (ρ)}.sub.4=û(τ.sub.1+τ.sub.2){circumflex over (ρ)}.sub.3û*(τ.sub.1+τ.sub.2)

5. After applying a 90-degree hard pulse along the y direction, the density operator becomes:

[00002]${\hat{\rho }}_5={\hat{u}}_y\left(\frac{\pi }{2}\right){\hat{\rho }}_4{\hat{u}}_y^{*}\left(\frac{\pi }{2}\right)$

6. After the evolution through time of τ.sub.3, the density operator becomes:

{circumflex over (ρ)}.sub.6=û(τ.sub.3){circumflex over (p)}.sub.5û*(τ.sub.3)

[0048] After the evolution under the operation operators of Hamilton and the hard pulse, the target operator of the single state is:

[00003]$\hat{Q}=\frac{1}{\sqrt{3}}\left({\hat{I}}_{1x}{\hat{I}}_{3x}+{\hat{I}}_{1y}{\hat{I}}_{3y}+{\hat{I}}_{1z}{\hat{I}}_{3z}\right)+\frac{1}{\sqrt{3}}\left({\hat{I}}_{2x}{\hat{I}}_{3x}+{\hat{I}}_{2y}{\hat{I}}_{3y}+{\hat{I}}_{2z}{\hat{I}}_{3z}\right)$

[0049] The trace of the product of {circumflex over (ρ)}.sub.6 and the target operator {circumflex over (Q)} yields the singlet conversion efficiency, S:

S=trace({circumflex over (Q)}{circumflex over (ρ)}.sub.6)

[0050] To maximize the singlet conversion efficiency, τ.sub.1, τ.sub.2, and τ.sub.3 need to be optimized.

7. After the evolution through time of τ.sub.3, the gradient field pulse, g.sub.1, is applied to remove all other signals except the spin singlets in the sample. The spin singlets are not affected by the gradient field pulse, g.sub.1.
8. After the gradient field pulse, g.sub.1, a decoupling pulse with a period of τ.sub.m is applied to preserve the spin singlets in the sample.
9. After the decoupling pulse, the gradient field pulse, g.sub.2, is applied to further suppress all other signals except the spin singlets in the sample.

[0051] Because the spin singlets cannot be directly detected in the NMR instrument, it needs to be converted into the observable signals. The signal evolution of the system after the gradient field pulse, g.sub.2, is as follows:

10. After the evolution through time of τ.sub.4, the density operator becomes:

{circumflex over (ρ)}.sub.7=û(τ.sub.4){circumflex over (Q)}û*(τ.sub.4)

11. After applying a 90-degree hard pulse along the y direction, the density operator becomes:

[00004]${\hat{\rho }}_8={\hat{u}}_y\left(\frac{\pi }{2}\right){\hat{\rho }}_7{\hat{u}}_y^{*}\left(\frac{\pi }{2}\right)$

12. After the evolution through time of τ.sub.5, the density operator becomes:

{circumflex over (ρ)}.sub.9=û(τ.sub.5){circumflex over (Q)}û*(τ.sub.5)

[0052] If using the quantum state after the singlet preparation as the initial state, to have the best optimized singlet preparation efficiency, the target quantum state is:

{circumflex over (P)}=½(Î.sub.1xÎ.sub.3z+Î.sub.1zÎ.sub.3x)+½(Î.sub.2xÎ.sub.3z+Î.sub.3xÎ.sub.2z)

[0053] The trace of the product of {circumflex over (ρ)}.sub.9 and the target quantum state {circumflex over (P)} yields the projection of the final state on the target quantum state:

R=trace({circumflex over (P)}{circumflex over (ρ)}.sub.9)

[0054] wherein, R represents the observable signal transformed from the spin singlets. In order to achieve the maximum value of R, τ.sub.4 and τ.sub.5 need to be optimized.

[0055] Step C: Use the pulse sequence designed in Step B to selectively detect the signals of dopamine. The experimental parameters need to be optimized. The optimized experimental parameters in this experiment are as follows: τ.sub.1=30.9 ms, τ.sub.2=6.8 ms, τ.sub.3=3.4 ms, τ.sub.4=6.8 ms, τ.sub.5=3.3 ms. The gradient field pulse, g.sub.1, is applied along the z direction, the pulse length is 1 ms and the amplitude is 5 Gauss/cm. The direction, length, and power of the gradient pulse, g.sub.2, is the same as g.sub.1. A CW pulse is used as the decoupling pulse. The decoupling time is 100 ms, and the decoupling amplitude is 3 watts. The NMR spectrum in FIG. 6 is obtained by using the above parameters in the experiment. The sample is a dopamine deuterium aqueous solution (i.e. dopamine dissolved in D.sub.2O).

[0056] This embodiment comprises the following steps:

[0057] The instrument used in this embodiment is Bruker AVANCE III 500 NMR spectroscopy.

[0058] The specific steps from the singlet preparation of dopamine to the realization of the filtering of dopamine signals are as follows:

[0059] Step 1: This embodiment is a system in which dopamine is dissolved in D.sub.2O. Preparing an aqueous solution of dopamine with a mass fraction of 3%. In the NMR experiment, a 90-degree hard pulse is applied, obtaining the signal peaks of the dopamine .sup.1H spectrum shown in FIG. 1. Among them, the peaks of three protons on the dopamine benzene ring are at the far left, and the water peak is at the middle of the spectrum, the signal on the right is the ethyl signal of the dopamine molecule. Selecting the three .sup.1H peaks of the benzene ring in the dopamine structure as the target peaks (The labeled peaks 1, 2, and 3 on the .sup.1H spectrum is shown in FIG. 2). After setting the center frequency between peaks 1 and 3 to zero, the J coupling values between the protons of the system and their corresponding chemical shifts read as:

ω.sub.1=36.5 Hz,ω.sub.2=−7.8 Hz,ω.sub.3=−36.5 Hz;

J.sub.12=0 HZ,J.sub.13=8.14 Hz,J.sub.23=2.18 Hz.

[0060] Step 2: Firstly, the Pauli operator is used to construct the base operator of the three-spin system. Then according to the chemical shifts and the J coupling values of the three .sup.1H spin system of dopamine obtained in Step 1, the spin Hamiltonians and the evolution operators of the three-spin system are as follows:

Ĥ=ω.sub.1Î.sub.1z+ω.sub.3Î.sub.2z+ω.sub.3Î.sub.3z+2πJ.sub.12Î.sub.1zÎ.sub.2z+2πJ.sub.13Î.sub.1zÎ.sub.3z+2πJ.sub.23Î.sub.2zÎ.sub.3z

û(τ)=e.sup.−iĤτ

û*(τ)=e.sup.iĤτ

[0061] Wherein, ω.sub.1, ω.sub.2, ω.sub.3 represent the values of the chemical shift of the three 41 spins on the benzene ring of dopamine respectively. Î.sub.1z,Î.sub.2z,Î.sub.3z are the base operators along the z axis of the three-spin system; T is the evolution time, û(τ) is the evolution operator after τ time, û*(τ) is the transposed complex conjugate of û(t), and Ĥ is the Hamiltonian of the spin system.

[0062] Shown as an example, the operator of the 90-degree hard pulse along the x direction is as follows:

[00005]${\hat{u}}_x\left(\frac{\pi }{2}\right)=e^{-i*\left({\hat{I}}_{1x}{\hat{I}}_{2x}+{\hat{I}}_{3x}\right)*\frac{\pi }{2}}$${\hat{u}}_x^{*}\left(\frac{\pi }{2}\right)=e^{i*\left({\hat{I}}_{1x}{\hat{I}}_{2x}+{\hat{I}}_{3x}\right)*\frac{\pi }{2}}$

[0063] Wherein, Î.sub.1x,Î.sub.2x,Î.sub.3x are the base operators along the x axis of the three-spin system, and

[00006]${\hat{u}}_x^{*}\left(\frac{\pi }{2}\right)$

is the transposed complex conjugate of

[00007]${\hat{u}}_x\left(\frac{\pi }{2}\right).$

[0064] The three-spin system is in the thermal equilibrium at room temperature, and thus the density operator is as follows:

{circumflex over (ρ)}.sub.0=Î.sub.1z+Î.sub.2z+Î.sub.3z

[0065] The pulse form of the two-spin system singlet prepared in the weakly coupled system is as shown in FIG. 3b

1. After applying a 90-degree hard pulse along the x direction, the density operator can be written as:

[00008]${\hat{\rho }}_1={\hat{u}}_x\left(\frac{\pi }{2}\right){\hat{\rho }}_0{\hat{u}}_x^{*}\left(\frac{\pi }{2}\right)$

2. After the evolution through time of τ.sub.1, the density operator becomes:

{circumflex over (ρ)}.sub.2=û(τ.sub.1){circumflex over (ρ)}.sub.1û*(τ.sub.1)

3. After applying a 180-degree hard pulse along the x direction, the density operator becomes:

{circumflex over (ρ)}.sub.3=û.sub.x(π){circumflex over (ρ)}.sub.2û*.sub.x(π)

4. After the evolution through time of τ.sub.1+τ.sub.2, the density operator becomes:

{circumflex over (ρ)}.sub.4=û(τ.sub.1+τ.sub.2){circumflex over (p)}.sub.3û*(τ.sub.1+τ.sub.2)

5. After applying a 90-degree hard pulse along the y direction, the density operator becomes:

[00009]${\hat{\rho }}_5={\hat{u}}_y\left(\frac{\pi }{2}\right){\hat{\rho }}_4{\hat{u}}_y^{*}\left(\frac{\pi }{2}\right)$

6. After the evolution through time of τ.sub.3, the density operator becomes:

{circumflex over (ρ)}.sub.6={circumflex over (u)}(τ.sub.3){circumflex over (ρ)}.sub.5{circumflex over (u)}*(τ.sub.3)

Wherein,

[0066] [00010]${\tau }_3=\frac{{\tau }_2}{2}$

[0067] The function of this pulse is to prepare the singlet state of the three-spin system consisting of three .sup.1H on the benzene ring of dopamine. It can be named as the preparation pulse for short.

[0068] After the evolution under the operation operators of Hamilton and hard pulse, the targeted operator of the singlet state is:

[00011]$\hat{Q}=\frac{1}{\sqrt{3}}\left({\hat{I}}_{1x}{\hat{I}}_{3x}+{\hat{I}}_{1y}{\hat{I}}_{3y}+{\hat{I}}_{1z}{\hat{I}}_{3z}\right)+\frac{1}{\sqrt{3}}\left({\hat{I}}_{2x}{\hat{I}}_{3x}+{\hat{I}}_{2y}{\hat{I}}_{3y}+{\hat{I}}_{2z}{\hat{I}}_{3z}\right)$

[0069] After the evolution under the operation operators of hard pulse and Hamilton, the trace of the product of {circumflex over (ρ)}.sub.6 and the target operator {circumflex over (Q)} yields the projection on the singlet state:

S=trace({circumflex over (Q)}{circumflex over (ρ)}.sub.6)

[0070] Wherein, S represents the conversion efficiency of the singlet state. In order to maximize the singlet state conversion efficiency, it is needed to use MATLAB to perform a simulation to obtain the appropriate parameters. Firstly, writing the above formulas into MATLAB scripts respectively; then continuously optimize the values of τ.sub.1, τ.sub.2 through programming, so that the singlet state conversion efficiency (the absolute value of S) reaches the optimal value. After optimizing, τ.sub.1=30.9 ms, τ.sub.2=6.8 ms, and the singlet efficiency can reach the optimal value.

[0071] Because the single state is not a single quantum signal, and the spin singlet cannot be directly detected in the NMR instrument. It is necessary to apply the detection pulse in order to detect the singlet state signals. The detection pulse form is as follows:

1. After the evolution through time of τ.sub.4, the density operator becomes:

{circumflex over (ρ)}.sub.7=û(τ.sub.4){circumflex over (Q)}û*(τ.sub.4)

2. After applying a 90-degree hard pulse along the y direction, the density operator becomes:

[00012]${\hat{\rho }}_8={\hat{u}}_y\left(\frac{\pi }{2}\right){\hat{\rho }}_7{\hat{u}}_y^{*}\left(\frac{\pi }{2}\right)$

3. After the evolution through time of τ.sub.5, the density operator becomes:

{circumflex over (ρ)}.sub.9=û(τ.sub.5){circumflex over (Q)}û*(τ.sub.5)

[0072] If using the quantum state after the singlet preparation as the initial state, to have the best-optimized singlet preparation efficiency, the target quantum state is:

{circumflex over (P)}=½({circumflex over (I)}.sub.1x{circumflex over (I)}.sub.3z+{circumflex over (I)}.sub.1z{circumflex over (I)}.sub.3x)+½({circumflex over (I)}.sub.2x{circumflex over (I)}.sub.3z+Î.sub.3x{circumflex over (I)}.sub.2z)

[0073] Similarly, after the detection pulse, the trace of the product of {circumflex over (ρ)}.sub.9 and the target quantum state {circumflex over (P)} yields the projection of the final state on the target quantum state:

R=trace({circumflex over (P)}{circumflex over (ρ)}.sub.9)

[0074] Wherein, R represents the conversion efficiency of the target quantum state. In order to achieve the maximum value of R, τ.sub.4 and τ.sub.5 need to be continuously optimized by MATLAB to maximize the absolute value of R. After optimizing, τ.sub.4=3.3 ms,τ.sub.5=6.8 ms.

[0075] After calculation, the complete pulse sequence to prepare and detect the three-spin system of dopamine is obtained by combining the pulse form of the two-spin system singlet under the weakly coupled system and the calculated pulse parameters.

[0076] Step 3: The complete pulse obtained in Step 2 is written into the computer by NMR instrument language; then, a D.sub.2O aqueous solution of dopamine is put into the NMR instrument, and then the field-locking, field-shimming, matching, and tuning are performed; finally, the radio frequency center of the transmitter is set to the three .sup.1H on the benzene ring of dopamine, and the complete pulse written into the computer is applied to prepare and detect the singlet states of dopamine. The detected NMR spectrum of dopamine is shown in FIG. 4. The three signals on the left are from the three protons on the dopamine benzene ring, and the right is the water signal in deuterium water. It can be seen that after the singlet preparation and detection, the signal intensity of water in deuterium water is still far greater than the signal of dopamine, and selective filtering has not been achieved. Therefore, it is needed to be improved on the basis of preparing single states.

[0077] The pulse sequence for preparing the spin singlets of dopamine can be designed according to actual needs, and different preparation singlet pulse sequences can be used to realize the singlet preparation of dopamine.

[0078] Step 4: Based on the preparation and detection of singlet state in Step 3, a new high power CW pulse and two gradient field pulses with different amplitudes are applied between the pulses for the singlet preparation and the pulses for the signal detection, these pulses form are shown in FIG. 5. The pulses before b are the preparation module of the dopamine singlet. The CW pulse and two gradient field pulses with different amplitudes are applied in the period between b and c. The pulses after c are the module to detect the singlet of dopamine. Wherein, the duration of the CW pulse is 100 ms, the amplitude is 3 watts. Two gradient field pulses with different amplitudes require a duration of 1 ms, an amplitude of 5 Gauss/cm along the z direction. In the process, the duration and amplitude of CW are needed to be optimized in order to achieve the best filtering effect, that is, filtering out non-singlet signal components while keeping the three .sup.1H signals from the benzene ring of dopamine. The experimental result is shown in FIG. 6. The three signals on the left are from the three .sup.1H on the dopamine ring, and the signal on the right is the water signal in deuterium water. The signal intensity of water in deuterium water is much smaller than the signal of dopamine molecules after applying the CW decoupling pulse and two gradient field pulses with different amplitudes between the preparation module of the spin singlets and the module of spin singlet detection. This achieves selectively filtering.

Example 2

[0079] NMR Instrument: Bruker AVANCE III 500 NMR spectroscopy.
The example was carried out as follows:

[0080] Step A: Prepare an aqueous solution of dopamine, creatine, inositol, and glutamine with a mass fraction of 1.5% for each solute. D.sub.2O instead of H.sub.2O is used to prepare the solution. The .sup.1H NMR spectrum of the sample is shown in FIG. 7. Analyze the .sup.1H NMR spectrum and make the signal assignment. The leftmost signals are from the three protons on the benzene ring of dopamine. The middle signal at 4.7 ppm is assigned to the water (HDO). The .sup.1H signals in the right side of the spectrum are from the ethyl group of dopamine, and the groups of creatine, inositol, and glutamine. If setting the center frequency between peaks 1 and 3 to zero (see FIG. 7), the resonance frequencies and J coupling values of peaks 1, 2 and 3 are:

ω.sub.1=36.5 Hz,ω.sub.2=−7.8 Hz,ω.sub.3=−36.5 Hz;

J.sub.12=0 Hz,J.sub.13=8.14 Hz,J.sub.23=2.18 Hz.

[0081] Step B: Because the resonance frequencies and J coupling values of peaks 1, 2 and 3 have not been changed in the complex system, and the spin Hamiltonian thus remains unchanged too. The pulse sequence in FIG. 5 was directly used in the experiments in this example.

[0082] Step C: This step is similar to Step C in Example 1. The parameters used in the experiments are the same as those used in Example 1. The spectrum in FIG. 8 was obtained.

[0083] This embodiment comprises the following steps:

[0084] The instrument used in this embodiment is Bruker AVANCE III 500 NMR spectroscopy.

[0085] The specific steps from the preparation of the spin singlet of dopamine to the realization of the signal filtering are as follows:

[0086] Steps 1: Prepare an aqueous solution of dopamine, creatine, inositol, and glutamine with a mass fraction of 1.5% for each solute. Put the solution sample into an NMR instrument. After applying a 90-degree hard pulse, the signals of dopamine as shown in the .sup.1H spectrum in FIG. 7 can be obtained. In the spectrum, the leftmost signals are the three .sup.1H on the benzene ring of dopamine. The middle signal at 4.7 ppm is assigned to the water (HDO). The .sup.1H signals in the right side of the spectrum are from the ethyl of dopamine, and the groups of creatine, inositol and glutamine. The J coupling values between the protons of the system and the corresponding chemical shifts read as:

ω.sub.1=36.5 Hz,ω.sub.2=−7.8 Hz,ω.sub.3=−36.5 Hz;

J.sub.12=0 Hz,J.sub.13=8.14 Hz,J.sub.23=2.18 Hz

Step 2: Because the chemical shifts and J coupling values of peaks 1, 2 and 3 have not been changed in the complex system, and the spin Hamiltonian thus remains unchanged too. The complete pulse sequence of Step 2 of Example 1 can be directly used, which is composed of the pulse sequence for preparing the singlet state of the two-spin system and the calculated pulse parameters under the weakly coupled system.

[0087] Step 3: Similar to Step 3 in Example 1, applying the pulse, which has been written into the computer, to prepare and detect the spin singlet of dopamine. The detected dopamine NMR spectrum is shown in FIG. 9. The intensity of the signals of the three protons from the benzene ring of the dopamine molecule at the far left barely changes. But the phase of the signals changes. The signal at 4.7 ppm is the water signal, which has been suppressed greatly. The .sup.1H signals in the right side of the spectrum are from the ethyl of dopamine, the groups of creatine, inositol and glutamine. The intensity of these signals changed. But compared with that of dopamine, these signals are still higher than those of dopamine. In general, compared with the .sup.1H NMR spectrum (FIG. 7), the intensity of the other non-singlet signals in the spectrum is suppressed to a certain extent, but the effect is not obvious except for the water signal.

[0088] Step 4: Based on the preparation and detection of singlet state in Step 3, a new high power CW pulse and two gradient field pulses with different amplitudes are added between the preparation module of the spin singlets and the module of spin singlet detection. These pulses are shown in FIG. 5. The pulses before b belong to the preparation module of the spin singlets of dopamine. The pulses between b and c are the newly applied CW pulse and the two gradient field pulses with different amplitudes. The pulses after c are the module of spin singlet detection. Wherein, the duration of the CW pulse is 100 ms, the amplitude is 3 watts. Two gradient field pulses with different amplitudes require a duration time of 1 ms, with an amplitude of 5 Gauss/cm along the z direction. The results are shown in FIG. 8. The .sup.1H signals which are from the benzene ring of dopamine at the far left of the spectrum show a clear decay, and meanwhile, the phases of the signals change. The signal at 4.7 ppm is the water signal and the signal intensity is almost 0. Calculation shows that compared with the water signal in the .sup.1H spectrum, the height of the residual water signal is less than 0.1%. The .sup.1H signals in the right side of the spectrum are from the ethyl of dopamine, the groups of creatine, inositol and glutamine. The intensities of these signals are obviously changed. Comparing with the .sup.1H signals which are from the benzene ring of dopamine, their intensities are almost negligible. In summary, the other non-singlet signals in the spectrum are almost completely suppressed in the signal intensity compared with the .sup.1H spectrum (FIG. 9). In particular, the signal of water is almost attenuated to 0 in intensity, and only the signals from the singlets of dopamine are retained, realizing the selective detection of the dopamine signals in this mixed system.

[0089] The present invention discloses a method for selectively detecting dopamine molecules based on the magnetic resonance nuclear spin singlet. The method utilizes the nuclear spin singlet of the three .sup.1H spins on the benzene ring of dopamine to realize the selective detection of dopamine signals in complex systems. The present invention detecting dopamine by the magnetic resonance technique, has good accuracy, sensitivity and selectivity, and can accurately detect dopamine signals from a system having complex components, and meanwhile can well eliminate the interference of signals from the other substances. At the same time, the present invention has the advantages of simple operation and non-intervention, and thus can be used to monitor the amount and distribution of the dopamine molecules in a living body. Thus this invention will have important applications in the fields of biology and medicine.

[0090] The content described above is only the preferred embodiments of the present invention. The skilled persons in the field may modify or equivalently make some substitution within the spirit of the present invention. All changes made in accordance with the spirit of the present invention should fall within the scope of protection claimed by the present invention.