Abstract
PROteolysis TArgeting Chimeras (PROTACs) are hetero-bifunctional small molecules that can simultaneously recruit target proteins and E3 ligases to form a ternary complex, promoting target protein ubiquitination and degradation via the Ubiquitin-Proteasome System (UPS). PROTACs have gained increasing attention in recent years due to certain advantages over traditional therapeutic modalities and enabling targeting of previously "undruggable" proteins. To better understand the mechanism of PROTAC-induced Target Protein Degradation (TPD), several computational approaches have recently been developed to study and predict ternary complex formation. However, mounting evidence suggests that ubiquitination can also be a rate-limiting step in PROTAC-induced TPD. Here, we propose a structure-based computational approach to predict target protein ubiquitination induced by cereblon (CRBN)-based PROTACs by leveraging available structural information of the CRL4A ligase complex (CRBN/DDB1/CUL4A/Rbx1/NEDD8/E2/Ub). We generated ternary complex ensembles with Rosetta, modeled multiple CRL4A ligase complex conformations, and predicted ubiquitination efficiency by separating the ternary ensemble into productive and unproductive complexes based on the proximity of the ubiquitin to accessible lysines on the target protein. We validated our CRL4A ligase complex models with published ternary complex structures and additionally employed our modeling workflow to predict ubiquitination efficiencies and sites of a series of cyclin-dependent kinases (CDKs) after treatment with TL12-186, a pan-kinase PROTAC. Our predictions are consistent with CDK ubiquitination and site-directed mutagenesis of specific CDK lysine residues as measured using a NanoBRET ubiquitination assay in HEK293 cells. This work structurally links PROTAC-induced ternary formation and ubiquitination, representing an important step toward prediction of target "degradability."