StreaMD: the toolkit for high-throughput molecular dynamics simulations

Pavel Polishchuk,Aleksandra Ivanova,Olena Mokshyna

DOI: https://doi.org/10.26434/chemrxiv-2024-2rjqz

2024-06-13

Abstract: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 module that streamlines all phases of molecular dynamics simulations, encompassing preparation, execution, and analysis. This module 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.

Chemistry

What problem does this paper attempt to address?

The paper primarily aims to address several key issues in molecular dynamics (MD) simulations and binding free energy calculations:

1. **Simplifying the MD simulation process**: Traditional MD simulation setup and execution require specialized knowledge and are prone to errors, especially when dealing with a large number of ligands and complexes. The paper proposes a Python module—the StreaMD toolkit, which aims to automate all stages of MD simulations, including preparation, execution, and analysis, thereby lowering the expertise threshold required for users.

2. **Supporting distributed computing environments**: Existing automation tools have limitations when simulating large compound sets, particularly in their ability to run on distributed server clusters. StreaMD leverages the Dask library to efficiently run high-throughput MD simulations on a single server or distributed systems without the need for a dedicated scheduler.

3. **Binding free energy calculations and interaction fingerprint analysis**: StreaMD not only automates MD trajectory simulations but also integrates binding free energy calculations (using MM-GBSA/PBSA methods) and interaction fingerprint analysis between protein-ligand complexes, making these computations more straightforward.

4. **Supporting simulations of special chemical components**: The tool specifically mentions support for systems containing special chemical components such as cofactors and boron atoms, a feature lacking in many existing tools, which helps improve the accuracy of simulation results.

In summary, StreaMD aims to simplify the MD simulation workflow by providing a fully automated, user-friendly tool that supports high-throughput computational needs and can handle complex chemical systems, thereby advancing research in structural biology and drug discovery.