Quantitative Structure–Activity Relationship (QSAR) Model for the Severity Prediction of Drug-Induced Rhabdomyolysis by Using Random Forest

Drug-induced rhabdomyolysis (DIR) is a rare and potentially life-threatening muscle injury that is characterized by low incidence and high risk. To our best knowledge, the performance of the current predictive models for the early detection of DIR is suboptimal because of the scarcity and dispersion of DIR cases. Therefore, on the basis of the curated drug information from the Drug-Induced Rhabdomyolysis Atlas (DIRA) database, we proposed a random forest (RF) model to predict the DIR severity of the marketed drugs. Compared with the state-of-art methods, our proposed model outperformed extreme gradient boosting, support vector machine, and logistic regression in distinguishing the Most-DIR concern drugs from the No-DIR concern drugs (Matthews correlation coefficient (MCC) and recall rate of our model were 0.46 and 0.81, respectively). Our model was subsequently applied to predicting the potentially serious DIR for 1402 drugs, which were reported to cause DIR by the postmarketing DIR surveillance data in the FDA Spontaneous Adverse Events Reporting System (FAERS). As a result, 62.7% (94) of drugs ranked in the top 150 drugs with the Most-DIR concerns in FAERS can be identified by our model. The top four drugs (odds ratio >30) including acepromazine, rapacuronium, oxyphenbutazone, and naringenin were correctly predicted by our model. In conclusion, the RF model can well predict the Most-DIR concern drug only based on the chemical structure information and can be a facilitated tool for early DIR detection.