PySTARC: GPU-accelerated Brownian dynamics for bimolecular association rate constants

Abstract

Drug-target association and dissociation rates often determine in vivo efficacy more than affinity alone. The association rate kon, however, is computationally challenging to predict since the productive encounter is a rare event in a large translational and orientational search space, which lies beyond the reach of conventional atomistic simulations. We present PySTARC (Python Simulation Toolkit for Association Rate Constants), a GPU-accelerated Brownian dynamics engine for estimating bimolecular association rate constants. PySTARC is a Python reimplementation of the BrownDye engine that converges even small reaction probabilities on a single GPU, at a throughput that would otherwise require a large CPU cluster. The current framework resolves reactions between integration steps with a closed-form Brownian bridge, models the internal flexibility of the solute through a coarse-grained bead chain, distributes trajectories across multiple GPUs, checkpoints long runs, monitors convergence, and automates the entire workflow from input structures to rate estimates. PySTARC is validated against protein-ligand and protein-protein complexes spanning five orders of magnitude in kon, with most estimates within one order of magnitude of experiment.