AI-Assisted Identification of a Putative Allosteric Ligand Targeting the CDK4/Cyclin D1 Protein-Protein Interface

Background/Objectives: First-generation CDK4/6 inhibitors (palbociclib, ribociclib, abemaciclib) target the conserved ATP-binding pocket of CDK4 and, despite clinical success, are limited by acquired resistance and insufficient exploration of alternative regulatory sites. This study aimed to identify a putative allosteric small-molecule candidate at the CDK4 αE-helix-Cyclin D1 α1-helix protein-protein interaction (PPI) interface within the CDK4/Cyclin D1/p21 ternary complex using RapidFunnel-AI, a decision-interpretable virtual-screening pipeline. Methods: Starting from 50,000 ChEMBL 33 molecules, the pipeline sequentially applied a Q-Fold/RapidFunnel topological Tanimoto scan based on clinical CDK4/6 inhibitor motifs, fragment-level electronic-property enrichment, ADMET/PAINS filtering, dry Vina-GPU docking, hydration-mediated AutoDock-GPU (Version 1.6) docking, explicit-solvent molecular dynamics, contact-retention analysis, and MM-GBSA energy decomposition. The Q-Fold Thermo-Core surrogate model provided fragment-level enrichment, predicting the HOMO-LUMO gap (R² = 0.93) and isotropic polarizability (R² = 0.98) on QM9. Candidate selection did not rely on the lowest docking or MM-GBSA score alone, but on pose persistence, contact continuity, and energy-component consistency. Results: The workflow reduced the initial library to 43 topologically prioritized candidates, 25 ADMET/PAINS-filtered ligands, and 9 docking-derived complexes for MD validation. Ligand_020 emerged as the only candidate that preserved a persistent binding mode at Site 2 during a 500 ns simulation-an interface engagement reproduced across three independent 500 ns replicates with no full dissociation in any replicate-with a protein Cα RMSD of 2.88 ± 0.32 Å, a ligand heavy-atom RMSD of 3.56 ± 0.28 Å, and a van der Waals-dominated MM-GBSA profile (ΔGbind = -28.23 ± 3.57 kcal/mol). In contrast, palbociclib and ribociclib, forcibly placed at Site 2 as negative controls, lost most initial contacts within 5 ns and tended to detach despite more favorable MM-GBSA values. Conclusions: These results suggest that single-score docking or MM-GBSA ranking can generate false positives at shallow PPI interfaces. By integrating AI-assisted prioritization, multipocket docking, explicit-solvent MD, contact-retention analysis, and energy-component consistency, RapidFunnel-AI nominated Ligand_020 as an experimentally testable putative allosteric hit targeting the CDK4/Cyclin D1 interface, offering a reusable platform for PPI-focused oncological drug discovery.