{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,9,4]],"date-time":"2026-09-04T05:16:26Z","timestamp":1788498986652,"version":"build-2803163510"},"reference-count":53,"publisher":"Frontiers Media SA","license":[{"start":{"date-parts":[[2026,9,4]],"date-time":"2026-09-04T00:00:00Z","timestamp":1788480000000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"content-domain":{"domain":["frontiersin.org"],"crossmark-restriction":true},"short-container-title":["Front. Physiol."],"abstract":"<jats:p>Prehospital pain management is challenging and is expected to become even more complex in future large-scale combat operations, where mass-casualty events and evacuation delay may be inevitable. An improved understanding of how analgesic drugs affect the physiological response of trauma patients can enhance treatment efficacy of combat casualties. We previously developed and validated a cardio-respiratory (CR) model for humans that accounts for vital-sign responses to hemorrhagic injuries, resuscitation with six fluid types, airway obstruction, and ketamine analgesia. Here, we extended the model to include the effect of fentanyl on vital signs by integrating existing fentanyl pharmacokinetic-pharmacodynamic models with the neuronal controller of the model. We calibrated and validated the extended model using experimental data from eight studies involving intravenous fentanyl administration (0.71\u201350.00 \u03bcg\/kg) to healthy humans and swine with hemorrhagic injury. The model predictions reasonably captured the trend of the experimental data, with root mean square errors (RMSEs) between model predictions and measured data of 0.83 L\/min for minute ventilation (MV), 1.20 mmHg for end-tidal carbon dioxide, 0.09 L for tidal volume, and 1.61 mmHg for mean arterial pressure, all of which were within 3\u201312% of their baseline values. For plasma fentanyl concentration, we obtained RMSEs of 0.79 \u03bcg\/L in humans and 25.83 \u03bcg\/L in swine. In simulations, we observed that as hemorrhage increased from 0 to 40% of blood volume, the fentanyl-induced decrease in MV increased from 19 to 34% of its baseline value prior to administration due to reduced fentanyl clearance. Similarly, in simulations of airway obstruction, the fentanyl-induced decrease in MV was 26% of its baseline value after a 100% obstruction compared to only 19% for a no-obstruction condition. Given that most combat casualties receive either fentanyl or ketamine for pain management, the ability to predict and quantify the physiological effects of these drugs will allow us to generate relevant synthetic datasets of diverse battlefield scenarios.<\/jats:p>","DOI":"10.3389\/fphys.2026.1907767","type":"journal-article","created":{"date-parts":[[2026,9,4]],"date-time":"2026-09-04T04:27:08Z","timestamp":1788496028000},"update-policy":"https:\/\/doi.org\/10.3389\/crossmark-policy","source":"Crossref","is-referenced-by-count":0,"title":["A model to simulate cardio-respiratory responses to fentanyl analgesia after traumatic injury"],"prefix":"10.3389","volume":"17","author":[{"given":"Varghese","family":"Kurian","sequence":"first","affiliation":[{"name":"Department of Defense Biotechnology High Performance Computing Software Applications Institute, Defense Health Agency Research & Development, Medical Research and Development Command","place":["Fort Detrick, MD, United States"]},{"name":"The Henry M. Jackson Foundation for the Advancement of Military Medicine, Inc.","place":["Bethesda, MD, United States"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Xin","family":"Jin","sequence":"additional","affiliation":[{"name":"Department of Defense Biotechnology High Performance Computing Software Applications Institute, Defense Health Agency Research & Development, Medical Research and Development Command","place":["Fort Detrick, MD, United States"]},{"name":"The Henry M. Jackson Foundation for the Advancement of Military Medicine, Inc.","place":["Bethesda, MD, United States"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Anders","family":"Wallqvist","sequence":"additional","affiliation":[{"name":"Department of Defense Biotechnology High Performance Computing Software Applications Institute, Defense Health Agency Research & Development, Medical Research and Development Command","place":["Fort Detrick, MD, United States"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Jaques","family":"Reifman","sequence":"additional","affiliation":[{"name":"Department of Defense Biotechnology High Performance Computing Software Applications Institute, Defense Health Agency Research & Development, Medical Research and Development Command","place":["Fort Detrick, MD, United States"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Sridevi","family":"Nagaraja","sequence":"additional","affiliation":[{"name":"Department of Defense Biotechnology High Performance Computing Software Applications Institute, Defense Health Agency Research & Development, Medical Research and Development Command","place":["Fort Detrick, MD, United States"]},{"name":"The Henry M. Jackson Foundation for the Advancement of Military Medicine, Inc.","place":["Bethesda, MD, United States"]}],"role":[{"vocabulary":"crossref","role":"author"}]}],"member":"1965","published-online":{"date-parts":[[2026,9,4]]},"reference":[{"key":"B1","doi-asserted-by":"publisher","first-page":"637","DOI":"10.1002\/cpt.2027","article-title":"Tolerance to opioid\u2010induced respiratory depression in chronic high\u2010dose opioid users: a model\u2010based comparison with opioid\u2010naive individuals","volume":"109","author":"Algera","year":"2021","journal-title":"Clin. Pharmacol. Ther."},{"key":"B2","doi-asserted-by":"publisher","first-page":"8","DOI":"10.1213\/00000539-199001000-00003","article-title":"Differences in magnitude and duration of opioid-induced respiratory depression and analgesia with fentanyl and sufentanil","volume":"70","author":"Bailey","year":"1990","journal-title":"Anesth. Analgesia"},{"key":"B3","doi-asserted-by":"publisher","first-page":"114","DOI":"10.1038\/s43856-024-00536-5","article-title":"Whole body physiology model to simulate respiratory depression of fentanyl and associated naloxone reversal","volume":"4","author":"Baird","year":"2024","journal-title":"Commun. Med."},{"key":"B4","doi-asserted-by":"publisher","first-page":"583","DOI":"10.1007\/s10928-020-09712-1","article-title":"The influence of cardiac output on propofol and fentanyl pharmacokinetics and pharmacodynamics in patients undergoing abdominal aortic surgery","volume":"47","author":"Bienert","year":"2020","journal-title":"J. Pharmacokinet. Pharmacodyn."},{"key":"B5","doi-asserted-by":"publisher","first-page":"912","DOI":"10.1097\/ALN.0000000000002160","article-title":"Influence of cardiac output on the pharmacokinetics of sufentanil in anesthetized pigs","volume":"128","author":"Birkholz","year":"2018","journal-title":"Anesthesiology"},{"key":"B6","doi-asserted-by":"publisher","first-page":"1337","DOI":"10.1152\/japplphysiol.01048.2020","article-title":"Acute overventilation does not cause lung damage in moderately hemorrhaged swine","volume":"130","author":"Blackburn","year":"2021","journal-title":"J. Appl. Physiol."},{"key":"B7","doi-asserted-by":"publisher","DOI":"10.14814\/phy2.14246","article-title":"Testing of novel spectral device sensor in swine model of airway obstruction","volume":"7","author":"Blackburn","year":"2019","journal-title":"Physiol. Rep."},{"key":"B8","doi-asserted-by":"publisher","DOI":"10.1136\/tsaco-2026-002259","article-title":"Ghost in the Machine: leveraging artificial intelligence in the trauma bay","volume":"11","author":"Burke","year":"2026","journal-title":"Trauma Surg. Acute Care Open"},{"key":"B9","doi-asserted-by":"publisher","first-page":"443","DOI":"10.30802\/AALAS-CM-18-000149","article-title":"Ethical and IACUC considerations regarding analgesia and pain management in laboratory rodents","volume":"69","author":"Carbone","year":"2019","journal-title":"Comp. Med."},{"key":"B11","first-page":"1467","article-title":"Pharmacokinetic and pharmacodynamic modeling in BioGears","author":"Clipp","year":"2016"},{"key":"B10","volume-title":"Tactical Combat Casualty Care (TCCC) Guidelines","year":"2024"},{"key":"B12","doi-asserted-by":"publisher","DOI":"10.1155\/2021\/6640638","article-title":"Seven mathematical models of hemorrhagic shock","volume":"2021","author":"Curcio","year":"2021","journal-title":"Comput. Math. Methods Med."},{"key":"B13","doi-asserted-by":"publisher","first-page":"81","DOI":"10.1097\/00000542-199901000-00013","article-title":"Volume kinetics of Ringer's solution in hypovolemic volunteers","volume":"90","author":"Drobin","year":"1999","journal-title":"Anesthesiology"},{"key":"B14","doi-asserted-by":"publisher","first-page":"S431","DOI":"10.1097\/TA.0b013e3182755dcc","article-title":"Death on the battlefield, (2001\u20132011): implications for the future of combat casualty care","volume":"73","author":"Eastridge","year":"2012","journal-title":"J. Trauma Acute Care Surg."},{"key":"B15","doi-asserted-by":"publisher","first-page":"156","DOI":"10.1097\/00000542-199907000-00024","article-title":"Fentanyl pharmacokinetics in hemorrhagic shock: a porcine model","volume":"91","author":"Egan","year":"1999","journal-title":"Anesthesiology"},{"key":"B16","doi-asserted-by":"publisher","first-page":"1380","DOI":"10.1345\/aph.1L168","article-title":"Pharmacokinetics and pharmacodynamics of intranasal versus intravenous fentanyl in patients with pain after oral surgery","volume":"42","author":"Foster","year":"2008","journal-title":"Ann. Pharmacother."},{"key":"B17","doi-asserted-by":"publisher","first-page":"160","DOI":"10.4103\/0970-9185.182085","article-title":"Intravenous sub-anesthetic ketamine for perioperative analgesia","volume":"32","author":"Gorlin","year":"2016","journal-title":"J. Anaesthesiology Clin. Pharmacol."},{"key":"B18","doi-asserted-by":"crossref","first-page":"19","DOI":"10.1007\/978-0-85729-148-6_2","article-title":"Euler\u2019s method","volume-title":"Numerical Methods for Ordinary Differential Equations: Initial Value Problems","author":"Griffiths","year":"2010"},{"key":"B19","doi-asserted-by":"publisher","first-page":"674","DOI":"10.1177\/0091270007299756","article-title":"Fentanyl clearance and volume of distribution are increased in patients with major burns","volume":"47","author":"Han","year":"2007","journal-title":"J. Clin. Pharmacol."},{"key":"B20","doi-asserted-by":"publisher","first-page":"S150","DOI":"10.1097\/TA.0000000000004722","article-title":"Assessment and training of Ukrainian trauma and combat casualty care via international symposia","volume":"99","author":"Holcomb","year":"2025","journal-title":"J. Trauma Acute Care Surg."},{"key":"B21","doi-asserted-by":"publisher","first-page":"108","DOI":"10.1093\/milmed\/usac211","article-title":"Patterns of palliation: a review of casualties that received pain management before reaching role 2 in Afghanistan","volume":"188","author":"Hudson","year":"2023","journal-title":"Military Med."},{"key":"B22","doi-asserted-by":"publisher","DOI":"10.3389\/fphys.2024.1327948","article-title":"AI algorithm for personalized resource allocation and treatment of hemorrhage casualties","volume":"15","author":"Jin","year":"2024","journal-title":"Front. Physiol."},{"key":"B23","doi-asserted-by":"publisher","DOI":"10.3389\/fphys.2025.1699315","article-title":"A model to simulate human cardio-respiratory responses to airway obstruction","volume":"16","author":"Jin","year":"2025","journal-title":"Front. Physiol."},{"key":"B24","doi-asserted-by":"publisher","DOI":"10.1002\/cnm.3662","article-title":"Development and validation of a mathematical model to simulate human cardiovascular and respiratory responses to battlefield trauma","volume":"39","author":"Jin","year":"2023","journal-title":"Int. J. Numer. Methods Biomed. Eng."},{"key":"B25","doi-asserted-by":"publisher","DOI":"10.3389\/fphys.2025.1613874","article-title":"A model to simulate human cardio-respiratory responses to different fluid resuscitation treatments after hemorrhagic injury","volume":"16","author":"Kurian","year":"2025","journal-title":"Front. Physiol."},{"key":"B26","doi-asserted-by":"publisher","DOI":"10.3389\/fphys.2025.1654768","article-title":"A model to simulate human cardio-respiratory responses to ketamine analgesia","volume":"16","author":"Kurian","year":"2025","journal-title":"Front. Physiol."},{"key":"B27","doi-asserted-by":"publisher","DOI":"10.1371\/journal.pone.0165461","article-title":"The effects of low-dose ketamine on acute pain in an emergency setting: a systematic review and meta-analysis","volume":"11","author":"Lee","year":"2016","journal-title":"PloS One"},{"key":"B28","doi-asserted-by":"publisher","first-page":"1","DOI":"10.1186\/s13049-025-01434-y","article-title":"Human and animal models for studying hemorrhagic shock","volume":"33","author":"Lie","year":"2025","journal-title":"Scandinavian J. Trauma Resuscitation Emergency Med."},{"key":"B29","doi-asserted-by":"publisher","first-page":"1115","DOI":"10.1109\/TBME.2004.827344","article-title":"Opioid-induced respiratory depression: a mathematical model for fentanyl","volume":"51","author":"Magosso","year":"2004","journal-title":"IEEE Trans. Biomed. Eng."},{"key":"B30","doi-asserted-by":"publisher","first-page":"1020","DOI":"10.1002\/cpt.2696","article-title":"Development of a translational model to assess the impact of opioid overdose and naloxone dosing on respiratory depression and cardiac arrest","volume":"112","author":"Mann","year":"2022","journal-title":"Clin. Pharmacol. Ther."},{"key":"B31","doi-asserted-by":"publisher","first-page":"215","DOI":"10.1016\/1350-4533(95)95712-j","article-title":"Validation of a computer model of haemorrhage and transcapillary refill","volume":"17","author":"Mardel","year":"1995","journal-title":"Med. Eng. Phys."},{"key":"B32","doi-asserted-by":"publisher","first-page":"37","DOI":"10.1046\/j.1365-2125.1998.00035.x","article-title":"Pulmonary administration of aerosolised fentanyl: pharmacokinetic analysis of systemic delivery","volume":"46","author":"Mather","year":"1998","journal-title":"Br. J. Clin. Pharmacol."},{"key":"B33","doi-asserted-by":"publisher","first-page":"965","DOI":"10.1046\/j.1365-2044.1998.00519.x","article-title":"The effect of low\u2010dose ketamine on fentanyl\u2010induced respiratory depression","volume":"53","author":"Mildh","year":"1998","journal-title":"Anaesthesia"},{"key":"B34","doi-asserted-by":"publisher","first-page":"939","DOI":"10.1097\/00000539-200110000-00028","article-title":"The concentration-effect relationship of the respiratory depressant effects of alfentanil and fentanyl","volume":"93","author":"Mildh","year":"2001","journal-title":"Anesth. Analgesia"},{"key":"B35","doi-asserted-by":"publisher","first-page":"283","DOI":"10.1038\/nrn2606","article-title":"Animal models of pain: progress and challenges","volume":"10","author":"Mogil","year":"2009","journal-title":"Nat. Rev. Neurosci."},{"key":"B36","doi-asserted-by":"publisher","DOI":"10.1172\/jci.insight.156973","article-title":"Modeling buprenorphine reduction of fentanyl-induced respiratory depression","volume":"7","author":"Olofsen","year":"2022","journal-title":"JCI Insight"},{"key":"B37","doi-asserted-by":"publisher","first-page":"598","DOI":"10.1093\/bja\/81.4.598","article-title":"The CYP 3A4 inhibitor itraconazole has no effect on the pharmacokinetics of iv fentanyl","volume":"81","author":"Palkama","year":"1998","journal-title":"Br. J. Anaesthesia"},{"key":"B38","doi-asserted-by":"publisher","first-page":"52","DOI":"10.1177\/0310057X211070292","article-title":"Clinical implications of opioid-induced ventilatory impairment","volume":"50","author":"Pattullo","year":"2022","journal-title":"Anaesthesia Intensive Care"},{"key":"B39","doi-asserted-by":"publisher","DOI":"10.1016\/j.imu.2020.100301","article-title":"Validation of new method for predicting human skin temperatures during cold exposure: The Cold Weather Ensemble Decision Aid (CoWEDA)","volume":"18","author":"Potter","year":"2020","journal-title":"Inf. Med. Unlocked"},{"key":"B40","doi-asserted-by":"publisher","first-page":"17","DOI":"10.1093\/milmed\/usad253","article-title":"Prehospital lessons from the war in Ukraine: damage control resuscitation and surgery experiences from point of injury to role 2","volume":"189","author":"Quinn","year":"2024","journal-title":"Military Med."},{"key":"B41","doi-asserted-by":"publisher","first-page":"744","DOI":"10.1080\/10903127.2017.1335814","article-title":"Multicenter, prospective study of prehospital administration of analgesia in the US combat theater of Afghanistan","volume":"21","author":"Schauer","year":"2017","journal-title":"Prehospital Emergency Care"},{"key":"B42","doi-asserted-by":"publisher","DOI":"10.1038\/s41598-026-57039-7","article-title":"A 3-D whole-body human thermoregulatory model to simulate extreme cold stress","author":"Singh","year":"2026","journal-title":"Sci. Rep"},{"key":"B43","doi-asserted-by":"publisher","first-page":"44534","DOI":"10.1038\/s41598-025-28021-6","article-title":"A reinforcement learning algorithm to optimize resource utilization in combat casualty care","volume":"15","author":"Subramaniyan","year":"2025","journal-title":"Sci. Rep."},{"key":"B44","doi-asserted-by":"publisher","first-page":"212","DOI":"10.1093\/bja\/aem165","article-title":"Crystalloid infusion rate during fluid resuscitation from acute haemorrhage","volume":"99","author":"Tatara","year":"2007","journal-title":"Br. J. Anaesthesia"},{"key":"B45","unstructured":"ACTIQ Prescribing Information\n          \n          2023"},{"key":"B46","doi-asserted-by":"publisher","first-page":"341","DOI":"10.1016\/j.jpain.2005.01.351","article-title":"Establishing a link between heart rate and pain in healthy subjects: a gender effect","volume":"6","author":"Tousignant-Laflamme","year":"2005","journal-title":"J. Pain"},{"key":"B47","doi-asserted-by":"publisher","DOI":"10.1186\/1471-2253-5-2","article-title":"Investigation of the potentiation of the analgesic effects of fentanyl by ketamine in humans: a double-blinded, randomised, placebo controlled, crossover study of experimental pain [ISRCTN83088383","volume":"5","author":"Tucker","year":"2005","journal-title":"BMC Anesthesiology"},{"key":"B48","doi-asserted-by":"publisher","first-page":"1171","DOI":"10.1097\/ALN.0000000000005710","article-title":"Fentanyl-induced ventilatory depression: population pharmacokinetic\u2013pharmacodynamic framework for evaluation of opioid-induced ventilatory depression","volume":"143","author":"van Lemmen","year":"2025","journal-title":"Anesthesiology"},{"key":"B49","doi-asserted-by":"publisher","DOI":"10.1093\/sleep\/zsad080","article-title":"2B-Alert App 2.0: personalized caffeine recommendations for optimal alertness","volume":"46","author":"Vital-Lopez","year":"2023","journal-title":"Sleep"},{"key":"B50","doi-asserted-by":"publisher","first-page":"600","DOI":"10.1080\/10903127.2023.2172493","article-title":"Comparing the effects of low-dose ketamine, fentanyl, and morphine on hemorrhagic tolerance and analgesia in humans","volume":"27","author":"Watso","year":"2023","journal-title":"Prehospital Emergency Care"},{"key":"B51","doi-asserted-by":"publisher","first-page":"S93","DOI":"10.1097\/TA.0000000000002604","article-title":"Fentanyl impairs but ketamine preserves the microcirculatory response to hemorrhage","volume":"89","author":"Xiang","year":"2020","journal-title":"J. Trauma Acute Care Surg."},{"key":"B52","doi-asserted-by":"publisher","first-page":"50","DOI":"10.1038\/sj.clpt.6100025","article-title":"Mechanism\u2010based PK\/PD modeling of the respiratory depressant effect of buprenorphine and fentanyl in healthy volunteers","volume":"81","author":"Yassen","year":"2007","journal-title":"Clin. Pharmacol. Ther."},{"key":"B53","article-title":"Propaq M operator's guide","year":"2016"}],"container-title":["Frontiers in Physiology"],"original-title":[],"link":[{"URL":"https:\/\/www.frontiersin.org\/articles\/10.3389\/fphys.2026.1907767\/full","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2026,9,4]],"date-time":"2026-09-04T04:27:10Z","timestamp":1788496030000},"score":1,"resource":{"primary":{"URL":"https:\/\/www.frontiersin.org\/articles\/10.3389\/fphys.2026.1907767\/full"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2026,9,4]]},"references-count":53,"alternative-id":["10.3389\/fphys.2026.1907767"],"URL":"https:\/\/doi.org\/10.3389\/fphys.2026.1907767","relation":{},"ISSN":["1664-042X"],"issn-type":[{"value":"1664-042X","type":"electronic"}],"subject":[],"published":{"date-parts":[[2026,9,4]]},"article-number":"1907767"}}