Quadratic unconstrained binary optimization and constraint programming approaches for lattice-based cyclic peptide
J Kyle Brubaker1, Kyle E C Booth2, Akihiko Arakawa3
1Amazon Advanced Solutions Lab, Seattle, WA, 98170, USA.
Scientific Reports
|July 1, 2025
Summary
Quadratic unconstrained binary optimization (QUBO) can model peptide-protein docking, but struggles with larger problems. Constraint programming offers a more scalable solution for this structural biology challenge.
Area of Science:
- Structural Biology
- Computational Biology
- Drug Design
Background:
- Peptide-protein interactions are crucial for biological processes and drug development.
- Accurate modeling of these interactions is essential for rational drug design.
- The peptide-protein docking problem presents significant computational challenges.
Purpose of the Study:
- To explore the application of quantum-amenable Quadratic Unconstrained Binary Optimization (QUBO) for peptide-protein docking.
- To incorporate peptide cyclization and docking constraints within a QUBO framework.
- To compare the performance of a novel QUBO approach against a constraint programming (CP) method.
Main Methods:
- Developed a resource-efficient QUBO encoding for peptide-protein docking on a tetrahedral lattice.
- Implemented an end-to-end framework for evaluating QUBO and CP methods using Protein Data Bank (PDB) instances.
- Utilized a classical simulated annealing solver for the QUBO approach.
- Developed a novel constraint programming (CP) approach for benchmarking.
Main Results:
- The QUBO approach successfully modeled feasible conformations for small peptide-protein docking problems (up to 6 peptide residues).
- QUBO performance degraded significantly with increasing problem size (e.g., PDB 3WNE, 5LSO).
- The CP approach demonstrated superior scalability, solving larger instances (up to 11 peptide residues, PDB 2F58).
Conclusions:
- QUBO can be applied to peptide-protein docking but exhibits scaling limitations.
- Constraint programming shows stronger performance and scalability for larger peptide-protein docking instances.
- CP may be a more suitable method than QUBO for complex peptide-protein docking problems in drug design.
More Related Videos
Related Concept Videos
Lattice Centering and Coordination Number
10.0K
The structure of a crystalline solid, whether a metal or not, is best described by considering its simplest repeating unit, which is referred to as its unit cell. The unit cell consists of lattice points that represent the locations of atoms or ions. The entire structure then consists of this unit cell repeating in three dimensions. The three different types of unit cells present in the cubic lattice are illustrated in Figure 1.
Types of Unit Cells
Imagine taking a large number of identical...
Types of Unit Cells
Imagine taking a large number of identical...
10.0K
Ligand Binding and Linkage
4.9K
Allosteric proteins have more than one ligand binding site; the binding of a ligand to any of these sites influences the binding of ligands to the other sites. When a protein is allosteric, its binding sites are called coupled or linked. In the case of enzymes, the site that binds to the substrate is known as the active site and the other site is known as the regulatory site. When a ligand binds to the regulatory site, this leads to conformational changes in the protein that can influence...
4.9K
Conserved Binding Sites
4.4K
Many proteins’ biological role depends on their interactions with their ligands, small molecules that bind to specific locations on the protein known as ligand-binding sites. Ligand-binding sites are often conserved among homologous proteins as these sites are critical for protein function.
Binding sites are often located in large pockets, and if their location on a protein’s surface is unknown, it can be predicted using various approaches. The energetic method computationally...
Binding sites are often located in large pockets, and if their location on a protein’s surface is unknown, it can be predicted using various approaches. The energetic method computationally...
4.4K
Ligand Binding Sites
13.4K
Proteins are dynamic macromolecules that carry out a wide variety of essential processes; however, the activities of most proteins depend on their interactions with other molecules or ions, known as ligands.
Protein-ligand interactions are quite specific; even though numerous potential ligands surround a cellular protein at any given time, only a particular ligand can bind to that protein. Moreover, a ligand binds only to a dedicated area on the surface of the protein, known as the...
Protein-ligand interactions are quite specific; even though numerous potential ligands surround a cellular protein at any given time, only a particular ligand can bind to that protein. Moreover, a ligand binds only to a dedicated area on the surface of the protein, known as the...
13.4K
Protein-protein Interfaces
13.5K
Many proteins form complexes to carry out their functions, making protein-protein interactions (PPIs) essential for an organism's survival. Most PPIs are stabilized by numerous weak noncovalent chemical forces. The physical shape of the interfaces determines the way two proteins interact. Many globular proteins have closely-matching shapes on their surfaces, which form a large number of weak bonds. Additionally, many PPIs occur between two helices or between a surface cleft and a...
13.5K
Conformations of Cyclohexane
13.2K
Cyclohexane does not exist in a planar form due to the high angle and torsional strain it would experience in the planar structure. Instead, it adopts non-planar chair and boat conformations.
The chair form is the most stable and derives its name from its resemblance to the “easy chair.” In the chair conformation, two carbon atoms are arranged out-of-plane — one above and one below, minimizing the torsional strain. In the chair form, the bond angle is very close to the ideal...
The chair form is the most stable and derives its name from its resemblance to the “easy chair.” In the chair conformation, two carbon atoms are arranged out-of-plane — one above and one below, minimizing the torsional strain. In the chair form, the bond angle is very close to the ideal...
13.2K


