Investigating myotoxicity following Australian red-bellied black snake (Pseudechis porphyriacus) envenomation

• Published in: PloS one 2021-09, Vol.16 (9), p.e0256653
• Main Authors: Sanhajariya, Suchaya, Duffull, Stephen B, Isbister, Geoffrey K
• Format: Article
• Language: English
• Subjects: Antivenom; Biology and Life Sciences; Bites and stings; Confidence intervals; Creatine; Creatine kinase; Ethics; Half-life; Health services; Kinases; Mathematical models; Medicine and Health Sciences; Parameter estimation; Patients; Pharmacodynamics; Pharmacokinetics; Pharmacology; Pharmacy; Physical Sciences; Prevention; Pseudechis porphyriacus; Snake bites; Snakes; Toxicology; Toxins; Venom
• ISSN: 1932-6203; 1932-6203
• DOI: 10.1371/journal.pone.0256653

Background Myotoxicity is one of the common clinical manifestations of red-bellied black snake (Pseudechis porphyriacus) envenomation characterised by elevated creatine kinase (CK) concentrations of greater than 1000 U/L. This study aimed to investigate the occurrence of myotoxicity in patients following envenomation. Methods/Principal findings Patient characteristics and serial blood samples (timed venom concentrations and CK concentrations, pre- and post- antivenom) from 114 patients (median age 41, 2-90y; 80 male) were extracted from the Australian Snakebite Project database. Patients were categorised into three groups based on peak CK concentrations [no myotoxicity (10,000 U/L)]. The odds of (mild or severe) myotoxicity was lower in patients that received early antivenom (within 6 hours post-bite) compared to those that received late or no antivenom (odd ratio was 0.186; 95% confidence interval, 0.052-0.664). A population pharmacokinetic-pharmacodynamic (PKPD) model was developed to describe the relationship between the time course of venom (a mixture of toxins) and effect (elevated CK). In addition, a kinetic-pharmacodynamic (KPD) model was developed to describe the relationship between time course of a theoretical toxin and effect. Model development and parameter estimation was performed using NONMEM v7.3. No single set of parameter values from either the PKPD or KPD models were found that could accurately describe the time course of different levels of severity of myotoxicity. The predicted theoretical toxin half-life from the KPD model was 11 ± 3.9 hours compared to the half-life of venom of 5.3 ± 0.36 hours. This indicates that the putative causative toxin's concentration-time profile does not parallel that of venom. Conclusion Early antivenom administration reduces the incidence of myotoxicity. The venom concentration profile does not appear to be the driver for myotoxicity following envenomation. Additional factors that affect the sensitivity of the patient to snake venom/toxins must be explored to understand the relationship with myotoxicity.