![]()
|
|
||||||||
J. Biol. Chem., Vol. 279, Issue 24, 25598-25604, June 11, 2004
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
1 Subunit Reduces Membrane Expression and Impairs Gating by Agonists*








From the
Division of Pharmacology, Department of Neuroscience, School of Medicine, University of Naples Federico II, 80131 Naples, Italy, the Departments of
Child Neuropsychiatry and ¶Pediatrics, Second University of Naples, 80138 Naples, Italy, the ||Department of Anesthesia, Beckman Program for Molecular and Genetic Medicine, Stanford University, Stanford, California 943055117, and the **Departments of Anesthesiology and Pharmacology, Weill Medical College, Cornell University, New York, New York 10021
Received for publication, October 7, 2003 , and in revised form, March 5, 2004.
| ABSTRACT |
|---|
|
|
|---|
1 subunit of the glycine receptor (GlyRA1) causing hyperekplexia, an inherited neurological channelopathy. In HEK-293 cells, the agonist EC50s for glycine-activated Cl- currents were increased from 26 µM in wtGlyRA1, to 5747, 135, and 129 µM in R218Q, V260M, and Q266H GlyRA1 channels, respectively. Cl- currents elicited by
-alanine and taurine, which behave as agonists at wtGlyRA1, were decreased in V260M and Q266H mutant receptors and virtually abolished in GlyRA1 R218Q receptors. Gly-gated Cl- currents were similarly antagonized by low concentrations of strychnine in both wild-type (wt) and R218Q GlyRA1 channels, suggesting that the Arg-218 residue plays a crucial role in GlyRA1 channel gating, with only minor effects on the agonist/antagonist binding site, a hypothesis supported by our molecular model of the GlyRA1 subunit. The R218Q mutation, but not the V260M or the Q266H mutation, caused a marked decrease of receptor subunit expression both in total cell lysates and in isolated plasma membrane proteins. This decreased expression does not seem to explain the reduced agonist sensitivity of GlyRA1 R218Q channels since no difference in the apparent sensitivity to glycine or taurine was observed when wtGlyRA1 receptors were expressed at levels comparable with those of R218Q mutant receptors. In conclusion, multiple mechanisms may explain the dramatic decrease in GlyR function caused by the R218Q mutation, possibly providing the molecular basis for its association with a more severe clinical phenotype. | INTRODUCTION |
|---|
|
|
|---|
-aminobutyric acid (GABA)1 and Gly, which, by activating Cl--selective ligand-gated ion channels (LGICs), hyperpolarize the membrane of the postsynaptic neuron and inhibit firing (1). Glycine receptors (GlyRs) are heteropentameric proteins made up of the assembly of two types of homologous subunits:
subunits, which carry the molecular determinants for agonist and antagonist binding, and
subunits, which determine the subsynaptic localization of GlyRs (2). In humans, only one gene encoding for
subunits has been identified (3); this gene is widely expressed in the central nervous system at all developmental stages of vertebrates. On the other hand, four genes encoding distinct
subunits (
14) are known, each showing a specific regional distribution and ontogenetic profile (4). In particular,
1 subunits seem to be the predominant isoform expressed postnatally in the spinal cord and brainstem; GlyRs of mature vertebrates are thought to have a 3
12
stoichiometry. The topological organization of the GlyR subunits is believed to be similar to that of other members of the LGIC super-family (5), with both the N and the C termini located extracellularly and four transmembrane segments; the second transmembrane segment (M2) determines permeation properties such as anion selectivity and conductance.
GlyRs play a central role in the pathophysiology of hyperekplexia (HE), a rare inherited human neurological disease characterized by muscle stiffness, hyperreflexia, and an exaggerated startle response to sensory stimuli (6). The muscle stiffness and hyperreflexia, usually present at birth, decline during the first year of life, although an exaggerated startle response persists during the adult life. The gene associated with hyperekplexia was linked to chromosomal locus 5q3335 in 1992 (7); within this region, mutation analysis in HE-affected families proved that the GlyRA1 gene, encoding the
1 subunit of GlyR, was the defective gene in HE (8). Several HE-causing mutations are known in the GlyRA1 gene; these include both dominant and recessive inherited mutations, as well as de novo mutations found in sporadic cases. Further support for the hypothesis that abnormal GlyR function causes neurological disturbances derives from the observation that the phenotypes of spasmoid (9), spastic (10), and oscillator (11) mice, all of which carry genetically determined defects in GlyR-encoding genes, bear a striking resemblance to some of the sensorimotor manifestations of human HE.
In the present study, we have characterized the functional properties of three GlyRA1 mutations associated with HE. The R218Q mutation, recently described by members of our group (12), affects a residue that is located at the boundary between the extracellular N-terminal domain and M1 and is highly conserved among subunits forming both anion-selective (GlyR and GABAA,CR) and cation-selective (nicotinic
17 and
9,
14,
,
,
; 5HT-3A and 5HT-3B; some purinergic receptors) LGICs (Supplemental Fig. 1). The V260M (13) and the Q266H (14) mutations affect residues located in the M2 domain of GlyRA1 subunits but are not believed to participate in permeation (5, 15); they are less well conserved among LGIC subunits. We have used pharmacological, biochemical, and molecular modeling techniques to investigate the molecular basis of these three human "channelopathies" and found evidence that all three residues may cause a variable degree of impairment in GlyR gating, with the most dramatic phenotype associated with the R218Q mutation.
| EXPERIMENTAL PROCEDURES |
|---|
|
|
|---|
Heterologous Expression of GlyRA1 SubunitsWild-type (wt) and mutant GlyRA1 subunits were transiently expressed via transfection of human embryonic kidney (HEK-293) cells. HEK-293 cells were grown in minimal essential medium containing 10% fetal bovine serum, nonessential amino acids (0.1 mM), penicillin (50 units/ml), and streptomycin (50 µg/ml), in an humidified atmosphere at 37 °C with 5% CO2 in 100-mm plastic Petri dishes. For electrophysiological experiments, the cells were seeded on glass coverslips (Carolina Biological Supply Company, Burlington, NC). One day after plating, the cells were transfected with 2 µg of the appropriate cDNAs using LipofectAMINE 2000 (Invitrogen), according to the manufacturer's protocol, together with 2 µg of a plasmid encoding the enhanced green fluorescent protein (Clontech) used as a transfection marker. When lower amounts of GlyRA1 cDNA were used for transfection (see Fig. 5), the amount of enhanced green fluorescent protein cDNA was proportionally increased. All of the experiments were performed 12 days after transfection.
|
Cell Surface Biotinylation and Western BlottingPlasma membrane expression of GlyRA1 subunits in HEK-293 cells was investigated by surface biotinylation of membrane proteins in intact cells 2 days after transfection using a cell membrane-impermeable reagent (sulfosuccinimidyl-6-(biotinamido) hexanoate; Pierce) according to the manufacturer's protocol (18). Following biotinylation, the cells were lysed, and biotinylated proteins were isolated by reacting cell lysates with ImmunoPure streptavidin beads (Pierce). GlyR
subunits in streptavidin precipitates and total lysates were revealed in Western blots using an anti-GlyRA monoclonal antibody (mAb4a; Alexis Biochemicals, Lausen, Switzerland; dilution 1:4000); an alkaline phosphatase-conjugated sheep anti-mouse IgG (Amersham Biosciences; dilution 1:2000) was used as secondary antibody. Reactive bands were detected by chemiluminescence (SuperSignal, Pierce), using a ChemiDoc station (Bio-Rad). To confirm that intracellular proteins were not labeled by the biotinylation reagent, blots were stripped and reprobed with an anti-
-tubulin antibody (Sigma; dilution 1:5000).
Molecular ModelingA molecular model of a GlyRA1 subunit was built using a previously described model of the GABAAR
1 subunit as a template (1922). We built loops to fill in the gaps between the GlyRA1 sequence and the template sequence and refined the side chains in the model using the auto-rotomer feature in the Biopolymer Module of Insight II (v. 2000.1, Accelrys, San Diego, CA) as described previously (19, 22). The backbone atoms (C, C
, N) of each GlyR residue were tethered to the coordinates of corresponding residues in the template with a force constant of 100 kcal/A2, and the structure was optimized with the Discover_3 module of Insight II. Coulombic interactions were calculated with a constant dielectric of one and a cutoff at 9.5 A. The structure was relaxed by performing 5000 2-fs steps of molecular dynamics at 298 K and was then reoptimized with the Discover_3 module to a derivative of 1 kcal/A.
| RESULTS |
|---|
|
|
|---|
1 subunit, highlighting the Arg-218, the Val-260, and the Gln-266 residues, is shown in Fig. 1A. To assess the functional consequences of the HE-causing GlyRA1 mutations affecting these residues, we transfected HEK-293 cells with the cDNAs encoding wtGlyRA1 or R218Q, V260M, or Q266H GlyRA1 mutants and recorded Cl- currents carried by the expressed channels upon perfusion with the agonist Gly. In wtGlyRA1-expressing cells, perfusion with concentrations of Gly above 3 µM (0.0130 mM) activated robust inward Cl- currents when the cell was held at -60 mV (Fig. 1B); the reversal potential of these currents approximated the Cl- equilibrium potential (data not shown), confirming anionic selectivity. The peak amplitude of these Gly-gated Cl- currents increased with Gly concentration up to 300 µM; larger Gly concentrations (330 mM) failed to increase Clcurrent peak size, although they accelerated the time course of Gly-gated current desensitization. Peak amplitudes of Gly-gated currents were normalized to the maximum, expressed as a function of Gly concentrations, and fitted to a standard form of the Hill equation (y = ymax ([Gly]nH/[Gly]nH + EC50nH). The EC50 for Gly was 26 ± 4 µM, with a Hill coefficient (nH) of 2.0 ± 0.1 (Fig. 1C). These values concur with those reported previously (23, 24).
|
The two other HE-associated GlyRA1 mutations, V260M and Q266H, caused a more modest 5-fold increase in the EC50 for the agonist Gly (the EC50 was 135 ± 3 µM for V260M and 129 ± 8 µM for Q266H mutants) (Fig. 1C). Also in this case, no significant differences were detected in the Hill coefficients of V260M (nH = 1.7 ± 0.2) and Q266H (nH = 1.6 ± 0.2) GlyRA1 receptors when compared with wtGlyRA1 receptors.
To gain further insight into the molecular defect prompted by the R218Q mutation, we measured the rates of channel activation, desensitization, and deactivation by fitting with a single exponential function the early phases of agonist (30 mM Gly)-induced current activation, the time course of agonist (30 mM Gly)-induced current desensitization, and the current decay upon agonist (3 mM) removal in both wtGlyRA1 and R218Q channels. Comparison of these kinetic properties showed that wtGlyRA1 receptors displayed an increased activation rate and a faster desensitization rate during agonist exposure when compared with GlyRA1 R218Q channels (Fig. 1, D and E), without significant changes in the deactivation rate upon removal of the agonist (Fig. 1E).
Taurine and
-Alanine Sensitivity of HE-causing GlyRA1 Mutants
-amino acids such as
-alanine (
-ALA) and taurine (TAU) behave as agonists at wtGlyR, but they appear to be of lower efficacy than Gly in some mutant receptors, acting as partial agonists (25, 26). As shown in Fig. 2A, the peak Cl- currents elicited by exposure to 3 mM
-ALA or TAU in wtGlyRA1 receptors were similar in amplitude to those elicited by an identical Gly concentration (Fig. 2, B and C). By contrast, the currents elicited by both
-ALA and TAU in GlyRA1 carrying the V260M or the Q266H mutation were significantly smaller than those triggered by an identical concentration of Gly (3 mM) in the same cells (Fig. 2, AC). In cells expressing GlyRA1 subunits carrying the R218Q mutation, 3 mM TAU or
-ALA application failed to elicit detectable Cl- currents (Fig. 2A). In the large majority of the R218Q-transfected cells, perfusion with a higher
-ALA concentration (30 mM) also failed to elicit measurable Cl- currents. In the only three cells in which these currents were clearly above background, the I
-ALA/IGly ratio was 0.13 ± 0.03 (p < 0.05 versus controls) (Fig. 2B); in R218Q-transfected cells, 30 mM TAU failed to activate measurable Cl- currents (Fig. 2C).
|
|
In cells expressing R218Q GlyRA1 subunits, preincubation with 400 nM STR markedly reduced the response to 6 mM Gly (always in the presence of 400 nM STR) (Fig. 3C); the peak amplitude of the Cl- current elicited by 6 mM Gly in the presence of 400 nM STR was 8.0 ± 1.1% that of control (n = 4). This result was very similar to that obtained in cells expressing wtGlyRA1 subunits using the EC50 concentration of 30 µM (Fig. 3C), in which the peak amplitude of the Cl- current elicited by 30 µM Gly in the presence of 400 nM STR was 5.3 ± 1.7% (n = 5) of those recorded in the absence of STR. Measurements of the true equilibrium dissociation constant for STR by the method of Schild (27, 28) were precluded in GlyRA1 R218Q channels due to the very low potency of Gly, requiring the use of concentrations higher than 10 mM to elicit even small currents.
Membrane Localization of GlyRA1 Subunits Carrying HE MutationsTo determine whether the HE-causing mutations modified membrane insertion of GlyRA1 subunits, in addition to causing changes in the functional properties of correctly processed channels, Western blots were performed in total cell lysates and streptavidin-isolated plasma membrane proteins from HEK-293 cells transfected with the cDNAs encoding wtGlyRA1 or mutant GlyRs. The mAb4a monoclonal antibody was used, which recognizes an extracellular epitope (amino acids 96105) shared by all GlyRA subunits; GlyRA1 subunits were revealed by the appearance of an immunoreactive band of an approximate molecular mass of 48 kDa (29). This band was not detected in extracts from untransfected HEK-293 cells (data not shown).
The intensity of this 48-kDa band did not show significant differences between cells expressing wtGlyRA1 subunits or GlyRA1 subunits carrying the V260M or the Q266H mutation, both in total lysates (Fig. 4A) and in streptavidin-isolated plasma membrane proteins (Fig. 4B). By contrast, cells expressing GlyRA1 R218Q subunits showed a significant decrease in the intensity of this 48-kDa immunoreactive band, both in total lysates (Fig. 4A) and in isolated plasma membrane fractions (Fig. 4B) when compared with wtGlyRA1-expressing cells. In fact, densitometric analysis revealed that the intensity ratios between the 48-kDa band in total lysates and the band corresponding to
-tubulin (
55 kDa) in total lysates were 2.61 ± 0.94, 1.65 ± 0.45, 1.88 ± 0.64, and 0.46 ± 0.17 for cells transfected with wt, V260M, Q266H, and R218Q GlyRA1 cDNAs, respectively (n = 4; p < 0.05 between R218Q GlyRA1- and wtGlyRA1-transfected cells). Furthermore, the ratios between the intensity of the GlyRA1-specific band in streptavidin-isolated plasma membrane proteins and the band corresponding to
-tubulin in total lysates were 3.67 ± 1.66, 1.96 ± 0.39, 1.74 ± 0.33, and 0.25 ± 0.16 for cells transfected with wt, V260M, Q266H, and R218Q GlyRA1 cDNAs, respectively (n = 5; p < 0.05 between R218Q GlyRA1- and wtGlyRA1-transfected cells).
|
| DISCUSSION |
|---|
|
|
|---|
|
-amino acids (TAU and
-ALA) suggests that the coupling mechanism between agonist binding and channel gating is rendered less efficient in these GlyRA1 mutants, implying that these two M2 residues may participate in the molecular motions that occur during the gating process (24). However, these changes are relatively modest when compared with those produced by mutations in other critical gating residues, most notably R271Q, the first HE mutation to be identified (8, 26).
Our results with the R218Q substitution (12) were remarkable. The R218Q mutation caused a 200-fold decrease in the apparent affinity of Gly at its receptor. In these channels, maximal concentrations of
-ALA and TAU were virtually unable to activate Cl- currents. Once again, these drastic changes could reflect alterations in the agonist binding site, but the observation that Gly-gated currents were similarly antagonized by low concentrations of the competitive antagonist STR in both wt and R218Q GlyRA1 channels suggests that any changes at the agonist/antagonist binding site must be subtle. The efficacy of
-ALA and TAU was reduced close to zero in R218Q GlyRA1 channels, suggesting a very deleterious change in the coupling machinery by the R218Q mutation in a residue close to the M1 segment of the GlyRA1 subunit. However, it should be emphasized that, because of technical limitations (limited solubility and osmolarity changes), agonist concentrations up to 30 mM were used in these experiments; thus, we cannot exclude the possibility that higher (>100 mM)
-ALA and TAU concentrations might have elicited partial agonist responses in R218Q GlyRA1 channels.
Interestingly, at the highest Gly concentration (30 mM), both activation and desensitization rates of wtGlyRA1 channels were faster than those observed in GlyRA1 R218Q channels, with no change in deactivation rates. These data are consistent with the hypothesis of a deleterious effect of the R218Q mutation on channel gating since at these high agonist concentrations, the activation rate approximates the channel opening rate. Given the positive coupling between channel desensitization and opening demonstrated for nicotinic acetylcholine receptors (32), it seems plausible that an enhancement in the channel opening rate also increases the proportion of receptors entering the desensitized state from the open state. An interference with the gating process has been also proposed to explain the results recently obtained upon mutation of an arginine residue (Arg-222) located immediately past the residue corresponding to Arg-218 in GlyRA1 receptors in another member of the LGIC family, the 5HT3A receptor (33). The possible alteration of GlyRA1 gating by the R218Q mutation is also supported by our molecular model (Fig. 6A), in which Arg-218 is located a long distance from the agonist binding site (which is close to Phe-159) (34). In fact, the distance between the C
atoms of Phe-159 and Arg-218 is 29 Å in our model; however, Arg-218 lies relatively close to other residues that have been implicated in the gating process, notably Arg-271 (8, 26), Lys-276, and Asp-148 (24). The function of Arg-218 cannot easily be surmised from this model, but several negatively charged residues (Glu-53, Asp-57, Asp-148) are located close to Arg-218 and might interact with it during the gating process (Fig. 6B), as has been proposed for other residues within the corresponding region of the closely related GABAA receptor (20). It is also striking that this residue is strictly conserved among subunits forming both cation-selective (nicotinic, serotoninergic, some purinergic) and anion-selective (GABAA,C, Gly) receptors within the LGIC super-family, indicating that it is likely to play a vitally important structural or functional role.
In addition to those mutations that cause functional changes in channels that have been properly inserted into the plasma membrane, disease-causing human mutations may affect channel function by interfering with the complex mechanism of protein folding, maturation, assembly, membrane trafficking, and retrieval (35, 36). Our results show that the steady-state levels of R218Q GlyRA1 subunit protein were markedly decreased when compared with wtGlyRA1 subunits, whereas the levels of the V260M or Q266H subunit protein were normal, as reported previously for another HE-causing GlyRA1 mutant (P250T) located in the intracellular linker between the transmembrane regions M1 and M2 (29). Thus, in addition to profoundly affecting gating, the R218Q mutation reduces GlyR function by drastically decreasing the number of receptors incorporated into the plasma-membrane, possibly because of a decreased synthesis and/or increased degradation of the receptor protein (37, 38). Interestingly, wtGlyRs can be internalized and cleaved proteolytically after ubiquitination at the plasma membrane (39). It will be of interest to verify whether GlyRA1 subunits carrying the R218Q mutation undergo an enhanced protein turnover as a consequence of an increased plasma membrane ubiquitination prompted by the mutation.
It should be underlined that, although in Xenopus oocytes expression level affects GlyR agonist sensitivity (30, 31), the reduced sensitivity of R218Q mutant receptors to Gly seems not to be a consequence of the decreased expression of mutant subunits in HEK-293 cells. In fact, agonist sensitivity of wtGlyRA1s was unchanged when these receptors were expressed at levels comparable with those of GlyRA1s carrying the R218Q mutation. Interestingly, the observation that the maximal macroscopic currents recorded from cells expressing GlyRA1 R218Q mutant subunits are considerably smaller than those recorded from cells expressing a comparable amount of wtGlyRA1 subunits seems to suggest that the R218Q mutation might affect single channel conductance or maximal opening probability, although this possibility remains to be verified.
When compared with other GlyRA1 mutations found among HE-affected patients, the R218Q mutation represents the most dramatic end of the spectrum in terms of its effects on the receptor function and expression. Interestingly, the clinical evolution of the disease in the child carrying this mutation also lies at one extreme end of an ample spectrum of phenotypic severity; in fact, in addition to the classical HE symptoms, the affected child also showed a considerable degree of neurocognitive impairment, requiring special education in school at 6 years of age (12). This result is of considerable interest in view of the fact that neurocognitive development is usually reported to be normal in HE patients (6), and no perinatal environmental factor could be identified that may explain the neurocognitive impairment of our patient (12). In a similar case of a 6-year-old child affected by a severe recessive form of HE (40), a progressive decrease of the classical HE symptoms occurred with time, but mental retardation developed during childhood. The GlyRA1 mutation described in this case was a complete deletion of exons 16 of the GlyRA1 gene, which would result in a complete loss of protein. This is in close analogy with our results in the R218Q mutant, in which a 200-fold decrease in Gly potency, regardless of the molecular nature of the defect, can be considered as a functional GlyRA1 knockout. Thus, GlyRA1 subunits, in addition to their role in brainstem reflex modulation and spinal control of muscle tone, might also participate in neurocognitive development in humans. GlyRs have been characterized in septal cholinergic neurons, which play a major role in learning and memory by means of their connections with the hippocampus and the cerebral cortex (41), and their function can be modulated by several psychotropic drugs (4). However, compensatory mechanisms for the lack of GlyRA1 subunits seem to operate more efficiently in humans than in mice. The frameshift mutation found in the oscillator mouse (11) produces a complete loss of GlyR
l polypeptide in the mouse central nervous system (42), and homozygous oscillator mice all die around the 10th postnatal day.
In conclusion, our analysis of the molecular defects associated with three HE-causing mutations in GlyRA1 subunits confirm a modular hypothesis for LGIC subunits, in which distinct regions participate in agonist and antagonist binding, gating, and post-translational stability (4), and reveals a potential novel role for the Arg-218 residue in GlyR gating. This mutation has a dramatic effect on several aspects of GlyR function and is associated with a more severe clinical phenotype, possibly highlighting a previously under-recognized role of GlyRA1 subunits in higher cognitive functions in humans.
| FOOTNOTES |
|---|
The on-line version of this article (available at http://www.jbc.org) contains a supplementary figure showing the alignment of the predicted primary sequence of LGIC subunits. ![]()

To whom correspondence should be addressed: Div. of Pharmacology, Dept. of Neuroscience, University of Naples Federico II, Via Pansini 5, 80131 Naples, Italy. Tel.: 39-081-7463310; Fax: 39-081-7463323; E-mail: mtaglial{at}unina.it.
1 The abbreviations used are: GABA,
-aminobutyric acid; GlyR, glycine receptor; LGIC, ligand-gated ion channel; HE, hyperekplexia; TAU, taurine;
-ALA,
-alanine; STR, strychnine; wt, wild-type. ![]()
| ACKNOWLEDGMENTS |
|---|
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
D. K. Crawford, D. I. Perkins, J. R. Trudell, E. J. Bertaccini, D. L. Davies, and R. L. Alkana Roles for Loop 2 Residues of {alpha}1 Glycine Receptors in Agonist Activation J. Biol. Chem., October 10, 2008; 283(41): 27698 - 27706. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. Purohit and A. Auerbach Acetylcholine Receptor Gating at Extracellular Transmembrane Domain Interface: the "Pre-M1" Linker J. Gen. Physiol., November 26, 2007; 130(6): 559 - 568. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Wang, H. A. Lester, and D. A. Dougherty Establishing an Ion Pair Interaction in the Homomeric {rho}1 {gamma}-Aminobutyric Acid Type A Receptor That Contributes to the Gating Pathway J. Biol. Chem., September 7, 2007; 282(36): 26210 - 26216. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. L. Price, K. S. Millen, and S. C. R. Lummis Transducing Agonist Binding to Channel Gating Involves Different Interactions in 5-HT3 and GABAC Receptors J. Biol. Chem., August 31, 2007; 282(35): 25623 - 25630. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. Bellini, F. Miceli, S. Mangano, E. M. del Giudice, G. Coppola, A. Barbagallo, M. Taglialatela, and A. Pascotto HYPEREKPLEXIA CAUSED BY DOMINANT-NEGATIVE SUPPRESSION OF GLYRA1 FUNCTION Neurology, May 29, 2007; 68(22): 1947 - 1949. [Full Text] [PDF] |
||||
![]() |
|