Reproducible Computational Workflow for Drug Discovery to Standardize Network Pharmacology and Molecular Docking...

Network pharmacology and molecular docking are widely applied in drug discovery, yet fragmented workflows and inconsistent operations frequently undermine result reproducibility. Here, a standardized protocol integrating these approaches into a reproducible framework for drug screening and mechanism exploration is described, with the workflow organized into three sequential phases: data preparation, computational analysis, and validation. In the preparation phase, compound libraries from public databases are filtered via Absorption, Distribution, Metabolism, Excretion, and Toxicity (ADMET) criteria, including oral bioavailability, drug-likeness, and toxicity prediction, while potential therapeutic targets are obtained through target prediction and integration of disease-related databases to comprehensively identify drug-disease interaction candidates. In the computational analysis phase, overlapping targets undergo Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses and protein-protein interaction network topological analysis to identify core targets; molecular docking is configured with two standardized optional strategies with distinct advantages. The two-step progressive strategy employs AutoDock Vina for high-throughput preliminary screening of the compound library, followed by precise re-docking with YASARA, which eliminates false positives from high-throughput screening and generates protein-ligand complexes natively compatible with subsequent YASARA molecular dynamics (MD) simulations to avoid structural deviations caused by cross-software format conversion. The one-step strategy completes the full docking process via YASARA alone, which simplifies the operation workflow, improves experimental efficiency, and is fully applicable for specific research goals. In the validation phase, standardized MD simulations evaluate ligand-protein complex stability via core metrics of root mean square deviation (RMSD) and root mean square fluctuation (RMSF). This unified, reproducible pipeline enhances the reliability of network pharmacology and docking studies and facilitates cross-study comparisons in computational drug discovery.