Optical Spectroscopic Analysis of Muscle Spasticity for Low-Level Laser Therapy (LLLT)
Optical Spectroscopic Analysis of Muscle Spasticity for Low-Level Laser Therapy (LLLT)
Journal of the Optical Society of Korea. 2011. Dec, 15(4): 373-379
Copyright ©2011, Optical Society of Korea
  • Received : September 09, 2011
  • Accepted : November 30, 2011
  • Published : December 25, 2011
Export by style
Cited by
About the Authors
Yeon Ui, Lee
Department of Biomedical Engineering, College of Medical Science, Catholic University of Daegu, Gyeongsan 712-702, Korea
Sangkwan, Lee
Department of Internal Medicine and Neuroscience, College of Oriental Medicine, WonKwang University, Iksan 570-749, Korea
Jong-In, Youn
Department of Biomedical Engineering, College of Medical Science, Catholic University of Daegu, Gyeongsan 712-702, Korea

Current therapeutic methods for suppressing muscle spasticity are intensive functional training, surgery, or pharmacological interventions. However, these methods have not been fully supported by confirmed efficacy due to the aggravation of the muscle spasticity in some patients. In this study, a combined system was developed to treat with a low-level laser and to monitor the region of the treatment using an optical spectroscopic probe that measures oxygen saturation and deoxygenation during low-level laser therapy (LLLT). The evaluation of the wavelength dependence for LLLT was performed using a Monte Carlo simulation and the results showed that the greatest amount of heat generation was seen in the deep tissue at λ = 830 nm. In the oxy- and deoxygenation measurements during and after the treatment, oxygen-Hb concentration was significantly increased in the laser-irradiated group when compared to the control group. These findings suggest that LLLT using λ = 830 nm may be of benefit in accelerating recovery of muscle spasticity. The combined system that we have developed can monitor the physiological condition of muscle spasticity during the laser treatment in real time and may also be applied to various myotonia conditions such as muscle fatigue, back-pain treatment/monitoring, and ulcer due to paralysis.
Cerebral palsy (CP) is a disorder of movement and posture resulting from a non-progressive injury to an immature brain. The most common type in CP is skeletal muscle spasticity that occurs secondary to upper motor neuron lesions and can result in serious complications for affected individuals [1] . It is known that the function of individuals with spastic muscles is severely compromised due to decreased range of motion, decreased voluntary strength, and increased muscle stiffness, but the basic mechanisms underlying the functional deficits after the development of spasticity are not clearly understood [2] . Since CP affects the muscular movements of its sufferers, many medications for CP target the muscles for relaxation [3] . The muscle relaxants used for treating CP work by stopping the muscles from contracting, but several side effects have been reported such as diarrhea, dizziness, drowsiness, fatigue, vague feeling of illness, and muscle weakness [3 - 5] . Physiotherapy such as massage has also thought to increase tissue blood circulation, thereby decreasing hypertonicity and enhancing recovery from muscle disorder [6] . Drust et al. demonstrated that there was an increase in skin and intramuscular temperature following physiotherapy, but it was not enough to increase deep intramuscular regions [7] .
Recently low-level laser therapy (LLLT) has been proposed as a therapeutic tool for treating CP. LLLT (<500 mW, <35 J/cm², typical wavelengths of 600-1300 nm) was introduced for the management of pain, wound healing and soft tissue injury, was developed in the 1960s [8] . LLLT is a form of phototherapy that involves the application of low-power, monochromatic, coherent light to injuries and lesions [9] . LLLT has been found to accelerate wound healing and reduce pain, possibly by stimulating oxidative phosphorylation and reducing inflammatory responses, thus producing several beneficial effects upon inflammation and healing [10 - 17] . In skeletal muscle, LLLT has been effective at promoting improvements in electrical activity, preventing fatigue induced by tetanic contractions and improving healing after traumatic injury [11] . At the molecular level, LLLT directly stimulates the activity of cytochrome oxidase in muscle fibers, triggering a series of biochemical cascades that improve cellular functions [12] .
Asagai et al. reported on the use of GaAlAs laser treatment ( λ = 810 nm) in a group of 1000 patients with CP. By observing infrared image and doppler ultrasonography they found that the laser reduces muscle spasm and increases the mobility of the muscles [13] . Leal Junior et al. demonstrated that LLLT with λ = 655 nm may delay the perceived onset of muscle fatigue and exhaustion through a reduction in oxidative stress or a decrease in the production of reactive oxygen species [14] . Walker et al. investigated the effects on improvement of blood flow through the autonomic nerve and the spinal reflex pathway after LLLT with λ = 632.5 nm and suggested that the mechanism of low-level laser therapy in spastic paralysis was an efferent inhibitory action and an effect on the autonomic nerve function [15] . Although many investigations have been performed for the evaluation of the effect of LLLT in muscle recovery, various laser wavelengths were treated and different results were presented. Since few studies have been done to compare photothermal effects depending on the wavelengths, Monte Carlo simulation was performed to compare different wavelengths for the effective LLLT in this study.
Recent studies have suggested that a diffuse optical signal measured over the near-infrared spectral range provides sufficient information of hemoglobin concentration and oxygen saturation in biological tissues non-invasively [18] . The advantages of near-infrared spectroscopic techniques (NIRS) include the possibility of obtaining information about the heterogeneity of local muscle metabolism with high temporal resolution [19] . Major chromophores that absorb light in living tissues are oxygenated and deoxygenated hemoglobin molecules (HbO 2 and Hb, respectively) with distinct optical-absorption spectra in the visible and nearinfrared ranges. The ability of NIRS to quantify changes in oxygenation is based on the difference in absorption spectra for oxygenated and deoxygenated heme groups [20] . A number of different studies have evaluated the influence of increased activity as well as decreased activity on muscle function using NIRS [20 - 22] . Costes et al. examined exercise induced adaptations in muscles determined by NIRS and showed that training did not change the pattern of muscle oxygenation, but a significant relationship was found between blood lactate and muscle oxygenation after exercise [21] . Tiziano Binzoni et al. measured local oxygen consumption in human muscle during rhythmically contractions of exercising skeletal muscle. They found that the saturation of the blood during the contraction phase was lower than that during the relaxation phase [22] .
For laser therapy applications, it is important to know the penetration and heat generation of laser in tissue in order to select optimal criteria such as wavelength, fluence, and beam size to afford the desired penetration for stimulating tissue healing in deep seated pathologies [23 , 24] . LLLT has demonstrated a number of biological and clinical efficacious results for muscle treatment, but this therapy has not been fully supported by scientific evidence due to the absence of quantitative assessments and treatment feedback data in real time. In this study, a NIRS system combined with laser treatment source was developed to monitor muscle tissue conditions such as oxygen saturation and deoxygenation during laser treatment using an optical spectroscopic probe.
- 2.1. Laser Wavelength Dependence for LLLT using Monte Carlo Simulation
A Monte Carlo simulation technique, based on the statistical nature of radiation interactions, has been widely applied in laser radiation transport studies [25 - 28] . Monte Carlo simulation methods trace the propagation of individual photons through a turbid medium until the photon experiences complete absorption by or escape from the medium [28] . In this study, the heat generation in skin was investigated by a computer simulation of appropriately weighted random absorption and scattering interactions. This mathematical modeling helps to determine the effective wavelength for deep tissue penetration in LLLT. The optical parameters of skin corresponding to the wavelengths of 660 nm, 750 nm, 830 nm, 980 nm, and 1064 nm are presented in Table 1 . The anisotropy factor for skin tissue, g, defined as the mean cosine of the deflection angle due to a scattering event, has a typical value of 0.9 [29] . These optical
Optical properties of skin for Monte Carlo simulation variablesμa: Absorption coefficient, μs: Scattering coefficient g : Anisotropic factor
Lager Image
Optical properties of skin for Monte Carlo simulation variables μa : Absorption coefficient, μs : Scattering coefficient g : Anisotropic factor
properties of the tissues for simulation were determined on the basis of literature values [29] .
- 2.2. The Combined Therapeutic Laser and NIRS System
Near infrared light (700 nm~1500 nm) can transmit through tissues about a few centimeters, depending on the separation between the light source and detector [30 , 31] . Some of the injected NIR light is lost as a result of scattering and absorption due to different structures and optical properties in the tissue. The attenuation of light between the source and detector can be formulated by using modified Beer-Lambert’s law [30 , 31] . It is mathematically expressed as
Lager Image
where, OD is the tissue’s optical density, I 0 is the incident light intensity, I is transmitted light intensity, ε is the extinction coefficient of the chromophores (mM -1 •cm -1 ), c is the concentration of the chromophores (mM -1 ), and L is the path length of light through tissue.
By detecting the light intensity changes at two or more wavelengths in the NIR region, the changes in concentrations of deoxy-hemoglobin ([Hb]) oxy-hemoglobin ([HbO 2 ]) at deep tissues can be estimated as follows
Lager Image
here OD λ is the optical density or absorbance at wave-length λ, ε Hb λ and ε HbO2 λ are the extinction coefficients at wavelength λ for molar concentrations of [Hb] and [HbO 2 ], respectively.
By employing two wavelengths, it follows that changes in [Hb] and [HbO 2 ] can be consequently given as
Lager Image
As seen in the Figure 1 , the system consists of a diode laser source (DL8142-201, Thorlabs, Newton, NJ, USA) for LLLT and a spectrometer (USB 4000, Ocean Optics, Dunedin, FL, USA) module to monitor dynamic changes of HbO 2 under the skin. The laser source is connected with an optical fiber at the center of the designed probe and the fibers for the light source and spectroscopic detection were connected at the sides of the probe. Because the thickness of the skin and adipose tissue layer in humans is approximately 3 mm the source-detector separation used in this study was 2 cm, which allowed for sufficient light penetration into the muscle [32] . The incident light source was generated by a tungsten-halogen lamp (HL-2000, Ocean Optics, Dunedin, FL, USA) and the acquired signal reflected from tissue was transmitted to the spectrometer with the spectral range from 300-1000 nm. A diode laser (Power =
Lager Image
The schematic diagram of the combined probe in the developed system.
Lager Image
The schematic diagram of the developed system that is combined with a therapeutic laser and spectroscopic system.
80 mW, λ= 830 nm, spot size = 3 mm) was used for the treatment of muscle spasticity ( Fig. 2 ). This wavelength provided high contrast and deep penetration in tissue that delivers enough laser fluence for stimulating chromophores in muscle based on Monte Carlo simulation.
- 2.3. Experimental Procedure
Twenty male and female subjects between 23 and 25 years of age agreed to participate in the study through a written informed consent. All subjects were healthy and physically active students with no recent experience with isokinetic eccentric exercise. Subjects agreed not to perform physical exercise or use any other therapeutic modalities during the data collection sessions. The subjects were randomly divided into two groups: the control group and the laserirradiated group. Both groups were instructed to perform repeated elbow flexion-extension from full elbow extension to 90° of flexion for 10 minutes. Prior to the experimental protocol, the relative changes of [HbO 2 ] and [Hb] concentration in normal state with the developed system were measured for all subjects. Immediately after flexion-extension exercise, the probe was placed on the trigger point in the brachioradial muscle of the left forearm, 4 cm distal to the elbow. A trigger point is defined as a point that elicits referred pain on deep palpation and demonstrates a lowered skin resistance in comparison with that of the surrounding tissue [33] . The probe was in contact with the skin at a right angle and the low-level laser at λ = 830 nm irradiated it with a power of 80 mW. During exercise and LLLT procedure, the relative changes in oxygen and deoxygenation concentrations in hemoglobin were measured simultaneously for the subjects to monitor muscle recovery through the developed system.”
Figure 3 shows the relative rate of heat generation versus optical penetration depth in skin with different laser wavelengths (λ= 660 nm, 750 nm, 830 nm, 980 nm, and 1064 nm). In Monte Carlo simulation with λ= 660 nm and 1064 nm, the greatest amount of heat is generated at
Lager Image
The relative heat generation vs. penetration depth in skin depending on the wavelengths for LLLT using Monte Carlo simulation (left), and the log scale plot (right).
the surface, but the heat is dramatically reduced as the penetration depth is increased. On the other hand, in λ= 750 nm, 830 nm, and 980 nm, the rate of heat generation is less at the surface, but higher in deep tissue. The highest heat generation in deep tissue was achieved at λ= 830 nm below a depth of 5 mm. According to this result, the maximum heat generation and deep penetration for skin was achieved at λ= 830 nm for effective LLLT.
Prior to the measurements of LLLT’s effect by a spectroscopic analysis, the feasibility test of the developed system was performed to see the hemoglobin concentration changes in muscle without low-level laser irradiation. The changes of oxy- ([HbO 2 ]) and deoxy-hemoglobin ([Hb]) concentrations were monitored by applying the occlusion of a cuff pressure around the upper arm. Figure 4 shows that the rapid increase of the cuff pressure induced rapid and substantial decrease in [HbO 2 ]. At the same time, [Hb] gradually increased as a result of continuous oxygen consumption in muscle and undelivered oxygen throughout the arterial occlusion by the cuff pressure. The typical reactive hyperemia following the release of cuff pressure was much higher than the normal state of [HbO 2 ]. This variation was primarily due to delays between venous and arterial occlusion during both pressure cuff and release cuff which leads to trapping of blood in the capillaries.
To evaluate the effect of LLLT, the subjects were divided into two groups, the control group and the laser irradiated group, and the relative changes in [Hb] and [HbO 2 ] concentrations in muscle were measured after exercise for both groups. Figure 5 shows the results of [Hb] for a single, typical member of each group. [Hb] had a significant difference between the two groups. In the laser-irradiated group, [Hb] concentration rapidly and dramatically decreased after the exercise. Conversely in the control group, [Hb] concentration decreased slowly over the time. Figure 6 shows data for a single, typical member of each group. The figure shows that [HbO 2 ] increased slowly to baseline level before exercise in the control group, whereas in the
Lager Image
The relative changes in [HbO2] and [Hb] concentrations from the application of a pressure cuff: before, during and after the blood flow occlusion on a forearm.
Lager Image
The relative changes in [Hb] concentrations: post exercise for one typical member of the control group (□), and one typical member of the laser irradiated group (X).
Lager Image
The relative changes in [HbO2] concentrations: post exercise for one typical member of the control group (□), and one typical member of the laser irradiated group (X).
laser-irradiated group, [HbO 2 ] increased rapidly toward their baseline level. All the participants showed similar tendencies within each group. The average relative changes in [Hb] before and after the treatment at 300 seconds were 0.014±0.005, showing the rapid decrease in deoxygenation of hemoglobin when compared to the control group. On the other hand, the average relative changes in [HbO 2 ] were 0.003±0.002, which represents the rapid recovery of hemoglobin oxygenation to the baseline after laser irradiation when compared to the control group.
In figure 7 , the relative changes in [Hb] and [HbO 2 ] concentration for the one of the participants showed that [Hb] concentration increased after the exercise began while [HbO 2 ] concentration decreased as expected. Once LLLT was irradiated after the exercise, the concentrations of [Hb] and [HbO 2 ] recovered rapidly to their baseline levels. It is expected that the intrinsic auto-regulation ability in muscle vasculature increases blood flow to meet the increased need for delivery and removal [34] . Since muscle contraction by intensive exercise or in patients with muscular disorders
Lager Image
The relative changes in [HbO2] and [Hb] concentrations before, during and after exercise with LLLT application.
have common symptoms such as exercise intolerance, fatigue and stiffness by insufficient blood circulation, the reduction of oxygen supply to muscles increases the rate of muscle fatigue and inadequate removal of metabolic waste products [35] . The increased oxygen delivery after low-level laser irradiation could result in microcirculation increase, lactic acid concentration reduction and muscle fatigue recovery [13] . LLLT has been proved to simulate tissue chromophores and an increase in blood flow, inducing electronically excited states that modify their redox potential and catalytic activity [12] . E. C. Leal Junior et al. also found that LLLT seemed to attenuate skeletal muscle fatigue and to reduce the muscle damage caused by titanic contractions induced by electrical stimulation [36] . Passarella et al. found that laser irradiation generates an extra electrochemical potential and an increase of ATP synthesis inside the isolated mitochondria when compared with a non-irradiated control group [37] .
The selected wavelength at λ = 830 nm based on Monte Carlo simulation in this study is especially attributed to a relatively oxidized form of cytochrome c oxidase that mediated the transfer of electrons from cytochrome c to molecular oxygen [38] . Kudoh et al. also observed a charge of Na-K ATPase in the rat saphenous nerve after treatment with a GaAlAs diode laser ( λ = 830 nm) and reported that this phenomenon may be considered very important for the attenuation of pain [39] .
From the results of this study in the control group without low-level laser irradiation, muscle fatigue condition lasted for minutes after exercise. However, in the laserirradiated group, muscle fatigue recovered quickly to normal muscle condition after exercise. The developed combined system that contains low-level laser source and NIRS system in this study evaluated the effect of therapeutic modality and monitored the physiological condition during the low-level laser treatment in real time. This study may lend insight into the heat generation event occurring inside tissue layers depending on the wavelengths and help effective phototherapy, LLLT.
In this study, the wavelength dependence of penetration depth in tissue was evaluated using the Monte Carlo simulation and a combined system that contains low-level laser treatment source and non-invasive optical monitoring system for the effective treatment of muscle spasticity was developed. The Monte Carlo simulation results show that the greatest amount of heat generation was observed in deep tissue at λ=830 nm. Using the developed system, healthy subjects were used for measurements of tissue condition like [HbO 2 ] and [Hb] concentration changes during LLLT after exercise. The [HbO 2 ] concentration was significantly increased in the laser-irradiated group and these findings suggest that LLLT using λ=830 nm may be of benefit in accelerating recovery of muscle spasticity. The developed system monitored the physiological condition of muscle spasm during the LLLT in real time and this system may also be applied to back-pain treatment/ monitoring, ulcer due to paralysis, and various myotonic conditions.
This research was supported by Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education, Science and Technology (2009-0065758)
Vaz D. V. , Mancini M. C. (2006) “Muscle stiffness and strength and their relation to hand function in children with hemiplegic cerebral palsy” Developmental Medicine & Child Neurology 48 728 - 733    DOI : 10.1017/S0012162206001563
Lieber R. L. , Steinman S. , Barash I. A. , Chambers H. (2004) “Structural and functional changes in spastic skeletal muscle” Muscle Nerve 29 615 - 27    DOI : 10.1002/mus.20059
Deyo R. A. , Bergman J. , Phillips W. R. (1996) “Drug therapy for back pain: Which drugs help which patients?,” Spine 21 2840 - 2850    DOI : 10.1097/00007632-199612150-00007
Chou R. , Peterson K. , Helfand M. (2004) “Comparative efficacy and safety of skeletal muscle relaxants for spasticity and musculoskeletal condition: a systematic review” J. Pain Symptom 28 140 - 175    DOI : 10.1016/j.jpainsymman.2004.05.002
Hulme A. , MacLennan W. J. , Ritchie R. T. (1985) “Baclofen in the elderly stroke patient its side effects and pharmacokinetics” Eur. J. Clin. Pharmacol 29 467 - 469    DOI : 10.1007/BF00613463
Mori H. , Ohsawa H. , Tanaka T. H. , Taniwaki E. , Leisman G. , Nishijo K. (2004) “Effect of massage on blood flow and muscle fatigue following isometric lumbar exercise” Med. Sci. Monit. 10 173 - 178
Drust B. , Atkinson G. , Gregson W. , French D. , Binningsley D. (2003) “The effects of massage on intra muscular temperature in the vastus lateralis in humans” Int. J. Sports Med. 24 395 - 399    DOI : 10.1055/s-2003-41182
Huang Y. Y , Chen A. C. , Carroll J. D. , Hamblin M. R. (2009) “Biphasic dose response in low level light therapy” Dose Response 7 358 - 383    DOI : 10.2203/dose-response.09-027.Hamblin
Sussai D. A. , Carvalho P. T. , Dourado D. M. , Belchior A. C. , dos Reis F. A , Pereira D. M. (2009) “Low-level laser therapy attenuates creatine kinase levels and apoptosis during forced swimming in rats” Lasers Med. Sci. 25 115 - 120
Medrado A. R. , Pugliese L. S. , Reis S. R. , Andrade Z. A. (2003) “Influence of low level laser therapy on wound healing and its biological action upon myofibroblasts” Lasers Surg. Med. 32 239 - 244    DOI : 10.1002/lsm.10126
Hayworth C. R. , Rojas J. C. , Padilla E. , Holmes G. M. , Sheridan E. C. , Gonzalez-Lima F. (2010) “In vivo low-level light therapy increases cytochrome oxidase in skeletal muscle” Photochem. Photobiol. 86 673 - 680    DOI : 10.1111/j.1751-1097.2010.00732.x
Karu T. (1999) “Primary and secondary mechanisms of action of visible to near-IR radiation on cells” J. Photochem. Photobiol. B: Biol. 49 1 - 17
Asagai Y. , Imakiire A. , Ohshiro T. (2000) “Thermographic effects of laser therapy in patients with cerebral palsy” Laser Therapy 12 12 - 15    DOI : 10.5978/islsm.12.12
Leal Junior E. C. , Lopes-Martins R. A. , Dalan F. , Ferrari M. , Sbabo F. M. , Generosi R. A. , Baroni B. M. , Penna S. C. , Iversen V. V. , Bjordal J. M. (2008) “Effect of 655-nm low-level laser therapy on exerciseinduced skeletal muscle fatigue in humans” Photomed. Laser Surg. 26 419 - 424    DOI : 10.1089/pho.2007.2160
Walker J. B. (1985) “Temporary suppression of clonus in humans by brief photostimulation. Brain research” Brain Res. 340 109 - 113    DOI : 10.1016/0006-8993(85)90779-6
MacDougall J. D. , Hicks A. L. , MacDonald J. R. , Mckelvie R. S. , Green H. J. , Smith K. M. (1998) “Muscle performance and enzymatic adaptations to sprint interval training” J. Appl. Physiol. 84 2138 - 2142    DOI : 10.1063/1.368275
Ahmaidi S. , Granier P. , Taoutaou Z. , Mercier J. , Dubouchaud H. , Prefaut C. (1996) “Effects of active recovery on plasma lactate and anaerobic power following repeated intensive exercise,” Med. Sci. Sports Exerc. 28 450 - 456    DOI : 10.1097/00005768-199604000-00009
McBride T. O. , Pogue B. W. , Gerety E. D. , Poplack S. B. , Osterberg U. L. , Paulsen K. D. (1999) “Spectroscopic diffuse optical tomography for the quantitative assessment of hemoglobin concentration and oxygen saturation in breast tissue” Appl. Opt. 38 5480 - 5490    DOI : 10.1364/AO.38.005480
Volianitis S. , Krustrup P. , Dawson E. , Secher N. H. (2003) “Arm blood flow and oxygenation on the transition from arm to combined arm and leg exercise in humans” J. Physiol. 547 641 - 648    DOI : 10.1113/jphysiol.2002.034496
Johns M. , Giller C. A. , Liu H. (2001) “Determination of hemoglobin oxygen saturation from turbid media using reflectance spectroscopy with small source-detector separations” Applied Spectroscopy 55 1686 - 1694    DOI : 10.1366/0003702011954026
Costes F. , Prieur F. , Reasson L. , Geyssant A. , Barthelemy J. C. , Denis C. (2001) “Influence of trainging on NIRS muscle oxygen saturation during submaximal exercise” Med. Sci. Sports Exercise 33 1484 - 1489    DOI : 10.1097/00005768-200109000-00010
Binzoni T. , Ngo L. , Hiltbrand E. , Springett R. , Delpy D. (2002) “Comsumption-temperature curves during rest and isometric exercise in human skeletal muscle” Comp. Biochem. Physiol. A 132 27 - 32
Karu T. (1999) “Primary and secondary mechanisms of action of visible to near-IR radiation on cells” J. Photochem. Photobiol. B 49 1 - 17
Carvalho C. M. , de Lacerda J. A. , dos Santos Neto F. P. , Cangussu M. C. T. , Marques A. M. C. , Pinheiro A. L. B. (2010) “Wavelength effect in temporomandibular joint pain : a clinical experience” Lasers Med. Sci. 25 229 - 232    DOI : 10.1007/s10103-009-0695-y
Dam J. S. , Andersen P. E. , Dalgaard T. , Fabricius P. E. (1998) “Determination of tissue optical properties from diffuse reflectance profiles by multivariate calibration” Appl. Opt. 37 772 - 778    DOI : 10.1364/AO.37.000772
Youn J. I. (2009) “A comparison of wavelength dependence for laser-assisted lipolysis effect using Monte Carlo simulation” J. Opt. Soc. Korea 13 267 - 271    DOI : 10.3807/JOSK.2009.13.2.267
Lee H. , Jeong Y. U. , Chan K. F. (2009) “The advent of laser therapies in dermatology and urology: underlying mechanisms, recent trends, and future” J. Opt. Soc. Korea 13 321 - 329    DOI : 10.3807/JOSK.2009.13.3.321
Wang L. , Jacques S. L. , Zheng L. (1995) “MCML-Monte Carlo modeling of light transport in multi-layered tissues” Computer Methods and Programs in Biomedicine 47 131 - 146    DOI : 10.1016/0169-2607(95)01640-F
Enwemeka C. S. (2001) “Attenuation and penetration depth of red 632.8 nm and invisible infrared 904 nm light in soft tissues” J. Laser Ther. 13 95 - 101
Bozkurt A. , Rosen A. , Rosen H. , Onaral B. (2005) “A portable near infrared spectroscopy system for bedside monitoring of newborn brain” BioMedical Engineering OnLine 4 1 - 11    DOI : 10.1186/1475-925X-4-1
Subudhi A. W. , Dimmen A. C. , Roach R. C. (2007) “Effects of acute hypoxia on cerebral and muscle oxygenation during incremental exercise” J. Appl. Physiol. 103 177 - 183    DOI : 10.1152/japplphysiol.01460.2006
Yu G. , Durduran T. , Zhou C. , Yodh A. G. (2005) “Timedependent blood flow and oxygenation in human skeletal muscles measured with noninvasive near-infrared diffuse optical spectroscopies” J. Biomed. Opt. 10 024027 -    DOI : 10.1117/1.1884603
Snyder M. L. , Bork C. , Bourbon B. , Trumbore D. (1986) “Effect of helium-neon laser on musculoskeletal trigger points” Phys. Ther. 66 1087 - 1090
Spector S. A. , Lemmer J. T. , Koffman B. M. , Fleisher T. A. , Feuerstein I. M. , Hurley B. F. , Dalakas M. C. (1997) “Safety and efficacy of strength training in patients with sporadic inclusion body myositis” Muscle Nerve. 20 1242 - 1248    DOI : 10.1002/(SICI)1097-4598(199710)20:10<1242::AID-MUS6>3.0.CO;2-C
Hepple R. T. , Hagen J. L. , Krause D. J. , Jackson C. C. (2003) “Aerobic power declines with aging in rat skeletal muscles perfused at matched convective O2delivery” J. Appl. Physiol. 94 744 - 751
Leal Junior E. C. , Lopes-Martins R. A , Vanin A. A. , Baroni B. M. , Grosselli D. , De Marchi T. , Iversen V. V. , Bjordal J. M. (2009) “Effect of 830 nm low-level laser therapy in exercise-induced skeletal muscle fatigue in humans” Lasers Med. Sci. 24 425 - 431    DOI : 10.1007/s10103-008-0592-9
Passarella S. , Casamassima E. , Molinari S. , Pastore D. , Quagliariello E. , Catalano I. M. , Cingolani A. (1984) “Increase of proton electrochemical potential and ATP synthesis in rat liver mitochondria irradiated in vitro by helium-neon laser” FEBS Letters 175 95 - 99    DOI : 10.1016/0014-5793(84)80577-3
Karu T. I. (2008) “Mitochondrial signaling in mammalian cells activated by red and near-IR radiation” Photochem. Photobiol. 84 1091 - 1099    DOI : 10.1016/0014-5793(84)80577-3
Kudoh C. , Inomata K. , Okajima K. , Motegi M. , Ohshiro T. (1988) “Low level laser therapy : pain attenuation mechanism” Laser ther. 10 3 - 6