Importance of Three-Body Problems and Protein-Protein Interactions in Proteolysis-Targeting Chimera Modeling: Insights from Molecular Dynamics Simulations

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 K-d values. The modeling of a unique property, namely, the a 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.

Automated summary

This study explores the use of molecular dynamics simulation and molecular mechanics combined with the generalized Born and surface area continuum solvation (MM/GBSA) method to solve the three-body problem in PROTAC modeling. The results validate the accuracy of the approach and demonstrate the qualitative and quantitative importance of PROTAC-induced protein-protein interactions in its modeling.