Atomistic Lattice Structure Optimization with Quantum Annealing

Xiaojun Wang, Zhenghuan Wang, Bowen Ni

Mark Fingerhuth, Tomáš Babej, Christopher E. Ing

Gate-based universal quantum computers form a rapidly evolving field of quantum computing hardware technology. In previous work, we presented a quantum algorithm for lattice protein folding on a cubic lattice, tailored for quantum annealers. In this paper, we introduce a novel approach for solving the lattice protein folding problem on universal gate-based quantum computing architectures. Lattice protein models are coarse-grained representations of proteins that have been used extensively over the past thirty years to examine the principles of protein folding and design.These models can be used to explore a vast number of possible protein conformations and to infer structural properties of more complex atomistic protein structures. We formulate the problem as a quantum alternating operator ansatz, a member of the wider class of variational quantum/classical hybrid algorithms. To increase the probability of sampling the ground state, we propose splitting the optimization problem into hard and soft constraints. This enables us to use a previously under-utilised component of the variational algorithm to constrain the search to the subspace of solutions that satisfy the hard constraints.

E. Arvidsson

Tomáš Babej, Christopher E. Ing, Mark Fingerhuth

Lattice models have been used extensively over the past thirty years to examine the principles of protein folding and design. These models can be used to determine the conformation of the lowest energy fold out of a large number of possible conformations. However, due to the size of the conformational space, new algorithms are required for folding longer proteins sequences. Preliminary work was performed by Babbush et al. (2012) to fold a small peptide on a planar lattice using a quantum annealing device. We extend this work by providing improved Ising-type Hamiltonian encodings for the problem of finding the lowest energy conformation of a lattice protein. We demonstrate a decrease in quantum circuit complexity from quadratic to quasilinear in certain cases. Additionally, we generalize to three spatial dimensions in order to obtain results with higher correlation to the actual atomistic 3D structure of the protein and outline our heuristic approach for splitting large problem instances into smaller subproblems that can be directly solved with the current D-Wave 2000Q architecture. To the best of our knowledge, this work sets a new record for lattice protein folding on a quantum annealer by folding Chignolin (10 residues) on a planar lattice and Trp-Cage (8 residues) on a cubic lattice.

Sei Suzuki

Anders Irbäck, Lucas Knuthson, Sandipan Mohanty +1 more

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.