60/210,040, filed June 7,2000, the contents of which are hereby incorporated herein by reference.

FIELD OF THE INVENTION A method of diagnosing a patient using a noninvasive method that allows for localization specification of treatment and following monitoring of brain thalamocortical dysrythmia.

BACKGROUND OF THE INVENTION Neuronal rhythmicity, ensemble neuronal oscillation, and resonance are closely related to the emergence of brain functions (Llinas, Science, 242: 1654-1664,1998; Llinas et al., Temporal Coding in the Brain, Springer, Berlin, 251-272,1994). Moreover, prior studies have shown that there is a link between high frequency neuronal oscillations and sensorimotor and cognitive functions. See, for example, Ahissar and Vaadia, Proc. Natl. Acad. Sci. USA, 87: 8935- cognitive functions. See, for example, Ahissar and Vaadia, Proc. Natl. Acad Sci. USA, 87: 8935- 8939, 1990; Eckhorn et al., Biol. Cybern., 60: 121-130,1988; Llinas and Parc, Neurosci., 44: 521- 535,1991; and Llinas, et al., Philos. Trans. R. Soc. London B, 353: 1841-1849.

Classical neurology has suggested that damage to the cerebral cortex leads to a variety of conditions and syndromes. The development of neurological conditions and syndromes is dependant upon the localization of the cortical damage. For example, damage to specific locations in the visual cortex may lead to different types of blindness depending on the insult location. Similar findings are also associated with other cortical structures, including the somatosensory, motor, and pre-motor corticies. Combined these studies have suggested that the majority of cognitive functions arise from coherent electrical actions that occur at the cortical level of the brain. An opposing view considers is that cognitive function is due to coherent events in thalamocortical re-entry dynamics and that neuropsychiatric conditions and syndromes may be due to altered electrical activity in these regions.

Treatment of patients with neurological disorders (such as Parkinson's disease, neurogenic pain, tinnitus, and depression) have historically been pharmacological in nature and guided by empirical findings. Even under such treatments, many patients showed drug resistance. Surgical methods have required cortical ablation. Presently, stereotactic cortical or thalamic lesion techniques have been used.

Patients with Parkinson's disease required a stereotactic medial thalamotomy.

Briefly, a stereotactic frame, which was compatible with a magnetic resonance image (MRI), allowed localization of the medial thalamic targets. An electrode was inserted and positioned according to a computer-calculated target. Lesions were then made to the medial thalamus region. This method also may be used to diagnose and classify patients with neurological disorders into a specific disease state by analyzing the MRI data that is obtained and correlating any observed altered brain activity to brain activity associated with a specific disease state.

Despite the current use of this treatment method for a variety of neurological disorders, a medial thalamotomy requires surgery and insertion of an electrode into the brain.

Minor movements of the implanted electrode could affect brain regions other than intended.

Furthermore, the use of this method to diagnose an individual requires a procedure that is risky and highly invasive in nature.

SUMMARY OF THE INVENTION The present invention describes a non-invasive method of diagnosing patients with a neuropsychiatric disease. The method includes placing a patient's head in a whole-head magnetoenchalography recording system; collecting spontaneous brain activity data, for 2 minutes to 1 hour depending on the disease state, from said patient, whose eyes are closed; analyzing the spontaneous brain activity data; and correlating the brain activity of said patient to the brain activity associated with a neuropsychiatric disorder to diagnose the patient.

A specific embodiment of the invention describes the method of diagnosing patients with a specific neuropsychiatric diseases, including neurogenic pain, obsessive compulsive disorder, depression, panic disorder, Parkinson's disease, schizophrenia, rigidity, dystonia, tinnitus, and epilepsy.

The present invention also describes a non-invasive method for the planning of treating patients with a neuropsychiatric disease. The method includes placing a patient's head in a whole-head magnetoenchalography recording system; collecting spontaneous brain activity data, for 10 minutes, from said patient, whose eyes are closed; analyzing the spontaneous brain activity data; and correlating the brain activity of said patient to the brain activity associated with a neuropsychiatric disorder to diagnose the patient; and specifying treatment procedures for the patient with the appropriate pharmacological or electrical treatments for the diagnosed neuropsychiatric disorder.

A specific embodiment of the invention describes the method of diagnosing, specifying treatments for, and monitoring follow-ups of patients with a specific neuropsychiatric disease, including neurogenic pain, obsessive compulsive disorder, depression, panic disorder, Parkinson's disease, schizophrenia, rigidity, dystonia, tinnitus, and epilepsy.

BRIEF DESCRIPTION OF THE DRAWINGS Figs. 1 a-c are diagrams of the cortical and thalamic anatomy; Fig. 2 is a diagram of the corticothalamic pathway interactions; Fig. 3 is a block diagram of the Thalamocortical Dysrhythmia Treatment process; Figs. 4a-e are graphs of the Power Spectra of Control Subjects and Patients; Figs. 5a-c are graphs of the Power Spectrum Versus Power Ratio of Control Subjects and Patients; Fig. 6 is a graph illustrative of the Power Spectrum Correlation Regions; Fig. 7 is a graph of the Power Spectrum Correlation for a Control Subject; Fig. 8 is a graph of the Power Spectrum Correlation for a Patient with Psychosis; Fig. 9 is a graph of the Power Spectrum Correlation for a Patient with OCD; Fig. 10 is a graph of the Power Spectrum Correlation for a Patient with Depression; Fig. 11 is a graph of the Power Spectrum Correlation for a Patient with Neuropathic Pain; Fig. 12 is a graph of the Power Spectrum Correlation for a Patient with Parkinson's ; Fig. 13 is a graph of the Power Spectrum Correlation for a Patient with Tinnitus; Fig. 14 is a graph of the Power Spectrum of a Patient Pre and Post Treatment ; and Figs. 15a-d are graphs of the Power Spectrum Correlation of a Patient Pre-and Post-Treatment.

DETAILED DESCRIPTION OF THE INVENTION The present invention advantageously uses the analysis of brain activity that is associated with specific neuropsychiatric disease states to diagnose and treat patients more effectively.

The diagnosing and treating of a patient for a neuropsychiatric disorder is obtained using a noninvasive method. This may be done by using any machine that can collect spontaneous brain activity which may then be transformed into data that can be analyzed for specific brain activity patterns. Preferably, the machine is a whole-head magnetoenchepalography recording system that can scan brain activity over several channels.

More preferably, the machine is a whole-head, 148-channel MEG system Magnes 260 WH (4D-Neuroimaging, San Diego, CA).

The collection of brain activity may occur for about 0-60 minutes with the patient's eyes open or closed. The time course of collection of brain activity will be dependent upon the disease state is being diagnosed in the patient. The patient may subjected to specific stimuli or perform particular cognitive functions during the collection of the brain activity.

Preferably, the collection of brain activity of the patient occurs with the patient's eyes closed.

More preferably, the collection of brain activity of the patient occurs for about 10 minutes with the patient's eyes closed.

The brain activity is then transferred to a computer system for analysis. This computer system may be a laptop, handheld, or desktop computer. The computer system may be used to analyze the information or be used as a storage medium. If the computer is a storage medium, then the data is stored until it can be transferred to a computer system that will analyze the information. Preferably the computer system that will analyze the data will be a desktop computer system. Most preferably, the computer system that will analyze the data will be a LINUX cluster computer system.

The computer system that will analyze the data may contain any computer, statistical, or mathematical modeling program that will sufficiently convert the collected brain activity to a format that can be used to correlate brain activity to specific neuropsychiatric diseases. Preferably, the computer system will analyze the brain activity to determine the spectral estimation of the brain activity. Preferably, the computer system will analyze the data with the computational methods that are described in the preferred embodiments.

The correlation of altered brain activity in specific brain regions to a particular disease state may be completed by a computer system or a person. The localization of the altered brain activity may define the specific disease state of the patient. Disease states that may be diagnosed and/or treated with this invention include, but are not limited to, neurogenic pain, obsessive compulsive disorder, depression, panic disorder, Parkinson's disease, schizophrenia, rigidity, dystonia, tinnitus, and epilepsy.

Preferably, altered brain activity in the cingulate cortex will be correlated to the diagnosis of neuogenic pain, obsessive compulsive disorder, depression, and panic disorders.

Preferably, altered brain activity in the medial dorsal cortex may be correlated to the diagnosis of schizophrenia.

Preferably, altered brain activity in areas 4 and 6 of the motor cortex may be associated with Parkinson's disease, rigidity, and dystonia.

Preferably, altered brain activity in the auditory cortex and medial geniculate nucleus may be associated with tinnitus.

In treatment of a patient, once any altered brain activity is associated with a neurological disease state the patient may be given any one of several treatment options that are appropriate for that disease state. These treatments include, but are not limited to, drug therapy, surgery, and electrical modification of the brain cells with the altered activity.

Drug therapy includes any pharmacological or toxicological compound that may ameliorate the symptoms associated with the disease state or any compound that may cure the disease state. Drugs may be administered as a tablet, capsule, caplet, suppository, injection, or intravenous fluid. The drug may be administered as many times per day as necessary to produce a therapeutically effective level of the drug in the patient's blood stream.

Surgery includes, but is not limited to, the removal of the cells that are producing the altered brain activity or addition of new cells to compensate for the defective cells.

Electrical modification includes, but is not limited to, implantation of electrodes into specific brain regions and stimulation or inhibition of specific brain regions to modify the brain activity. The electrodes may be implanted for any time period. The oscillations and the frequency of the implanted probe will vary depending on the diagnosed disease state, the severity of the disease, and The system allows the monitoring of a patient during or following treatment with the appropriate pharmacological agents or using surgical procedures using non-invasive techniques, preferably MEG recordings of brain activity.

As used herein, the term"noninvasive"refers to a method that does not involve penetration, by surgery or needle, of the outer layer of an intact organism. For example, a noninvasive method would not penetrate the skin of a mammal, preferably a human.

PREFERRED EMBODIMENTS Subjects Spontaneous brain activity was obtained from nine healthy control subjects and nine test patients with chronic, severe, and therapy resistant neurological or neuropsychiatric disorders, including Parkinson's disease, tinnitus, neurogenic pain, and major depression.

Control subjects ranged between 24 and 45 years old, test subjects ranged between 28 and 73 years old.

Magnetoencephalography Recordings and Analysis Magnetoencephalography (MEG) recordings were obtained with a whole-head, 148-channel MEG system Magnes 2500 WH (4D-Neuroimaging, San Diego, CA). Spontaneous brain activity was continuously recorded for 10 minutes (bandpass: 0.1-100 Hz; sample rate: 508 Hz) while the subject rested with eyes closed. Auxiliary EKG measurements were recorded simultaneously. Continuous MEG raw data were analyzed on a LINUX xluster computer system using developed software and commercial MATLAB analysis software packages.

Computational Methods Bias and variance in spectral estimation was resolved by averaging over a set of orthogonal basis functions, such as the Slepian sequences (Slepian and Pollak, Biophys., 40: 43- 63,1961). These sequences [w* (t)] form a sequence of orthogonal functions, defined on the time interval where t = 1 to T. The functions are parameterized by a bandwidth parameter W, so that there are K = [2WT] basis functions with spectra that are confined to a frequency, f, band of [f- W, f + W] around the frequency of interest.

The basic quantity involved in the spectral estimation, for a given data sequence, is the tapered Fourier transformation: Using these quantities, a direct estimation of the spectrum was determined by the equation: The cross-correlation between spectral amplitudes at different frequencies were obtained by computation of the multipaper spectrogram, SMT (F, t) and using a slow moving analysis window. The correlation coefficient of two time series, log (SMT (fol, t)) and log (SMTCf2, t)), with their means removed and wheref, andf2 represent points in a two dimensional grid in"2 space. Using this computational method, a two-dimensional image of the spectral correlation can be generated.

Noise from cardiac and distant sources of magnetic fields was removed from the raw data prior to performing any spectral analysis calculations. To remove cardiac magnetic fields from the raw data generated by the MEG, time-domain modeling of the cardiac spikes on a channel-by-channel basis were performed. Using the auxiliary EKG measurement, where spike times were extracted by using a matched filter to the spike shapes, spike times were used to extract interpolated spike shapes from the channels. The shapes were used to create smooth templates that were then subtracted from the raw MEG data. Reference channels were subtracted in the frequency domain by using a linear model for the contribution of the reference channels to the MEG data. The coefficients of the linear model could be estimated by using a multipaper transfer function estimate. For example, if y (t) was a data channel and x (t) was a reference channel, the subtracted spectral coefficient to be used in the spectral estimate could be defined as: k (/') = (/)- (/) (/). where Results Plots of the overall frequency content for the rostral and caudal halves of the brains showed of control and patient subjects showed that the peak frequencies. (Figures 1A and 1B). The average power spectra from all the patient subjects and control subjects are shown in (Figures 2A and 2B). The aggregate of all the frequency channels of the control and test subjects are shown in Figure 3. The total power in the 5 to 15 Hz band against the power ratio between 5 to 10 Hz band and between the 10 to 15 Hz band are shown in Figure 4. Correlation plots of power spectra over a period of 10 minutes for a control subject and a test subject with Parkinson's disease are shown in Figures 5A and 5B. Correlation plots of the average power spectra over a period of 10 minutes for control and test subjects are shown in Figures 5C and 5D.

Discussion The overall frequency data in Figures 4a-e show that in test patients, there was a shift from the normal a rhythm to a low-frequency 0 rhythmicity. The shift was less prominent in the rostral pole, compared to the caudal pole, which is in agreement with prior studies that indicate that in ceratin circumstances 0 rhytmicity is present in normal individuals (Sasaki, et al., Cognit. Brain Res., 5: 165-174,1996). In addition, the average power spectra of Figures 4a-e also show a decrease in the a power and an increase in the lower-frequency activity in the 0 range. Furthermore, these was an increase in the global power, which is indicative of an overall increased coherence in the test patients. These figures indicate that in those subjects with neuropsychiatric disorders there are two distinct characteristics that distinguish the low- frequency thalamocortical activity in test patients from the normal 0 rhythmicity in control subjects: (i) the persistent low-frequency thalamocortical resonance present during the awake state, and (ii) the wide coherence over the recorded channels.

Figures 5a-c, which show the relationship of the total power in the 5 to 15 Hz band to the power ratio between 5 to 10 Hz band and between the 10 to 15 Hz band supports the above statement. The figure shows that control subjects tend to be clustered in spaces with high- frequency and lower global power. Comparatively, the test patients were clustered in regions with lower-frequency and higher global power.

Figures 6-13 show individual and average correlation plots of power spectra over a period of 10 minutes for control subject and test patients. The plots indicate that the increase in O activity in the test patients was in accordance with the presence of low-threshold spike bursting activity. Furthermore, the plots indicated a correlation with harmonics at the y frequency. This indicates the development of the edge effect. Such an effect occurs when certain cortical structures are forced to generate y frequencies in a continuous stereotyped manner. This leads to the productions of cognitive and motor behaviors in the absence of external context or intention.

It is proposed that the edge effect is responsible for the positive symptoms that are reported by the test subjects.

* * * The present invention is not to be limited in scope by the specific embodiments described herein. Indeed, various modifications of the invention in addition to those described herein will become apparent to those skilled in the art from the foregoing description and the accompanying figures. Such modifications are intended to fall within the scope of the appended claims.

Patents, patent applications, and publications are cited throughout this application, the disclosures of which are incorporated herein by reference in their entireties.