![]()
|
|
||||||||
J. Biol. Chem., Vol. 279, Issue 18, 18614-18622, April 30, 2004
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||


¶

**
From the
Center for Neurologic Diseases, Brigham and Women's Hospital and ||Department of Endocrinology, Beth Israel Deaconess Medical Center, Harvard Medical School, Boston, Massachusetts 02115 and ¶Institute for Human Genetics, University Clinic Charité, D-13353 Berlin, Germany
Received for publication, February 2, 2004 , and in revised form, February 19, 2004.
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
Exposure to selective neurotoxins, including paraquat and 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP), has been linked to either increased risk of PD or chemically induced parkinsonism (8, 9). These compounds have been shown experimentally to decrease mitochondrial function and selectively inhibit the activity of complex I (10). In vitro chemical inhibition of complex I results in reduced oxidative phosphorylation and increased mitochondrial generation of ROS, similar to what was observed in mitochondria from PD patients (1113). Pathological examinations of PD brains have revealed increases in protein and lipid byproducts of ROS, including protein carbonyls and 4-hydroxynonenal (4HNE) (14, 15). Furthermore, 4HNE forms adducts with and inhibits the activities of the D1 dopamine (DA) receptor and the DA transporter (DAT), suggesting that lipid peroxides may contribute to the disruption of DA signaling (16, 17). Cultured dopaminergic neurons have been shown to exhibit enhanced sensitivity to paraquat and MPTP (1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine) as well as ROS (18). These findings suggest that mitochondrial dysfunction and accompanying ROS generation could be a common mechanism for the selective loss of substantia nigra neurons and the nigrostriatal DA signal in PD (19).
In addition to the more prevalent, idiopathic form, a subset of PD patients exhibits familial inheritance patterns. Large numbers and varieties of autosomal recessively inherited mutations in parkin are the predominant cause of familial PD (20). Initially described as juvenile-onset, atypical parkinsonism lacking Lewy bodies, subsequently identified cases are often clinically and pathologically indistinguishable from early onset familial or sporadic PD, including the presence of Lewy bodies in a single case (2123). We have recently reported that the loss of parkin function in mice results in nigrostriatal dysfunction, as evidenced by increased extracellular dopamine concentration in the striatum, reduced synaptic excitability in the striatal neurons, and behavioral deficits in paradigms that are sensitive to alterations in the nigrostriatal pathway (24). Despite measurable differences in nigrostriatal function in parkin-/- mice (24), no reduction in the number of dopaminergic neurons was observed in two independently generated parkin-/- mice (24, 25).
Parkin has been reported as an E3 ubiquitin-protein ligase (26). Previous reports described several substrates for parkin-mediated ubiquitinylation (27). It has been suggested that the loss of parkin function results in aberrant accumulation of substrate proteins including PAEL receptor, synphlin-1, and CDC-rel1. Accumulation of one or more of these proteins has been postulated to confer toxicity upon dopaminergic neurons in the substantia nigra (28). However, steady-state levels of these substrates are unchanged in parkin-/- brains (24).2 Recent evidence has also suggested a role for parkin in the protection of monoaminergic neurons against proteasomal dysfunction,
-synuclein overexpression-mediated cell death (29), and kainic acid-induced toxicity (30). It was shown that parkin is localized in mitochondria and inhibits mitochondria-dependent cell death (31). Other studies demonstrate that overexpression of mutant parkin elevates cellular markers of oxidative stress, whereas overexpression of wild-type parkin results in reduced levels of these markers (32). These observations are consistent with findings from parkin-null flies, which exhibit altered mitochondrial morphology and degeneration of tissues such as wing flight muscles and spermatocytes (33). These results raised the possibility that parkin may be involved in mitochondrial function. Based on these observations, we hypothesized that lack of parkin function may cause impairment of mitochondrial function in parkin-/- mice.
To determine whether a lack of parkin causes changes in protein abundance and/or modification, we conducted a nonbiased proteomic analysis of the ventral midbrain of parkin-/- and wild-type mice. Using a well established method for two-dimensional analysis of brain lysates (34), we were able to detect
8000 discrete protein spots from extracts of the ventral midbrain of parkin-/- and wild-type mice. Comparative analysis of 10 pairs of wild-type and parkin-/- brain samples revealed reproducible, quantitative changes of fifteen protein spots by silver staining. Subsequent mass spectrometric (MS) analysis revealed that these 15 spots represented 14 distinct proteins, 13 of which exhibited decreases in abundance in brains of parkin-/- mice and 1 additional protein which exhibited altered electrophoretic mobility, consistent with differential post-translational modification. Eight of these proteins were involved in either oxidative phosphorylation or antioxidant activities. Consistent with these findings, parkin-/- mice exhibited decreases in oxidative phosphorylation, weight gain, and antioxidant capacity as well as increased ROS-mediated tissue damage, suggesting an essential role for parkin in regulating normal respiratory function of mitochondria as well as in the protection of cells from oxidative stress.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Two-dimensional Gel Electrophoresis and Mass SpectrometryProtein samples for two-dimensional gel electrophoresis were prepared from the dissected ventral midbrain (including the substantia nigra) of each of the 10 pairs of parkin-/- and wild-type mice as previously described (35) with the following modifications. The solutions used for extractions were100 mM phosphate buffer, pH 7.1 (0.2 M KCl, 20% w/v glycerol, and 4% w/v 3-[(3-chloramidopropyl) dimethylammonio]-1-propanesulfonate) (A), protease inhibitor solution I (1 CompleteTM tablet (Roche Applied Science) dissolved in 2 ml of buffer A) (B), and protease inhibitor solution II (1.4 µM pepstatin A and 1 mM phenylmethylsulfonyl fluoride in ethanol) (C). The frozen tissue was transferred into a mortar placed in a liquid nitrogen bath. An aliquot of 1.25 parts v/w of A, 0.08 parts v/w of protease inhibitor I, and 0.02 parts v/w of protease inhibitor II were added to the tissue and ground to a fine powder. The resulting powder was filled into a 2-ml microtube, quickly thawed, supplied with 0.034 parts of glass beads, and then sonicated in an ice-cold water bath 6 times for 10 s with intervals of 1 min 50 s. The homogenate was stirred 30 min in the presence of 0.025 parts v/w of Benzonase (Merck). 6.5 M urea, 2 M thiourea, and 70 mM dithiothreitol solution were added, and stirring was continued for an additional 30 min. The protein extract was supplied with 0.1 parts v/w of ampholyte mixture Servalyte pH 24 (Serva, Heidelberg, Germany) and stored at -80 °C or analyzed immediately.
Proteins were separated by large two-dimensional gels as described previously (34, 35). Briefly, the gel format was 40 cm (isoelectric focusing, prepared with carrier ampholyte mixture covering pH 310) x 30 cm (SDS-PAGE, 15%) x 0.75 mm. The amount of the protein sample applied to the gel was 5 µl (60 µg/µl). For sample comparisons brain extracts from each pair of parkin-/- and control mice were run and stained in parallel. High sensitivity silver staining was used to visualize proteins (35). Two-dimensional gels were evaluated visually pairwise, and changes of spots were considered with respect to variation in the presence or absence, quantity, and spot position. Protein spots found to be reproducibly altered in parkin-/- patterns versus wild type were evaluated with the Proteomweaver imaging software Version 2.1 (Definiens, Munich, Germany). Although the mice we used are in a homogenous genetic background (129/Sv inbred strain), we still observed individual variations. Protein alterations confirmed in more than six pairs of mice were scored. All 10 parkin-/- mice investigated were affected at least in 7 of 14 proteins, and 5 mice were affected in more than 12 proteins. Data were analyzed by Student's t test.
For protein identification using MS, 18-µl (60 µg/µl) samples were electrophoresed on 1.5-mm gels and stained with MS-compatible silver stain or colloidal Coomassie Brilliant Blue G-250. Protein spots of interest were excised from gels and subjected to in-gel trypsin digestion without reduction or alkylation. Tryptic fragments were analyzed by a combination of matrix-assisted laser desorption ionization time-of-flight and liquid chromatography/electrospray ionization ion trap MS. The mass spectra were analyzed using Protein Prospector (MS-Fit) and Sequest Version 3.1 software.
Mitochondrial RespirationMice were euthanized by CO2 inhalation, and tissues were rapidly dissected on ice. Brains were removed, and striata were isolated as described previously (24). Striata from 2 mice of each genotype were pooled for mitochondrial isolation. Tissue samples were homogenized in 10 ml of buffer A (320 mM sucrose, 5 mM Tris, 2 mM EGTA, pH 7.4, at 4 °C) with 5 strokes of a Teflon Dounce. Samples were centrifuged for 3 min at 2000 x g to remove nuclei and tissue particles. Supernatants were collected and centrifuged for 10 min at 12,000 x g to pellet mitochondria and synaptosomes. The crude pellet was resuspended in 10 ml of buffer A with the addition of 0.02% w/v of digitonin to disrupt synaptosomal membranes and release trapped mitochondria. The resuspended pellet was centrifuged for 10 min at 12,000 x g to pellet mitochondria, which was resuspended in 100 µl of buffer A, and protein content was determined by BCA assay (Pierce). Mitochondria were resuspended at a final concentration of 0.4 mg/ml protein in 0.5 ml of buffer B (120 mM KCl, 3 mM HEPES, 1 mM EGTA, 5 mM KH2PO4, pH 7.2) with 1% w/v of bovine serum albumin and assayed for respiration using an excess of 8 mM glutamate, 8 mM malate (complex I), 4 mM succinate (complex II), or 0.4 mM N,N,N',N'-tetramethylphenylenediamine (TMPD)/1 mM ascorbate (complex III/IV) as electron donors. ADP was added in limiting amounts (14 µM), and state 3 respiration was measured. After depletion of ADP, state 4 respiration was measured. After determination of coupled respiration, 400 nM carbonyl cyanide 4-(trifluoromethoxy)phenylhydrazone (FCCP) was added to the reaction chamber, and respiration was measured in the absence of a proton gradient. Mitochondrial respiration was determined using a platinum electrode with a 1-ml buffer chamber (DM-10, Rank Bros., Ltd., UK). Due to intra-day variations in respiration rates, oxygen consumption values are represented as a fraction of the wild-type state 3 respiration for succinate. Data were analyzed by unpaired Student's t test.
Electron MicroscopyWild-type and parkin-/- mice were euthanized by CO2 inhalation and transcardially perfused with 20 ml of phosphate-buffered saline followed by 10 ml of fresh 2.5% glutaraldehyde plus 2.5% formaldehyde in 100 mM cacodylate buffer. Brains were removed and post-fixed in the above fixative for an additional 18 h at 4 °C. Brains were washed in phosphate-buffered saline and cut into 1-mm thick coronal sections, and striata were dissected and processed for standard electron microscopy by epon embedding with osmium tetroxide fixation and uranyl acetate counterstaining. Mitochondrial number and morphology were determined in images from 25 different fields from 2 mice per genotype by an investigator blind to the genotype.
Body Weight MeasurementsBody weights were measured in parkin-/- and wild-type mice at regular intervals beginning 10 days after weaning. Weights were analyzed by two-way analysis of variance followed by Bonferroni post-hoc analysis. Data for adult male mice were collected in conjunction with behavioral analysis and were analyzed by unpaired Student's t test.
Serum Antioxidant Capacity AssayMice were fasted overnight to minimize variability due to dietary uptake. Mice were euthanized by CO2 inhalation, blood was collected in heparinized tubes, and serum was isolated by centrifugation at 1200 rpm for 10 min. Serum total antioxidant capacity was measured by the conversion of Cu2+ to Cu1+ using a colorimetric assay and is represented as µM uric acid equivalents and carried out as per manufacturer's instructions (Total Antioxidant Potential, Oxis Research). Data were analyzed by unpaired Student's t test.
Protein Carbonyl AssayBrains were homogenized in 50 mM Tris, 150 mM NaCl, and 1% v/v Triton X-100, pH 7.5, and insoluble material was removed by centrifugation. Supernatants were assayed for protein content (BCA, Pierce) and 20 µg of protein was assayed for protein carbonyls as per the manufacturer's instructions (OxyBlot, Chemicon). Briefly, proteins were diluted into a final concentration of 6% SDS and reacted with 2,4-dinitrophenylhydrazine for 15 min. After the reaction, samples were neutralized, electrophoresed on NuPage gels (Invitrogen), transferred to nitrocellulose (Protran BA-83, Schleicher & Schuell), and Western-blotted using an antibody specific to the dinitrophenylhydrazone-derivatized residues on oxidatively damaged proteins. Blots were stripped and subsequently reprobed for actin (AC-15, Abcam, Cambridge, MA) to confirm equivalent protein loading.
ImmunohistochemistryBrains were removed, fixed in neutral-buffered formalin for 2 h at room temperature, paraffin-embedded, and sectioned as previously described (24). Tissue sections were blocked with 10% goat serum, incubated with an antibody to Michael adducts of 4HNE (#393207, 1:200, Calbiochem) and then with an avidin-conjugated secondary. Immunoreactivity was visualized with diaminobenzidine.
| RESULTS |
|---|
|
|
|---|
8 months of age. Proteins were resolved in the first dimension by their isoelectric point (pI) on a 40-cm tube gel using carrier ampholytes and subsequently resolved in the second dimension by their molecular weight on 40 x 30-cm SDS-PAGE gels. After silver staining, we detected reproducible, specific changes in 15 of
8000 discrete spots between the genotypes (Fig. 1). Contrary to our expectations, the staining intensity of all but one of these 15 protein spots was decreased in parkin-/- mice. Isolation of protein spots from these gels followed by trypsin digestion and subsequent matrix-assisted laser desorption ionization and electrospray ionization MS provided the identification of these proteins (Table I). Fourteen of the 15 spots represented distinct proteins, whereas one protein was detected twice in 2 adjacent spots (A5 and A6) with varying pI, suggesting a post-translational modification. Changes in protein spot intensity that were confirmed in more than six pairs of parkin-/- and wild-type mice were scored (Table I).
|
|
-crystallin). parkin-/- Mice Exhibit Reduced Mitochondrial RespirationPrevious reports suggest that PD may be linked to mitochondrial dysfunction (39). The biochemical changes detected in mitochondrial markers in the ventral midbrain of parkin-/- mice (Fig. 1) suggested possible mitochondrial deficiency in these mice. To address this issue directly, we examined the metabolic capacity of mitochondria isolated from the striatum of parkin-/- mice and wild-type controls by assaying for state 3 and state 4 respiration using substrates for complex I (glutamate/malate), complex II (succinate), and complex III/IV (TMPD/ascorbate). State 3 respiration measures the capacity of mitochondria to metabolize oxygen and the selected substrate in the presence of a limiting quantity of ADP, which is a substrate for complex V, the ATP synthase. In contrast, state 4 respiration occurs when all ADP is exhausted and measures respiration associated with proton leak across the inner mitochondrial membrane and represents a "basal-coupled" rate of respiration for each metabolic substrate.
Consistent with reductions in several units of complexes I and IV, we observed reduced rates of State 3 respiration (wild type versus parkin-/-: 75.6 ± 13.1 versus 62.8 ± 11.2, 147.7 ± 17.3 versus 110.3 ± 11.1, and 366.6 ± 72.2 versus 282.2 ± 59.7 nmol of O2/min/mg of protein for complex I, II, and III/IV, respectively; data are the mean ± S.E.), indicating a reduced capacity of mitochondrial metabolism. Representative traces for succinate-mediated respiration are shown in Fig. 2, A and B. Due to intra-day variations in absolute respiratory rate, data were analyzed as a fraction of the wild-type state 3 respiration for succinate (Fig. 2C). Consistent, statistically significant reductions were seen in mitochondria from parkin-/- mice using substrates that enter the electron transport chain at complex I (glutamate/malate), complex II (succinate), and complex III/IV (TMPD/ascorbate), suggesting an overall reduction in respiratory capacity. Measurement of state 4 respiration using either succinate or TMPD/ascorbate similarly exhibited a significant decrease in capacity in parkin-/- mice (Fig. 2D). The comparable magnitude of reduction in both state 3 and 4 respiration further supported the indication of reduced mitochondrial capacity for electron transport in parkin-/- mice.
|
A reduction in electron transport chain capacity is further confirmed by the significant reduction in respiratory rate for both succinate and TMPD/ascorbate after treatment with FCCP (Fig. 2F). FCCP treatment results in a collapse of the proton gradient across the inner mitochondrial membrane and allows mitochondrial respiration to proceed at the maximal capacity for the components of the electron transport chain without regard for either the capacity of complex V or proton leak. Thus, the decrease in mitochondrial respiration can be linked to a decrease in the capacity of the electron transport chain rather than a defect in ATP synthase capacity or function.
Mitochondrial Dysfunction Is Not Coupled with Alterations in Mitochondrial MorphologyThe decreased abundance of several mitochondrial proteins and the altered mitochondrial respiratory chain activity prompted us to examine the morphology of striatal mitochondria from parkin-/- mice and wild-type controls. Previous studies indicate that parkin-null flies exhibit morphological changes in mitochondria of the wing muscle and spermatocytes (33). Electron microscopy of striatal sections revealed no apparent change in the total number or size of mitochondria present in parkin-/- mice (data not shown). Examination of 100 µm2 of striatal area from two mice per genotype showed no gross alterations in mitochondrial morphology in parkin-/- mice (Fig. 3, A and B). Further examination at higher magnification (40,000x) showed that both parkin-/- and wild-type mitochondria had well formed cristae with no apparent swelling or separation of the inner and outer membranes (Fig. 3, C and D). These results suggest that reduced mitochondrial respiratory chain activity was not associated with detectable morphological changes in the mitochondria of parkin-/- mice.
|
20 days of age were monitored for weight gain starting at
30 days of age. Animals were weighed regularly for an additional 3 months. Data were pooled in either 10- or 20-day bins (±2 days). Over the course of the measurements, both male and female parkin-/- mice exhibited a significant decrease in the rate and total amount of weight gain (Fig. 4, A and B). Adult (612-month-old) male parkin-/- mice also exhibited significantly lower body weights (Fig. 4C).
|
|
| DISCUSSION |
|---|
|
|
|---|
-synuclein, in which missense and triplication mutations have been reported (4143). The linkage of
-synuclein mutations to PD and the accumulation of insoluble
-synuclein in Lewy bodies provided a strong link between protein aggregation and PD pathogenesis. Similarly, mutations in parkin and UCHL-1, two proteins linked to the ubiquitin-proteasome pathway, supported a possible involvement of aberrant accumulation of insoluble protein substrate(s) in PD pathogenesis. However, most parkin-linked PD cases lack Lewy bodies, suggesting that parkin mutations may cause PD by mechanisms distinct from protein aggregation and Lewy body formation (22, 23). The association of parkin expression with cytoprotection in cultured cells (29) and flies (44), and its amelioration of mitochondrial dysfunction-induced cell death (31), provided tantalizing clues to an alternative mode of action, namely a regulatory function for parkin in mitochondrial activity.
In the present study we compared the proteome of the ventral midbrain of parkin-/- and wild-type mice and found specific and consistent decreases in the steady-state abundance of 13 proteins and altered electrophoretic mobility of 1 additional protein. Eight of these proteins are linked to mitochondrial respiration or detoxification of byproducts of mitochondrial respiration. Four proteins that are specifically reduced, i.e. the E1
subunit of pyruvate dehydrogenase, the 24- and 30-kDa subunits of complex I and subunit Vb of complex IV, are all directly involved in mitochondrial oxidative phosphorylation activity. An additional one, lactolylglutathione lyase, is part of a two-step enzymatic cascade to metabolize methylglyoxal to lactate (45, 46). Methylglyoxal is the non-enzymatic breakdown product of glycerol 3-phosphate and is linked to formation of advanced glycation end products. Accumulations of advanced glycation end product-modified proteins are present in Lewy bodies of PD patients (47). Three additional proteins showing selective and reproducible reductions, PRDX 1, 2, and 6, are small, thiol-dependent peroxidases. Members of the PRDX family have been linked to protection from ROS stresses, and mouse models of loss of PRDX proteins result in enhanced sensitivity to oxidative stress (4850).
The decreased steady-state levels of proteins essential for mitochondrial respiratory chain activity are accompanied by measurable physiological changes in mitochondria from brains of parkin-/- mice. These alterations resulted in reduced electron transport chain capacity in mitochondria from parkin-/- mice and were analogous to functional alterations described in mitochondria isolated from PD patients (51). Furthermore, studies have shown that overexpression of
-synuclein induces decreases in mitochondrial respiration and increases sensitivity to mitochondrial stress (5254), providing further support for a common mitochondrial mechanism in the pathogenesis of PD. Coupled to these mitochondrial changes were decreased levels of three antioxidant proteins and one protein linked to detoxifying mitochondrial glycation byproducts along with an overall reduction in serum antioxidant capacity. Others have reported lower levels of reduced glutathione in another parkin-/- mouse model (25), analogous to changes seen in PD patients (55). In similar experiments, overexpression of mutant parkin results in increased markers of oxidative stress (32). These decreases in antioxidant capacity are likely to render parkin-/- mice more susceptible to ROS-mediated damage.
In addition to decreases in several proteins linked to mitochondria and ROS, levels of several other proteins with possible connections to PD pathology were decreased in parkin-/- mice. The small chaperone protein,
-crystallin, regulates the stability of and inhibits aggregation of proteins in many tissues, notably the lens of the eye (56). Vacuolar protein sorting-29 is a mammalian homolog of a yeast protein linked to protein sorting and trafficking between the endoplasmic reticulum/Golgi and the endosome (57) and, thus, may be involved in the synaptic transmission changes we have previously observed in the parkin-/- mice. Lasp1 is a molecular scaffold protein that has been associated with process extension (58), and may be involved in the regulation of synapse maintenance. Similarly, profilin II is a small, actin-binding protein that has been shown to negatively regulate neurite formation, suggesting a role in the formation and maintenance of synapses (59). Finally, heterogeneous nuclear ribonucleoprotein A1 is a component of an RNA binding complex that regulates RNA splicing (60) and may be functionally related to the p38 subunit of tRNA synthase, a substrate for parkin-mediated ubiquitinylation (61).
Consistent with increased susceptibility to ROS damage, we detected elevated levels of oxidized proteins and lipids in the brains of parkin-/- mice. Increases in 4HNE immunoreactive species are present in Lewy bodies of PD brains (40). Furthermore, in vitro data indicates that 4HNE forms adducts with, and inhibits the function of the D1 receptor and DAT (16, 17), providing a functional correlation between ROS damage and decreases in DA signaling in PD. This may suggest a connection between the mitochondrial dysfunction and our findings of increased extracellular DA in the striatum of parkin-/- mice. Elevated extracellular DA has been shown to create an oxidizing environment both in vitro and in vivo (62, 63). This environment results in ROS damage and mitochondrial dysfunction analogous to our findings (64). Oxidative damage to lipids may result in increased lipid peroxides, including 4HNE. These labile lipid peroxides may damage DAT, resulting in reduced DAT activity (17), leading to elevated extracellular DA, thereby exacerbating the oxidative stress. A recent paper shows that parkin is inactivated by ROS (65), suggesting that increased ROS mediated damage may result in further decreases of parkin activity in PD brains. Clinical reports indicate that heterozygous parkin mutations may represent a risk factor for late-onset PD in a subset of cohorts (66). This presents the hypothetical scenario of a vicious cycle where decreased levels of parkin result in mitochondrial dysfunction, which in turn lead to increased ROS formation and further inactivation of parkin and so forth. Our data provide compelling evidence for mitochondrial dysfunction and oxidative stress in the absence of nigral degeneration. This suggests that these events may be proximal in the cascade of the pathogenesis of both idiopathic and parkin-linked familial PD.
| FOOTNOTES |
|---|
These authors contributed equally to this work. ![]()
** To whom correspondence should be addressed: Center for Neurologic Diseases, Harvard Medical School, New Research Bldg. 636E, 77 Avenue Louis Pasteur, Boston, MA 02115. Tel.: 617-525-5561; Fax: 617-525-5522; E-mail: jshen{at}rics.bwh.harvard.edu.
1 The abbreviations used are: PD, Parkinson's disease; ROS, reactive oxygen species; 4HNE, 4-hydroxynonenal; DA, dopamine; DAT, dopamine transporter; MS, mass spectrometry; pI, isoelectric point; TMPD, N,N,N',N'-tetramethylphenylenediamine; FCCP, carbonyl cyanide 4-(trifluoromethoxy)phenylhydrazine; PRDX, peroxiredoxin. ![]()
2 J. J. Palacino and J. Shen, unpublished results. ![]()
| ACKNOWLEDGMENTS |
|---|
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |