Classical and current quantum approaches to protein structure prediction (QPSP) face limitations, notably massive qubit requirements restricting near-term models to simplistic on-lattice simulations. We propose a logarithmic-scale variational quantum eigensolver (VQE) that reduces qubit requirements for N torsional degrees of freedom to O(log2N), enabling off-lattice, all-atom simulations. Our architecture extracts molecular torsions from relative phases in statevector simulations. On quantum hardware, a decoder maps the empirical cumulative distribution function (CDF) from basis-state probabilities to bounded torsional variables. These feed a classical algorithm to build heavy-atom coordinates. We use an EfficientSU2 ansatz and multi-stage relaxation to mitigate barren plateaus. Structures are evaluated via a custom hybrid quantum-classical Hamiltonian, alongside Rosetta and OpenMM benchmarks.
Evaluation on chignolin and Trp-cage yielded native-like conformations. Chignolin reached a 0.623 Å Cα RMSD in retained snapshots and 1.199 Å in final models; Trp-cage achieved a 2.501 Å RMSD among snapshots (3.512 Å in final models). Execution on IBM processors (ibm_cleveland, ibm_miami) successfully recovered native-like structures with a best RMSD of 1.758 Å. The custom energy function performed best overall, though energy-ranking imbalances persisted across sampled landscapes for all functions.
This introduces the first all-atom, continuous-space quantum algorithm for QPSP. By converting physical qubit constraints into circuit depth constraints, it proves high-resolution prediction is feasible with exponentially fewer qubits. Despite current limits like computational overhead and energy function sensitivity, it establishes a scalable foundation for hybrid quantum biophysics.
Classical and current quantum approaches to protein structure prediction (QPSP) face limitations, notably massive qubit requirements restricting near-term models to simplistic on-lattice simulations. We propose a logarithmic-scale variational quantum eigensolver (VQE) that reduces qubit requirements for N torsional degrees of freedom to O(log2N), enabling off-lattice, all-atom simulations. Our architecture extracts molecular torsions from relative phases in statevector simulations. On quantum hardware, a decoder maps the empirical cumulative distribution function (CDF) from basis-state probabilities to bounded torsional variables. These feed a classical algorithm to build heavy-atom coordinates. We use an EfficientSU2 ansatz and multi-stage relaxation to mitigate barren plateaus. Structures are evaluated via a custom hybrid quantum-classical Hamiltonian, alongside Rosetta and OpenMM benchmarks.
Evaluation on chignolin and Trp-cage yielded native-like conformations. Chignolin reached a 0.623 Å Cα RMSD in retained snapshots and 1.199 Å in final models; Trp-cage achieved a 2.501 Å RMSD among snapshots (3.512 Å in final models). Execution on IBM processors (ibm_cleveland, ibm_miami) successfully recovered native-like structures with a best RMSD of 1.758 Å. The custom energy function performed best overall, though energy-ranking imbalances persisted across sampled landscapes for all functions.
This introduces the first all-atom, continuous-space quantum algorithm for QPSP. By converting physical qubit constraints into circuit depth constraints, it proves high-resolution prediction is feasible with exponentially fewer qubits. Despite current limits like computational overhead and energy function sensitivity, it establishes a scalable foundation for hybrid quantum biophysics.
Variational quantum algorithms provide a direct, physics-based approach to protein structure prediction, but their accuracy is limited by the coarse resolution of the energy landscapes generated on current noisy devices. We propose a hybrid framework that combines quantum computation with deep learning, formulating structure prediction as a problem of energy fusion. Candidate conformations are obtained through the Variational Quantum Eigensolver (VQE) executed on IBM’s 127 qubit superconducting processor, which defines a global yet low resolution quantum energy surface. To refine these basins, secondary structure probabilities and dihedral angle distributions predicted by the NSP3 neural network are incorporated as statistical potentials. These additional terms sharpen the valleys of the quantum landscape, resulting in a fused energy function that enhances effective resolution and better distinguishes native like structures. Evaluation on 375 conformations from 75 protein fragments shows that, on this short fragment-only benchmark (10–14 residues without surrounding chain context), the hybrid pipeline achieves a mean RMSD of 4.9 Å with statistical significance ( p < 0.001), outperforming both quantum-only re-rankings and fragment-level runs of AlphaFold3 and ColabFold. We caution that this conclusion is restricted to the fragment-scale regime accessible to current quantum hardware and does not extend to the full-chain regime in which AlphaFold3 and ColabFold are usually evaluated. The findings demonstrate that energy fusion offers a systematic method for combining data driven models with quantum algorithms, improving the practical applicability of near term quantum computing to molecular and structural biology.
Protein structure prediction remains central to computational biology and drug discovery. While deep learning models such as AlphaFold have achieved remarkable accuracy for known fold families, they struggle with novel folds, intrinsically disordered regions, and conformational dynamics, and they do not explicitly optimize physics-based energy functions. We present a quantum-assisted approach to protein conformation optimization that uses parameterized quantum circuits within a variational framework to explore the energy landscape of small proteins and fragments from first-principles or coarse-grained physics. Our method encodes backbone torsion angles (ϕ, ψ) and selected side-chain degrees of freedom into qubit registers and employs a variational quantum eigensolver (VQE)-inspired optimization to minimize an approximate Hamiltonian incorporating bonded interactions, sterics, hydrogen bonding, and hydrophobic contacts. We introduce problem-inspired ansätze that respect chain connectivity and secondary-structure motifs (helices, beta-sheets) and propose hierarchical circuit growth strategies for scaling from peptides to larger domains. Through classical simulations of small peptides and secondary-structure motifs, we demonstrate that quantum variational optimization can identify low-energy conformations consistent with experimentally determined structures. We provide explicit resource estimates for NISQ and early fault-tolerant quantum computers and discuss hybrid classical-quantum strategies for refining quantum solutions with classical molecular dynamics and machinelearning priors. While full protein-folding quantum advantage remains distant, we identify concrete subproblems where nearterm quantum devices offer complementary physics-based insight beyond data-driven approaches, opening pathways toward practical quantum-classical protein structure and drug discovery pipelines.
Accurate prediction of protein active-site structures remains a central challenge in structural biology, especially for short and flexible peptide fragments where conventional and simulation-based methods often fail. Here, we present a quantum computing framework designed for utility-level quantum processors to address this problem. Starting from an amino acid sequence, we cast structure prediction as a ground-state energy minimization task using the Variational Quantum Eigensolver (VQE). Amino acid connectivity is represented on a tetrahedral lattice, and steric, geometric, and chirality constraints are encoded into a problem-specific Hamiltonian expressed as sparse Pauli operators. A two-stage architecture separates energy estimation from measurement decoding, enabling noise mitigation under realistic device conditions. We evaluate the method on 23 real protein fragments from the PDBbind dataset and 7 fragments from therapeutically relevant proteins, executing all experiments on the IBM-Cleveland Clinic quantum processor. Structural predictions are benchmarked against AlphaFold3 (AF3) and classical simulation-based approaches using identical postprocessing and docking procedures. Our quantum framework outperforms both AF3 and classical baselines in Root-Mean-Square Deviation (RMSD) and docking efficacy, demonstrating a practical end-to-end pipeline for biologically relevant structure prediction on real quantum hardware and highlighting its engineering feasibility for near-term quantum devices.
Abstract Near-term quantum devices remain constrained by limited qubit counts, noise, and shallow circuit depth, restricting their direct application to peptide folding. This work addresses this limitation through a hybrid quantum–classical framework that separates global contact-space exploration from local configuration-space refinement. The pipeline was evaluated on the p53 transactivation peptide. The peptide was encoded using discrete torsional states and Miyazawa–Jernigan contact potentials. Global conformational search was performed using simulated annealing, followed by fragment-level CVaR-based variational quantum eigensolver refinement. Overlapping fragments were stitched and reconstructed into three-dimensional backbone coordinates for structural comparison. Three independent simulated annealing restarts produced negative-energy conformers. The workflow identified a best global conformer at E = −3.70 MJ using simulated annealing. Fragment-level CVaR-VQE refinement generated 30 candidate conformers across three overlapping fragments, the best assembled conformer reaching E = −2.58 MJ after DP stitching and post-stitch refinement. In this trial, the classical SA solution outperformed the VQE-assembled structure and was selected as the final output. The final 65-atom model was successfully reconstructed, but full-backbone alignment against 1YCR yielded an RMSD of 4.65 Å, improving to 2.90 Å after trimmed alignment. Our work demonstrates a practical quantum-enabled peptide folding workflow for short, flexible peptides, but the present structural accuracy is insufficient for claims of native fold recovery. Its main value is as an auditable framework for testing NISQ-compatible folding strategies.
- A study published in the international journal npj Quantum Information (Nature Portfolio) presents a resource-efficient quantum algorithm framework based on the Variational Quantum Eigensolver (VQE) for peptide and protein folding. - By translating protein folding into a tetrahedral lattice model and combining Conditional Value-at-Risk (CVaR) VQE with a classical evolutionary algorithm, the number of required qubits was drastically reduced. - It successfully folded a 7-amino acid neuropeptide using only 9 qubits on actual IBM 20-qubit quantum hardware, demonstrating its practical potential for next-generation drug discovery in the NISQ (Noisy Intermediate-Scale Quantum) era. [Quantum Biology Society] The protein folding problem, which involves predicting the three-dimensional structure of proteins and peptides from their primary amino acid sequence, plays a central role in various applications such as chemistry, biology, and medicine. However, as the sequence lengthens, the number of possible structures to explore explodes exponentially, making it difficult for classical computers to handle its intrinsic NP-hard complexity. Recently, a new quantum algorithm that can efficiently solve this long-standing challenge even on noisy intermediate-scale quantum (NISQ) computers has been developed, drawing the attention of the academic community. The paper titled "Resource-efficient quantum algorithm for protein folding," published in the international journal npj Quantum Information, presents an innovative computing framework for the discovery of next-generation peptide drugs. The research team of Anton Robert and Ivano Tavernelli from IBM Quantum, IBM Research - Zurich in Switzerland successfully led this groundbreaking study. ■ The Meeting of Lattice Models and Qubit Optimization Until now, attempts to introduce quantum computing to protein folding have required an unattainably high number of qubits that current, near-term quantum hardware cannot handle. To solve this, the research team converted protein folding into a tetrahedral lattice model. This model can efficiently depict the growth of a polymer chain while maintaining chemically plausible bond angles. In particular, they adopted a method of adding penalty terms to the Hamiltonian to prevent unphysical amino acid overlaps (e.g., beads occupying the same lattice site). Through this, they succeeded in drastically reducing the number of qubits required for calculation while suppressing the complexity of the algorithm to the level of O(N⁴). ■ Demonstration of Folding on Actual Quantum Hardware This algorithm maximized efficiency by combining CVaR-VQE, a variation of the Variational Quantum Eigensolver (VQE), with a classical Differential Evolution optimizer. Applying this technique, the researchers simulated the folding of the 10-amino acid Angiotensin peptide using only 22 qubits, successfully finding the most stable energy state (native fold). Furthermore, they successfully demonstrated the folding of a 7-amino acid neuropeptide (APRLRFY) using only 9 qubits on IBM's 20-qubit quantum computer hardware. This study is a groundbreaking achievement proving that the protein folding problem, which previously required massive computing resources, can be successfully solved even in a noisy real quantum processor environment. It is expected to bring tremendous innovation to next-generation /de novo/ drug design and disease treatment research utilizing quantum computing in the future.
The variational quantum eigensolver (VQE) is a hybrid quantum classical algorithm designed for current and near-term quantum devices. Despite its initial success, there is a lack of understanding involving several of its key aspects. There are problems with VQE that forbid a favorable scaling towards quantum advantage. In order to alleviate the problems, we propose and extensively test a quantum annealing inspired heuristic that supplements VQE. The improved VQE enables an efficient initial state-preparation mechanism, in a recursive manner, for a quasidynamical unitary evolution. We conduct an in-depth scaling analysis of finding the ground-state energies with increasing lattice sizes of the Heisenberg model, employing simulations of up to 40 qubits that manipulate the complete state vector. In addition to systematically finding the ground-state energy, we observe that it avoids barren plateaus, escapes local minima, and works with low-depth circuits. For the current devices, we further propose a benchmarking toolkit using a mean-field model and test it on IBM Q devices. Realistic gate execution times estimate a longer computational time to complete the same computation on a fully functional error-free quantum computer than on a quantum computer emulator implemented on a classical computer. However, our proposal can be expected to help accurate estimations of the ground-state energies beyond 50 qubits when the complete state vector can no longer be stored on a classical computer, thus enabling quantum advantage.