A Scalable Heuristic for Molecular Docking on Neutral-Atom Quantum Processors

Molecular docking is a critical computational method in drug discovery used to predict the binding conformation and orientation of a ligand within a protein's binding site. Mapping this challenge onto a graph-based problem, specifically the Maximum Weighted Independent Set (MWIS) problem, allows it to be addressed by specialized hardware such as neutral-atom quantum processors. However, a significant bottleneck has been the size mismatch between biologically relevant molecular systems and the limited capacity of near-term quantum devices. In this work, we overcome this scaling limitation by the use of a novel divide-and-conquer heuristic introduced in Cazals et al. This algorithm enables the solution of large-scale MWIS problems by decomposing a single, intractable graph instance into smaller sub-problems that can be solved sequentially on a neutral-atom quantum emulator, incurring only a linear computational overhead. We demonstrate the power of this approach by solving a 540-node MWIS problem representing the docking of an inhibitor to the Tumor necrosis factor-$α$ Converting Enzyme--thiol-containing Aryl Sulfonamide (TACE-AS) complex. Our work enables the application of quantum methods to more complex and physically realistic molecular systems than previously possible, paving the way for tackling large-scale docking challenges on near-term quantum hardware.