Biomolecular Multiscale Simulation (BMS25) Dataset to Train Neural Network Potentials for QM/MM Settings with Electrostatic Embedding

Abstract

Neural network potentials (NNPs) can provide insight into biological processes at atomic resolution. Training these NNPs requires large and diverse datasets of molecules, conformations, and configurations. However, so far little attention has been paid to the description of solvation, despite its importance for biomolecular systems. This work lays the foundation for NNPs where solvation is an integral part of the model. Following a quantum-mechanics/molecular-mechanics (QM/MM) formalism with an electrostatic embedding scheme, systems are decomposed into a QM zone with the solute(s), which is electrostatically coupled to the point charges from surrounding solvent molecules (MM zone). Using an accelerated sampling approach, we generate the biomolecular multiscale simulation (BMS25) dataset with over 50,000 topologies and more than 1.5 million unique conformations of peptides and miniproteins as well as small molecules and transition states from chemical reactions. The dataset includes energies, gradients, and multipoles of solute molecules as well as gradients on solvent molecules at the ω B97M-D4/ma-def2-TZVPP level of theory, enabling the development of multiscale NNPs for simulating large biomolecular systems.