Structural Characterization of Degrader-Induced Ternary Complexes Using Hydrogen-Deuterium Exchange Mass Spectrometry and Computational Modeling: Implications for Structure-Based Design

Examining DNA Structures with In-droplet Hydrogen/Deuterium Exchange Mass Spectrometry

Capillary vibrating sharp-edge spray ionization (cVSSI) combined with hydrogen/deuterium exchange-mass spectrometry (HDX-MS) has been utilized to characterize different solution-phase DNA conformers including DNA G-quadruplex topologies as well as triplex DNA and duplex DNA. In general, G-quadruplex DNA shows a wide range of protection of hydrogens extending from ~12% to ~21% deuterium incorporation. Additionally, the DNA sequences selected to represent parallel, antiparallel, and hybrid G-quadruplex topologies exhibit slight differences in deuterium uptake levels which appear to loosely relate to overall conformer stability. Notably, the exchange level for one of the hybrid sequence sub topologies of G-quadruplex DNA (24 TTG) is significantly different (compared with the others studied here) despite the DNA sequences being highly comparable. For the quadruplex-forming sequences, correlation analysis suggests protection of base hydrogens involved in tetrad hydrogen bonding. For duplex DNA ~19% deuterium incorporation is observed while only ~16% is observed for triplex DNA. This increased protection of hydrogens may be due to the added backbone scaffolding and Hoogsteen base pairing of the latter species. These experiments lay the groundwork for future studies aimed at determining the structural source of this protection as well as the applicability of the approach for ascertaining different oligonucleotide folds, co-existing conformations, and/or overall conformer flexibility.

Structural Basis of Conformational Dynamics in the PROTAC-Induced Protein Degradation

Pronounced conformational dynamics is unveiled upon analyzing multiple crystal structures of the same proteins recruited to the same E3 ligases by PROTACs, and yet, is largely permissive for targeted protein degradation due to the intrinsic mobility of E3 assemblies creating a large ubiquitylation zone. Mathematical modelling of ternary dynamics on ubiquitylation probability confirms the experimental finding that ternary complex rigidification need not correlate with enhanced protein degradation. Salt bridges are found to prevail in the PROTAC-induced ternary complexes, and may contribute to a positive cooperativity and prolonged half-life. The analysis highlights the importance of presenting lysines close to the active site of the E2 enzyme while constraining ternary dynamics in PROTAC design to achieve high degradation efficiency.

Proteomic approaches advancing targeted protein degradation

Targeted protein degradation (TPD) is an emerging modality for research and therapeutics. Most TPD approaches harness cellular ubiquitin-dependent proteolytic pathways. Proteolysis-targeting chimeras (PROTACs) and molecular glue (MG) degraders (MGDs) represent the most advanced TPD approaches, with some already used in clinical settings. Despite these advances, TPD still faces many challenges, pertaining to both the development of effective, selective, and tissue-penetrant degraders and understanding their mode of action. In this review, we focus on progress made in addressing these challenges. In particular, we discuss the utility and application of recent proteomic approaches as indispensable tools to enable insights into degrader development, including target engagement, degradation selectivity, efficacy, safety, and mode of action.

The 17th EFMC Short Course on Medicinal Chemistry on Small Molecule Protein Degraders

The 17th EFMC Short Course on Medicinal Chemistry took place April 23-26, 2023 in Oegstgeest, near Leiden in the Netherlands. It covered for the first time the exciting topic of Targeted Protein Degradation (full title: Small Molecule Protein Degraders: A New Opportunity for Drug Design and Development). The course was oversubscribed, with 35 attendees and 6 instructors mainly from Europe but also from the US and South Africa, and representing both industry and academia. This report summarizes the successful event, key lectures given and topics discussed.

The potential of cross-linking mass spectrometry in the development of protein-protein interaction modulators

Cross-linking mass spectrometry (XL-MS) can provide a wealth of information on endogenous protein-protein interaction (PPI) networks and protein binding interfaces. These features make XL-MS an attractive tool to support the development of PPI-targeting drugs. Though not yet widely used, applications of XL-MS to drug characterization are beginning to emerge. Here, we compare XL-MS to established structural proteomics methods in drug research, discuss the current state and remaining challenges of XL-MS technology, and provide a perspective on the future role XL-MS can play in drug development, with a particular emphasis on PPI modulators.

VHL-Modified PROteolysis TArgeting Chimeras (PROTACs) as a Strategy to Evade Metabolic Degradation in In Vitro Applications

PROteolysis TArgeting Chimeras (PROTACs) are tripartite molecules consisting of a linker connecting a ligand for a protein of interest to an E3 ligase recruiter, whose rationale relies on proteasome-based protein degradation. PROTACs have expanded as a therapeutic strategy to open new avenues for unmet medical needs. Leveraging our expertise, we undertook a series of in vitro experiments aimed at elucidating PROTAC metabolism. In particular, we focused on PROTACs recruiting the von Hippel-Lindau (VHL) E3 ligase. After high-resolution mass spectrometry measurements, a characteristic metabolite with mass reduction of 200 units was detected and successively confirmed as a product deriving from the cleavage of the VHL ligand moiety. Subsequently, we identified hepatic and extrahepatic prolyl endopeptidases as the main putative metabolic enzymes involved. Finally, we designed and synthesized analogs of the VHL ligands that we further exploited for the synthesis of novel VHL-directed PROTACs with an improved metabolic stability in in vitro applications.

Targeted Protein Degradation: Advances, Challenges, and Prospects for Computational Methods

The therapeutic approach of targeted protein degradation (TPD) is gaining momentum due to its potentially superior effects compared with protein inhibition. Recent advancements in the biotech and pharmaceutical sectors have led to the development of compounds that are currently in human trials, with some showing promising clinical results. However, the use of computational tools in TPD is still limited, as it has distinct characteristics compared with traditional computational drug design methods. TPD involves creating a ternary structure (protein-degrader-ligase) responsible for the biological function, such as ubiquitination and subsequent proteasomal degradation, which depends on the spatial orientation of the protein of interest (POI) relative to E2-loaded ubiquitin. Modeling this structure necessitates a unique blend of tools initially developed for small molecules (e.g., docking) and biologics (e.g., protein-protein interaction modeling). Additionally, degrader molecules, particularly heterobifunctional degraders, are generally larger than conventional small molecule drugs, leading to challenges in determining drug-like properties like solubility and permeability. Furthermore, the catalytic nature of TPD makes occupancy-based modeling insufficient. TPD consists of multiple interconnected yet distinct steps, such as POI binding, E3 ligase binding, ternary structure interactions, ubiquitination, and degradation, along with traditional small molecule properties. A comprehensive set of tools is needed to address the dynamic nature of the induced proximity ternary complex and its implications for ubiquitination. In this Perspective, we discuss the current state of computational tools for TPD. We start by describing the series of steps involved in the degradation process and the experimental methods used to characterize them. Then, we delve into a detailed analysis of the computational tools employed in TPD. We also present an integrative approach that has proven successful for degrader design and its impact on project decisions. Finally, we examine the future prospects of computational methods in TPD and the areas with the greatest potential for impact.

The role of loop dynamics in the prediction of ligand-protein binding enthalpy

The enthalpic and entropic components of ligand-protein binding free energy reflect the interactions and dynamics between ligand and protein. Despite decades of study, our understanding and hence our ability to predict these individual components remains poor. In recent years, there has been substantial effort and success in the prediction of relative and absolute binding free energies, but the prediction of the enthalpic (and entropic) contributions in biomolecular systems remains challenging. Indeed, it is not even clear what kind of performance in terms of accuracy could currently be obtained for such systems. It is, however, relatively straight-forward to compute the enthalpy of binding. We thus evaluated the performance of absolute enthalpy of binding calculations using molecular dynamics simulation for ten inhibitors against a member of the bromodomain family, BRD4-1, against isothermal titration calorimetry data. Initial calculations, with the AMBER force-field showed good agreement with experiment (R² = 0.60) and surprisingly good accuracy with an average of root-mean-square error (RMSE) = 2.49 kcal mol^-1. Of the ten predictions, three were obvious outliers that were all over-predicted compared to experiment. Analysis of various simulation factors, including parameterization, buffer concentration and conformational dynamics, revealed that the behaviour of a loop (the ZA loop on the periphery of the binding site) strongly dictates the enthalpic prediction. Consistent with previous observations, the loop exists in two distinct conformational states and by considering one or the other or both states, the prediction for the three outliers can be improved dramatically to the point where the R² = 0.95 and the accuracy in terms of RMSE improves to 0.90 kcal mol^-1. However, performance across force-fields is not consistent: if OPLS and CHARMM are used, different outliers are observed and the correlation with the ZA loop behaviour is not recapitulated, likely reflecting parameterization as a confounding problem. The results provide a benchmark standard for future study and comparison.

On Ternary Complex Stability in Protein Degradation: In Silico Molecular Glue Binding Affinity Calculations

Molecular glues are small molecules that simultaneously bind to two proteins, creating a chemically induced protein-protein interface. CELMoDs (cereblon E3 ligase modulators) are a class of molecular glues that promote recruitment of neosubstrate proteins to the E3 ubiquitin ligase cereblon (CRBN) for poly-Lys48-ubiquitination and proteasomal degradation. Ternary complex structures of clinical CELMoDs CC-885 and CC-90009 bound to CRBN and neosubstrate G1 to S phase transition protein 1 (GSPT1) have been experimentally determined. Although cellular degradation is a downstream event, dependent not only on the affinity of the glue CELMoD in the ternary complex, we test the applicability of established structure-based drug design principles to predict binding affinity of CELMoDs to the protein-protein neointerface and correlation to measured cellular degradation for the neosubstrates GSPT1 and zinc finger Aiolos (IKZF3). For a congeneric series of CELMoDs, which have a similar sequence of binding events and resultant binding modes, we conclude that well-established structure-based methods that measure in silico ternary complex stabilities can predict relative degradation potency by CELMoDs.

Breaking free from the crystal lattice: Structural biology in solution to study protein degraders

Structural biology offers a versatile arsenal of techniques and methods to investigate the structure and conformational dynamics of proteins and their assemblies. The growing field of targeted protein degradation centres on the premise of developing small molecules, termed degraders, to induce proximity between an E3 ligase and a protein of interest to be signalled for degradation. This new drug modality brings with it new opportunities and challenges to structural biologists. Here we discuss how several structural biology techniques, including nuclear magnetic resonance, cryo-electron microscopy, structural mass spectrometry and small angle scattering, have been explored to complement X-ray crystallography in studying degraders and their ternary complexes. Together the studies covered in this review make a case for the invaluable perspectives that integrative structural biology techniques in solution can bring to understanding ternary complexes and designing degraders.

A deep dive into degrader-induced protein-protein interfaces

Targeted protein degradation (TPD) relies on a comprehensive understanding of interfaces between hijacked E3 ligases and their substrates. In vitro techniques often do not capture the interaction dynamics. Recently, Hanzl et al. introduced deep mutational scanning (DMS) in combination with structural and biochemical approaches to identify residues crucial for degrader activity.

Functional E3 ligase hotspots and resistance mechanisms to small-molecule degraders

Targeted protein degradation is a novel pharmacology established by drugs that recruit target proteins to E3 ubiquitin ligases. Based on the structure of the degrader and the target, different E3 interfaces are critically involved, thus forming defined 'functional hotspots'. Understanding disruptive mutations in functional hotspots informs on the architecture of the assembly, and highlights residues susceptible to acquire resistance phenotypes. Here we employ haploid genetics to show that hotspot mutations cluster in substrate receptors of hijacked ligases, where mutation type and frequency correlate with gene essentiality. Intersection with deep mutational scanning revealed hotspots that are conserved or specific for chemically distinct degraders and targets. Biophysical and structural validation suggests that hotspot mutations frequently converge on altered ternary complex assembly. Moreover, we validated hotspots mutated in patients that relapse from degrader treatment. In sum, we present a fast and widely accessible methodology to characterize small-molecule degraders and associated resistance mechanisms.

Crystallization of VHL-based PROTAC-induced ternary complexes

X-ray crystal structures of PROTAC-induced ternary complexes provide invaluable insights into the critical species underpinning PROTAC mode of action, explain protein degradation selectivity profiles, and can guide rational degrader design. Nevertheless, crystallization of the ternary complexes formed by PROTACs remains an important bottleneck in employing this method. This is mainly due to the potential flexibility and heterogeneity that is inherent to a non-native protein-protein complex mediated by a small molecule, which together can hamper crystallization of the desired species. To overcome this limitation, selecting PROTAC compounds that enable the formation of stable, high-affinity and preferably cooperative ternary complexes in stoichiometric amount is, in our experience, critical to the success of co-crystallization studies. In this chapter, examples of stable PROTAC-mediated ternary complexes are illustrated. Learnings from biophysical & biochemical data are used as a guideline in achieving the highest "crystallizability" of ternary complexes. A case study of VHL-based SMARCA2 PROTAC degrader ternary complex crystallization is described. The procedure includes over-expression and purification of the E3 ligase and target protein, forming (and sometimes isolating) the ternary complex, and crystallizing it. The protocols can be applied for other combinations of E3 ligase, PROTAC and target protein.

Exploring PROTAC Cooperativity with Coarse-Grained Alchemical Methods

Proteolysis targeting chimera (PROTAC) is a novel drug modality that facilitates the degradation of a target protein by inducing proximity with an E3 ligase. In this work, we present a new computational framework to model the cooperativity between PROTAC-E3 binding and PROTAC-target binding principally through protein-protein interactions (PPIs) induced by the PROTAC. Due to the scarcity and low resolution of experimental measurements, the physical and chemical drivers of these non-native PPIs remain to be elucidated. We develop a coarse-grained (CG) approach to model interactions in the target-PROTAC-E3 complexes, which enables converged thermodynamic estimations using alchemical free energy calculation methods despite an unconventional scale of perturbations. With minimal parametrization, we successfully capture fundamental principles of cooperativity, including the optimality of intermediate PROTAC linker lengths that originates from configurational entropy. We qualitatively characterize the dependency of cooperativity on PROTAC linker lengths and protein charges and shapes. Minimal inclusion of sequence- and conformation-specific features in our current force field, however, limits quantitative modeling to reproduce experimental measurements, but further development of the CG model may allow for efficient computational screening to optimize PROTAC cooperativity.

Investigation into the Use of Encorafenib to Develop Potential PROTACs Directed against BRAFV600E Protein

BRAF is a serine/threonine kinase frequently mutated in human cancers. BRAF^V600E mutated protein is targeted through the use of kinase inhibitors which are approved for the treatment of melanoma; however, their long-term efficacy is hampered by resistance mechanisms. The PROTAC-induced degradation of BRAF^V600E has been proposed as an alternative strategy to avoid the onset of resistance. In this study, we designed a series of compounds where the BRAF kinase inhibitor encorafenib was conjugated to pomalidomide through different linkers. The synthesized compounds maintained their ability to inhibit the kinase activity of mutated BRAF with IC₅₀ values in the 40-88 nM range. Selected compounds inhibited BRAF^V600E signaling and cellular proliferation of A375 and Colo205 tumor cell lines. Compounds 10 and 11, the most active of the series, were not able to induce degradation of mutated BRAF. Docking and molecular dynamic studies, conducted in comparison with the efficient BRAF degrader P5B, suggest that a different orientation of the linker bearing the pomalidomide substructure, together with a decreased mobility of the solvent-exposed part of the conjugates, could explain this behavior.

Investigating how intrinsically disordered regions contribute to protein function using HDX-MS

A large amount of the human proteome is composed of highly dynamic regions that do not adopt a single static conformation. These regions are defined as intrinsically disordered, and they are found in a third of all eukaryotic proteins. They play instrumental roles in many aspects of protein signaling, but can be challenging to characterize by biophysical methods. Intriguingly, many of these regions can adopt stable secondary structure upon interaction with a variety of binding partners, including proteins, lipids, and ligands. This review will discuss the application of Hydrogen-deuterium exchange mass spectrometry (HDX-MS) as a powerful biophysical tool that is particularly well suited for structural and functional characterization of intrinsically disordered regions in proteins. A focus will be on the theory of hydrogen exchange, and its practical application to identify disordered regions, as well as characterize how they participate in protein–protein and protein–membrane interfaces. A particular emphasis will be on how HDX-MS data can be presented specifically tailored for analysis of intrinsically disordered regions, as well as the technical aspects that are critical to consider when designing HDX-MS experiments for proteins containing intrinsically disordered regions.

Chasing molecular glue degraders: screening approaches

Protein-protein interactions (PPIs) govern all biological processes. Some small molecules modulate PPIs through induced protein proximity. In particular, molecular glue degraders are monovalent compounds that orchestrate interactions between a target protein and an E3 ubiquitin ligase, prompting the proteasomal degradation of the former. This and other pharmacological strategies of targeted protein degradation (e.g. proteolysis-targeting chimeras - PROTACs) overcome some limitations of traditional occupancy-based therapeutics. Here, we provide an overview of the "molecular glue" concept, with a special focus on natural and synthetic inducers of proximity to E3s. We then briefly highlight the serendipitous discoveries of some clinical and preclinical molecular glue degraders, and discuss the first examples of intentional discoveries. Specifically, we outline the different screening strategies reported in this rapidly evolving arena and our thoughts on future perspectives. By mastering the ability to influence PPIs, molecular glue degraders can induce the degradation of unligandable proteins, thus providing an exciting path forward to broaden the targetable proteome.

E3 ligase ligand chemistries: from building blocks to protein degraders

In recent years, proteolysis-targeting chimeras (PROTACs), capable of achieving targeted protein degradation, have proven their great therapeutic potential and usefulness as molecular biology tools. These heterobifunctional compounds are comprised of a protein-targeting ligand, an appropriate linker, and a ligand binding to the E3 ligase of choice. A successful PROTAC induces the formation of a ternary complex, leading to the E3 ligase-mediated ubiquitination of the targeted protein and its proteasomal degradation. In over 20 years since the concept was first demonstrated, the field has grown substantially, mainly due to the advancements in the discovery of non-peptidic E3 ligase ligands. Development of small-molecule E3 binders with favourable physicochemical profiles aided the design of PROTACs, which are known for breaking the rules of established guidelines for discovering small molecules. Synthetic accessibility of the ligands and numerous successful applications led to the prevalent use of cereblon and von Hippel-Lindau as the hijacked E3 ligase. However, the pool of over 600 human E3 ligases is full of untapped potential, which is why expanding the artillery of E3 ligands could contribute to broadening the scope of targeted protein degradation. In this comprehensive review, we focus on the chemistry aspect of the PROTAC design process by providing an overview of liganded E3 ligases, their chemistries, appropriate derivatisation, and synthetic approaches towards their incorporation into heterobifunctional degraders. By covering syntheses of both established and underexploited E3 ligases, this review can serve as a chemistry blueprint for PROTAC researchers during their future ventures into the complex field of targeted protein degradation.

Modeling the CRL4A ligase complex to predict target protein ubiquitination induced by cereblon-recruiting PROTACs

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.”

Computational Modeling of Molecular Structures Guided by Hydrogen-Exchange Data

Data produced by hydrogen-exchange monitoring experiments have been used in structural studies of molecules for several decades. Despite uncertainties about the structural determinants of hydrogen exchange itself, such data have successfully helped guide the structural modeling of challenging molecular systems, such as membrane proteins or large macromolecular complexes. As hydrogen-exchange monitoring provides information on the dynamics of molecules in solution, it can complement other experimental techniques in so-called integrative modeling approaches. However, hydrogen-exchange data have often only been used to qualitatively assess molecular structures produced by computational modeling tools. In this paper, we look beyond qualitative approaches and survey the various paradigms under which hydrogen-exchange data have been used to quantitatively guide the computational modeling of molecular structures. Although numerous prediction models have been proposed to link molecular structure and hydrogen exchange, none of them has been widely accepted by the structural biology community. Here, we present as many hydrogen-exchange prediction models as we could find in the literature, with the aim of providing the first exhaustive list of its kind. From purely structure-based models to so-called fractional-population models or knowledge-based models, the field is quite vast. We aspire for this paper to become a resource for practitioners to gain a broader perspective on the field and guide research toward the definition of better prediction models. This will eventually improve synergies between hydrogen-exchange monitoring and molecular modeling.