Abstract
The rational design of novel drug candidates presents a formidable challenge in modern drug discovery. Proteolysis-targeting chimeras (PROTACs) drug design is particularly demanding due to their limited crystal structure availability and design of a viable small molecule to bridge the protein of interest (POI) and ubiquitin-protein ligase (E3). An integrated approach that combines superimposition techniques and deep neural networks is demonstrated in this study to leverage the power of deep learning and structural biology to generate structurally diverse molecules with enhanced binding affinities. The superimposition technique ensures the congruence of initial and new protein-ligand pairs, which are evaluated via subsequent comprehensive screening using the root-mean-square deviation (RMSD), binding free energy (BFE), and buried solvent-accessible surface area (SASA). The final candidates are subjected to the incorporation of molecular dynamics (MD) and free energy perturbation (FEP) simulations to provide a quantitative evaluation of relative binding energies, reinforcing the efficacy and reliability of the generated molecules. The outcomes of the generated novel PROTACs molecules exhibit comparable structural attributes while demonstrating superior binding affinities within the binding pockets when contrasted with those of the established cocrystal ternary complexes. To enhance the generalizability of the workflow, we chose the ternary structure of the cellular inhibitor of apoptosis protein 1 (cIAP1) and Bruton's Tyrosine Kinase (BTK) for validating the chemical properties generated from the processes. The new linker molecules additionally showed superior affinity from the simulations. In summary, this methodology serves as an effective workflow to align computational predictions with current limitations, thereby introducing a novel paradigm in AI-driven drug design.