Q-Score is introduced, encoding GNN-predicted orbital donor-acceptor energies into a weighted graph and scoring binding by solving a maximum-weight vertex clique problem via Digitized-Counterdiabatic QAOA, enriching for strong orbital interactions at twice the random rate.
Abstract
Molecular docking predicts how a small molecule binds to a protein and is a key bottleneck in drug discovery. Classical scoring functions sum empirical pairwise contacts, blind to quantum-mechanical effects like orbital charge transfer that govern binding specificity. We introduce Q-Score, encoding GNN-predicted orbital donor-acceptor energies into a weighted graph and scoring binding by solving a maximum-weight vertex clique problem via Digitized-Counterdiabatic QAOA. Each interaction anchor maps to one qubit and compatibility constraints become edges. Across 11 protein targets, DC-QAOA recovers the exact optimum on 8 at 10 qubits. On 1000 AI-generated molecules, Q-Score is orthogonal to classical scoring with Spearman rho of 0.05, driven by orbital quality with rho of 0.90, and free of molecular-weight bias, enriching for strong orbital interactions at twice the random rate. DC-QAOA achieves a mean approximation ratio of 0.94 with 52 percent exact. Execution of 1000 circuits on IBM Eagle confirms 6-qubit solvability on NISQ hardware.
Covalent virtual screening requires ranking compounds according to both noncovalent recognition and their ability to adopt a reaction-competent geometry with an appropriately reactive warhead. Here, we introduce BCover, a reaction-aware scoring suite that combines pre-reactive docking with quantum-chemistry-derived ligand reactivity descriptors, the electrostatic properties of the protein pocket. Quantum-chemical descriptors are aggregated using a nonlinear tree-based scoring model. BCover was evaluated retrospectively on the COValid benchmark, comprising nine targets and ten reactive sites, and compared with AutoDock, DOCK6, DOCKovalent, AlphaFold3 with Rosetta rescoring, and AlphaFold3 confidence-based ranking. BCover achieved an average adjusted LogAUC of 31% (57% max), an average ROC-AUC of 0.88 (0.96 max), an average EF1 of 18 (33 max). Its average LogAUC exceeded those of the classical docking methods and AlphaFold3-Rosetta, although AlphaFold3-mPAE provided the strongest overall enrichment. At an average runtime of about 10~s per ligand, BCover was approximately 25-fold faster than the evaluated AlphaFold3 workflows and achieved the highest average time-adjusted virtual-screening productivity index. Redocking experiments further showed that the method recovered near-native ligand conformations. These results demonstrate that combining docking-derived geometry with ligand local electronic reactivity and pocket electrostatics provides an efficient and interpretable strategy for covalent ligand prioritization. BCover is intended as a high-throughput screening method that complements more computationally demanding QM/MM and free-energy calculations during subsequent lead optimization.
Accurate identification of near-native ligand binding poses is a central challenge in structure-based drug design. From a physical point of view, the successful construction of a protein–ligand complex structure is dependent on whether protein and ligand can form enough atomically pairwise interactions that result in a global energy minimum. In this work, we report a machine learning scoring strategy for protein–ligand screening which explicitly considers the Native Contact Ratio (NCR), a topology inspired metric that quantifies the preservation of protein–ligand interfacial contacts as well as interaction energy. This physics-awared supervision strategy provides a simple but efficient gradient field that faithfully reflects the complicated protein energy landscape than conventional 3D coordinate-based objectives. Building on this principle, we present DeepNCR, an energy-informed Transformer framework that encodes approximate Coulombic and dispersive interaction potentials across the protein–ligand binding interface. Furthermore, we introduce a feature pruning step that compresses the interaction tensor from 1470 to 868 dimensions, further improving signal-to-noise ratio and directing model attention toward the interaction motifs critical for binding specificity. The model optimizes topological objectives and at inference drives pose refinement through a differentiable hybrid gradient field integrating predicted NCR and AutoDock Vina energetics. Extensive evaluation on the CASF-2016 benchmark and the 3D-DISCO cross-docking data set demonstrates consistently high performance: a Top-1 docking success rate of 94.7%, a 1% Enrichment Factor of 21.21 in virtual screening, and a Top-1 cross-docking success rate of 34.8%. Mechanistic analysis reveals that NCR-guided optimization enables decoy escaping from local energy minima and drives the recovery of disrupted native interactions, confirming that NCR captures the physical determinants of binding rather than mere geometric proximity.
Zhenqiang Zhang, Zhihao Wang, Yang Liu et al.· Journal of Chemical Informat...· 0 citations
A hybrid quantum-classical approach for molecular docking is proposed leveraging the MVWCP formalism with a variational full-basis encoding (FBE) strategy, which enables efficient encoding of classical binary variables with Bloch sphere vectors and proves that a global minimizer of the FBE objective can always be chosen to be a pure product state.
Tianqi Chen, A. Mak, Jianguo Li et al.· 0 citations
Results indicate that the topology-aligned inductive bias is the active ingredient driving parameter efficiency at QM9 scale, with implications for matched-baseline benchmarking in quantum machine learning.
Electrostatic embedding improved every accuracy and correlation metric for TYK2 but performed comparably to the classical and mechanical-embedding baselines for CDK2, thrombin, p38 and JNK1, and standard single-molecule energy and charge benchmarks were not good predictors of this target-dependent outcome.