How to Use Quantum Computers for Biomolecular Free Energies

Free energy calculations are at the heart of physics-based analyses of biochemical processes. They allow us to quantify molecular recognition mechanisms, which determine a wide range of biological phenomena, from how cells send and receive signals to how pharmaceutical compounds can be used to treat diseases. Quantitative and predictive free energy calculations require computational models that accurately capture both the varied and intricate electronic interactions between molecules as well as the entropic contributions from the motions of these molecules and their aqueous environment. However, accurate quantum-mechanical energies and forces can be obtained only for small atomistic models and not for large biomacromolecules. Here, we demonstrate how to consistently link accurate quantum-mechanical data obtained for substructures to the overall potential energy of biomolecular complexes using machine learning in an integrated algorithm. We do so using a two-fold quantum embedding strategy where the innermost quantum cores are treated at a very high level of accuracy. We demonstrate the viability of this approach for the molecular recognition of a ruthenium-based anticancer drug by its protein target by applying traditional quantum chemical methods. As such methods scale unfavorably with system size, we analyze the requirements for quantum computers to provide highly accurate energies that affect the resulting free energies. Once the requirements are met, our computational pipeline, FreeQuantum, is able to make efficient use of the quantum-computed energies, thereby enabling quantum computing-enhanced modeling of biochemical processes. This approach combines the exponential speedups of quantum computers for simulating interacting electrons with modern classical simulation techniques that incorporate machine learning to model large molecules.