Abstract
The SARS-CoV-2 main protease (Mpro) is a validated therapeutic target for inhibiting viral replication. Few compounds have advanced clinically, underscoring the difficulty in optimizing both target affinity and drug-like properties. To address this challenge, we integrated machine learning (ML), molecular docking, and molecular dynamics (MD) simulations to investigate the balance between pharmacodynamic (PD) and pharmacokinetic (PK) properties in Mpro inhibitor design. We developed ML models to classify Mpro inhibitors based on experimental IC(50) data, combining molecular descriptors with structural insights from MD simulations. Our Support Vector Machine (SVM) model achieved strong performance (training accuracy = 0.84, ROC AUC = 0.91; test accuracy = 0.79, ROC AUC = 0.86), while our Logistic Regression model (training accuracy = 0.78, ROC AUC = 0.85; test accuracy = 0.76, ROC AUC = 0.83). Notably, PK descriptors often exhibited opposing trends to binding affinity: hydrophilic features enhanced binding affinity but compromised PK properties, whereas hydrogen bonding, hydrophobic, and π-π interactions in Mpro subsites S2 and S3/S4 are fundamental for binding affinity. Our findings highlight the need for a balanced approach in Mpro inhibitor design, strategically targeting these subsites may balance PD and PK properties. For the first time, we demonstrate antagonistic trends between pharmacokinetic (PK) and pharmacodynamic (PD) features through the integrated application of ML/MD. This study provides a computational framework for rational Mpro inhibitors, combining ML and MD to investigate the complex interplay between enzyme inhibition and drug likeness. These insights may guide the hit-to-lead optimization of the novel next-generation Mpro inhibitors of SARS-CoV-2 with preclinical and clinical potential.