Reuse & Permissions

It is not necessary to obtain permission to reuse this article or its components as it is available under the terms of the Creative Commons Attribution 4.0 International license. This license permits unrestricted use, distribution, and reproduction in any medium, provided attribution to the author(s) and the published article's title, journal citation, and DOI are maintained. Please note that some figures may have been included with permission from other third parties. It is your responsibility to obtain the proper permission from the rights holder directly for these figures.

Export citation

Export citation

Choose format for download:

Download Citation
  • Open Access

Practical protein-pocket hydration-site prediction for drug discovery on a quantum computer

Daniele Loco*

Kisa Barkemeyer* and Andre R. R. Carvalho

Jean-Philip Piquemal†

  • Q-CTRL, Berlin, Germany and Q-CTRL, Sydney, New South Wales, Australia

  • *These authors contributed equally to this work.
  • †Contact author: jean-philip.piquemal@sorbonne-universite.fr

Phys. Rev. Research 8, 033344 – Published 21 September, 2026

DOI: https://doi.org/10.1103/gyqr-mv1h

Abstract

Demonstrating the practical utility of noisy intermediate-scale quantum hardware for recurrent tasks in computer-aided drug discovery is of paramount importance. We tackle this challenge by performing three-dimensional hydration-site prediction in protein pockets on a quantum computer. Formulating the water placement problem as a quadratic unconstrained binary optimization, we use a hybrid approach coupling a classical three-dimensional reference-interaction site model to an efficient quantum optimization solver to run various hardware experiments up to 123 qubits. Matching the precision of classical approaches, our results reproduced experimental predictions on real-life protein-ligand complexes. Furthermore, through a detailed resource estimation analysis, we show that accuracy can be systematically improved with increasing number of qubits, indicating that full quantum utility is within reach as devices scale in the near term. The method has potential for assisting simulations of protein-ligand complexes for drug lead optimization and setup of docking calculations.

View figure in article

Physics Subject Headings (PhySH)

Article Text

Supplemental Material

References (80)

  1. A. B. Gurung, M. A. Ali, J. Lee, M. A. Farah, and K. M. Al-Anazi, An updated review of computer-aided drug design and its application to COVID-19, Biomed Res. Int. 2021, 8853056 (2021).
  2. D. Vemula, P. Jayasurya, V. Sushmitha, Y. N. Kumar, and V. Bhandari, CADD, AI and ML in drug discovery: A comprehensive review, Eur. J. Pharm. Sci. 181, 106324 (2023).
  3. Y. Chang, B. A. Hawkins, J. J. Du, P. W. Groundwater, D. E. Hibbs, and F. Lai, A guide to in silico drug design, Pharmaceutics 15, 49 (2023).
  4. D. Bucher, P. Stouten, and N. Triballeau, Shedding light on important waters for drug design: Simulations versus grid-based methods, J. Chem. Inf. Model. 58, 692 (2018).
  5. M. Blazhynska, N. Ansari, L. Lagardère, and J.-P. Piquemal, From water networks to binding affinities: Resolving solvation dynamics for accurate protein-ligand predictions, bioRxiv (2026), doi:10.1101/2025.10.17.683050.
  6. G. A. Ross, G. M. Morris, and P. C. Biggin, Rapid and accurate prediction and scoring of water molecules in protein binding sites, PLoS One 7, e32036 (2012).
  7. Y. Lu, R. Wang, C.-Y. Yang, and S. Wang, Analysis of ligand-bound water molecules in high-resolution crystal structures of protein–ligand complexes, J. Chem. Inf. Model. 47, 668 (2007).
  8. J. Michel, J. Tirado-Rives, and W. L. Jorgensen, Prediction of the water content in protein binding sites, J. Phys. Chem. B 113, 13337 (2009).
  9. S. M. Ivanov, Calculated hydration free energies become less accurate with increases in molecular weight, PLoS One 19, e0309996 (2024).
  10. M. L. Samways, R. D. Taylor, H. E. B. Macdonald, and J. W. Essex, Water molecules at protein–drug interfaces: Computational prediction and analysis methods, Chem. Soc. Rev. 50, 9104 (2021).
  11. L. Lagardère, L. Maurin, O. Adjoua, K. El Hage, P. Monmarché, J.-P. Piquemal, and J. Hénin, Lambda-ABF: Simplified, portable, accurate, and cost-effective alchemical free-energy computation, J. Chem. Theory Comput. 20, 4481 (2024).
  12. N. Ansari, Z. F. Jing, A. Gagelin, F. Hédin, F. Aviat, J. Hénin, J.-P. Piquemal, and L. Lagardère, Lambda-ABF-OPES: Faster convergence with high accuracy in alchemical free energy calculations, J. Phys. Chem. Lett. 16, 4626 (2025).
  13. W. L. Jorgensen, Efficient drug lead discovery and optimization, Acc. Chem. Res. 42, 724 (2009).
  14. S. A. Adcock and J. A. McCammon, Molecular dynamics: Survey of methods for simulating the activity of proteins, Chem. Rev. 106, 1589 (2006).
  15. J. Hénin, T. Lelièvre, M. R. Shirts, O. Valsson, and L. Delemotte, Enhanced sampling methods for molecular dynamics simulations [article v1.0], Living J. Comput. Mol. Sci. 4, 1583 (2022).
  16. M. S. Bodnarchuk, R. Viner, J. Michel, and J. W. Essex, Strategies to calculate water binding free energies in protein–ligand complexes, J. Chem. Inf. Model. 54, 1623 (2014).
  17. S. de Beer, N. P. Vermeulen, and C. Oostenbrink, The role of water molecules in computational drug design, CTMC 10, 55 (2010).
  18. J. Setiadi, F. Biedermann, W. M. Nau, and M. K. Gilson, Thermodynamics of water displacement from binding sites and its contributions to supramolecular and biomolecular affinity, Angew. Chem. Int. Ed. 64, e202505713 (2025).
  19. Y. Ge, O. J. Melling, W. Dong, J. W. Essex, and D. L. Mobley, Enhancing sampling of water rehydration upon ligand binding using variants of grand canonical Monte Carlo, J. Comput.-Aided Mol. Des. 36, 767 (2022).
  20. T. J. Inizan, F. Célerse, O. Adjoua, D. El Ahdab, L.-H. Jolly, C. Liu, P. Ren, M. Montes, N. Lagarde, L. Lagardère, P. Monmarché, and J.-P. Piquemal, High-resolution mining of the SARS-CoV-2 main protease conformational space: Supercomputer-driven unsupervised adaptive sampling, Chem. Sci. 12, 4889 (2021).
  21. D. El Ahdab, L. Lagardère, T. J. Inizan, F. Célerse, C. Liu, O. Adjoua, L.-H. Jolly, N. Gresh, Z. Hobaika, P. Ren, R. G. Maroun, and J.-P. Piquemal, Interfacial water many-body effects drive structural dynamics and allosteric interactions in SARS-CoV-2 main protease dimerization interface, J. Phys. Chem. Lett. 12, 6218 (2021).
  22. A. Schahl, L. Lagardère, B. Walker, P. Ren, H. Wioland, M. Ballet, A. Jégou, M. Chavent, and J.-P. Piquemal, Histidine 73 methylation coordinates β-actin plasticity in response to key environmental factors, Nat. Commun. 16, 2304 (2025).
  23. L. El Khoury, Z. Jing, A. Cuzzolin, A. Deplano, D. Loco, B. Sattarov, F. Hédin, S. Wendeborn, C. Ho, D. El Ahdab, et al., Computationally driven discovery of SARS-CoV-2 Mpro inhibitors: From design to experimental validation, Chem. Sci. 13, 3674 (2022).
  24. N. Ansari, C. Liu, F. Hédin, J. Hénin, J. W. Ponder, P. Ren, J.-P. Piquemal, L. Lagardère, and K. El Hage, Targeting RNA with small molecules using state-of-the-art methods provides highly predictive affinities of riboswitch inhibitors, Commun. Biol. 8, 1405 (2025).
  25. L. Lagardère, L.-H. Jolly, F. Lipparini, F. Aviat, B. Stamm, Z. F. Jing, M. Harger, H. Torabifard, G. A. Cisneros, M. J. Schnieders, N. Gresh, Y. Maday, P. Y. Ren, J. W. Ponder, and J.-P. Piquemal, Tinker-HP: A massively parallel molecular dynamics package for multiscale simulations of large complex systems with advanced point dipole polarizable force fields, Chem. Sci. 9, 956 (2018).
  26. O. Adjoua, L. Lagardère, L.-H. Jolly, A. Durocher, T. Very, I. Dupays, Z. Wang, T. J. Inizan, F. Célerse, P. Ren, J. W. Ponder, and J.-P. Piquemal, Tinker-HP: Accelerating molecular dynamics simulations of large complex systems with advanced point dipole polarizable force fields using GPUs and multi-GPU systems, J. Chem. Theory Comput. 17, 2034 (2021).
  27. B. Kohnke, C. Kutzner, and H. Grubmüller, A GPU-accelerated fast multipole method for GROMACS: Performance and accuracy, J. Chem. Theory Comput. 16, 6938 (2020).
  28. J. C. Phillips, D. J. Hardy, J. D. Maia, J. E. Stone, J. V. Ribeiro, R. C. Bernardi, R. Buch, G. Fiorin, J. Hénin, W. Jiang, et al., Scalable molecular dynamics on CPU and GPU architectures with NAMD, J. Chem. Phys. 153, 044130 (2020).
  29. D. A. Case, D. S. Cerutti, V. W. D. Cruzeiro, T. A. Darden, R. E. Duke, M. Ghazimirsaeed, G. M. Giambaşu, T. J. Giese, A. W. Götz, J. A. Harris, et al., Recent developments in Amber biomolecular simulations, J. Chem. Inf. Model. 65, 7835 (2025).
  30. F. Bigi, S. Chong, A. Kristiadi, and M. Ceriotti, FlashMD: Long-stride, universal prediction of molecular dynamics, arXiv:2505.19350.
  31. M. Plainer, H. Wu, L. Klein, S. Günnemann, and F. Noé, Consistent sampling and simulation: Molecular dynamics with energy-based diffusion models, arXiv:2506.17139.
  32. F. Noé, A. Tkatchenko, K.-R. Müller, and C. Clementi, Machine learning for molecular simulation, Annu. Rev. Phys. Chem. 71, 361 (2020).
  33. A. Agathangelou, D. Manawadu, and I. Tavernelli, Quantum algorithm for protein side-chain optimisation: Comparing quantum to classical methods, arXiv:2507.19383.
  34. B. Claudon, J.-P. Piquemal, and P. Monmarché, Quantum speedup for nonreversible Markov chains, Nat. Commun. 16, 10732 (2025).
  35. B. Claudon, P. Rodenas-Ruiz, J.-P. Piquemal, and P. Monmarché, Quantum circuits for the Metropolis-Hastings algorithm, J. Phys. A: Math. Theor. 59, 305304 (2026).
  36. A. Robert, P. K. Barkoutsos, S. Woerner, and I. Tavernelli, Resource-efficient quantum algorithm for protein folding, npj Quantum Inf. 7, 38 (2021).
  37. I. Liliopoulos, G. D. Varsamis, T. Karamanidou, C. Papalitsas, G. Koulouras, V. Pantazopoulos, T. G. Stavropoulos, and I. G. Karafyllidis, Quantum algorithm for protein-ligand docking sites identification in the interaction space, J. Comput.-Aided Mol. Des. 39, 40 (2025).
  38. Q.-M. Ding, Y.-M. Huang, and X. Yuan, Molecular docking via quantum approximate optimization algorithm, Phys. Rev. Appl. 21, 034036 (2024).
  39. W. Li, Z. Yin, X. Li, D. Ma, S. Yi, Z. Zhang, C. Zou, K. Bu, M. Dai, J. Yue, Y. Chen, X. Zhang, and S. Zhang, A hybrid quantum computing pipeline for real world drug discovery, Sci. Rep. 14, 16942 (2024).
  40. P. Das, A. Ray, S. Bhattacharyya, J. Platos, V. Snasel, L. Mrsic, T. Huang, and I. Zelinka, A brief review on quantum computing based drug design, WIREs Data Min. Knowl. Discov. 14, e1553 (2024).
  41. M. D'Arcangelo, L.-P. Henry, L. Henriet, D. Loco, N. Gouraud, S. Angebault, J. Sueiro, J. Forêt, P. Monmarché, and J.-P. Piquemal, Leveraging analog quantum computing with neutral atoms for solvent configuration prediction in drug discovery, Phys. Rev. Res. 6, 043020 (2024).
  42. P. S. Mundada, A. Barbosa, S. Maity, Y. Wang, T. Merkh, T. Stace, F. Nielson, A. R. R. Carvalho, M. Hush, M. J. Biercuk, and Y. Baum, Experimental benchmarking of an automated deterministic error-suppression workflow for quantum algorithms, Phys. Rev. Appl. 20, 024034 (2023).
  43. N. Sachdeva, G. S. Hartnett, S. Maity, S. Marsh, Y. Wang, A. Winick, R. Dougherty, D. Canuto, Y. Q. Chong, M. Hush, P. S. Mundada, C. D. B. Bentley, M. J. Biercuk, and Y. Baum, Quantum optimization using a 127-qubit gate-model IBM Quantum computer can outperform quantum annealers for nontrivial binary optimization problems, New J. Phys. 28, 064506 (2026).
  44. M. L. Samways, H. E. B. Macdonald, R. D. Taylor, and J. W. Essex, Water networks in complexes between proteins and FDA-approved drugs, J. Chem. Inf. Model. 63, 387 (2023).
  45. https://www.rcsb.org/.
  46. See Supplemental Material at http://link.aps.org/supplemental/10.1103/gyqr-mv1h for additional details on parameter choices in Sec. S2, cost gaps, metrics, and method descriptions, which include Refs. [10, 49, 80].
  47. https://github.com/dwavesystems/dwave-neal.
  48. IBM Quantum computing|Technology and roadmap, https://www.ibm.com/quantum/technology.
  49. E. Nittinger, F. Flachsenberg, S. Bietz, G. Lange, R. Klein, and M. Rarey, Placement of water molecules in protein structures: From large-scale evaluations to single-case examples, J. Chem. Inf. Model. 58, 1625 (2018).
  50. D. Roy and A. Kovalenko, Performance of 3D-RISM-KH in predicting hydration free energy: Effect of solute parameters, J. Phys. Chem. A 123, 4087 (2019).
  51. A. Zamanos, G. Ioannakis, and I. Z. Emiris, Hydraprot: A new deep learning tool for fast and accurate prediction of water molecule positions for protein structures, J. Chem. Inf. Model. 64, 2594 (2024).
  52. D. J. Sindhikara, N. Yoshida, and F. Hirata, Placevent: An algorithm for prediction of explicit solvent atom distribution-application to HIV-1 protease and F-ATP synthase, J. Comput. Chem. 33, 1536 (2012).
  53. H.-H. Bui, A. J. Schiewe, and I. S. Haworth, Watgen: An algorithm for modeling water networks at protein–protein interfaces, J. Comput. Chem. 28, 2241 (2007).
  54. A. Morozenko and A. A. Stuchebrukhov, Dowser++, a new method of hydrating protein structures, Proteins: Struct. Funct. Bioinf. 84, 1347 (2016).
  55. Avoiding an unexpected roadblock in quantum computing—Compilation|Q-CTRL, https://q-ctrl.com/topics/avoiding-an-unexpected-roadblock-in-quantum-computing-compilation.
  56. D. S. Johnson, C. R. Aragon, L. A. McGeoch, and C. Schevon, Optimization by simulated annealing: An experimental evaluation; Part I, graph partitioning, Oper. Res. 37, 865 (1989.
  57. D. S. Johnson, C. R. Aragon, L. A. McGeoch, and C. Schevon, Optimization by simulated annealing: An experimental evaluation; Part II, graph coloring and number partitioning, Oper. Res. 39, 378 (1991.
  58. D. Henderson, S. H. Jacobson, and A. W. Johnson, The theory and practice of simulated annealing, in Handbook of Metaheuristics, edited by F. Glover and G. A. Kochenberger (Springer US, Boston, MA, 2003), pp. 287–319.
  59. M. Streif and M. Leib, Comparison of QAOA with quantum and simulated annealing, arXiv:1901.01903.
  60. L. Zhou, S.-T. Wang, S. Choi, H. Pichler, and M. D. Lukin, Quantum approximate optimization algorithm: Performance, mechanism, and implementation on near-term devices, Phys. Rev. X 10, 021067 (2020).
  61. T. Luchko, S. Gusarov, D. R. Roe, C. Simmerling, D. A. Case, J. Tuszynski, and A. Kovalenko, Three-dimensional molecular theory of solvation coupled with molecular dynamics in Amber, J. Chem. Theory Comput. 6, 607 (2010).
  62. D. A. Case, H. M. Aktulga, K. Belfon, D. S. Cerutti, G. A. Cisneros, V. W. D. Cruzeiro, N. Forouzesh, T. J. Giese, A. W. Götz, H. Gohlke, et al., AmberTools, J. Chem. Inf. Model. 63, 6183 (2023).
  63. J. A. Maier, C. Martinez, K. Kasavajhala, L. Wickstrom, K. E. Hauser, and C. Simmerling, ff14SB: Improving the accuracy of protein side chain and backbone parameters from ff99SB, J. Chem. Theory Comput. 11, 3696 (2015).
  64. H. J. C. Berendsen, J. R. Grigera, and T. P. Straatsma, The missing term in effective pair potentials, J. Phys. Chem. 91, 6269 (1987).
  65. https://ambermd.org/doc12/Amber24.pdf.
  66. K. Osaki, T. Ekimoto, T. Yamane, and M. Ikeguchi, 3D-RISM-AI: A machine learning approach to predict protein–ligand binding affinity using 3D-RISM, J. Phys. Chem. B 126, 6148 (2022).
  67. D. J. Fowles and D. S. Palmer, Solvation entropy, enthalpy and free energy prediction using a multi-task deep learning functional in 1D-RISM, Phys. Chem. Chem. Phys. 25, 6944 (2023).
  68. L. Fusani, I. Wall, D. Palmer, and A. Cortes, Optimal water networks in protein cavities with GAsol and 3D-RISM, Bioinformatics 34, 1947 (2018).
  69. E. Farhi, J. Goldstone, and S. Gutmann, A quantum approximate optimization algorithm, arXiv:1411.4028.
  70. R. Shaydulin, C. Li, S. Chakrabarti, M. DeCross, D. Herman, N. Kumar, J. Larson, D. Lykov, P. Minssen, Y. Sun, et al., Evidence of scaling advantage for the quantum approximate optimization algorithm on a classically intractable problem, Sci. Adv. 10, eadm6761 (2024).
  71. S. Boulebnane and A. Montanaro, Solving Boolean satisfiability problems with the quantum approximate optimization algorithm, PRX Quantum 5, 030348 (2024).
  72. S. Omanakuttan, Z. He, Z. Zhang, T. Hao, A. Babakhani, S. Boulebnane, S. Chakrabarti, D. Herman, J. Sullivan, M. A. Perlin, R. Shaydulin, and M. Pistoia, Threshold for fault-tolerant quantum advantage with the quantum approximate optimization algorithm, arXiv:2504.01897.
  73. S. Wang, E. Fontana, M. Cerezo, K. Sharma, A. Sone, L. Cincio, and P. J. Coles, Noise-induced barren plateaus in variational quantum algorithms, Nat. Commun. 12, 6961 (2021).
  74. J. Weidenfeller, L. C. Valor, J. Gacon, C. Tornow, L. Bello, S. Woerner, and D. J. Egger, Scaling of the quantum approximate optimization algorithm on superconducting qubit based hardware, Quantum 6, 870 (2022).
  75. E. Pelofske, A. Bärtschi, and S. Eidenbenz, Short-depth QAOA circuits and quantum annealing on higher-order Ising models, npj Quantum Inf. 10, 30 (2024).
  76. Accelerating the schedule for quantum-enhanced rail|Q-CTRL, Section: Case Studies, 2025, https://q-ctrl.com/case-study/accelerating-the-schedule-for-quantum-enhanced-rail.
  77. Predicting protein-pocket hydration sites using Fire Opal|Q-CTRL, accessed 13 April 2026, https://docs.q-ctrl.com/fire-opal/apply/predicting-protein-pocket-hydration-sites-using-fire-opal.
  78. QUbit Pharmaceuticals, Protein pocket hydration QUBO data. In Practical protein-pocket hydration-site prediction for drug discovery on a quantum computer, Zenodo, 2026, https://doi.org/10.5281/zenodo.19912753.
  79. Both Q-CTRL and IBM offer limited free access to their platforms.
  80. J. Eberhardt and S. Forli, WaterKit: Thermodynamic profiling of protein hydration sites, J. Chem. Theory Comput. 19, 2535 (2023).

Outline

Information

Sign In to Your Journals Account

Filter

Filter

Article Lookup

Enter a citation