|
|
||||||||
SPECIAL SECTION ON SYSTEMS BIOLOGY OF THE MITOCHONDRION
1Biotechnology and Bioengineering Center and Department of Physiology, Medical College of Wisconsin, Milwaukee, Wisconsin; and 2Biomedical NMR Laboratory, Department of Biomedical Engineering, Eindhoven University of Technology, Eindhoven, The Netherlands
Submitted 3 May 2006 ; accepted in final form 7 July 2006
| ABSTRACT |
|---|
|
|
|---|
computational model; mitochondria; cellular energetics; oxidative phosphorylation; 31P-NMR spectroscopy
Here the question of whether this same mechanism can explain the observed data on the control of oxidative metabolism in skeletal muscle is investigated. Previous analyses of 31P-NMR spectroscopy (31P-MRS) data on energy balance in exercising skeletal muscle have mainly focused on testing ADP feedback control of mitochondrial ATP synthesis using black box descriptions of the mitochondrial ATP synthetic pathway (8, 1416, 28), Pi acceptor control (7), and thermodynamic control involving quasi-linear relations between cytoplasmic Gibbs free energy of ATP hydrolysis and mitochondrial ATP synthesis flux (13, 18, 31). Yet, to date, these 31P-MRS data have not been adequately explained based on a detailed mechanistic model of oxidative phosphorylation and cellular energetics.
To analyze and interpret data from skeletal muscle, our previously published model of oxidative ATP synthesis and metabolism in cardiomyocytes (5) is adapted to skeletal muscle by setting intracellular concentration pools of creatine and phosphate to appropriate values (based on measured data) and appropriately adjusting the cellular mitochondrial content to match the available morphometric data. Data on cytoplasmic ADP and Pi concentrations as a function of work rate in human forearm flexor muscle from 20 untrained healthy subjects (13), 6 subjects with complex I deficiency in skeletal muscle (11, 22), and a single patient with deficiency in adenine nucleotide translocase (ANT) (2, 3) were analyzed based on the computational model.
Results of this analysis indicate that the existing computational model of the kinetics of mitochondrial oxidative phosphorylation accurately captures the in vivo kinetics of oxidative ATP synthesis and transport of phosphate metabolites between mitochondria and cytoplasm in skeletal muscle. The mechanism of control of oxidative phosphorylation is through feedback of substrates for ATP synthesis. No additional regulatory mechanisms, such as feed-forward control of certain enzymes via cytosolic calcium levels (9) or functional coupling between mitochondrial creatine kinase and ANT (23, 26, 27), are necessary to explain the majority of the observed data. In addition, the impact of specific protein deficiencies on the relationship between oxidative phosphorylation and cytoplasmic ADP and Pi is successfully explained by making the appropriate modifications to the mitochondrial enzymes altered in the diseases.
| METHODS |
|---|
|
|
|---|
ATP utilization, cytoplasmic phosphoenergetic buffers, and oxidative ATP synthesis are simulated in a model of skeletal muscle energetics illustrated in Fig. 1. The cell is divided into cytoplasmic and mitochondrial compartments; the variables simulated within the compartments are listed in Table 1, with a brief description of the variables and units associated with each variable. The computational model for cellular energetics and oxidative phosphorylation is derived from recently published computational models developed for cardiac mitochondria (4) and cardiomyocytes (5). A complete description of the computational model is provided in the APPENDIX.
|
|
Model parameter descriptions and assigned values are listed in Table 2. With the exception of the total pool of exchangeable phosphate (TPP) in the cell, all parameter values in the model are fixed at values justified by previous studies. The total exchangeable phosphate pool is computed from the equation
![]() | (1) |
|
|
31P-NMR spectroscopic measurements of phosphate metabolites in the ulnar finger flexor muscle in the human forearm were acquired at 1.5 T at rest and during voluntary ramp exercise using 1H imaging-guided localization in all subjects and analyzed according to the methods described in detail elsewhere (12).
| RESULTS |
|---|
|
|
|---|
Model predictions of steady-state concentrations of cytoplasmic ADP and Pi as functions of the work rate (rate of cellular ATP consumption) are compared with experimental measures obtained from healthy subjects in Fig. 2. The model predicts that ADP concentration increases from
17 µM at rest to
110 µM at the maximum observed work rate and that Pi concentration increases from 0.3 mM at rest to
18 mM at the maximum observed work rate.
|
The agreement between the model simulations and the observed data is striking considering that a single adjustable parameter (TPP) was varied to match model simulations to the data. In fact, the model-predicted [ADP]c values are not sensitive to the value of TPP. A 10% increase in the value of this parameter results no significant change in the model predicted [ADP]c and a >100% increase in the mean-squared difference between model predictions of [Pi]c and the observed data. Thus the nature of the relationship between work rate and [ADP]c does not depend on the value of TPP. The predicted [ADP]c and [Pi]c as a functions of workload at the value of TPP = 40.5 mM (10% greater than the optimal value) is plotted as a dashed lines in Fig. 2. It is apparent that the higher value of TPP results in an improvement in the model fit to the [Pi]c data at low work rates, but an overall agreement between the model predictions and experimental data that is worse than for the optimal value.
To investigate the factors controlling this relationship, we analyzed the sensitivity of the model predictions of Fig. 2A to the parameter values used in the kinetic model for the ANT flux of Eq. A10. This expression invokes three parameters, XANT,
, and Km,ADP (assumed values for these parameters are listed in Table 2). To quantify the impact of variation in these parameters on the work-ADP relationship, we fit the predicted data to the function V = Vmax([ADP]c xo)/([ADP]c xo + Km), where V represents the ATPase flux and Vmax, xo, and Km are fitting parameters. The predictions plotted in Fig. 2A are well represented by this function for Vmax = 0.44 mmols1 (l cell)1, xo = 16.5 µM, and Km = 0.11 mM. Sensitivity analysis reveals that the apparent Km is highly sensitive to the value of
(sensitivity coefficient |
ln Km/
ln
|
5), and less sensitive to XANT (sensitivity coefficient |
ln Km /
ln XANT |
0.5). The Vmax is most sensitive to the value of XANT (sensitivity coefficient |
ln Vmax /
ln XANT|
0.6) and less sensitive to
(sensitivity coefficient |
ln Vmax/
ln
|
0.3). Neither Km nor Vmax is sensitive to the value of Km,ADP. Thus, in terms of the ANT transporter model, the apparent Km for the relationship between ADP and work is primarily controlled by the parameter
.
Analysis of Data from Patients with Complex I Deficiency
Model predictions and experimental data from six patients (three of which were first- and second-degree relatives) who were diagnosed with mitochondrial complex I deficiency at the isolated mitochondria level are plotted in Fig. 3. These data were collected under the same protocol as for the healthy subjects; data were first published in Ref. 11. The experimental data show steeper relationships between [ADP]c (or [Pi]c) and ATP hydrolysis rate than are observed in healthy subjects, resulting from an impaired capacity of the mitochondria to synthesize ATP and transport it to the cytoplasm as levels of ADP and Pi increase.
|
1/1,000 of the normal value to match the observed data. Thus, to obtain the model predictions illustrated in Fig. 3, the activity of complex I was reduced by approximate 1,000-fold compared with the normal case. Here the model shows excellent agreement to the observed data on both [ADP]c and [Pi]c. As is the case for the data from healthy subjects, the predictions of [Pi]c are sensitive to the value of TPP, whereas the predictions of [ADP]c are not. Thus the relationship between work rate and substrates for ATP synthesis is explained by a drastic reduction in the activity and a loss in proton pumping of mitochondrial complex I.
Analysis of Data from Patient with ANT Deficiency
Figure 4 shows data on [ADP]c and [Pi]c measured in a single patient characterized as having an ANT deficiency in muscle (2, 3). The ANT transporter exchanges mitochondrial ATP for cytoplasmic ADP. Thus impairment in the activity of ANT results in a reduction in the ability of the mitochondrion to deliver ATP to the cytoplasm. Data were collected under the protocol described for healthy and complex I-deficient subjects; a subset of these data was first published in Ref. 11.
|
Figure 4 shows that ADP was much higher at rest and increased more rapidly with exercise in this patient than in healthy subjects and in the complex I-deficient patients. The ADP concentration of 150 µM that was measured at the modest work rate of 0.081 mmol ATP consumed per second per liter cell was greater than any value measured in the healthy or complex I deficient subjects. Model predictions are similar to the measured data; the model predicts [ADP]c is 42 µM at rest (compared with the measured value of 64 µM) and increases sharply with work rate.
Maintenance of free energy of ATP hydrolysis.
Figure 5 shows model-predicted free energy of ATP hydrolysis vs. ATP hydrolysis rate, computed using the parameters obtained for the healthy subjects, complex I-deficient patients, and the ANT deficient patient. Since cellular pH values are 7.0 ± 0.1 for all the subjects, the free energy of ATP hydrolysis is computed by the relationship
![]() | (2) |
GATP is the standard free energy of ATP hydrolysis at pH 7.0, R is the universal gas constant, and T is temperature in degrees Kelvin (for parameter values, see Table 2). [fADP]c and [fATP]c denote magnesium unbound ADP and ATP concentration in the cytoplasm, respectively. For the normal subjects, Jeneson et al. (9) observed quasi-linear relationship between the free energy of ATP hydrolysis and power output; the model predictions verify this observation except at low work rates where the predicted values of [Pi]c tend to be lower than the observed values.
|
| DISCUSSION |
|---|
|
|
|---|
The analysis predicts that the rate of oxidative phosphorylation is primarily regulated through concentrations of the substrates for ATP synthesis (ADP and Pi), since no additional control mechanisms, such as feed-forward control of certain enzymes via cytosolic calcium levels (9) and functional coupling between mitochondrial creatine kinase and ANT (23, 26, 27) that have been proposed to operate in the heart, were incorporated into the model. The current analysis does not rule out the possibility that ancillary control mechanisms are active in skeletal muscle (16, 28); however, it shows that major contributions of such mechanisms to the overall regulation of the mitochondrial ATP synthetic pathway are not necessary to explain the thrust of the observed data.
Although the model predicts cytoplasmic ADP, Pi, and ATP (not shown) concentrations that agree well with observed data, the present model systematically underpredicts Pi concentration in the resting state. As illustrated in Fig. 2B, 31P MRS measurements indicate that inorganic phosphate concentrations are
3 mM at rest, whereas the model predicts resting [Pi]c to be only 0.3 mM. This may be due to in part to the fact that the mitochondrial model was constructed to match data obtained from mitochondria isolated from cardiomyocytes. Inorganic phosphate concentrations in heart are significantly lower than in mixed fiber type skeletal muscle such as human skeletal muscle (19) if not undetectable at low work rates (32).
We used the empirical fitting function V = Vmax([ADP]c xo)/([ADP]c xo + Km) to capture key trends in the model prediction for the control case and to characterize the sensitivity of the model predictions to key model parameters in terms of sensitivity coefficients for the fitting parameters. Sensitivity analysis revealed that the parameter
for ANT flux in the model has a significant impact on the apparent affinity for ADP in the relationship between cytoplasmic ADP concentration and rate of oxidative metabolism. Future studies including a full-scale sensitivity and metabolic control analysis of the current model as well as next-generation models incorporating ancillary biochemical detail will be necessary to further improve agreement of predicted data with empirical knowledge.
The control of oxidative phosphorylation by substrate concentrations allows the mitochondria to maintain a free energy of ATP hydrolysis of less than 55 kJ/mol over the observed range of work rates in human forearm muscle of healthy subjects. However, as shown in Fig. 5, the magnitude of
GATP drops more quickly with increasing work in the ANT-deficient patients than in normal subjects, and the predicted magnitude of
GATP in the complex I-deficient patients at rest is significantly lower than that of the other two sets of subjects. These abnormalities in the complex I- and ANT-deficient subjects result in reduced capacity to do work.
Data on complex I-deficient patients are fit by reducing the activity of complex I compared with the normal case, and assuming that complex I pumps no protons in the complex I deficient patients. However, while the closest fit to the observed data on ADP and Pi is obtained by reducing the complex I activity by a factor of
1/1,000 compared with the normal case, it is necessary to be cautious in interpreting the scaling factor of
1/1,000. Since the complex I activity for the normal case was not identified with significant sensitivity (4), the ratio between activities in normal and complex I deficient patients also cannot be determined sensitively. Yet, regardless of the sensitivity of the estimate, it was clear from oxygen polarographic studies on isolated mitochondria from patient leg muscle biopsies (H. R. Scholte, unpublished observations) that the activity of complex I is significantly diminished in the complex I deficient patients.
To analyze data from the patient with a deficiency of ANT, it was not necessary to introduce an arbitrary scaling factor to fit the measured data. Since the level of ANT expressed in mitochondria of the patient was directly assayed and found to be 25% of that in healthy subjects, it was possible to incorporate this measurement directly in the model by scaling the healthy ANT activity by 0.25. The fact that differences observed in cytoplasmic phosphoenergetic compounds between healthy subjects and this patient are explained based on this single adjustment to the model provides independent validation of the mitochondrial model and the conclusions drawn from its behavior.
| APPENDIX |
|---|
|
|
|---|
The computational model for skeletal muscle energetics and oxidative phosphorylation is derived from a previously published model applied to cardiac tissue (5). For the current application, the cardiac model has been modified in two ways to adapt the model to analyze data from human skeletal muscle. First, the oxygen transport component of the previous model has been removed. It is assumed that the skeletal muscle remains normoxic during the experiments and the cellular oxygen concentration, [O2], which is a variable in the model of Beard (5), appears as a fixed parameter in the current model. Second, the mitochondrial volume of the muscle cell, which is
30% of cell volume in cardiomyocytes (29), is set to Vmito = 0.056 (l mito) (l cell) 1, a value measured from biopsies of human vastus lateralis muscle (30).
The model is expressed in terms of the following set of differential equations:
![]() | (A1) |
In the above set of equations, the subscripts "x", "i", and "c", denote mitochondrial matrix, intermembrane space, and cytoplasm, respectively. All of the variables in this set of equations are defined in Table 1.
In addition to the state variables treated in Eq. (A1), the concentrations of several species are computed
![]() | (A2) |
Parameters that appear in the above equations are described in detail below. The fluxes that appear on the right-hand side of the governing equations are tabulated in Table 2. For mitochondrial species, the governing equations follow from Ref. 4. For cytoplasmic species, the reactions modeled are ATP consumption, creatine kinase reaction, adenylate kinase reaction, and transport between the cytoplasm and the mitochondrial intermembrane space.
Mathematical Expressions for Mitochondrial Fluxes
The expressions for the mitochondrial fluxes in the model are described in detail in Ref. 4, and are listed here without detailed explanations. Definitions of the variables and parameters that appear in the following expressions are listed in Tables 2 and 3.
![]() | (A3) |
![]() | (A4) |
GH = F
+ RT ln([H+]c/[H+]x).
![]() | (A5) |
![]() | (A6) |
![]() | (A7) |
![]() | (A8) |
![]() | (A9) |
![]() | (A10) |
is an empirical parameter with value set to 0.35.
The phosphate-hydrogen cotransporter flux
![]() | (A11) |
![]() | (A12) |
Mitochondrial adenylate kinase flux
![]() | (A13) |
![]() | (A14) |
![]() | (A15) |
Mathematical expressions for cytoplasmic reaction fluxes
Four biochemical processes are modeled in the cytoplasm-the adenylate kinase reaction, the creatine kinase reaction, ATP hydrolysis, and binding of magnesium ions to ADP and ATP.
The binding of magnesium to ATP and ADP in the cytoplasm takes the same form as the binding fluxes in the mitochondria
![]() | (A16) |
![]() | (A17) |
In Eq. A18, KAK is the equilibrium constant for the reaction 2ADP
ATP + AMP, and XAK is the enzyme activity, which is set to a large enough value so that the reaction is effectively maintained in equilibrium.
The creatine kinase flux is modeled using the expression
![]() | (A18) |
The flux JAtC is defined as the flux through the reaction ATP
ADP + Pi. Mathematical models for the ATP consumption flux are considered in the Results section.
Modification of the Model for Patients with Complex I Deficiency
Equation A1 assumes that 4 protons are pumped from the matrix to the intermembrane space for each pair of electrons transferred via the reaction at complex I. For patients with complex I deficiency, it is assumed that no protons are pumped at complex I, and the equations for [H+]x and 
are modified as follows
![]() | (A19) |
| GRANTS |
|---|
|
|
|---|
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
| REFERENCES |
|---|
|
|
|---|
2. Bakker HD, Scholte HR, Jeneson JA. Vitamin E in a mitochondrial myopathy with proliferating mitochondria. Lancet 342: 175176, 1993.[ISI][Medline]
3. Bakker HD, Scholte HR, Van den Bogert C, Ruitenbeek W, Jeneson JA, Wanders RJ, Abeling NG, Dorland B, Sengers RC, Van Gennip AH. Deficiency of the adenine nucleotide translocator in muscle of a patient with myopathy and lactic acidosis: a new mitochondrial defect. Pediatr Res 33: 412417, 1993.[ISI][Medline]
4. Beard DA. A biophysical model of the mitochondrial respiratory system and oxidative phosphorylation. PLoS Comput Biol 1: e36, 2005.[CrossRef][Medline]
5. Beard DA. Integrated computational modeling of oxygen transport and cellular energetics explains observations on in vivo cardiac oxygen consumption and energy metabolites. PLoS Comput Biol. In press.
6. Blei ML, Conley KE, Kushmerick MJ. Separate measures of ATP utilization and recovery in human skeletal muscle. J Physiol 465: 203222, 1993.
7. Chance B, Eleff S, Bank W, Leigh JS Jr, Warnell R. 31P NMR studies of control of mitochondrial function in phosphofructokinase-deficient human skeletal muscle. Proc Natl Acad Sci USA 79: 77147718, 1982.
8. Chance B, Leigh JS Jr, Clark BJ, Maris J, Kent J, Nioka S, Smith D. Control of oxidative metabolism and oxygen delivery in human skeletal muscle: a steady-state analysis of the work/energy cost transfer function. Proc Natl Acad Sci USA 82: 83848388, 1985.
9. Cortassa S, Aon MA, Marban E, Winslow RL, O'Rourke B. An integrated model of cardiac mitochondrial energy metabolism and calcium dynamics. Biophys J 84: 27342755, 2003.
10. Gentet LJ, Stuart GJ, Clements JD. Direct measurement of specific membrane capacitance in neurons. Biophys J 79: 314320, 2000.
11. Jeneson JA. In vivo 31P NMR Studies of Cellular Bioenergetics in Healthy and Diseased Human Skeletal Muscle (PhD thesis). Utrecht, The Netherlands: Utrecht University, 1992.
12. Jeneson JA, van Dobbenburgh JO, van Echteld CJ, Lekkerkerk C, Janssen WJ, Dorland L, Berger R, Brown TR. Experimental design of 31P MRS assessment of human forearm muscle function: restrictions imposed by functional anatomy. Magn Reson Med 30: 634640, 1993.[ISI][Medline]
13. Jeneson JA, Westerhoff HV, Brown TR, Van Echteld CJ, Berger R. Quasi-linear relationship between Gibbs free energy of ATP hydrolysis and power output in human forearm muscle. Am J Physiol Cell Physiol 268: C1474C1484, 1995.
14. Jeneson JA, Wiseman RW, Westerhoff HV, Kushmerick MJ. The signal transduction function for oxidative phosphorylation is at least second order in ADP. J Biol Chem 271: 2799527998, 1996.
15. Kushmerick MJ. Energy balance in muscle activity: simulations of ATPase coupled to oxidative phosphorylation and to creatine kinase. Comp Biochem Physiol B Biochem Mol Biol 120: 109123, 1998.[CrossRef][Medline]
16. Kushmerick MJ, Meyer RA, Brown TR. Regulation of oxygen consumption in fast- and slow-twitch muscle. Am J Physiol Cell Physiol 263: C598C606, 1992.
17. Lee AC, Zizi M, Colombini M.
-NADH decreases the permeability of the mitochondrial outer membrane to ADP by a factor of 6. J Biol Chem 269: 3097430980, 1994.
18. Meyer RA. A linear model of muscle respiration explains monoexponential phosphocreatine changes. Am J Physiol Cell Physiol 254: C548C553, 1988.
19. Mizuno M, Horn A, Secher NH, Quistorff B. Exercise-induced 31P-NMR metabolic response of human wrist flexor muscles during partial neuromuscular blockade. Am J Physiol Regul Integr Comp Physiol 267: R408R414, 1994.
20. Munoz DR, de Almeida M, Lopes EA, Iwamura ES. Potential definition of the time of death from autolytic myocardial cells: a morphometric study. Forensic Sci Int 104: 8189, 1999.[CrossRef][ISI][Medline]
21. Pybus J, Tregear RT. The relationship of adenosine triphosphatase activity to tension and power output of insect flight muscle. J Physiol 247: 7189, 1975.
22. Roef MJ, Reijngoud DJ, Jeneson JA, Berger R, de Meer K. Resting oxygen consumption and in vivo ADP are increased in myopathy due to complex I deficiency. Neurology 58: 10881093, 2002.
23. Saks VA, Kongas O, Vendelin M, Kay L. Role of the creatine/phosphocreatine system in the regulation of mitochondrial respiration. Acta Physiol Scand 168: 635641, 2000.[CrossRef][ISI][Medline]
24. Tomashek JJ, Brusilow WS. Stoichiometry of energy coupling by proton-translocating ATPases: a history of variability. J Bioenerg Biomembr 32: 493500, 2000.[CrossRef][ISI][Medline]
25. Veech RL, Lawson JW, Cornell NW, Krebs HA. Cytosolic phosphorylation potential. J Biol Chem 254: 65386547, 1979.
26. Vendelin M, Kongas O, Saks V. Regulation of mitochondrial respiration in heart cells analyzed by reaction-diffusion model of energy transfer. Am J Physiol Cell Physiol 278: C747C764, 2000.
27. Vendelin M, Lemba M, Saks VA. Analysis of functional coupling: mitochondrial creatine kinase and adenine nucleotide translocase. Biophys J 87: 696713, 2004.
28. Vicini P, Kushmerick MJ. Cellular energetics analysis by a mathematical model of energy balance: estimation of parameters in human skeletal muscle. Am J Physiol Cell Physiol 279: C213C224, 2000.
29. Vinnakota KC, Bassingthwaighte JB. Myocardial density and composition: a basis for calculating intracellular metabolite concentrations. Am J Physiol Heart Circ Physiol 286: H1742H1749, 2004.
30. Vogt M, Puntschart A, Geiser J, Zuleger C, Billeter R, Hoppeler H. Molecular adaptations in human skeletal muscle to endurance training under simulated hypoxic conditions. J Appl Physiol 91: 173182, 2001.
31. Westerhoff HV, van Echteld CJ, Jeneson JA. On the expected relationship between Gibbs energy of ATP hydrolysis and muscle performance. Biophys Chem 54: 137142, 1995.[CrossRef][ISI][Medline]
32. Zhang J, Ugurbil K, From AH, Bache RJ. Myocardial oxygenation and high-energy phosphate levels during graded coronary hypoperfusion. Am J Physiol Heart Circ Physiol 280: H318H326, 2001.
This article has been cited by other articles:
![]() |
R. K. Dash and D. A. Beard Analysis of cardiac mitochondrial Na+-Ca2+ exchanger kinetics with a biophysical model of mitochondrial Ca2+ handing suggests a 3: 1 stoichiometry J. Physiol., July 1, 2008; 586(13): 3267 - 3285. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. Laurent Reply to Kemp: a clarification on the interpretation of muscular ATP synthase flux data obtained by 31P saturation transfer Am J Physiol Endocrinol Metab, March 1, 2008; 294(3): E643 - E644. [Full Text] [PDF] |
||||
![]() |
G. J. Kemp The interpretation of abnormal 31P magnetic resonance saturation transfer measurements of Pi/ATP exchange in insulin-resistant skeletal muscle Am J Physiol Endocrinol Metab, March 1, 2008; 294(3): E640 - E642. [Full Text] [PDF] |
||||
![]() |
N. Lai, G. M. Saidel, B. Grassi, L. B. Gladden, and M. E. Cabrera Model of oxygen transport and metabolism predicts effect of hyperoxia on canine muscle oxygen uptake dynamics J Appl Physiol, October 1, 2007; 103(4): 1366 - 1378. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. P. Smith, E. J. Crampin, S. A. Niederer, J. B. Bassingthwaighte, and D. A. Beard Computational biology of cardiac myocytes: proposed standards for the physiome J. Exp. Biol., May 1, 2007; 210(9): 1576 - 1583. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |