Recent advancements in quantum computing have created new opportunities for tackling long-standing, complex combinatorial optimization problems that are intractable for classical computers. Predicting the secondary structure of mRNA is one such notoriously difficult problem that can benefit from the rapid advancements in quantum computing technology. Accurate prediction of mRNA secondary structure is crucial for designing RNA-based therapeutics, as it dictates various stages of the mRNA life cycle, including transcription, translation, and decay. The current generation of quantum computers has reached a utility scale, enabling us to explore relatively large problem sizes. In this paper, we examine the feasibility of solving mRNA secondary structures on a quantum computer for sequence lengths up to 60 nucleotides, representing problems in the qubit range of 10 to 80. We use the Conditional Value at Risk (CVaR)-based VQE algorithm to solve optimization problems derived from mRNA secondary structure prediction on the IBM Eagle and Heron quantum processors. Encouragingly, even with minimal error mitigation and fixed-depth circuits, our hardware runs yield accurate predictions of minimum free energy (MFE) structures that match the results of the classical solver CPLEX. Our results provide substantial evidence for the viability of solving mRNA structure prediction problems on a quantum computer and motivate continued research in this direction.
Recent advancements in quantum computing have created new opportunities for tackling long-standing, complex combinatorial optimization problems that are intractable for classical computers. Predicting the secondary structure of mRNA is one such notoriously difficult problem that can benefit from the rapid advancements in quantum computing technology. Accurate prediction of mRNA secondary structure is crucial for designing RNA-based therapeutics, as it dictates various stages of the mRNA life cycle, including transcription, translation, and decay. The current generation of quantum computers has reached a utility scale, enabling us to explore relatively large problem sizes. In this paper, we examine the feasibility of solving mRNA secondary structures on a quantum computer for sequence lengths up to 60 nucleotides, representing problems in the qubit range of 10 to 80. We use the Conditional Value at Risk (CVaR)-based VQE algorithm to solve optimization problems derived from mRNA secondary structure prediction on the IBM Eagle and Heron quantum processors. Encouragingly, even with minimal error mitigation and fixed-depth circuits, our hardware runs yield accurate predictions of minimum free energy (MFE) structures that match the results of the classical solver CPLEX. Our results provide substantial evidence for the viability of solving mRNA structure prediction problems on a quantum computer and motivate continued research in this direction.
Accurate prediction of mRNA secondary structure is critical for understanding gene expression, translation efficiency, and advancing mRNA-based therapeutics. However, the combinatorial complexity of possible foldings, especially in long sequences, poses significant computational challenges for classical algorithms. In this work, we propose a scalable, quantum-centric optimization framework that integrates quantum sampling with classical post-processing to tackle this problem. Building on a Quadratic Unconstrained Binary Optimization (QUBO) formulation of the mRNA folding task, we develop two complementary workflows: a Conditional Value at Risk (CVaR)-based variational quantum algorithm enhanced with gauge transformations and local search, and an Instantaneous Quantum Polynomial (IQP) circuit-based scheme where training is done classically and sampling is delegated to quantum hardware. We demonstrate the effectiveness of these approaches using IBM quantum processors, solving problem instances with up to 156 qubits and circuits containing up to 950 nonlocal gates, corresponding to mRNA sequences of up to 60 nucleotides. Additionally, we validate scalability of the CVaR algorithm on a tensor network simulator, reaching up to 354 qubits in noiseless settings. These results demonstrate the growing practical capabilities of hybrid quantum-classical methods for tackling large-scale biological optimization problems.
Co-optimizing mRNA sequences for both codon optimality and secondary structure is crucial for producing stable and efficacious mRNA therapeutics. Codon optimization, which adjusts nucleotide sequences to enhance translational efficiency, inherently influences mRNA secondary structure - a key determinant of molecular stability both in-vial and in-cell. Because both properties are governed by the same underlying sequence, optimizing one directly impacts the other. To address this interdependence, we introduce a novel variational framework that simultaneously optimizes codon usage and secondary structure. Our method employs a dual-objective function that balances the codon adaptation index (CAI) and minimum free energy (MFE), incorporating variational parameters for codon selection. Leveraging a hybrid quantum-classical computational strategy and building on prior advancements in quantum algorithms for secondary structure prediction, we effectively navigate this complex optimization space. We demonstrate the feasibility of executing this end-to-end workflow on real quantum hardware, using IBM's 127-qubit Eagle processor. We validate our approach through both simulations and quantum hardware experiments on sequences of up to 30 nucleotides. These results highlight the potential of our framework to accelerate the design of optimal mRNA constructs for therapeutic and research applications.
It is the first step for understanding how RNA structure folds from base sequences that to know how its secondary structure is formed. Traditional energy-based algorithms are short of precision, particularly for non-nested sequences, while learning-based algorithms face challenges in obtaining high-quality training data. Recently, quantum annealer has rapidly predicted the folding of the secondary structure, highlighting that quantum computing is a promising solution to this problem. However, gate model algorithms for universal quantum computing are not available. In this paper, gate-based quantum algorithms will be presented, which are highly flexible and can be applied to various physical devices. Mapped all possible secondary structure to the state of a quadratic Hamiltonian, the whole folding process is described as a quadratic unconstrained binary optimization model. Then the model can be solved through quantum approximation optimization algorithm. We demonstrate the performance with both numerical simulation and experimental realization. Throughout our benchmark dataset, simulation results suggest that our quantum approach is comparable in accuracy to classical methods. For non-nested sequences, our quantum approach outperforms classical energy-based methods. Experimental results also indicate our method is robust in current noisy devices. It is the first instance of universal quantum algorithms being employed to tackle RNA folding problems, and our work provides a valuable model for utilizing universal quantum computers in solving RNA folding problems.
Protein folding has been one of the most difficult problems for over a half-century,The random thermal motions causing conformational changes that lead energetically downhill towards the native structure, a principle captured in funnel-shaped energy landscapes.Unfolded polypeptides have a wide range of possible conformations. The search problem becomes intractable for classical computers due to the exponential growth of potential conformations with chain length. So far, there is theoretical and experimental evidence that solving such optimization problems using Quantum Computing approaches such as Quantum Annealing, VQE, and QAOA has an advantage. Although Google's DeepMind-AlphaFold has accomplished much,But we can go even further with the quantum approach.Here we show how to predict structure of protein as well as RNA folding using the Variational Quantum Eigensolver with Conditional Value at Risk (CVaR) expectation values for the solution of the problem and for finding the minimum configuration energy and our task is to identify a protein's minimal energy structure. The protein's structure is optimized to reduce energy. Also making sure that all physical constraints are met and encoding the protein folding problem into a qubit operator.
Optimizing the mRNA codon has an essential impact on gene expression for a specific target protein. It is an NP-hard problem; thus, exact solutions to such optimization problems become computationally intractable for realistic problem sizes on both classical and quantum computers. However, approximate solutions via heuristics can substantially impact the application they enable. Quantum approximate optimization is an alternative computation paradigm promising for tackling such problems. Recently, there has been some research in quantum algorithms for bioinformatics, specifically for mRNA codon optimization. This research presents a denser way to encode codons for implementing mRNA codon optimization via the variational quantum eigensolver algorithms on a gate-based quantum computer. This reduces the qubit requirement by half compared to the existing quantum approach, thus allowing longer sequences to be executed on existing quantum processors. The performance of the proposed algorithm is evaluated by comparing its results to exact solutions, showing well-matching results.
Pauli Correlation Encoding (PCE) compresses m binary variables onto n=O(m1/k) qubits by mapping them to commuting Pauli correlators, but its continuous expectation values must be decoded into feasible binary solutions, a challenge for dense-constraint problems. We apply PCE to mRNA secondary-structure prediction, formulated as a densely constrained QUBO, and train with a QUBO-space sigmoid loss thatpreserves the QUBO penalty structure. For decoding, we introduce the Problem-Aware Guided Decoder (PAGD), which scores candidate variable commitments by combining marginal QUBO energy reduction with a trained expectation-value prior and constraint-aware feasibility pruning. On six benchmark mRNA sequences (30-60 nt, 50-240 variables, 7-14 qubits), PAGD with 100 restarts achieves 75-100 percent near-optimal recovery, defined as P(gap<1%), for sequences up to 152 variables, compared with 0-30 percent for a sign-rounding plus local-search baseline. On the 240-variable instance, trained PAGD reaches 50 percent P(gap<1%) at 200 restarts, outperforming untrained-circuit and random-expectation-value controls. Hardware-scale tests extend the pipeline to three 102-105 nt instances (694-745 variables, 172,000-193,000 pair constraints, 23 qubits) on IBM Heron processors. The circuits transpile SWAP-free into 480 native two-qubit gates at depth 256, and PAGD decoded gaps on QPU runs match or beat simulator means for all three instances, including exact CPLEX-optimum recovery for one sequence. These results show that PCE-trained priors can survive deployment to noisy superconducting hardware at biologically relevant scale.