LAUSR

Proteolysis-targeting chimeras (PROTACs) are a class of bifunctional molecules that can induce the ubiquitin degradation of its target protein by hijacking the E3 ligase to form a target protein-PROTAC-E3 ligase ternary complex. Its underlying principle has inspired the development of a wide range of protein degraders that are similar to or beyond PROTACs in recent years. The formation of the ternary complexes is the key to the success of PROTAC-induced protein degradation. Nevertheless, the lack of effective ternary complex modeling techniques has limited the application of computer-aided drug discovery tools to this emerging and fast developing new land in drug industry. Thus, in this study, we explored the application of the more physically sound molecular dynamics simulation and the molecular mechanics combined with the generalized Born and surface area continuum solvation (MM/GBSA) method to solve the underlying three-body problem in PROTAC modeling. We first verified the accuracy of our approach using a series of known Brd4 BD2 degraders. The calculated binding energy showed a good correlation with the experimental Kd values. The modeling of a unique property, namely, the α value, for PROTACs was also first and accurately performed to our best knowledge. The results also demonstrated the importance of PROTAC-induced protein-protein interactions in its modeling, either qualitatively or quantitatively. Finally, by standing on the success of earlier docking-based approaches, our protocol was also applied as a rescoring function in pose prediction. The results showed a notable improvement in reranking the initial poses generated from a modified Rosetta method, which was reportedly one of the best among a handful of PROTAC modeling approaches available in this field. We hope this work could provide a practical protocol and more insights to study the binding and the design of PROTACs and other protein degraders.