M: Musculo-skeletal system

ATC group M for modelers: four unrelated PD paradigms under one letter, with notes held at M01, M03, M04 and M05.
Modified

September 22, 2026

Drugs for joint, bone and muscle disease, plus the neuromuscular blockers given under anesthesia [1]. Published models here span a subjective pain score [2] and an evoked muscle twitch [1].

In 1938 Richard Gill returned from the Amazon with 11 kg of curare, from which Squibb isolated a standardized extract; Harold Griffith and Enid Johnson gave it in the operating room in 1942 [3,4].

Table 1: ATC M level 2 codes
Level 2 code Description In plain terms Example drug
M01 Antiinflammatory and antirheumatic products Ibuprofen (M01AE01) / naproxen (M01AE02)
M02 Topical products for joint and muscular pain Topical diclofenac (M02AA15) / topical ketoprofen (M02AA10)
M03 Muscle relaxants Rocuronium bromide (M03AC09) for neuromuscular block / baclofen (M03BX01) for spasticity
M04 Antigout preparations Allopurinol (M04AA01) / febuxostat (M04AA03)
M05 Drugs for treatment of bone diseases Alendronic acid (M05BA04) / denosumab (M05BX04)
M09 Other drugs for disorders of the musculo-skeletal system leftover bone, joint and muscle drugs Residual catch-all group, but it now carries real drug-level standards for rare neuromuscular disease: nusinersen (M09AX07) and risdiplam (M09AX10) for spinal muscular atrophy

Modeling notes

Endpoint. Each level 2 code brings its own: WOMAC pain in M01 [5], serum urate in M04 [6], bone markers and density in M05 [7], evoked twitch in M03 [8].

Easy to overlook. For M02 the relevant exposure is local, in skin and synovial fluid, and dermal PBPK is what links measured plasma concentration to that site, so bioequivalence is argued at the site of action alongside, not instead of, a plasma study [9].

Mechanistic extrapolation. A multiscale model of calcium homeostasis and bone remodeling [10] was extended to predict lumbar spine density under denosumab, including the decline after discontinuation [11]. For M02, dermal PBPK replaced a comparative clinical endpoint study in an FDA generic approval of topical diclofenac gel [9].

References

[1]
Sheiner LB, Stanski DR, Vozeh S, Miller RD, Ham J. Simultaneous modeling of pharmacokinetics and pharmacodynamics: Application to d-tubocurarine. Clinical Pharmacology & Therapeutics 1979;25:358–71. https://doi.org/10.1002/cpt1979253358.
[2]
Mandema JW, Stanski DR. Population pharmacodynamic model for ketorolac analgesia. Clinical Pharmacology & Therapeutics 1996;60:619–35. https://doi.org/10.1016/S0009-9236(96)90210-6.
[3]
Elsherbini N, Backman SB. The connection between Dr. Harold Griffith and Richard Gill: New insights into the history of curare use in anesthesia from previously unknown correspondence. Canadian Journal of Anesthesia 2024;71:1664–71. https://doi.org/10.1007/s12630-024-02750-0.
[4]
Griffith HR, Johnson GE. The use of curare in general anesthesia. Anesthesiology 1942;3:418–20. https://doi.org/10.1097/00000542-194207000-00006.
[5]
Wen X, Luo J, Mai Y, Li Y, Cao Y, Li Z, et al. Placebo response to oral administration in osteoarthritis clinical trials and its associated factors: A model-based meta-analysis. JAMA Network Open 2022;5:e2235060. https://doi.org/10.1001/jamanetworkopen.2022.35060.
[6]
Muraki S, Moriki K, Shigematsu S, Fukae M, Kakara M, Yamashita D, et al. Population pharmacodynamic analysis of uric acid-lowering effects of febuxostat based on electronic medical records in two hospitals. The Journal of Clinical Pharmacology 2018;58:304–13. https://doi.org/10.1002/jcph.1023.
[7]
Wu J, Wang C, Li GF, Tang ET, Zheng Q. Quantitative prediction of bone mineral density by using bone turnover markers in response to antiresorptive agents in postmenopausal osteoporosis: A model-based meta-analysis. British Journal of Clinical Pharmacology 2021;87:1175–86. https://doi.org/10.1111/bcp.14487.
[8]
Ji SH, Huh KY, Oh J, Jeong HJ, Jang YE, Kim EH, et al. Reversal of rocuronium-induced intense neuromuscular blockade by sugammadex in Korean children: A pharmacokinetic and pharmacodynamic analysis. Clinical and Translational Science 2023;16:92–103. https://doi.org/10.1111/cts.13429.
[9]
Tsakalozou E, Babiskin A, Zhao L. Physiologically-based pharmacokinetic modeling to support bioequivalence and approval of generic products: A case for diclofenac sodium topical gel, 1%. CPT: Pharmacometrics & Systems Pharmacology 2021;10:399–411. https://doi.org/10.1002/psp4.12600.
[10]
Peterson MC, Riggs MM. A physiologically based mathematical model of integrated calcium homeostasis and bone remodeling. Bone 2010;46:49–63. https://doi.org/10.1016/j.bone.2009.08.053.
[11]
Peterson MC, Riggs MM. Predicting nonlinear changes in bone mineral density over time using a multiscale systems pharmacology model. CPT: Pharmacometrics & Systems Pharmacology 2012;1:e14. https://doi.org/10.1038/psp.2012.15.