Classification of Herbal Drug Effects by Discriminant Analysis of Quantitative Human EEG Data ()
1. Introduction
Widespread clinical use of behavior-modifying drugs having defined biochemical actions has provided observations and generated hypotheses that relate neurochemical findings to normal and abnormal behaviors in man [1] . Since many drugs act on the central nervous system by interference with neurochemical processes leading to changes of electric activity, measurement of the electric activity of the brain is a very promising approach to characterize the effectiveness of drugs acting on the central nervous system [2] . To characterize and classify pharmacological effects of herbal drugs on the brain remains a major challenge. Within clinical studies many different parameters can be measured like filling out questionnaires, psychometric testing or collecting neurophysiological parameters like ECG and EEG.
In day-night converted rats recording of field potentials from the depth of the brain revealed significant changes in the presence of drugs in comparison to placebo. Quantitative analysis of these data by frequency analysis using Fast Fourier Transformation (FFT) documented time and dose dependent changes of spectral power in four brain areas (frontal brain, hippocampus, striatum and reticular formation). Results from frequency analysis were documented according to 6 specifically defined ranges (frequency bands). In order to discriminate effects of 40 drugs with known clinical indications from each other, 24 parameters (spectral power from 4 brain areas × 6 frequency bands) were fed into linear discriminant analysis. Surprisingly, drugs with similar indication clustered together (could not be discriminated well from each other) and were clearly separated from drugs with other indications [3] . The question now arose, if also quantitative human EEG data could serve to separate drug effects from each other. Therefore data from several clinical studies, which had been performed at NeuroCode AG in the past, were re-analyzed. In this case 102 parameters (spectral power from 17 brain areas × 6 frequency bands) were fed into discriminant analysis. As will be shown below, also human quantitative EEG data allow separation of effects from each other and comparison of the effects of drugs acting on the central nervous system in an analog manner i.e. drugs with a similar indication group together whereas others are well separated from each other. Thus, discriminant analysis allows for documentation of their unique pharmacological action also in humans. This leads to a direct comparison and a new method for classification of herbal drugs in humans.
2. Material and Method
Data from 10 clinical studies were reanalyzed. Quantitative EEG data were obtained by Fast Fourier Transformation of analog data [4] . Data were calculated as “source density” corresponding to magnet-encephalographic analysis [5] [6] . The resulting frequency spectra were cut into six frequency bands: delta: 1.25 - 4.5 Hz; theta: 4.75 - 6.75 Hz; alpha1: 7.0 - 9.5 Hz; alpha2: 9.75 - 12.5 Hz; beta1: 12.75 - 18.5 Hz; beta2: 18.75 - 35.0 Hz. Spectral power obtained in a relaxed state with open eyes in the presence of drugs was compared to placebo at 17 electrode positions corresponding to 17 brain areas. Technology CATEEM® from MEWICON CATEEM-Tec GmbH, Schwarzenberg, Austria was used.
All experiments were performed under identical conditions as published earlier [7] . In short: EEG was recorded for 5 minutes from subjects in a relaxed state before intake of the particular drug. Up to three hours later (as indicated individually for each drug) this procedure was repeated and the result compared to the baseline recording in terms of spectral power changes in 6 frequency ranges at all brain regions. Brain maps were constructed by transforming spectral power into spectral colours followed by additive colour mixture according to “RGB” [8] . Data were statistically compared to the spectral changes as observed under the corresponding placebo condition (representing the circadian rhythm) using the non-parametric Wilcoxon-Mann-Whitney Test. Finally, all 102 parameters were fed into discriminant analysis and documented as 6 dimensional display (results from the first 3 discriminant functions were displayed in space, result of next 3 functions into colour using the RGB mode of additive colour mixture). Detailed information on the drugs is given in Table 1 and Table 2.
3. Results
3.1 Circadian Rhythm
Electric activity of the human brain changes with time and follows a circadian rhythm [9] . Thus, during the 3 hours between baseline recording and the second recording spectral power has changed. This is documented as average spectral power in 148 subjects collected from 10 studies in Figure 1. Mainly alpha spectral power increases in nearly all brain areas. But also beta power increases in several brain areas. Therefore, spectral power changes induced in the presence of drugs have to be compared to this pattern statistically.
3.2. Calmvalera® Tablets
Calmvalera® tablets are prescribed for nervous restlessness and better coping with stress. This homeopathic remedy has been tested within a clinical study in young subjects suffering from examination anxiety. Clinical efficacy was demonstrated [10] . Statistically significant spectral power changes with respect to alpha waves in the presence of 6 tablets Calmvalera® were mainly confined to the left hemisphere as indicated by uneven numbers of electrode positions (F3, C3, P3, T3, O1). Increases of alpha and beta1 waves were statistically significantly different from spectral changes due to the circadian rhythm. Further details are given in Figure 2.
Table 1. Drugs with ingredients and given dosages.
Table 2. Clinical indications of the drugs investigated in the studies.
Figure 1. Spectral power changes in % of the baseline values (Ordinate) according to the circadian rhythm at all electrode positions (100% is indicated by a line). Red: spectral delta power; orange: theta power; yellow: alpha 1 power; green: alpha 2 power; turquoise: beta 1 power and blue: beta 2 power (for frequency band definitions see methods).
Figure 2. Spectral power changes in % of the baseline values in the presence of Calmvalera®. Statistical significance in comparison the circadian changes is documented by stars: *p < 0.1, **p < 0.05, ***p < 0.01. “+” underneath the bars indicates the direction of spectral changes (increase of spectral power).
3.3. Lasea®
Lasea® is an herbal drug prescribed to treat general anxiety [11] [12] . Its active ingredient Silexan (special extract of lavender oil) has been positively tested for cerebral bioavailability [13] . Spectral power changes in the presence of 2 capsules Lasea® consisted mainly in significant increases of alpha 1 or alpha2 power in most of the tested brain regions. Strongest increases were observed in central and temporal region as documented in Figure 3. This increase is in line with clinical relaxation. Psychometric testing did not reveal any negative results.
3.4. MemoLoges®
MemoLoges® has been clinically tested in subjects suffering from mild cognitive impairment [14] [15] . A dominant feature after intake of 2 capsules memoLoges® consisted in significant increases of alpha 1 spectral power in most of the brain areas except for the occipital lobe and parts of the frontal lobe (electrode positions F3 and F4). For details see Figure 4.
3.5. Neurapas®
Neurapas® balance has been developed to fight stress, anxiety and depression. The effect of this drug has been successfully characterized within a clinical study [16] . Intake of 6 tablets Neurapas® balance only led to a few statistically significant changes. In the parietal lobe delta and theta spectral power decreased somewhat, whereas alpha 1 power increased in the temporal lobe. Details are given in Figure 5.
3.6. Neuravena®
The next herbal drug, which has been analyzed, was Neuravena® (2500 mg).
Figure 3. Spectral power changes in % of the baseline values in the presence of Lasea®. Statistical significance in comparison the circadian changes is documented by stars: *p < 0.1, **p < 0.05, ***p < 0.01. “+” underneath the bars indicates the direction of spectral changes (increase of spectral power).
Figure 4. Spectral power changes in % of the baseline values in the presence of memoLoges®. Statistical significance in comparison the circadian changes is documented by stars: *p < 0.1, **p < 0.05, ***p < 0.01. “+” and “−” underneath the bars indicate the direction of spectral changes (increase or decrease of spectral power, respectively).
Figure 5. Spectral power changes in % of the baseline values in the presence of Neurapas® balance. Statistical significance in comparison the circadian changes is documented by stars: *p < 0.1, **p < 0.05. “−” underneath the bars indicates the direction of spectral changes (decrease of spectral power).
This drug has been tested and found to have positive effects on cognition [17] [18] . In comparison to the circadian rhythm increases of alpha1 power dominate this pattern of spectral changes, which, however, only reach statistical significance at electrode positions P4 and T4. Decreases of delta power in many brain regions did not reach statistical significance. All spectral power changes are displayed in Figure 6.
3.7. Neurexan®
Neurexan® is sold as a homeopathic drug for the treatment of nervousness. A double-blind, randomized clinical trial proved efficacy [19] . A special pattern of spectral changes was observed after intake of 4 tablets Neurexan® since only statistically significant decreases of power dominated the pattern of changes. Most of the changes were confined to alpha waves (Figure 7).
3.8. Nutrifin Relax®
Another relaxing drug is Nutrifin Relax®. EEG was recorded after intake of 2 lozenges. Results have been published earlier [20] . The pattern of alpha power changes is not much different from circadian changes except for left frontal and temporal area (electrode position F7 and T3). However, beta 1 power increased statistically significantly at central, frontal and occipital brain regions. Details are given in Figure 8.
3.9. Pascoflair®
Pascoflair® is used to treat anxiety, nervousness and generalized anxiety disorder. In preoperative surgery patients it reduced anxiety [21] . Its clinical efficacy
Figure 6. Spectral power changes in % of the baseline values in the presence of 2500 mg Neuravena®. Statistical significance in comparison to circadian changes is documented by stars: *p < 0.1, **p < 0.05, ***p < 0.01. “+” underneath the bars indicates the direction of spectral changes (increase of spectral power).
Figure 7. Spectral power changes in % of the baseline values in the presence of Neurexan®. Statistical significance in comparison the circadian changes is documented by stars: *p < 0.1, **p < 0.05, ***p < 0.01. “−” underneath the bars indicates the direction of spectral changes (decrease of spectral power).
was also shown by a combination of psychometry and quantitative EEG [22] . Spectral power changes in the presence of 3 tablets Pascoflair® reveal mainly some statistically significant decreases of delta and theta power within the parietal and occipital lobe as depicted in Figure 9. The pattern of alpha power did not differ significantly from the changes as observed according to circadian rhythm.
3.10. Suntheanine®
An herbal calming drug is Suntheanine®. Its active constituent is L-Theanine, a compound found in green tea [23] . The drug is known for its calming properties without producing sedation. In the presence of 250 mg of this preparation only increases of spectral power are observed. Strongest statistically significant effects are seen with respect to alpha power in nearly all brain regions. However, also increases of beta power are observed in some brain areas, which are statistically significant at electrode positions P3 and O1. For details see Figure 10.
Figure 8. Spectral power changes in % of the baseline values in the presence of 2 lozenges of Nutrifin Relax®. Statistical significance in comparison to circadian changes is documented by stars: *p < 0.1, **p < 0.05. “+” underneath the bars indicates the direction of spectral changes (increase of spectral power).
Figure 9. Spectral power changes in % of the baseline values in the presence of Pascoflair®. Statistical significance in comparison the circadian changes is documented by stars: *p < 0.1, **p < 0.05, ***p < 0.01. “−” underneath the bars indicates the direction of spectral changes (decrease of spectral power).
Figure 10. Spectral power changes in % of the baseline values in the presence of 250 mg Suntheanine®. Statistical significance in comparison the circadian changes is documented by stars: *p < 0.1, **p < 0.05, ***p < 0.01. “+” underneath the bars indicates the direction of spectral changes (increase of spectral power).
3.11. Valium®
Valium® is one of the best-known drugs with a strong sedative, anxiolytic and anticonvulsive action [24] . This non-herbal drug (developed in 1963) shows an interesting pattern of spectral changes much different from the circadian rhythm. Theta and alpha frequencies are suppressed in a highly significant manner, whereas spectral beta power has increased in most of the brain areas as documented in Figure 11.
3.12. Zembrin®
Zembrin® (a special extract of Sceletium tortuosum, which has been consumed by original inhabitants in South Africa for a long time) is a food supplement marketed also in the USA and Canada for improvement of cognitive performance [25] [26] . In the presence of 50 mg Zembrin® increases of alpha 1 spectral power dominated the pattern of changes in a statistically significant manner. Only within the central and parietal area decreases of delta and theta power emerged. For details see Figure 12.
Statistically significant changes of single frequencies in the presence of drugs at any of the brain regions are listed in Table 3. Despite the common indication of calming nervousness (except for Zembrin® and memoLoges®) there is no common change of frequencies indicating a similar mechanism of action for all these calming drugs. However, there is some similarity between Lasea® and Suntheanine® with regard to the fact, that spectral power of all frequencies increase in a statistically conspicuous or significant manner. This can be interpreted in the sense that single drugs might have additional effects like being able to treat anxiety and/or depression.
Figure 11. Spectral power changes in % of the baseline values in the presence of 5 mg Valium®. Statistical significance in comparison the circadian changes is documented by stars: *p < 0.1, **p < 0.05, ***p < 0.01. “+” and “−” underneath the bars indicate the direction of spectral changes (increase or decrease of spectral power, respectively).
Figure 12. Spectral power changes in % of the baseline values in the presence of Zembrin®. Statistical significance in comparison the circadian changes is documented by stars: *p < 0.1, **p < 0.05, ***p < 0.01. “+” underneath the bars indicates the direction of spectral changes (increase of spectral power).
Table 3. Listing of statistically significant frequency changes at any of the brain regions.
“+” or “−” indicate the direction of changes. Brackets indicate p < 0.10, otherwise change of the frequency is statistically significantly different from circadian rhythm at least at the level of p < 0.05.
3.13. Brain Maps
Another way to document the results from quantitative EEG recording is to construct a brain map. This approach allows an overview on the involvement of different brain regions during the action of drugs. However, there is an influence of the circadian rhythm on electric brain activity during relaxation, which is recognized when placebo is given (Figure 1). These changes are also visible in the brain map in Figure 13. Depending on the action of drugs with respect to different frequencies, colours indicate dominant changes of spectral power. For example, increases of beta 2 waves—like in the presence of Valium®—lead to a prevalence of blue colour as depicted in Figure 14. Quantitative EEG maps in the presence of individual drugs differ strongly with respect to change of spectral frequencies as exemplified for all drugs in Figure 13 and Figure 14 for direct comparison.
3.14. Discriminant Analysis
The pattern of frequency changes in the brain in the presence of drugs seems to be very specific and can be differentiated from each other by looking at the spectral power in 17 regions. Nevertheless, it would be desirable to find an approach to compress the information for individual characterization and classification of the drug effects. This can be achieved by using discriminant analysis. This type of analysis is based on all 102 parameters (spectral power in 17 brain regions) and separates the effects of all drugs from each other best (Figure 15). This means, that, if drugs are projected in rather close vicinity in space, their action is probably similar. Within the present analysis Nutrifin Relax®, Pascoflair® and Suntheanine® appear in rather close vicinity and all three drugs are prescribed for the same indication: relaxation. However, all calming herbal drugs are positioned
Figure 13. Documentation of drug effects by quantitative EEG in terms of brain mapping. Maps are constructed by transforming spectral power into spectral colours followed by additive colour mixture according to “RGB”. (No false colour coding, but showing direct result of FFT!). Data are interpolated in a non-linear manner from 17 regions to give results for 99 channels in order to produce a map corresponding mathematically to a 64 channel EEG.
Figure 14. Documentation of drug effects by quantitative EEG in terms of brain mapping. Maps are constructed by transforming spectral power into spectral colours followed by additive colour mixture according to “RGB”. (No false colour coding, but showing direct result of FFT!).
Figure 15. Result of discriminant analysis. Please note that the image is 6-dimensional! Result of first three discriminant functions is displayed in space (X; Y and Z axis). Result from next three discriminant functions is displayed according to additive colour mixture (so-called RGB). Neurapas® is positioned far in the foreground, Suntheanine® is positioned far in the background.
far away from the synthetic strongly sedative drug Valium®. Two drugs acting positively on cognition (Zembrin® and memoLoges®) are not very far from each other, but show a different colour due to the result of the 4th to 6th discriminant function. Furthermore, the unique efficacy of Lasea® and Neurapas® is documented by their isolated position as shown in Figure 15.
4. Discussion
Quantitative EEG data have been used extensively to characterize drug effects on the brain in the past under the term Pharmaco-EEG [2] . However, direct comparison of several drugs with each other is only documented for single parameters like a particular frequency within a brain region. Recording of the EEG has been used successfully to compare different sedative drugs [27] . A pattern of frequency changes consisting of decreases of delta, theta and alpha activity in combination of strong increases of beta activity in the presence of Diazepam has been known for a long time as reported earlier [28] . Frequency changes observed in the freely moving rat correspond very much to those observed in humans as exemplified for Diazepam (Valium®) in Figure 16.
Individual patterns of frequency changes in the presence of different drugs also give hints on their mechanism of action since changes of particular frequencies
Figure 16. Comparison field potential analysis in the rat with quantitative EEG in humans for Diazepam (Valium®).
have been related to neurotransmitter action. Extensive animal studies have revealed that delta waves are influenced by the activity of acetylcholine. Theta activity seems to be controlled by norepinephrine, whereas alpha 1 activity reflects serotonergic actions. Changes in dopaminergic activity are represented by alpha 2 spectral power [29] . Beta 1 waves are under the control of glutamate [30] . Finally, beta 2 waves change with the activity of GABA [31] . There is some indication that Pascoflair® is able to modulate neuronal GABA receptor mediated transmission [32] . However, according to our data, beta 2 waves are not augmented in a statistically significant manner in the presence of Pascoflair®. For the action of Valium® involvement of the GABA-ergic system has been proven [33] [34] . Accordingly, spectral power of beta 2 waves has increased.
A methodological approach is now presented, which allows evaluation of all experimental parameters at the same time: discriminant analysis [35] . Projection of the results from the first three discriminant functions into space already allows discrimination of drug action from each other. Regarding the result of next three mathematical discriminant functions (displayed by colour) allows even more detailed characterization of drug effects. If discriminant analysis is not able to separate the effects of different drugs from each other (being projected in close neighborhood), a similar clinical indication can be assumed. If the effect of a drug is projected in an isolated place, a unique indication might be suspected. Among the tested preparations this seems to be correct for Lasea® and Neurapas®.
After successful application in preclinical research [3] [36] , the same approach has now been performed with human EEG data with regard to 11 herbal drugs in comparison to the synthetic classic Valium®. As proven in the rat for analysis of quantitative field potentials, drugs acting on the human brain can now also be classified according to their overall effects on the brain.
5. Conclusion
In conclusion, discriminant analysis is not only used successfully in the pharmacological characterization of drug effects in rats, but proves to be a very convenient tool in describing the unique effects of herbal drugs in humans. This is a major progress and should be used in the future for characterization and classification of pharmacological effects of herbal drugs on the brain.