Back to Journals » Drug Design, Development and Therapy » Volume 19
Standards for Computational Methods in Drug Design and Discovery: Simplified Guidance for Authors and Reviewers
Authors Ibrahim TM, Hussain M
, Boeckler FM
Received 10 October 2025
Accepted for publication 20 October 2025
Published 27 October 2025 Volume 2025:19 Pages 9649—9652
DOI https://doi.org/10.2147/DDDT.S573358
Checked for plagiarism Yes
Editor who approved publication: Prof. Dr. Georgios Panos
Tamer M Ibrahim,1,2 Muzammal Hussain,3,4 Frank M Boeckler5,6
1Department of Pharmaceutical Chemistry, Faculty of Pharmacy, Kafrelsheikh University, Kafrelsheikh, 33516, Egypt; 2Center for Informatics Science (CIS), School of Information Technology and Computer Science (ITCS), Nile University, Giza, Egypt; 3Department of Biochemistry and Molecular Pharmacology, New York University Grossman School of Medicine, New York, NY, 10016, USA; 4Howard Hughes Medical Institute, New York University Grossman School of Medicine, New York, NY, 10016, USA; 5Department of Pharmacy and Biochemistry, Eberhard Karls Universität Tübingen, Laboratory for Molecular Design and Pharmaceutical Biophysics, Institute of Pharmaceutical Sciences, Tübingen, 72076, Germany; 6Interfaculty Institute for Biomedical Informatics (IBMI), Eberhard Karls Universität Tübingen, Tübingen, 72076, Germany
Correspondence: Tamer M Ibrahim, Department of Pharmaceutical Chemistry, Faculty of Pharmacy, Kafrelsheikh University, Kafrelsheikh, 33516, Egypt, Email [email protected] Frank M Boeckler, Institute of Pharmaceutical Sciences, Department of Pharmacy and Biochemistry, Eberhard Karls Universität Tübingen, Auf der Morgenstelle 8, Tübingen, 72076, Germany, Tel +49 7071 29 74567, Email [email protected]
Part I: Policy Statement
Computational methods are integral to modern drug discovery and development. However, manuscripts that apply these approaches superficially or without methodological rigor undermine scientific progress. Drug Design, Development and Therapy will therefore only consider computational studies that are transparent, reproducible, validated and biologically meaningful.
- 2D-QSAR studies will be rejected immediately.
- Predictions of bioactivities, ADMET, pharmacokinetics, or toxicity without experimental or external validation will be rejected.
- Black-box applications of computational servers or uncritical pipelines will not be accepted.
- Docking, virtual screening, molecular dynamics (MD), quantum mechanics/molecular mechanics (QM/MM), artificial intelligence and machine learning (AI/ML), and free energy approaches must be conducted and reported according to established best practices.
Submissions not meeting these minimal requirements will not be sent for review.
Part II: Detailed Guidelines
QSAR and Molecular Modeling
2D-QSAR studies are obsolete and unacceptable. 3D-QSAR or more advanced models are permissible if applied to sufficiently large datasets, with rigorous alignment justification, proper cross-validation, and interpretable mechanistic insights.1 Molecular modeling may generate structural hypotheses but must not be overstated or presented as equivalent to experimental methods. Authors must highlight the model character of such work. Misrepresentation through selective omission or manipulative depictions will not be tolerated.
Docking and Virtual Screening
Docking studies are welcome when rationalizing structure–activity relationships or providing mechanistic support for experimental findings.2,3 Docking scores must not be reported as absolute binding energies. Binding site preparation, protonation states, tautomerization (e.g., histidine residues), and docking parameters must be fully described.4,5 Benchmarking of docking protocols or scoring functions—including ML scoring functions—is encouraged if robust metrics and clear recommendations are provided.6 Virtual screening campaigns must include experimental validation of proposed hits; unvalidated hit lists will not be accepted.
Artificial Intelligence and Machine Learning
The journal encourages the responsible use of AI/ML. Generative models for de novo drug design must be supported by preliminary experimental testing or confirmatory in silico validation such as MD simulations. Development of new ML models or scoring functions is welcome if trained on curated datasets, validated on independent sets, and evaluated with accepted performance metrics. Black-box models without interpretability or benchmarking will not be considered.7
MD and MM-PBSA/GBSA
MD simulations can provide insights into protein flexibility and ligand binding, but submissions based on single, short, or poorly prepared trajectories will be rejected. Best practices include initiating runs from high-quality structures, assigning correct protonation states, including relevant cofactors, ensuring production timescales are adequate, and using multiple replicas – if possible – to assess reproducibility.8 MM-PBSA/GBSA calculations may be used only for relative comparisons within consistent protocols and must be based on equilibrated ensembles. Considering absolute free energies of MM-PBSA/GBSA or single post-docking snapshot calculations is unacceptable.9 It is essential to consider what types of movements or structural changes are realistically observable on the timescale of the dynamics simulation. When interpreting the results of MD simulations any kind of bias, particularly unconscious bias, must be avoided at all costs. Video files with a meaningful graphical representation of major findings should be submitted as supplementary data or as part of the actual manuscript.
QM-Methods (HF, Post-HF, DFT, Semiempirical Methods) and QM/MM
QM and QM/MM calculations are valuable when appropriately and thoroughly applied. Authors must justify the choice of theory level, functional, and basis set, ideally benchmarking against higher-level methods. Reports of meaningless QM-derived descriptors (e.g., isolated orbital energies) will be rejected. FMO analysis (including depictions of HOMO/LUMOs) will only be tolerated, if it properly illustrates reactivities or reaction mechanisms with a suitable experimental basis. Depictions of molecular electrostatic potentials (MEPs) can be used to illustrate specific electrostatic features of molecular interactions, particularly when comparing similar compounds. However, they are not acceptable without proper justification and discussion. We encourage the use of QM methods for acceptable applications, such as evaluating protonation states, pKa values, conformational energy analysis, solvation energies, and key interactions in ligand–target complexes.10
Free Energy Methods
Free energy perturbation (FEP) and related alchemical methods are increasingly considered best-in-class for predicting relative binding affinities. These approaches will be accepted when system preparation, sampling, convergence analysis, and benchmarking against experimental data are rigorously documented. Submissions without error analysis or correlation to experiment will not be considered.11
Prospective Considerations for “Beyond Rule of Five (bRo5)” Modalities
The bRo5 space (MW > 500 Da) primarily includes macrocycles, constrained peptides, bulky natural products (e.g., paclitaxel, cyclosporine, etc.), and emerging proximity-based therapeutics such as PROTACs and other hetero-bifunctional agents designed to form ternary complexes for gain- or loss-of-function outcomes. In certain cases, it also encompasses molecular glue-like natural products with large molecular footprints (e.g., rapamycin). Owing to their size, conformational flexibility, chameleonicity, and multi-protein binding modes, these molecules often challenge conventional in silico methods.
Key considerations include:
- Adopting bRo5-tailored descriptors and ML models as conventional QSARs is unreliable in this chemical regimen.12–14
- Avoid over-interpreting single docking poses or raw docking scores, as they can be particularly misleading for highly flexible chemotypes such as macrocyclic and hetero-bifunctional ligands.15–17
- Avoid off-the-shelf (default) scoring functions as they may often fail; instead, calibrate or adapt scoring functions to capture the complexity of ternary systems.15,18
- Employ enhanced sampling methods (e.g., metadynamics, TTMD, etc.) to adequately explore conformational ensembles and transition states.19
- Use free energy methods (MM-GBSA, MM-PBSA) with caution; they may be useful for relative rankings but require wise interpretation with bRo5 compounds.20,21
Part III: Reject vs Best Practices Checklist
Immediate Rejection if Manuscripts Contain:
- Predictions of bioactivities, ADMET, pharmacokinetics, or toxicity without validation.
- 2D-QSAR studies or unvalidated QSAR models.
- Black-box computational servers or uncritical pipelines.
- Docking studies reporting raw scores as binding energies or lacking methodological transparency.
- Virtual screening recommending hits without experimental validation.
- AI/ML models without proper interpretability, benchmarking, or independent test sets.
- MD studies with single, short trajectories, poor setup, or overinterpretation.
- MD studies selectively highlighting certain frames or expected outcomes.
- MD studies including videos/graphics as decorative add-ons with no scientific value.
- MM-PBSA/GBSA based on single snapshots or reported as absolute values.
- QM/MM studies presenting meaningless QM-derived properties.
- FMO analysis presented as standalone explanations of activity or potency without mechanistic or experimental backing.
- MEPs used as decorative or speculative with no interpretation.
- FEP studies lacking convergence checks, error estimates, or experimental correlation.
Best Practices Required for Consideration:
- Use 3D-QSAR/advanced models with rigorous validation and interpretable insights.
- Employ molecular modeling strictly as hypothesis generation or rationalization, with accurate 3D depictions.
- For docking/VS: disclose preparation parameters, benchmark against known ligands, and validate hits experimentally.
- For benchmarking: provide robust metrics and clear recommendations.
- For AI/ML: justify algorithm choice, train on curated datasets, validate independently, ensure proper interpretability.
- For MD: use high-quality starting structures, correct protonation/tautomer states, sufficient timescales, and multiple replicas if possible. Use objective, reproducible analyses and provide videos that illustrate major findings, not cosmetic motions.
- For MM-PBSA/GBSA: apply only to equilibrated ensembles, clearly stating snapshot protocols.
- For QM/QM-MM: justify method choice, benchmark where possible, focus on biologically relevant properties.
- For FMO analysis: link to specific reactivity hypotheses with supporting experimental or mechanistic discussion.
- For HOMO/LUMO illustrations: they must serve an explanatory purpose.
- For MEPs: illustrate well-defined electrostatic features rationalizing site-specific interactions in a ligand–target complex.
- For FEP: apply rigorous sampling, convergence analysis, and experimental benchmarking.
Disclosure
Frank M. Boeckler is the Editor-in-Chief of Drug Design, Development and Therapy. Tamer M. Ibrahim and Muzammal Hussain are Associate Editors of Drug Design, Development and Therapy. The authors declare that they have no other competing interests in this work.
References
1. Cherkasov A, Muratov EN, Fourches D, et al. QSAR modeling: where have you been? Where are you going to? J Med Chem. 2014;57(12):4977–5010. doi:10.1021/jm4004285
2. Bauer MR, Ibrahim TM, Vogel SM, Boeckler FM. Evaluation and optimization of virtual screening workflows with DEKOIS 2.0 – a public library of challenging Docking Benchmark sets. J Chem Inf Model. 2013;53(6):1447–1462. doi:10.1021/ci400115b
3. Ibrahim TM, Bauer MR, Dörr A, Veyisoglu E, Boeckler FM. pROC-chemotype plots enhance the interpretability of benchmarking results in structure-based virtual screening. J Chem Inf Model. 2015;55(11):2297–2307. doi:10.1021/acs.jcim.5b00475
4. Ibrahim TM, Bauer MR, Boeckler FM. Probing the impact of protein and ligand preparation procedures on chemotype enrichment in structure-based virtual screening using DEKOIS 2.0 benchmark sets. J Cheminf. 2014;6(1):P19. doi:10.1186/1758-2946-6-S1-P19
5. Ibrahim TM, Bauer MR, Boeckler FM. Applying DEKOIS 2.0 in structure-based virtual screening to probe the impact of preparation procedures and score normalization. J Cheminf. 2015;7(1):21. doi:10.1186/s13321-015-0074-6
6. Warren GL, Andrews CW, Capelli A-M, et al. A critical assessment of docking programs and scoring functions. J Med Chem. 2006;49(20):5912–5931. doi:10.1021/jm050362n
7. Vamathevan J, Clark D, Czodrowski P, et al. Applications of machine learning in drug discovery and development. Nat Rev Drug Discov. 2019;18(6):463–477. doi:10.1038/s41573-019-0024-5
8. Hollingsworth SA, Dror RO. Molecular dynamics simulation for all. Neuron. 2018;99(6):1129–1143. doi:10.1016/j.neuron.2018.08.011
9. Homeyer N, Gohlke H. Free energy calculations by the molecular mechanics Poisson-Boltzmann surface area method. Mol Inform. 2012;31(2):114–122. doi:10.1002/minf.201100135
10. Vaas S, Zimmermann MO, Schollmeyer D, et al. Principles and applications of CF2X moieties as unconventional halogen bond donors in medicinal chemistry, chemical biology, and drug discovery. J Med Chem. 2023;66(15):10202–10225. doi:10.1021/acs.jmedchem.3c00634
11. Chodera JD, Mobley DL, Shirts MR, Dixon RW, Branson K, Pande VS. Alchemical free energy methods for drug discovery: progress and challenges. Curr Opin Struct Biol. 2011;21(2):150–160. doi:10.1016/j.sbi.2011.01.011
12. Price E, Weinheimer M, Rivkin A, et al. Beyond rule of five and PROTACs in modern drug discovery: polarity reducers, chameleonicity, and the evolving physicochemical landscape. J Med Chem. 2024;67(7):5683–5698. doi:10.1021/acs.jmedchem.3c02332
13. DeGoey DA, Chen H-J, Cox PB, Wendt MD. Beyond the rule of 5: lessons learned from AbbVie’s drugs and compound collection. J Med Chem. 2018;61(7):2636–2651. doi:10.1021/acs.jmedchem.7b00717
14. Poongavanam V, Doak BC, Kihlberg J. Opportunities and guidelines for discovery of orally absorbed drugs in beyond rule of 5 space. Curr Opin Chem Biol. 2018;44:23–29. doi:10.1016/j.cbpa.2018.05.010
15. Pereira GP, Gouzien C, Souza PCT, Martin J. Challenges in predicting PROTAC-mediated protein–protein interfaces with AlphaFold reveal a general limitation on small interfaces. Bioinform Adv. 2025;5(1):vbaf056. doi:10.1093/bioadv/vbaf056
16. Liao J, Nie X, Unarta IC, Ericksen SS, Tang W. In silico modeling and scoring of PROTAC-mediated ternary complex poses. J Med Chem. 2022;65(8):6116–6132. doi:10.1021/acs.jmedchem.1c02155
17. Allen SE, Dokholyan NV, Bowers AA. Dynamic docking of conformationally constrained macrocycles: methods and applications. ACS Chem Biol. 2016;11(1):10–24. doi:10.1021/acschembio.5b00663
18. Drummond ML, Henry A, Li H, Williams CI. Improved Accuracy for Modeling PROTAC-mediated ternary complex formation and targeted protein degradation via new in silico methodologies. J Chem Inf Model. 2020;60(10):5234–5254. doi:10.1021/acs.jcim.0c00897
19. AS BG, Agrawal D, Kulkarni NM, Vetrivel R, Gurram K. PROTAC-Design-Evaluator (PRODE): an advanced method for in-silico PROTAC design. ACS Omega. 2024;9(11):12611–12621. doi:10.1021/acsomega.3c07318
20. Roux B, Chipot C. Editorial guidelines for computational studies of ligand binding using MM/PBSA and MM/GBSA approximations wisely. J Phys Chem A. 2024;128(49):12027–12029. doi:10.1021/acs.jpcb.4c06614
21. Genheden S, Ryde U. The MM/PBSA and MM/GBSA methods to estimate ligand-binding affinities. Expert Opin Drug Discov. 2015;10(5):449–461. doi:10.1517/17460441.2015.1032936
© 2025 The Author(s). This work is published and licensed by Dove Medical Press Limited. The
full terms of this license are available at https://www.dovepress.com/terms
and incorporate the Creative Commons Attribution
- Non Commercial (unported, 4.0) License.
By accessing the work you hereby accept the Terms. Non-commercial uses of the work are permitted
without any further permission from Dove Medical Press Limited, provided the work is properly
attributed. For permission for commercial use of this work, please see paragraphs 4.2 and 5 of our Terms.
