Quantum annealing is a promising approach for obtaining good approximate solutions to difficult optimization problems.Folding a protein sequence into its minimum-energy structure represents such a problem.For testing new algorithms and technologies for this task, the minimal lattice-based [hydrophobic (H) or polar (P) beads] HP model is well suited, as it represents a considerable challenge despite its simplicity.The HP model has favorable interactions between adjacent, not directly bound hydrophobic residues.Here, we develop a novel spin representation for lattice protein folding tailored for quantum annealing.With a distributed encoding onto the lattice, it differs from earlier attempts to fold lattice proteins on quantum annealers, which were based upon chain growth techniques.With our encoding, the Hamiltonian by design has the quadratic structure required for calculations on an Ising-type annealer, without having to introduce any auxiliary spin variables.This property greatly facilitates the study of long chains.The approach is robust to changes in the parameters required to constrain the spin system to chainlike configurations, and performs very well in terms of solution quality.The results are evaluated against existing exact results for HP chains with up to N = 30 beads with 100% hit rate, thereby also outperforming classical simulated annealing.In addition, the method allows us to recover the lowest known energies for N = 48 and N = 64 HP chains, with similar hit rates.These results are obtained by the commonly used hybrid quantum-classical approach.For pure quantum annealing, our method successfully folds an N = 14 HP chain.The calculations were performed on a D-Wave Advantage quantum annealer.
Quantum annealing is a promising approach for obtaining good approximate solutions to difficult optimization problems.Folding a protein sequence into its minimum-energy structure represents such a problem.For testing new algorithms and technologies for this task, the minimal lattice-based [hydrophobic (H) or polar (P) beads] HP model is well suited, as it represents a considerable challenge despite its simplicity.The HP model has favorable interactions between adjacent, not directly bound hydrophobic residues.Here, we develop a novel spin representation for lattice protein folding tailored for quantum annealing.With a distributed encoding onto the lattice, it differs from earlier attempts to fold lattice proteins on quantum annealers, which were based upon chain growth techniques.With our encoding, the Hamiltonian by design has the quadratic structure required for calculations on an Ising-type annealer, without having to introduce any auxiliary spin variables.This property greatly facilitates the study of long chains.The approach is robust to changes in the parameters required to constrain the spin system to chainlike configurations, and performs very well in terms of solution quality.The results are evaluated against existing exact results for HP chains with up to N = 30 beads with 100% hit rate, thereby also outperforming classical simulated annealing.In addition, the method allows us to recover the lowest known energies for N = 48 and N = 64 HP chains, with similar hit rates.These results are obtained by the commonly used hybrid quantum-classical approach.For pure quantum annealing, our method successfully folds an N = 14 HP chain.The calculations were performed on a D-Wave Advantage quantum annealer.
Quantum annealing has shown promise for finding solutions to difficult optimization problems, including protein folding. Recently, we used the D-Wave Advantage quantum annealer to explore the folding problem in a coarse-grained lattice model, the HP model, in which amino acids are classified into two broad groups: hydrophobic (H) and polar (P). Using a set of 22 HP sequences with up to 64 amino acids, we demonstrated the fast and consistent identification of the correct HP model ground states using the D-Wave hybrid quantum-classical solver. An equally relevant biophysical challenge, called the protein design problem, is the inverse of the above, where the task is to predict protein sequences that fold to a given structure. Here, we approach the design problem by a two-step procedure implemented and executed on a D-Wave machine. In the first step, we perform a pure sequence-space search by varying the type of amino acid at each sequence position, and seek sequences which minimize the HP-model energy of the target structure. After mapping this task onto an Ising spin-glass representation, we employ a hybrid quantum-classical solver to deliver energy-optimal sequences for structures with 30–64 amino acids, with a 100% success rate. In the second step, we filter the optimized sequences from the first step according to their ability to fold to the intended structure. In addition, we try solving the sequence optimization problem using only the quantum processing unit (QPU), which confines us to sizes ≤20 , due to exponentially decreasing success rates. To shed light on the pure QPU results, we investigate the effects of control errors caused by an imperfect implementation of the intended Hamiltonian on the QPU, by numerically analyzing the Schrödinger equation. We find that the simulated success rates in the presence of control noise semiquantitatively reproduce the modest pure QPU results for larger chains. Published by the American Physical Society 2024
Using a simple lattice model with a 2 letter amino acid alphabet (H: hydrophobic and P: polar), we have explored the important biophysical problems of protein folding and protein design. Our formulation of these problems features a simple algebraic form for the Hamiltonian irrespective of the system size and composition. Since exact results are available from for lattice HP chains up to a size of 30, we were able to thoroughly validate our approach. The D-Wave advantage quantum annealer successfully identifies the ground state of the HP model protein chain in 100% of cases. For a few longer protein chains with up to 64 amino acids, where exact enumerations were not available but extensive Monte Carlo studies exist, the D-Wave hybrid annealer found the correct ground states within minutes, once again with a 100% success rate. Applying the same technique to the protein design problem, we found novel HP protein sequences with the same ground state as the largest systems we studied for folding. In contrast to the sequences used for folding, some of our newly found sequences had unique rather than degenerate ground states.
Protein folding is a computationally challenging problem because the number of possible structures increases rapidly with sequence length. This paper presents a hybrid quantum-classical framework for protein folding using the Hydrophobic-Polar (HP) model in a three-dimensional lattice. Valid protein structures are generated using a self-avoiding walk, and their energy is evaluated based on hydrophobic interactions. Simulated annealing is used as the classical opti- mization method, while the Quantum Approximate Optimization Algorithm (QAOA) is used for quantum optimization. Experimental results show that simulated annealing achieves lower energy values (e.g., -5.0 eV) compared to QAOA (-4.45 eV). Although quantum optimization does not yet outperform classical methods, it produces close approximations and shows potential for future improvements.
Lattice protein folding models are a cornerstone of computational biophysics. Although these models are a coarse grained representation, they provide useful insight into the energy landscape of natural proteins. Finding low-energy threedimensional structures is an intractable problem even in the simplest model, the Hydrophobic-Polar (HP) model. Description of protein-like properties are more accurately described by generalized models, such as the one proposed by Miyazawa and Jernigan (MJ), which explicitly take into account the unique interactions among all 20 amino acids. There is theoretical and experimental evidence of the advantage of solving classical optimization problems using quantum annealing over its classical analogue (simulated annealing). In this report, we present a benchmark implementation of quantum annealing for lattice protein folding problems (six different experiments up to 81 superconducting quantum bits). This first implementation of a biophysical problem paves the way towards studying optimization problems in biophysics and statistical mechanics using quantum devices.
Quantum computing has progressed from a theoretical concept to a practical experimental technology. By processing information according to quantum mechanics rather than classical rules, quantum computers may offer advantages for certain problems. Although a few algorithms have proven speedups, it remains unclear where quantum computers will be practically useful, especially since current devices are noisy, small, and not yet fault tolerant. This thesis explores quantum computing for biophysically motivated combinatorial problems, with a focus on lattice protein models. In these simplified models, a protein is represented as a chain on a discrete lattice. Despite their simplicity, lattice proteins retain important features of real proteins, such as such as the folding of a chain with exponentially many possible shapes to a unique minimum-energy structure and steric constraints.To run lattice protein problems on quantum hardware, they must first be formulated using binary variables and objective functions. The latter typically contain soft constraints needed to ensure physical solutions. In Papers~\I--\II, we developed and extended a coordinate-based mapping for lattice protein folding, enabling D-Wave's hybrid quantum-classical annealer to fold both 2D hydrophobic-polar and 3D Miyazawa-Jernigan lattice proteins. In Paper~\III, we adapted the approach to the gate-based paradigm and showed how quadratic soft constraints could be removed by a mixer developed for the maximum-independent-set problem. Papers~\IV--\V addressed lattice protein design using quantum annealing, quantum approximate optimization algorithm, quantum alternating operator ansatz and hardware-efficient ansatz, while Paper~\VI extended the framework to multi-chain thermodynamics and the reconstruction of densities of states and heat-capacity features.Together, these papers show how folding, design, and multi-chain thermodynamics can be formulated, implemented, and analyzed using present-day quantum and hybrid hardware. A recurring limitation is the treatment of soft constraints. In quantum annealing, their high connectivity leads to error accumulation. On gate-based hardware, they increase circuit depth and sensitivity to noise. At the same time, temporary constraint violations can provide useful freedom during optimization. Further progress will therefore require improved formulations, algorithms and hardware that reduce the cost of enforcing physical solutions while retaining this flexibility.
In computational biology, protein folding continues to be a key challenge with important implications for drug design and biological function comprehension. Using both classical and quantum computing methods, this paper provides an in-depth analysis of various possible lattice configurations and how theyare used in protein folding. To assess the protein design, we investigate classic models including the more intricate Miyazawa-Jernigan (MJ) interaction matrix and the hydrophobic-hydrophilic (HP) lattice model. Furthermore, we use integer programming methods to protein folding problems on lattices, addressing their NP-hardness and proving their effectiveness in obtaining optimal solutions. Expanding upon these conventional techniques, we present sophisticated graph neural network (GNN) models, particularly the RG-MPNN, which integrate pharmacophore ideas to predict chemical properties by using task-related representations. In the latter section of our paper, we examine how protein folding simulations on lattice structures could be transformed using quantum algorithms. Compared to classical methods, quantum computational techniques promise to be more accurate and efficient at handling the enormous diversity and complexity of protein folding.