QUBO Problem Formulation of Fragment-Based Protein-Ligand Flexible Docking

• Published in: Entropy (Basel, Switzerland) Vol. 26; no. 5; p. 397
• Main Authors: Yanagisawa, Keisuke, Fujie, Takuya, Takabatake, Kazuki, Akiyama, Yutaka
• Format: Journal Article
• Language: English
• Published: Switzerland MDPI AG 01.05.2024
• Subjects: Aldose reductase; Algorithms; Binding; Chemical bonds; Chemical compounds; combinatorial optimization; compound fragment; Covalent bonds; Degrees of freedom; Docking; Drug discovery; Field programmable gate arrays; flexible docking; Fragments; Free energy; Ligands; Optimization; Proteins; protein–ligand docking; quantum annealing; R&D; Reductases; Research & development; simulated quantum annealer; Japan; combinatorial optimization; simulated quantum annealer; flexible docking; protein–ligand docking; quantum annealing; SQBM; compound fragment; quadratic unconstrained binary optimization (QUBO)
• ISSN: 1099-4300; 1099-4300
• DOI: 10.3390/e26050397

Protein-ligand docking plays a significant role in structure-based drug discovery. This methodology aims to estimate the binding mode and binding free energy between the drug-targeted protein and candidate chemical compounds, utilizing protein tertiary structure information. Reformulation of this docking as a quadratic unconstrained binary optimization (QUBO) problem to obtain solutions via quantum annealing has been attempted. However, previous studies did not consider the internal degrees of freedom of the compound that is mandatory and essential. In this study, we formulated fragment-based protein-ligand flexible docking, considering the internal degrees of freedom of the compound by focusing on fragments (rigid chemical substructures of compounds) as a QUBO problem. We introduced four factors essential for fragment-based docking in the Hamiltonian: (1) interaction energy between the target protein and each fragment, (2) clashes between fragments, (3) covalent bonds between fragments, and (4) the constraint that each fragment of the compound is selected for a single placement. We also implemented a proof-of-concept system and conducted redocking for the protein-compound complex structure of Aldose reductase (a drug target protein) using SQBM+, which is a simulated quantum annealer. The predicted binding pose reconstructed from the best solution was near-native (RMSD = 1.26 Å), which can be further improved (RMSD = 0.27 Å) using conventional energy minimization. The results indicate the validity of our QUBO problem formulation.