Scalable Decoupling Contact and Configuration Spaces in a Quantum-Enabled Protein Folding Workflow
- 1R V College of Engineering
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 angstrom, improving to 2.90 angstrom 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.
Related
- Scalable Contact-Gated Stochastic Operators and Universal Reduction Laws for Macromolecular Protein Dynamics, Allostery, and Kinetics
- A Scalable Workflow for Extracting Species Interactions From Ecological Literature Using Large Language Models
- Machine-Learned Compact Subspace Generation for Quantum Selected Configuration Interaction within Density Matrix Embedding Framework
- Evaluation and validation of a reproducible mNGS workflow for gut Microbiome analysis
- 93 Swiss beads calibration enables quantitative quality control for high-dimensional rare-cell mass cytometry in an NIH-supported clinical research workflow
BibTeX
@misc{uttarkar2026scalable,
title = {{Scalable Decoupling Contact and Configuration Spaces in a Quantum-Enabled Protein Folding Workflow}},
author = {A. Uttarkar and V. Niranjan},
howpublished = {bioRxiv},
year = {2026},
month = oct,
doi = {10.64898/2026.10.09.758014},
url = {https://www.biorxiv.org/content/10.64898/2026.10.09.758014v1},
note = {CC BY-NC-ND},
}
RIS
TY - UNPB
TI - Scalable Decoupling Contact and Configuration Spaces in a Quantum-Enabled Protein Folding Workflow
AU - A. Uttarkar
AU - V. Niranjan
PY - 2026
DA - 2026/10/09
PB - bioRxiv
DO - 10.64898/2026.10.09.758014
UR - https://www.biorxiv.org/content/10.64898/2026.10.09.758014v1
AB - 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 angstrom, improving to 2.90 angstrom 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.
ER -