Rational Prediction of PROTAC-Compatible Protein-Protein Interfaces by Molecular Docking

In silico modeling of targeted protein degradation

Targeted protein degradation (TPD) techniques, particularly proteolysis-targeting chimeras (PROTAC) and molecular glue degraders (MGD), have offered novel strategies in drug discovery. With rapid advancement of computer-aided drug design (CADD) and artificial intelligence-driven drug discovery (AIDD) in the biomedical field, a major focus has become how to effectively integrate these technologies into the TPD drug discovery pipeline to accelerate development, shorten timelines, and reduce costs. Currently, the main research directions for applying CADD and AIDD in TPD include: 1) ternary complex modeling; 2) linker generation; 3) strategies to predict degrader targets, activities and ADME/T properties; 4) In silico degrader design and discovery. Models developed in these areas play a crucial role in target identification, drug design, and optimization at various stages of the discovery process. However, the limited size and quality of datasets related to TPD present challenges, leaving room for further improvement in these models. TPD involves the complex ubiquitin-proteasome system, with numerous factors influencing outcomes. Most current models adopt a static perspective to interpret and predict relevant tasks. In the future, it may be necessary to shift toward dynamic approaches that better capture the intricate relationships among these components. Furthermore, incorporating new and diverse chemical spaces will enhance the precision design and application of TPD agents.

Ternary Complex Modeling, Induced Fit Docking and Molecular Dynamics Simulations as a Successful Approach for the Design of VHL-Mediated PROTACs Targeting the Kinase FLT3

Proteolysis targeting chimeras (PROTACs) have proven to be a novel approach for the degradation of disease-causing proteins in drug discovery. One of the E3 ligases for which efficient PROTACs have been described is the Von Hippel-Lindau factor (VHL). However, the development of PROTACs has so far often relied on a minimum of computational tools, so that it is mostly based on a trial-and-error process. Therefore, there is a great need for resource- and time-efficient structure-based or computational approaches to streamline PROTAC design. In this study, we present a combined computational approach that integrates static ternary complex formation, induced-fit docking, and molecular dynamics (MD) simulations. Our methodology was tested using four experimentally derived ternary complex structures of VHL PROTACs, reported for BRD4, SMARCA2, FAK, and WEE1. In addition, we applied the validated approach to model a recently in-house developed FLT3-targeted PROTAC (MA49). The results show that static ternary models generated with a protein-protein docking method implemented in the software MOE have a high predictive power for reproducing the experimental 3D structures. The induced-fit docking of different active PROTACs to their respective models showed the reliability of this model for the development of new VHL-mediated degraders. In particular, the induced-fit docking was sensitive to structural changes in the PROTACs, as evidenced by the failed binding modes of the PROTAC negative controls. Furthermore, MD simulations confirmed the stability of the generated complexes and emphasized the importance of dynamic studies for understanding the relationship between PROTAC structure and function.

PROTAC-induced protein structural dynamics in targeted protein degradation

PROteolysis TArgeting Chimeras (PROTACs) are small molecules that induce target protein degradation via the ubiquitin-proteasome system. PROTACs recruit the target protein and E3 ligase; a critical first step is forming a ternary complex. However, while the formation of a ternary complex is crucial, it may not always guarantee successful protein degradation. The dynamics of the PROTAC-induced degradation complex play a key role in ubiquitination and subsequent degradation. In this study, we computationally modelled protein complex structures and dynamics associated with a series of PROTACs featuring different linkers to investigate why these PROTACs, all of which formed ternary complexes with Cereblon (CRBN) E3 ligase and the target protein bromodomain-containing protein 4 (BRD4BD1), exhibited varying degrees of degradation potency. We constructed the degradation machinery complexes with Culling-Ring Ligase 4A (CRL4A) E3 ligase scaffolds. Through atomistic molecular dynamics simulations, we illustrated how PROTAC-dependent protein dynamics facilitating the arrangement of surface lysine residues of BRD4BD1 into the catalytic pocket of E2/ubiquitin cascade for ubiquitination. Despite featuring identical warheads in this PROTAC series, the linkers were found to affect the residue-interaction networks, and thus governing the essential motions of the entire degradation machine for ubiquitination. These findings offer a structural dynamic perspective on ligand-induced protein degradation, providing insights to guide future PROTAC design endeavors.

PROTAC-induced Protein Structural Dynamics in Targeted Protein Degradation

PROteolysis TArgeting Chimeras (PROTACs) are small molecules that induce target protein degradation via the ubiquitin-proteasome system. PROTACs recruit the target protein and E3 ligase; a critical first step is forming a ternary complex. However, while the formation a ternary complex is crucial, it may not always guarantee successful protein degradation. The dynamics of the PROTAC-induced degradation complex play a key role in ubiquitination and subsequent degradation. In this study, we computationally modelled protein complex structures and dynamics associated with a series of PROTACs featuring different linkers to investigate why these PROTACs, all of which formed ternary complexes with Cereblon (CRBN) E3 ligase and the target protein bromodomain-containing protein 4 (BRD4 BD1 ), exhibited varying degrees of degradation potency. We constructed the degradation machinery complexes with Culling-Ring Ligase 4A (CRL4A) E3 ligase scaffolds. Through atomistic molecular dynamics simulations, we illustrated how PROTAC-dependent protein dynamics facilitating the arrangement of surface lysine residues of BRD4 BD1 into the catalytic pocket of E2/ubiquitin cascade for ubiquitination. Despite featuring identical warheads in this PROTAC series, the linkers were found to affect the residue-interaction networks, and thus governing the essential motions of the entire degradation machine for ubiquitination. These findings offer a structural dynamic perspective on ligand-induced protein degradation, providing insights to guide future PROTAC design endeavors.

Development of PROTACs using computational approaches

Proteolysis-targeting chimeras (PROTACs) are drugs designed to degrade target proteins via the ubiquitin-proteasome system. With the application of computational biology/chemistry technique in drug design, numerous computer-aided drug design and artificial intelligence (AI)-driven drug design (CADD/AIDD) methods have recently emerged to facilitate the development of PROTAC drugs. We systematically review the role of in silico tools in PROTAC drug design, emphasizing how computational software can model PROTAC action and structure, predict activity, and assist in molecule design. We also discuss current challenges in the rational design of PROTACs from an in silico perspective, such as deviations from small-molecule druggability and the limited availability of training data. We provide an overview of recent discoveries and emerging research in this field, and discuss their potential impact on PROTAC design strategies.

Advances in designing ternary complexes: Integrating in-silico and biochemical methods for PROTAC optimisation in target protein degradation

Target protein degradation (TPD) is an emerging approach to mitigate disease-causing proteins. TPD contains several strategies, and one of the strategies that gained immersive importance in recent times is Proteolysis Targeting Chimeras (PROTACs); the PROTACs recruit small molecules to induce the poly-ubiquitination of disease-causing protein by hijacking the ubiquitin-proteasome system (UPS) by bringing the E3 ligase and protein of interest (POI) into appropriate proximity. The steps involved in designing and evaluating the PROTACs remain critical in optimising the PROTACs to degrade the POI. It is observed that using in-silico and biochemical methods to study the ternary complexes (TCs) of the POI-PROTAC-E3 ligase is essential to understanding the structural activity, cooperativity, and stability of formed TCs. A better understanding of the above-mentioned leads to an appropriate rationale for designing the PROTACs targeting the disease-causing proteins. In this review, we tried to summarise the approaches used to design the ternary complexes, i.e., in-silico and in-vitro methods, to understand the behaviour of the PROTAC-induced ternary complexes.

Characterizing the Cooperative Effect of PROTAC Systems with End-Point Binding Free Energy Calculation

Proteolytic targeting chimeras (PROTACs), as an emerging type of drug, function by proximity-based modalities that narrow the distance between a target protein and the E3 ubiquitin ligase to facilitate the ubiquitination labeling of the target protein for degradation. Although it is evidenced that the cooperativity of the PROTAC ternary interaction is one of the key factors affecting the degradation rate of a target protein, PROTAC design utilizing this indicator is still challenging because of the complicated/flexible interactions in a target-PROTAC-E3 ternary system. Therefore, developing reliable and practicable computational methods is of great interest for PROTAC design. Hence, in this study, we investigate the feasibility of using the end-point binding free energy calculation method, represented by molecular mechanics/Poisson-Boltzmann (generalized-Born) surface area (MM/PB(GB)SA), for characterizing cooperativity (including the stabilization and hook effects) of the PROTAC systems. The result shows that MM/GBSA is a good predictor in characterizing these effects under a relatively long molecular dynamics adjustment (50-100 ns) and low dielectric constant (εin = 1), with the Pearson correlation coefficient (rp) > 0.5 and 0.6 for the stabilization and hook effect, respectively. This study provides a feasible strategy for characterizing the cooperativity of the PROTAC systems, facilitating the rational design of PROTAC molecules.

Revolutionizing Drug Targeting Strategies: Integrating Artificial Intelligence and Structure-Based Methods in PROTAC Development

PROteolysis TArgeting Chimera (PROTAC) is an emerging technology in chemical biology and drug discovery. This technique facilitates the complete removal of the target proteins that are "undruggable" or challenging to target through chemical molecules via the Ubiquitin-Proteasome System (UPS). PROTACs have been widely explored and outperformed not only in cancer but also in other diseases. During the past few decades, several academic institutes and pharma companies have poured more efforts into PROTAC-related technologies, setting the stage for several major degrader trial readouts in clinical phases. Despite their promising results, the formation of robust ternary orientation, off-target activity, poor permeability, and binding affinity are some of the limitations that hinder their development. Recent advancements in computational technologies have facilitated progress in the development of PROTACs. Researchers have been able to utilize these technologies to explore a wider range of E3 ligases and optimize linkers, thereby gaining a better understanding of the effectiveness and safety of PROTACs in clinical settings. In this review, we briefly explore the computational strategies reported to date for the formation of PROTAC components and discuss the key challenges and opportunities for further research in this area.