REVIEW 3 major objections 4 minor 1 cited by
Designing lattice proteins with variational quantum algorithms
T0 review · 3 major / 4 minor · reviewed 2026-08-06 · deepseek-v4-flash
Pith's one-line read For lattice-protein sequence optimization, this paper shows that shallow problem-agnostic circuits outperform problem-informed QAOA circuits once device noise is included, because QAOA's circuit depth becomes prohibitive.
desk verdict A careful, reproducible QAOA-vs-HEA comparison on HP sequence optimization; the qualitative noise-tolerance ranking holds, but the hardware threshold N≤12 is one step beyond the data and the instance selection is narrow. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The carrying object is the QUBO reformulation of the HP contact energy. For a fixed target structure, the contact matrix $w_{ij}$ is known, so the only variables are the bead-type bits $s_i$; energy is a quadratic binary function with a Lagrange-penalty term $\lambda(\sum_i s_i - N_H)^2$ that fixes the number of hydrophobic beads. The paper contrasts two circuit families built on this cost: QAOA, whose phase separator and XY-mixers encode the problem and the composition constraint directly (with Dicke states preparing the feasible subspace), and HEA, which ignores problem structure and uses alternating parameterized $R_y/R_z$ rotations and CNOT entangling layers matched to the hardware topology. The metric that carries the comparison is the success rate—the fraction of runs that return the known unique optimal sequence, averaged over ten independent optimizations.
What would settle it
Run the same one-layer-HEA versus QAOA comparison on a randomized set of HP design instances that includes degenerately designable structures, multiple compositions per chain length, and targets without a known unique sequence; if QAOA matches or beats HEA on such instances, or if the N≈12 boundary shifts substantially, the paper's conclusion would not generalize. A second check is to repeat the Torino experiments with a noise model that includes temporal and correlated multi-qubit errors: if the simulation-hardware gap persists, the paper's stated explanation for the discrepancy is called into question.
Extended reading notes
Core claim
On the paper's own terms, the central finding is a comparison: problem-informed QAOA variants are not usable for HP sequence optimization on noisy intermediate-scale hardware, whereas a minimal one-layer HEA is, at least for short chains. The authors encode each bead type in one qubit and minimize $E(s)=-\sum_{i<j} w_{ij}s_i s_j + \lambda(\sum_i s_i - N_H)^2$, with $s_i=1$ for H and $0$ for P. Five QAOA variants—standard X-mixer, fully connected and ring XY-mixers, each with basis or Dicke initial states—yield acceptable or high success rates in noiseless simulation, but all five drop sharply when a hardware-derived noise model is applied; the fully connected XY-mixer with Dicke states, the best performer, has circuit depth above 2000 at N=16. One- and two-layer HEAs, whose circuits are much shallower, achieve comparable noiseless success rates and substantially better noisy success rates, with one-layer HEA the most noise-tolerant. Hardware runs on the Torino device confirm the HEA approach works for N≤11, with the paper concluding it can be used for N≤12, while noting that simulated success rates overestimate hardware performance, an effect the authors attribute to temporal and correlated multi-qubit errors absent from the noise model.
Load-bearing premise
The load-bearing premise is that the chosen test instances—one target structure and one composition value per chain length, each with a unique known solution that folds correctly—are representative of protein sequence optimization; the irregular size dependence of the results shows that instance-specific effects may dominate, so the ordering of algorithms could change on typical or degenerate problems.
Editorial extensions
If this is right
- One-layer HEA with warm-started parameters can be used to solve HP sequence optimization for short chains (N up to about 12) on current superconducting devices.
- QAOA, including constraint-preserving XY-mixer variants, is not competitive on this problem until circuit depth is reduced or error mitigation becomes effective.
- Hardware-tailored but problem-agnostic circuits can outperform problem-informed circuits under noise, so hardware compatibility is a key design criterion for NISQ-era variational optimization.
- Success rates from existing noise models are optimistic: ignoring temporal and correlated multi-qubit errors leads to systematic overestimation of hardware performance.
- Parameter donation across increasing chain lengths makes HEA optimization work for larger systems where random initialization fails.
Reading between the lines
- Extension: Because the depth penalty comes from the fully connected quadratic form, the QAOA-versus-HEA ordering may change under alternative encodings, such as domain-wall or parity encodings, that reduce entangling-layer depth.
- Extension: The success-rate metric, measured on uniquely solvable instances, is a good algorithm benchmark but not the criterion for real design, where near-optimal sequences that still fold to the target would be useful; a folding-aware evaluation could rank the algorithms differently.
- Extension: The same tradeoff likely applies to other fixed-backbone design problems with pairwise contact energies, not only HP lattice proteins, because the resource bottleneck is the quadratic all-to-all structure rather than protein-specific detail.
- Extension: The observed simulation-hardware gap makes hardware-aware noise models with temporal correlations, or noise-tailored variational training, a natural next test; if such a model closes the gap, it would support the paper's stated explanation.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper investigates the use of variational quantum algorithms (VQAs) for the first step of lattice-protein design: finding HP sequences that minimize the energy in a fixed target structure. The authors formulate the problem as a QUBO (Eq. 2) and compare five QAOA variants (Table I) with hardware-efficient ansatz (HEA) circuits, using noiseless simulations, simulations with the IBM Torino noise model, and hardware experiments on Torino. Problem instances are taken from Appendix A; for each chain length N (4 to 28) one target structure and one NH are selected so that the sequence-optimization problem has a unique known solution that folds to the target. The main findings are that QAOA variants perform well only in noiseless simulations and are defeated by circuit-depth-induced noise, whereas one-layer HEA has better noise tolerance in simulation but shows a sharp drop on hardware for N ≥ 12. The paper concludes that HEA could be used to solve the sequence optimization problem on Torino for short chains (N ≤ 12).
Significance. If the results are taken with appropriate qualifications, the paper provides a useful empirical comparison of problem-informed and problem-agnostic VQAs on a biologically motivated combinatorial optimization problem. Its strengths include the use of instances with externally verified unique ground states, statistics over 10 runs with standard errors, the combination of noiseless, noisy, and hardware measurements, and the public availability of the code. The qualitative conclusion that QAOA's depth makes it less noise-tolerant than shallow HEA is well supported. However, the quantitative hardware claim and its generalization to 'short chains' are not fully supported by the single-instance data, as detailed in the major comments.
major comments (3)
- [Sec. V versus Sec. III B / Fig. 8a] The concluding claim that HEA 'could be used to solve the sequence optimization problem on IBM's Torino device for short chains (N ≤ 12)' is one step beyond the reported hardware data. The text in Sec. III B states that one-layer HEA success rates are significant or high for N ≤ 11 but tiny for N ≥ 12, and Fig. 8a shows a sharp decline already at N = 12. The stated boundary should be N ≤ 11, or additional N = 12 hardware data (preferably on more than one instance) should be supplied.
- [Appendix A and Sec. III] The generality of the central claim is not established because only one hand-selected instance is used per chain length. For each N, Appendix A fixes one target structure and one NH, chosen so that the instance has a unique solution and known folding; these are the most designable structures, while only about 2% of HP sequences have a unique ground state. The irregular N-dependence visible in Figs. 3 and 7 shows that instance-specific effects likely dominate the observed success-rate behavior. Either test multiple instances per N (including instances with degenerate energy landscapes) or explicitly restrict the conclusion to the specific instances studied.
- [Sec. V and Sec. III B] The phrase 'could be used to solve' is not quantified against any classical baseline on the same instances. Because the instances are small and have known unique solutions, a comparison with a standard classical optimizer on the identical instances is needed to establish that the measured success rates represent a useful capability rather than trivial success on selected easy instances.
minor comments (4)
- [Sec. III A] The sentence 'Second, for a given initial state, the fully connected XY-mixer gives the results.' is missing a word; it should presumably read 'gives the best results.'
- [Sec. V] 'VQAs o ffers a promising approach' should be 'VQAs offer a promising approach', and the ligature should be removed.
- [Sec. II D] The phrase 'both angles initialized toπ' should have a space before π ('initialized to π').
- [Sec. III B and Fig. 8] The hardware results in Fig. 8 are shown for N up to 26, while the text discusses simulations up to N = 28; please clarify the range of the hardware experiments in the main text.
Circularity Check
No significant circularity: benchmarks are external exhaustive enumerations, and no prediction reduces to a fitted value.
full rationale
The paper's central claims are empirical success-rate measurements of QAOA and HEA variants against known ground states. Equation (2) is a standard QUBO formulation whose only hand-set constant, λ = 1.1, is declared and justified as robust once above a threshold, not fitted to the observed success rates. The ground-truth solutions and the guarantee that chosen instances have unique solutions that fold to the target structure are taken from exhaustive enumerations [33,34]. Although those references have author overlap with the present paper, they are parameter-free external results that supply the benchmark answers; they do not incorporate or depend on the VQA performance being reported, so the citation is not load-bearing in a circular way. Parameter donation is described as an initialization heuristic: optimized parameters from smaller instances seed larger ones, but the reported success rates are measured after optimizing on the target instance, so the measurement is not the fit itself. The paper's own caveats about irregular N-dependence, a single instance per chain length, and the overestimation of hardware success by noisy simulations are limitations on generality and model fidelity, not evidence that any claimed result is equivalent to its inputs by construction. No self-definitional step, fitted-input-as-prediction, or uniqueness-imported-from-authors pattern is present. Therefore the derivation chain is self-contained as an empirical benchmark study, and the appropriate circularity score is 0.
Assumptions & free parameters
free parameters (1)
- Lagrange multiplier λ =
1.1
assumptions (4)
- domain assumption The 2D HP model with contact energy EHP = -NHH is an adequate test bed for protein design.
- standard math Exhaustive enumerations in Refs. [33,34] correctly identify unique ground states and designability for N ≤ 30.
- standard math The QUBO objective in Eq. (2) with the λ penalty is equivalent to minimizing EHP subject to the composition constraint NH.
- domain assumption The IBM Torino noise model (gate errors, readout errors, thermal relaxation) captures the dominant error sources in noisy simulations.
Cite this review
Pith. "Pith review of Designing lattice proteins with variational quantum algorithms." pith.science (2026). https://pith.science/paper/M2VXDZSM
@misc{pith2026250802369,
author = {Pith},
title = {Pith review of: Designing lattice proteins with variational quantum algorithms},
year = {2026},
howpublished = {\url{https://pith.science/paper/M2VXDZSM}},
note = {Machine review of arXiv:2508.02369}
}
read the original abstract
Quantum heuristics have shown promise in solving various optimization problems, including lattice protein folding. Equally relevant is the inverse problem, protein design, where one seeks sequences that fold to a given target structure. The latter problem is often split into two steps: (i) searching for sequences that minimize the energy in the target structure, and (ii) testing whether the generated sequences fold to the desired structure. Here, we investigate the utility of variational quantum algorithms for the first of these two steps on today's noisy intermediate-scale quantum devices. We focus on the sequence optimization task, which is less resource-demanding than folding computations. We test the quantum approximate optimization algorithm and variants of it, with problem-informed quantum circuits, as well as the hardware-efficient ansatz, with problem-agnostic quantum circuits. While the former algorithms yield acceptable results in noiseless simulations, their performance drops under noise. With the problem-agnostic circuits, which are more compatible with hardware constraints, an improved performance is observed in both noisy and noiseless simulations. However, the results deteriorate when running on a real quantum device. We attribute this discrepancy to features not captured by the simulated noise model, such as the temporal aspect of the hardware noise.
Figures
Figures from the paper (7 more)
Forward citations
Cited by 1 Pith paper
-
Penalty-free quantum optimization applied to lattice protein folding
A QAOA variant without quadratic penalties, using independent sets in a conflict graph, is applied to lattice protein folding and validated on proteins up to length 14 via simulation and heuristic search.
Reference graph
Works this paper leans on
-
[1]
K. Wang, Z. Lu, C. Zhang, et al. , Demonstration of low-overhead quantum error correction codes (2025), arXiv:2505.09684 [quant-ph]
arXiv 2025
- [2]
- [3]
-
[4]
R. Acharya, D. A. Abanin, L. Aghababaie-Beni, et al., Nature 638, 920–926 (2024)
work page 2024
- [5]
- [6]
- [7]
- [8]
Show all 51 references
-
[9]
Preskill, Quantum 2, 79 (2018)
J. Preskill, Quantum 2, 79 (2018)
2018
-
[10]
J. R. McClean, S. Boixo, V . N. Smelyanskiy, et al., Nat. Com- mun. 9, 4812 (2018)
2018
-
[11]
Langfitt, J
Q. Langfitt, J. Falla, I. Safro, and Y . Alexeev, in2023 IEEE In- ternational Conference on Quantum Computing and Engineer- ing (QCE), V ol. 02 (2023) pp. 300–301
2023
-
[12]
Montañez-Barrera, D
A. Montañez-Barrera, D. Willsch, A. Maldonado-Romo, and K. Michielsen, Quantum Sci. Technol. 9, 025022 (2024)
2024
-
[13]
Galda, E
A. Galda, E. Gupta, J. Falla, et al., Front. Quantum Sci. Tech- nol. 2, 1200975 (2023)
2023
-
[14]
Lyngfelt and L
I. Lyngfelt and L. García-Álvarez, Phys. Rev. A 111, 022418 (2025)
2025
-
[15]
Kosen, H.-X
S. Kosen, H.-X. Li, M. Rommel, et al., Quantum Sci. Technol. 7, 035018 (2022)
2022
- [16]
-
[17]
Robert, P
A. Robert, P. K. Barkoutsos, S. Woerner, and I. Tavernelli, Npj Quantum Inf. 7, 38 (2021)
2021
-
[18]
Boulebnane, X
S. Boulebnane, X. Lucas, A. Meyder, et al., Npj Quantum Inf. 9, 70 (2023)
2023
-
[19]
Perdomo-Ortiz, N
A. Perdomo-Ortiz, N. Dickson, M. Drew-Brook, et al. , Sci. Rep. 2, 248 (2012)
2012
-
[20]
Outeiral, G
C. Outeiral, G. M. Morris, J. Shi,et al., New J. Phys.23, 103030 (2021)
2021
-
[21]
Irbäck, L
A. Irbäck, L. Knuthson, S. Mohanty, and C. Peterson, Phys. Rev. Res. 4, 043013 (2022)
2022
-
[22]
H. Linn, I. Brundin, L. García-Álvarez, and G. Johansson, Phys. Rev. Res. 6, 033112 (2024)
2024
-
[23]
Kuhlman, G
B. Kuhlman, G. Dantas, G. C. Ireton, et al., Science 302, 1364 (2003)
2003
-
[24]
Bhardwaj, V
G. Bhardwaj, V . K. Mulligan, C. D. Bahl, et al. , Nature 538, 329 (2016)
2016
-
[25]
K. K. Yang, Z. Wu, and F. H. Arnold, Nat. Methods 16, 687 (2019)
2019
-
[26]
Kuhlman and P
B. Kuhlman and P. Bradley, Nat. Rev. Mol. Cell Biol. 20, 681 (2019)
2019
-
[27]
L. Cao, I. Goreshnik, B. Coventry, et al. , Science 370, 426 (2020)
2020
-
[28]
V . K. Mulligan, H. Melo, H. I. Merritt, et al. , bioRxiv:10.1101/752485 (2020)
2020 doi
-
[29]
Irbäck, L
A. Irbäck, L. Knuthson, S. Mohanty, and C. Peterson, Phys. Rev. Res. 6 (2024)
2024
-
[30]
Panizza, P
V . Panizza, P. Hauke, C. Micheletti, and P. Faccioli, PRX Life 2, 043012 (2024)
2024
-
[31]
M. H. Khatami, U. C. Mendes, N. Wiebe, and P. M. Kim, PLOS Comput. Biol. 19, 1 (2023)
2023
-
[32]
K. F. Lau and K. A. Dill, Macromolecules 22, 3986 (1989)
1989
-
[33]
Irbäck and C
A. Irbäck and C. Troein, J. Biol. Phys. 28, 1 (2002)
2002
-
[34]
Holzgräfe, A
C. Holzgräfe, A. Irbäck, and C. Troein, J. Chem. Phys. 135, 195101 (2011)
2011
-
[35]
Irbäck, C
A. Irbäck, C. Peterson, F. Potthast, and E. Sandelin, Structure 7, 347 (1999)
1999
-
[36]
Aina and S
A. Aina and S. Wallin, J. Chem. Phys. 147, 095102 (2017)
2017
-
[37]
Nilsson and A
D. Nilsson and A. Irbäck, Phys. Rev. E 101, 022413 (2020)
2020
-
[38]
Statt, H
A. Statt, H. Casademunt, C. P. Brangwynne, and A. Z. Pana- giotopoulos, J. Chem. Phys. 152, 075101 (2020)
2020
-
[39]
Bornberg-Bauer and H
E. Bornberg-Bauer and H. S. Chan, Proc. Natl. Acad. Sci. USA 96, 10689 (1999)
1999
-
[40]
Aguirre, P
J. Aguirre, P. Catalán, J. A. Cuesta, and S. Manrubia, Open Biol. 8, 180069 (2018)
2018
-
[41]
Kadowaki and H
T. Kadowaki and H. Nishimori, Phys. Rev. E 58, 5355 (1998)
1998
-
[42]
Farhi, J
E. Farhi, J. Goldstone, S. Gutmann, et al. , Science 292, 472 (2001)
2001
-
[43]
Bärtschi and S
A. Bärtschi and S. Eidenbenz, in Fundamentals of Computation Theory, edited by L. A. G˛ asieniec, J. Jansson, and C. Levcopou- los (Springer, Cham, 2019) pp. 126–139
2019
-
[44]
Bittel and M
L. Bittel and M. Kliesch, Phys. Rev. Lett. 127, 120502 (2021)
2021
-
[45]
Zhou, S.-T
L. Zhou, S.-T. Wang, S. Choi, et al., Phys. Rev. X 10, 021067 (2020)
2020
-
[46]
J. A. Montañez-Barrera, D. Willsch, and K. Michielsen, Quan- tum Inf. Process. 24, 129 (2025)
2025
-
[47]
C. R. Harris, K. J. Millman, S. J. van der Walt, et al. , Nature 585, 357 (2020)
2020
-
[48]
Virtanen, R
P. Virtanen, R. Gommers, T. E. Oliphant, et al., Nat. Methods 17, 261 (2020)
2020
-
[49]
J. D. Hunter, Computing in Science & Engineering9, 90 (2007)
2007
-
[50]
Javadi-Abhari, M
A. Javadi-Abhari, M. Treinish, K. Krsulich, et al. , Quantum computing with Qiskit (2024), arXiv:2405.08810 [quant-ph]
2024 arXiv
-
[51]
M. J. D. Powell, A direct search optimization method that mod- els the objective and constraint functions by linear interpo- lation, in Advances in Optimization and Numerical Analysis , edited by S. Gomez and J.-P. Hennart (Springer Netherlands, Dordrecht, 1994) pp. 51–67
1994
Reviewed August 6, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.