Protein-Ligand Free Energy Perturbation on Quantum Hardware
arXiv QuantumArchived Apr 14, 2026✓ Full text saved
arXiv:2604.09857v1 Announce Type: new Abstract: The use of free energy perturbation (FEP) methods to study protein-ligand complexes is one of the most important tools in structure-based drug design. Because FEP methods typically rely on force fields, they may suffer from force field parameter-related issues. Herein, we present a quantum mechanics/molecular mechanics (QM/MM) hybrid method to overcome deficiencies in force-field models by using QM bookending approaches on both classical and quantu
Full text archived locally
✦ AI Summary· Claude Sonnet
Quantum Physics
[Submitted on 10 Apr 2026]
Protein-Ligand Free Energy Perturbation on Quantum Hardware
Zhen Li, Milana Bazayeva, Thaddeus Pellegrini, Mario Motta, Subhamoy Bhowmik, Susanta Das, Danil Kaliakin, Fangchun Liang, Akhil Shajan, Kenneth M. Merz Jr
The use of free energy perturbation (FEP) methods to study protein-ligand complexes is one of the most important tools in structure-based drug design. Because FEP methods typically rely on force fields, they may suffer from force field parameter-related issues. Herein, we present a quantum mechanics/molecular mechanics (QM/MM) hybrid method to overcome deficiencies in force-field models by using QM bookending approaches on both classical and quantum hardware. In the MM part of this QM/MM FEP method, AMBER is used to simulate the protein receptor and the unperturbed moiety of the ligand, with the ff19SB and GAFF2 force fields. In the QM part, QUICK was used to conduct Hartree-Fock (HF) calculations, followed by heat-bath configuration interaction (HCI) as a benchmark on classical devices. To enable the HCI function in QUICK, we developed a Python-based interface to execute HCI from IBM's qiskit-addon-dice-solver. Moreover, the same interface also enabled this work to execute QM/MM FEP calculations on quantum hardware using the Local Unitary Cluster Jastrow (LUCJ) ansatz, followed by sample-based diagonalization (SQD) and extended-SQD (extSQD) post-processing. Using a series of thermolysis inhibitors as an example, we find reasonable agreement with experiment between the classical HCI method and the LUCJ-SQD/extSQD method, with the latter yielding a result closer to the experimental value. The execution time between the HCI-based FEP method and the LUCJ-SQD/extSQD-based FEP method is also comparable, indicating a high potential for utility in the noisy intermediate-scale quantum (NISQ) era.
Subjects: Quantum Physics (quant-ph)
Cite as: arXiv:2604.09857 [quant-ph]
(or arXiv:2604.09857v1 [quant-ph] for this version)
https://doi.org/10.48550/arXiv.2604.09857
Focus to learn more
Submission history
From: Akhil Shajan [view email]
[v1] Fri, 10 Apr 2026 19:37:19 UTC (3,564 KB)
Access Paper:
HTML (experimental)
view license
Current browse context:
quant-ph
< prev | next >
new | recent | 2026-04
References & Citations
INSPIRE HEP
NASA ADS
Google Scholar
Semantic Scholar
Export BibTeX Citation
Bookmark
Bibliographic Tools
Bibliographic and Citation Tools
Bibliographic Explorer Toggle
Bibliographic Explorer (What is the Explorer?)
Connected Papers Toggle
Connected Papers (What is Connected Papers?)
Litmaps Toggle
Litmaps (What is Litmaps?)
scite.ai Toggle
scite Smart Citations (What are Smart Citations?)
Code, Data, Media
Demos
Related Papers
About arXivLabs
Which authors of this paper are endorsers? | Disable MathJax (What is MathJax?)