Molecular dynamics simulations serve as a prevalent approach for investigating the dynamic behaviour of proteins and protein-ligand complexes. Due to its versatility and speed, GROMACS stands out as a commonly utilized software platform for executing molecular dynamics simulations. However, its effective utilization requires substantial expertise in configuring, executing, and interpreting molecular dynamics trajectories. Existing automation tools are constrained in their capability to conduct simulations for large sets of compounds with minimal user intervention, or in their ability to distribute simulations across multiple servers. To address these challenges, we developed a Python-based tool that streamlines all phases of molecular dynamics simulations, encompassing preparation, execution, and analysis. This tool minimizes the required knowledge for users engaging in molecular dynamics simulations and can efficiently operate across multiple servers within a network or a cluster. Notably, the tool not only automates trajectory simulation but also facilitates the computation of free binding energies for protein-ligand complexes and generates interaction fingerprints across the trajectory. Our study demonstrated the applicability of this tool on several benchmark datasets. Additionally, we provided recommendations for end-users to effectively utilize the tool.Scientific contributionThe developed tool, StreaMD, is applicable to different systems (proteins, ligands and their complexes including co-factors) and requires a little user knowledge to setup and run molecular dynamics simulations. Other features of StreaMD are seamless integration with calculation of MM-GBSA/PBSA binding free energies and protein-ligand interaction fingerprints, and running of simulations within distributed environments. All these will facilitate routine and massive molecular dynamics simulations.

Institute of Molecular and Translational Medicine Faculty of Medicine and Dentistry Palacky University Hnevotinska 5 77900 Olomouc Czech Republic

Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences Flemingovo náměstí 542 2 160 00 Praha 6 Czech Republic

J Cheminform. 2024 Dec 23;16(1):144. doi: 10.1186/s13321-024-00942-w. PubMed

See more in PubMed

King E, Aitchison E, Li H, Luo R (2021) Recent developments in free energy calculations for drug discovery. Front Molecular Biosci. 10.3389/fmolb.2021.712085 PubMed PMC

De Vivo M, Masetti M, Bottegoni G, Cavalli A (2016) Role of molecular dynamics and related methods in drug discovery. J Med Chem 59(9):4035–4061 PubMed

Limongelli V (2020) Ligand binding free energy and kinetics calculation in 2020. WIREs Comput Mol Sci 10(4):e1455

Liao J, Nie X, Unarta IC, Ericksen SS, Tang W (2022) In silico modeling and scoring of PROTAC-mediated ternary complex poses. J Med Chem 65(8):6116–6132 PubMed

Åqvist J, Medina C, Samuelsson J-E (1994) A new method for predicting binding affinity in computer-aided drug design. Protein Eng Des Sel 7(3):385–391 PubMed

Wang E, Sun H, Wang J, Wang Z, Liu H, Zhang JZH, Hou T (2019) End-point binding free energy calculation with MM/PBSA and MM/GBSA: strategies and applications in drug design. Chem Rev 119(16):9478–9508 PubMed

Chai X, Sun H, Zhou W, Chen C, Shan L, Yang Y, He J, Pang J, Yang L, Wang X et al (2022) Discovery of N-(4-(Benzyloxy)-phenyl)-sulfonamide derivatives as novel antagonists of the human androgen receptor targeting the activation function 2. J Med Chem 65(3):2507–2521 PubMed

Eastman P, Swails J, Chodera JD, McGibbon RT, Zhao Y, Beauchamp KA, Wang L-P, Simmonett AC, Harrigan MP, Stern CD et al (2017) OpenMM 7: rapid development of high performance algorithms for molecular dynamics. PLoS Comput Biol 13(7):e1005659 PubMed PMC

Jo S, Kim T, Iyer VG, Im W (2008) CHARMM-GUI: A web-based graphical user interface for CHARMM. J Comput Chem 29(11):1859–1865 PubMed

Doerr S, Harvey MJ, Noé F, De Fabritiis G (2016) HTMD: high-throughput molecular dynamics for molecular discovery. J Chem Theory Comput 12(4):1845–1852 PubMed

Harvey MJ, Giupponi G, Fabritiis GD (2009) ACEMD: accelerating biomolecular dynamics in the microsecond time scale. J Chem Theory Comput 5(6):1632–1639 PubMed

Bray SA, Senapathi T, Barnett CB, Grüning BA (2020) Intuitive, reproducible high-throughput molecular dynamics in Galaxy: a tutorial. J Cheminf 12(1):54 PubMed PMC

Yang M, Bo Z, Xu T, Xu B, Wang D, Zheng H (2023) Uni-GBSA: an open-source and web-based automatic workflow to perform MM/GB(PB)SA calculations for virtual screening. Brief Bioinform 24:4 PubMed

Barreto Gomes DE, Galentino K, Sisquellas M, Monari L, Bouysset C, Cecchini M (2023) ChemFlow─From 2D chemical libraries to protein-ligand binding free energies. J Chem Inf Model 63(2):407–411 PubMed PMC

Minibaeva G, Ivanova A, Polishchuk P (2023) EasyDock: customizable and scalable docking tool. J Cheminf 15(1):102 PubMed PMC

Lindorff-Larsen K, Piana S, Palmo Kim, Maragakis Paul, Klepeis John L, Dror Ron O, Shaw David E (2010) Improved side-chain torsion potentials for the amber ff99SB protein force field. Proteins 78(8):1950–1958 PubMed PMC

Li P, Merz KM (2016) MCPB.py: a python based metal center parameter builder. J Chem Inf Model 56(4):599–604 PubMed

Valdés-Tresanco MS, Valdés-Tresanco ME, Valiente PA, Moreno E (2021) gmx_MMPBSA: a new tool to perform end-state free energy calculations with GROMACS. J Chem Theory Comput 17(10):6281–6291 PubMed

Duan L, Liu X, Zhang JZH (2016) Interaction entropy: a new paradigm for highly efficient and reliable computation of protein-ligand binding free energy. J Am Chem Soc 138(17):5722–5728 PubMed

Kohut G, Liwo A, Bősze S, Beke-Somfai T, Samsonov SA (2018) Protein-ligand interaction energy-based entropy calculations: fundamental challenges for flexible systems. J Phys Chem B 122(32):7821–7827 PubMed

Ekberg V, Ryde U (2021) On the use of interaction entropy and related methods to estimate binding entropies. J Chem Theory Comput 17(8):5379–5391 PubMed PMC

Bouysset C, Fiorucci S (2021) ProLIF: a library to encode molecular interactions as fingerprints. J Cheminf 13(1):72 PubMed PMC

Rehulka J, Subtelna I, Kryshchyshyn-Dylevych A, Cherniienko A, Ivanova A, Matveieva M, Polishchuk P, Gurska S, Hajduch M, Zagrijtschuk O et al (2022) Anticancer 5-arylidene-2-(4-hydroxyphenyl)aminothiazol-4(5H)-ones as tubulin inhibitors. Archiv der Pharmazie. 10.1002/ardp.202200419 PubMed

Jurášek M, Řehulka J, Hrubá L, Ivanová A, Gurská S, Mokshyna O, Trousil P, Huml L, Polishchuk P, Hajdúch M et al (2023) Triazole-based estradiol dimers prepared via CuAAC from 17α-ethinyl estradiol with five-atom linkers causing G2/M arrest and tubulin inhibition. Bioorg Chem 131:106334 PubMed

Řehulka J, Jurášek M, Dráber P, Ivanová A, Gurská S, Ječmeňová K, Mokshyna O, Hajdúch M, Polishchuk P, Drašar PB et al (2024) Click estradiol dimers with novel aromatic bridging units: synthesis and anticancer evaluation. J Enzyme Inhib Med Chem 39(1):2367139 PubMed PMC

Bahia MS, Kaspi O, Touitou M, Binayev I, Dhail S, Spiegel J, Khazanov N, Yosipof A, Senderowitz H (2023) A comparison between 2D and 3D descriptors in QSAR modeling based on bio-active conformations. Mol Inf 42(4):2200186 PubMed

Trott O, Olson AJ (2009) AutoDock Vina: Improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. J Comput Chem 31(2):455–461 PubMed PMC

McNutt AT, Francoeur P, Aggarwal R, Masuda T, Meli R, Ragoza M, Sunseri J, Koes DR (2021) GNINA 10: molecular docking with deep learning. J Cheminf 13(1):43 PubMed PMC

McGrath ME, Sprengeler PA, Hirschbein B, Somoza JR, Lehoux I, Janc JW, Gjerstad E, Graupe M, Estiarte A, Venkataramani C et al (2006) Structure-guided design of peptide-based tryptase inhibitors. Biochemistry 45(19):5964–5973 PubMed

Correction: StreaMD: the toolkit for high-throughput molecular dynamics simulations