In silico Modeling of PROTAC-Mediated Ternary Complexes

Feb 4, 2019 - In silico Modeling of PROTAC-Mediated Ternary Complexes: Validation and Application. Michael L Drummond and Chris Willaims. J. Chem...
2 downloads 0 Views 685KB Size
Subscriber access provided by Iowa State University | Library

Pharmaceutical Modeling

In silico Modeling of PROTAC-Mediated Ternary Complexes: Validation and Application Michael L Drummond, and Chris Willaims J. Chem. Inf. Model., Just Accepted Manuscript • DOI: 10.1021/acs.jcim.8b00872 • Publication Date (Web): 04 Feb 2019 Downloaded from http://pubs.acs.org on February 5, 2019

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

In silico Modeling of PROTAC-Mediated Ternary Complexes: Validation and Application Michael L. Drummond* and Christopher I. Williams Chemical Computing Group, Montreal, Quebec, H3A 2R7, Canada

ABSTRACT. In this work, four methods are described and validated for generating in silico ensembles of PROTAC-mediated ternary complexes. Filters based on characteristics of the proposed ternary complexes are developed to identify those that resemble known crystal structures. We then show how to use these modeling techniques a priori to discriminate the PROTACmediated degradation behavior of a mutant protein vs. its wild type, of three closely related targets, and among three different PROTAC molecules.

INTRODUCTION. Over the past decade, intense drug discovery research has focused on a class of bifunctional molecules commonly known as proteolysis-targeting chimeras, or PROTACs. These molecules, a subclass of the larger categories of molecular glues1 and chemical inducers of proximity,2 function by non-covalently linking a targeted protein to the substrate adaptor protein of an E3 ligase,3 which naturally functions to polyubiquitinate proteins in vivo; after being so tagged, the labeled target protein is then sent to the Ubiquitin Proteasome System (UPS), where it is degraded by the 26S proteasome.4 A PROTAC molecule consists of three parts: a moiety that

ACS Paragon Plus Environment

1

Journal of Chemical Information and Modeling 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 37

binds to the target protein, a (generally flexible) linker, and a second moiety that binds and recruits the E3 ligase. While originally PROTACs used peptidic E3 ligase substrates,5 in 2010 Itoh et al. reported6 PROTACs constructed entirely from small molecules. This finding reignited interest in PROTAC research, as can be seen from numerous reviews,7–13 as well as a recent issue14 of J. Med. Chem. dedicated entirely to protein degradation as a therapeutic tool orthogonal to classic small molecule-driven inhibition. The reported advantages of the PROTAC modality are many, with potentially game-changing significance. Foremost is the potential to target proteins traditionally regarded as “undruggable,” such those that can only be targeted with weak binders15 or that only possess non-enzymatic functions (e.g., scaffolds16,17). PROTACs that act substoichiometrically – i.e., capable of labeling multiple proteins for degradation in a catalytic fashion – have also been reported,18 as have PROTACs that successfully cause degradation despite the occurrence of mutations known to be deleterious to small molecule inhibition.19 PROTAC-mediated degradation has also been shown to be effective even when traditional inhibition is not.20 Moreover, the multifaceted nature of PROTAC molecules often imparts a selectivity not observed in the original target binding moiety,21 due to the structure of the ternary (three-body) target-PROTAC-ligase complexes they engender. Despite these many advantages, at the time of this writing no PROTACs have yet proceeded to the clinic. Although issues pertaining to PROTAC properties, such as their bioavailability and metabolic stability, are certainly substantial (due to their beyond-rule-of-five nature),22 the initial discovery phase for these molecules is in and of itself particularly challenging. Multiple reports have shown that each of the constituent parts of a PROTAC can impact its efficacy. That is, not only do the identity of the E3 and target binding moieties impact activity, but so too does the linker,

ACS Paragon Plus Environment

2

Page 3 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

in particular with respect to its length, its chemical nature, its flexibility, and how and where it is covalently bonded to the two binding moieties.23–28 Moreover, analysis of the first crystal structure29 of a target-ligase-PROTAC ternary complex (of only two total reports to date30) revealed the existence of all combinations of possible interactions: between the linker and both proteins, as well as between the target binding moiety and the E3 ligase protein and between the ligase binding moiety and the target protein. Thus, optimization of any individual PROTAC component may impact the behavior of the other PROTAC components, ultimately leading to complicated SAR.24 Although further experimental work will surely add clarity, computational modeling can also be used to provide much needed insight into these systems. Indeed, computational tools have already been applied to PROTAC design, albeit at varying levels of rigor. The most common practice to date is to simply dock a PROTAC into the pocket of one of its two cognate proteins, select a docking pose with an extended conformation, and then superpose the second protein onto its corresponding solvent-exposed binder.19,31,32 Other techniques used to date include those based on protein-protein docking,30,33 molecular dynamics,15,29 and an analysis of the conformational ensemble of the PROTAC absent its proteins.34 Perhaps the most extensive work to date was detailed by researchers at Pfizer,35 who developed a computational workflow based on evaluating the steric compatibility between a PROTAC conformational ensemble and post hoc appended proteins; importantly, this work also utilized a simple scoring function – a count of the number of non-clashing poses – to judge the ability of a PROTAC to facilitate productive ternary complexes. Although PROTAC design has certainly already benefitted from these computational approaches, none of the published methods thus far have been demonstrated to be applicable across multiple targets, using different E3 ligases, different binding moieties, or different linkers. In this

ACS Paragon Plus Environment

3

Journal of Chemical Information and Modeling 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 37

work, we describe a suite of in silico tools suitable for addressing the many varying facets inherent to PROTAC development, validated across multiple design scenarios. As will be shown, these tools not only successfully generate an ensemble of productive ternary complexes, but more importantly score the likelihood that a PROTAC can lead to effective protein degradation. MATERIALS AND METHODS. Due to the aforementioned tendency for the different components of a ternary complex to interact with one another,29 it is a priori unclear how a robust, diverse ensemble of putative ternary complexes should be generated. Conceptually, at one extreme, PROTAC conformations could be sampled entirely separately from their binding proteins; at the other extreme, all three pieces of a ternary complex would be present during sampling. In this work, rather than prejudge which approach is most appropriate, instead we present results for four different computational Methods covering these various scenarios. Full computational details are presented below, and are summarized graphically in Figure 1. Due to differences among the Methods, the details of scoring, i.e., the identification of potentially useful poses amidst the complete ensembles, differ across the four Methods, as will be discussed in the Results section. All Methods were written in the Scientific Vector Language (SVL) integrated into the MOE modeling package,36 and are freely available upon request.

ACS Paragon Plus Environment

4

Page 5 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

Figure 1. Graphical summary of the four Methods detailed in this paper. See text for full details. In brief, the four Methods can be summarized as follows: 

Method 1: The entire ternary complex is sampled all at once



Method 2: PROTAC conformations are sampled independently, followed by post

hoc addition of rigid-body proteins 

Method 3: The PROTAC is sampled in the context of one of the proteins, with the

second added afterwards 

Method 4: PROTAC conformations are sampled independently of the proteins (and,

indeed, a conformational database from Method 2 can be used, or vice versa), but possible ligase-target arrangements are provided via protein-protein docking. Method 1. For Method 1, it is first necessary to provide a starting orientation of each proteinbinding moiety in its respective protein, such as would be produced via docking or (ideally) a cocrystallized protein-binder complex; this binding moiety is referred to as an “anchor” below.

ACS Paragon Plus Environment

5

Journal of Chemical Information and Modeling 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 37

With this input data, a starting conformation of the PROTAC is constructed first by attaching each anchor group (fixed in its protein-bound conformation) to the initial, user-provided conformation of the PROTAC’s linker group. The proteins are then positioned around their respective anchor groups so as to match the initially provided protein-anchor group orientations, thereby forming an initial configuration of the entire ternary complex. This initial configuration often has significant protein-ligand and protein-protein clashes, due to the separate starting environments for the two protein/anchor group complexes and the initial linker conformation. In order to minimize the production of sterically clashing ternary complexes, the PROTAC conformation is automatically adjusted to adopt an extended conformation. Rotatable bonds of the PROTAC are modified to values that minimize volume overlap, and then restricted to these values subsequent to PROTAC conformational sampling. These restricted rotatable bonds are shown in black in Figure 1. This procedure is repeated for the next adjacent rotatable bond. If the search encounters a rotatable bond where there is no overlap between the dihedral terminal atom and the protein at any dihedral angle, the dihedral angle is set to 180°, so as to generate a fully extended ligand conformation. These non-restricted rotatable bonds are shown in green in Figure 1. After generating the extended conformation of the PROTAC, additional putative conformations of the ternary complex are generated by holding one protein and its anchor group fixed in space while sampling the unrestricted rotatable bonds of the linker (green bonds in Figure 1). Only one rotatable bond, randomly chosen, is sampled at each individual step of the search, and the second protein is moved in concert with the motion of the anchor group connected to the sampled linker rotatable bond. The generated conformation of the ternary complex is then scored using a variety of metrics (see below), and the conformation is either retained or rejected based on user-specified

ACS Paragon Plus Environment

6

Page 7 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

thresholds in ligand strain energy and protein-protein volume overlap. In Monte Carlo style, if a conformation is rejected, the search reverts to the last accepted conformation and then proceeds by sampling a new, randomly chosen rotatable bond. If instead a generated conformation is accepted, the sampling continues starting from this acceptable conformation. The search terminates when either the total number of sampling steps exceeds a user-defined maximum threshold, or if a new conformation is not produced after a user-specified maximum number of contiguous failures. Method 2. As in Method 1, it is necessary to provide two protein-ligand complexes, such as produced via X-ray crystallography or through docking, consisting of each protein with its binding moiety. Additionally, a core for each ligand is specified via SMILES (or through selection of the atoms using the MOE GUI, from which a SMILES string can be automatically extracted). The full PROTAC must also be provided via a MOE database (MDB), which can contain multiple PROTACs and can optionally contain pre-generated conformations for each PROTAC. If conformations are not provided, they are generated on-the-fly using any of the conformational generating methods in MOE (LowModeMD, Conformation Import, Stochastic, or Systematic search). The Potential settings (forcefield, implicit solvation model, dielectric settings, etc.) are taken from the main MOE settings. After a PROTAC conformational ensemble is either generated on-the-fly or read in, the two original protein/ligand complexes are then superposed onto their respective binding moieties for each conformation. During the conformational search, the userspecified core may deform from its cocrystallized position. A Core RMSD threshold is set (default: 2.5 Å, determined through visual inspection), and if either binding moiety in a free PROTAC conformation differs from the originally provided cocrystallized binder geometry by more than this threshold (after rigid body superposition), the conformation is rejected. Benchmarking data

ACS Paragon Plus Environment

7

Journal of Chemical Information and Modeling 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 37

(not shown) suggests that a core should be generally small and rigid, as larger cores result in more rejected poses and lower hit rates. Because the PROTAC sampling occurs in the absence of any proteins, there is often significant overlap after the proteins are reintroduced (if, say, a PROTAC conformation is especially collapsed, for example). In the interface, it is possible to pre-filter dramatically clashing conformations, based on all-atom protein-protein volume overlaps (default: 10000 Å3) and backbone-backbone volume overlaps (default: 2500 Å3). These settings are applied before the subsequent filter S_all described later in the main text of this article. Any generated ternary complexes that pass the specified filters are written to a MOE database, along with automatically calculated information that may prove useful for identifying competent ternary complexes (energetic criteria, PROTAC strain, patch-based descriptors, volume overlaps, etc.). Method 3. In Method 3, only one of the proteins is included in the conformational sampling phase; including the smaller of the two is more computationally efficient. Additionally, the entirety of the PROTAC molecule, with one end bound to the pocket of the included protein, is required. As in Method 2, the user defines two (generally rigid) cores. The core belonging to the binding moiety of the protein included in the sampling is tethered during the conformational sampling, so that it cannot dissociate from the binding pocket. The core that corresponds to the protein not included in the conformational search is instead defined as a rigid body, to prevent the deformation of the binding moiety during conformational sampling that sometimes occurred in Method 2. Additionally, protein residues near the first core are tethered during the conformational search, to afford some flexibility to the protein-ligand binding pocket. The remainder of the PROTAC molecule, i.e., the (generally flexible) linker and any non-rigid portions of either binding moiety, are unconstrained during the conformational search.

ACS Paragon Plus Environment

8

Page 9 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

The conformational search in Method 3 uses LowModeMD exclusively – settings for this algorithm are exposed to the user, but by default are set to be expansive, so as to robustly sample potential protein-PROTAC conformations. After the search completes, the second protein is superposed onto its respective binder, which was held rigid during the conformational search. As with Method 2, poses can be rejected at this phase based on protein-protein all atom or backbonebackbone overlaps. All poses that pass these filtering criteria are written to an output MDB, where patch-based descriptors are then calculated. Method 4. As with Method 2, the input required for Method 4 is two protein-ligand structures, where the ligand is the corresponding PROTAC binder moiety. Also similar to Method 2, cores for each binding moiety are defined by the user, and will be used as described later. As with the preceding Methods, the ideal starting inputs for Method 4 are two separate protein-binder crystal structures, although future work will demonstrate that useful predictions can still be generated if small-molecule docking is used to place the actual target binding moiety of the PROTAC into the protein pocket taken from a crystal structure of the target cocrystallized with a similar (but not necessarily identical) target binder. Once the input is provided, protein-protein docking of these two protein-ligand complexes is the first phase of Method 4. Protein sidechains are automatically repacked during this procedure, and putative arrangements of the two protein-ligand complexes are also subjected to rigid-body minimization. It is also possible to specify a database with pregenerated docking poses, such as can be produced with an earlier run of Method 4, or as a standalone MOE calculation, or as imported from a third-party software. It was decided early in the construction of Method 4 to dock two protein-ligand complexes (rather than two apo proteins), as it is known that certain PROTAC binding moieties, particularly the immunomodulatory imide drugs (IMiDs) such as lenalidomide

ACS Paragon Plus Environment

9

Journal of Chemical Information and Modeling 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 37

and pomalidomide,37,38 “fill in” gaps on the surface of cereblon, thereby expanding its degradation profile beyond its native substrate.1 As will be discussed in the Results section below, multiple protein-protein docking runs can optionally be performed, where each simulation is biased to match up a hydrophobic protein patch near (i.e., within 10 Å) the PROTAC binding pocket against a similar hydrophobic patch on the other protein, in an all-against-all fashion. A similar approach has been used to generate initial positions prior to generating trajectories, in a molecular dynamics context.15 If this biasing option is selected, all docking poses are collated, and then poses that are within 4 Å Cα RMSD of an already observed pose are discarded as duplicates. After the final protein-protein docking database is constructed, protein patch-based descriptors (see Results) are (optionally, but on by default and highly recommended) calculated. Only protein-protein docked poses with patch-patch contact surface areas above those specified by the user are considered for superposition against the PROTAC conformational ensemble (see below). The next phase involves generation of the conformational ensemble of the specified PROTAC molecule, which proceeds in a fashion identical to that described above for Method 2. (Indeed, the PROTAC conformational databases generated with either Method are interchangeable). After this database is generated, an initial check is performed to evaluate whether an individual proteinprotein docked pose is roughly compatible with an individual PROTAC conformation. The pocketpocket distance for all protein-protein docked poses is calculated, taken as the distance between the centroid of the two user-specified cores. Similarly, for each PROTAC conformation, the centroid-centroid distance, again based on the provided core definitions, is calculated. If the interpocket distance for a protein-protein pose is within 5 Å of the intraligand distance, then the two separate conformations are deemed roughly compatible, and they are combined. This pre-

ACS Paragon Plus Environment

10

Page 11 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

screen prevents grossly distinct conformations from being combined, a necessity given that, if unscreened all-against-all matching were attempted, ~10,000 PROTAC conformations would be combined with 100-1000 protein-protein docking poses, resulting in ~10,000,000 putative combinations, the vast majority of which would be incompatible. If a PROTAC conformation is deemed compatible with a protein-protein docking pose, then the two are combined in the following manner, graphically summarized in Figure 2. First, the PROTAC conformation is brought into the frame of reference of one of the cores in the proteinprotein docked pose by rigidly moving the entire PROTAC, based on superposing it onto its match in the protein-protein docked pose. Next, the other end of the PROTAC conformation, as defined by its user-specified core, is rigidly superposed onto its counterpart in the protein-protein docked pose, but the remainder of the PROTAC is left in its position defined by the first superposition. As a result, there is some strain in the PROTAC, which is then relieved by minimization (performed absent any proteins). To gauge the amount of strain, first the initial doubly-superposed conformation is stored. Next, three minimizations are performed. First, the matching core atoms are fixed, and the remaining atoms are strongly tethered, followed by minimization. Next, the core atoms are unfixed, all atoms are weakly tethered, followed by minimization. Finally, all tethers are removed, and the PROTAC is relaxed without constraints. This three-part minimization scheme reduces the effect of bonds unduly stretched by the double-core-superposition procedure described above. After the final, fully unconstrained minimization, the coordinates of the minimized PROTAC are compared to the coordinates after the double-superposition but before the three-part minimization scheme. This all-atom RMSD is termed Core RMSD in the main text. A low value (say, < 1 Å) indicates that the independently sampled PROTAC conformation was a close match

ACS Paragon Plus Environment

11

Journal of Chemical Information and Modeling 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 37

for the accompanying protein-protein docked pose, and there was minimal distortion as the PROTAC conformation was superposed onto the protein-protein docked pose.

Figure 2. Graphical summary of the procedure used to combine conformations from the sampling of the independent PROTAC with protein-protein docked poses of two protein-binding moiety complexes. a) The starting point is a single docked complex of the ligase-binder (cyan ribbons & atoms) and target-binder (green ribbons & atoms). From there, the proteins are removed, and a single conformation of the PROTAC is introduced (gold carbon atoms), leading to b). The userspecified Core 1 of the PROTAC conformation (gray atoms) is superposed onto the corresponding positions of Core 1 in the protein-protein docked position (black atoms). The remainder of the PROTAC is superposed as well, generating c). The user-specified Core 2 of the PROTAC conformation (magenta atoms) is superposed onto the corresponding positions of Core 2 in the protein-protein docked positions (purple atoms). Unlike the previous superposition, the remainder

ACS Paragon Plus Environment

12

Page 13 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

of the PROTAC is not rototranslated, resulting in d). The protein-protein docked atoms are then removed to form distorted structure in e), the isolated PROTAC with its two cores in their positions in the protein-protein docked complex. From there, three minimizations with differing levels of PROTAC restraints (see text) are applied, generating the final predicted PROTAC conformation in f). The Core RMSD metric used throughout this work is between the distorted PROTAC structure and e) and the final, smoothly minimized structure in f). If the fully formed ternary complex has passed both the patch-based threshold and the Core RMSD threshold, it is written to the final database. However, there may still be significant clashing between the proteins and the PROTAC, as they were sampled independently. (Note that, unlike Method 2, there is no appreciable protein-protein clashing, as a result of the protein-protein docking procedure). Optionally, to alleviate this clashing, a restrained minimization protocol may be performed on the final ternary complexes. The backbone atoms of each of the two proteins are defined as independent rigid bodies, while all sidechain atoms and the entirety of the PROTAC are unconstrained, followed by minimization. This protocol generally keeps the overall configuration of the ternary complex intact, but allows for movement of sidechains and/or the PROTAC to alleviate clashing. Occasionally, strands of atoms interpenetrate between the PROTAC and one or both of the proteins, and thus minimization is unable to relieve the clashing. These poses can readily be identified via a large, positive total interaction energies, or through PROTAC-protein overlap values, which are also calculated. It is likely that a subsequent repacking step would fully ameliorate these interpenetrating structures, but this improvement has not yet been implemented.

ACS Paragon Plus Environment

13

Journal of Chemical Information and Modeling 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 37

RESULTS. Validation: Generation and identification of crystal-like ternary complexes. The first crystal structure of a ligase-PROTAC-target ternary complex was reported in 2017,29 where the E3 ligase was the von Hippel-Landau (VHL); the target was the second bromodomain of Brd4, a member of the bromo- and extra-terminal (BET) family; and the PROTAC molecule was MZ1,21 built from JQ1,39 a pan-BET inhibitor, linked via three ethylene glycol monomers to VH032,40 a small molecule ligand of VHL. This crystal structure, deposited as PDB code 5T35 (resolution 2.7 Å), contains two additional chains, Elongin B & C, which are part of the larger E3 ligase structure and were manually removed prior to all modeling, to focus on the most parsimonious view of a ternary complex. The definition of a “crystal-like” pose used in this Validation phase is an in silico pose with a Cα RMSD of ≤10Å after rigid body superposition onto the coordinates of 5T35; this definition is commonly used as the upper limit for an “acceptable” pose in protein-protein docking contexts.41 Each of the four Methods were applied to generate ensembles of possible ternary complexes. The various Methods generate different numbers of possible conformations of ternary complexes: the sampling in Methods 1-3 is capped (by default) at 10,000 conformations, whereas the independent PROTAC and protein-protein samplings in Method 4, when combined, can generate hundreds of thousands of conformations. Ensembles of these magnitudes are impractical for downstream analysis, and thus it is critical to score these initial ensembles, to a priori filter out ternary complexes that are far from the crystallographically observed configuration. It may also be useful to cluster generated ternary complexes and use, e.g. cluster size and cluster diversity as additional scoring metrics. Efforts along these lines are currently underway, but will not be further discussed below in this current work.

ACS Paragon Plus Environment

14

Page 15 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

The filters proposed below depend on the Method and were determined empirically. In identifying possible filters, first each computationally generated ternary complex was superposed as a rigid-body onto a reference crystal structure, such that the target protein was perfectly superposed, followed by evaluation of the Cα RMSD for the ligase. Next, various properties were evaluated – a full listing is reported as Supporting Information, but in general includes interaction energies between the various parts of each ternary complex, properties of the independent PROTAC conformations, and metrics based on both volume and protein patch-overlap (see below). Correlation plots between Cα RMSD and all calculated properties were investigated, and any property that separated the 24 106 44 86 1467 1660 BD2 Brd3 > Brd3 MZ2 34 88 70 75 1034 1595 Brd2 Brd2BD2 MZ2 25 98 54 67 542 745 Agreement?** N N N N Y Y BD2 Brd4 MZ3 Brd4 > 37 133 62 93 2008 2427 Brd3 ~ Brd3BD2 MZ3 40 145 94 127 1418 2213 BD2 Brd2 Brd2 MZ3 34 125 67 82 631 950 Agreement?** N N N N ~ ~ *Qualitative assessment of the Western blots in Figure 2a of Ref. 21 (BRD4 long was used for the Brd4 data). **Qualitative assessment of whether the modeled results agree with the experimental data. A tilde indicates “somewhat.” As can be seen in Table 4, often the number of ternary complexes generated with Method 4 correlates quite well with the experimental ordering, particularly for the E filter set (with some ambiguity resulting due to a degree of subjectivity). However, more disappointing is the fact that, across the three PROTAC molecules, although there is a clear trend in the experimental data for MZ1 to be most effective at enabling protein degradation, followed by MZ2, and then MZ3, for the computational results only the A filter set consistently ranks MZ1 as the most effective PROTAC irrespective of the target BET. This discrepancy may be a consequence of differing properties (e.g., cell permeability, chemical stability, etc.) between MZ1, MZ2, and MZ3, a complication beyond the scope of any of the Methods in this work to account for. However, it should be noted that only two Method 4 simulations were performed for this analysis: one with a single protein-protein docking run, and a second with multiple protein-protein docking runs biased to favor matching hydrophobic patches across the two proteins. From there, filtering is a simple post-processing operation: corresponding filters were applied to the former simulation to generate data for 4A and 4C, and to the latter to generate the Biased datasets (4B, 4D, 4E, and 4F). Thus, while

ACS Paragon Plus Environment

24

Page 25 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

more simulations are needed to confirm or refine these findings, if it indeed holds true that the E set of filtering criteria is better at comparing different targets and the A criteria is better at comparing different PROTACs, only two variants of Method 4 need be run, which is not an undue computational burden. DISCUSSION AND CONCLUSIONS. In both of the publications to date that report29,30 crystal structures of PROTAC-mediated ternary complexes, the authors used the structure – along with computational modeling – to rationally propose new PROTAC derivatives with enhanced selectivity, thereby decisively demonstrating the power of this paired structure-based/modeling approach to PROTAC design. While further crystal structures will surely follow, we have demonstrated that the Methods in this work, particularly the protein-protein docking-based Method 4, are accurate enough to complement ongoing structural campaigns, offering reasonable, crystallike geometries for ternary complexes, which can serve as starting points for testing hypotheses and refining principles of PROTAC design. Moreover, after applying judicious filters – which were determined based on reproducing known X-ray crystal structures – computationally generated ensembles of ternary complexes can be scored via a simple count of acceptable poses. This count was shown to be able to discriminate the degradation tendencies for a mutant vs. its wild type, amongst closely related targets, and to a lesser degree among different PROTAC molecules as well. It is also noteworthy that the non-crystal-like “false positives” generated by Methods 1-4 are not necessarily wrong, for as noted by Nowak et al.,30 there is a certain amount of plasticity in the arrangements of PROTAC-mediated ternary complexes, i.e., multiple low-energy configurations may exist. In that sense, ternary complexes that differ from the observed crystal structure but otherwise pass the filtering criteria may instead be merely thus far unobserved.

ACS Paragon Plus Environment

25

Journal of Chemical Information and Modeling 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 26 of 37

However, even in the most successful cases, where the Methods of this paper produced filtered databases of ternary complexes with ~40% crystal-like poses, it is still difficult to a priori identify which 40% of the filtered output is crystal-like. Efforts to address this remaining challenge by applying various clustering algorithms and thresholds are currently underway in our lab, and of course further work by us and others to develop even more effective means to score or rank productive ternary complexes can continue. Particularly insightful will be the application of the Methods of this paper to additional systems, affording further refinement of the filtering criteria proposed in this initial work. In the meantime, it seems reasonable to state that if the ultimate goal is to use modeling techniques such as those detailed in this work as a surrogate for solving X-ray structures of ternary complexes, then this goal will likely be unmet for the foreseeable future – analogous to how computational protein-protein docking in a generic, non-PROTAC context cannot yet routinely generate crystal-like poses as the single best scoring answer,43 likely due to the typical inaccuracy of these scoring functions.44 Indeed, MOE’s scoring function for the docked target-binder/ligase-binder complexes was unable to reliably identify crystal-like poses when it was applied as a putative filter. Nonetheless, even with this limitation, there is still value in applying the ternary complex modeling algorithms of this work. Hypothetical structural arrangements of target-ligase complexes can be proposed in atomistic detail, which can then be confirmed or rejected via experimental techniques such as structure-guided site-directed mutatgenesis. Moreover, the data presented in this work indicates that using these tools in a relative context is informative. For example, the modeled QVK mutant leads to fewer “reasonable” ternary complexes vis-à-vis the wild type. Specific questions about whether the WT is 1.5x more active (based on the data for Method 4A [51 WT/34 QVK] from Table 3), or 2.02x more active (4B), etc. seem beyond the precision possible

ACS Paragon Plus Environment

26

Page 27 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

with Method 4. (Indeed, achieving such quantitative precision is also a challenge for the experimental, qualitative blotting techniques.) However, despite this imprecision, the overall result of Method 4 is still valuable, insofar as the QVK mutant is generally less capable of forming known (i.e., crystal-like) ternary complexes, and thus it is predicted to be less susceptible to degradation via a PROTAC-mediated mechanism – in line with experiment. With respect to the design of the PROTAC molecules themselves, in a recent paper, Qin et al. reported25 an exquisitely potent (single-digit picomolar) PROTAC, QCA570, targeting the same BET family targeted by the PROTACs discussed in this current work. When comparing QCA570 to the other BET-targeting PROTACs noted above, it is striking how much more rigid QCA570 is, particularly due to the use of two alkynyl spacers. QCA570 was developed through traditional (and very labor-intensive) medicinal chemistry optimization, but our hope is that tools such as Method 4 can suggest likely arrangements of ternary complexes. From there, the parent PROTAC linker can be modeled using well-known tools such as scaffold hopping or another recent tool, also implemented in MOE, that maps out the shortest path along receptor surfaces.45 When paired with these techniques, Method 4 should facilitate the targeting of novel or allosteric pockets – even on the E3 ligases themselves46 – and can be applied to new E3 ligases3,28 beyond those most commonly utilized (VHL and cereblon). Finally, there is another consideration to UPS-driven protein degradation that has not yet been addressed: the necessity for an accessible lysine on a putative ternary complex that can be ubiquitinated. General characteristics of these lysines have been elucidated using bioinformatics approaches,47 and it has also been demonstrated that multiple lysines can be ubiquitinated in a PROTAC-bridged ternary complex.28,29 A deeper analysis of the characteristics of these lysines can potentially add further precision to the modeling results detailed here, as of course will further

ACS Paragon Plus Environment

27

Journal of Chemical Information and Modeling 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 28 of 37

simulations using Method 4 to test and refine the filtering criteria described above. In the meantime, structure-based modeling of PROTAC-bridged ternary complexes, such as enabled via the Methods in this work, should facilitate rational PROTAC design, alongside new advancements in real-time kinetic profiling48 and continued efforts in experimental structural characterization. ASSOCIATED CONTENT Supporting

Information.

The

following

files

are

available

free

of

charge:

Calculated properties of ternary complexes investigated for their ability to identify crystal-like predictions from larger ensembles; Example of a calculated property of ternary complexes that can identify crystal-like poses; Full filtering possibilities for Method 4 for the VHL-MZ1-Brd4BD2 system; Pairwise Cα RMSD matrix for the crystal-like poses produced by Method 4B for the VHLMZ1-Brd4BD2 system. AUTHOR INFORMATION Corresponding Author *Correspondence to: [email protected] ACKNOWLEDGMENT The authors would like to acknowledge Yilin Meng and Ye Che of Pfizer for helpful discussions. ABBREVIATIONS PROTAC, Proteolysis-Targeting Chimera; UPS, Ubiquitin Proteasome System; SAR, Structure Activity Relationship; SVL, Scientific Vector Language; MOE, Molecular Operating Environment; MD, Molecular Dynamics; VHL, von Hippel-Landau; BET, Bromo- and ExtraTerminal; RMSD, Root-Mean-Square Deviation

ACS Paragon Plus Environment

28

Page 29 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

REFERENCES (1)

Che, Y.; Gilbert, A. M.; Shanmugasundaram, V.; Noe, M. C. Inducing Protein-Protein Interactions with Molecular Glues. Bioorganic Med. Chem. Lett. 2018, 28, 2585–2592.

(2)

Stanton, B. Z.; Chory, E. J.; Crabtree, G. R. Chemically Induced Proximity in Biology and Medicine. Science 2018, 359, eaao5902.

(3)

Bulatov, E.; Ciulli, A. Targeting Cullin–RING E3 Ubiquitin Ligases for Drug Discovery: Structure, Assembly and Small-Molecule Modulation. Biochem. J. 2015, 467, 365–386.

(4)

Buckley, D. L.; Crews, C. M. Small-Molecule Control of Intracellular Protein Levels through Modulation of the Ubiquitin Proteasome System. Angew. Chemie - Int. Ed. 2014, 53, 2312–2330.

(5)

Sakamoto, K. M.; Kim, K. B.; Kumagai, A.; Mercurio, F.; Crews, C. M.; Deshaies, R. J. Protacs: Chimeric Molecules That Target Proteins to the Skp1-Cullin-F Box Complex for Ubiquitination and Degradation. Proc. Natl. Acad. Sci. 2001, 98, 8554–8559.

(6)

Itoh, Y.; Ishikawa, M.; Naito, M.; Hashimoto, Y. Protein Knockdown Using Methyl Bestatin-Ligand Hybrid Molecules: Design and Synthesis of Inducers of UbiquitinationMediated Degradation of Cellular Retinoic Acid-Binding Proteins. J. Am. Chem. Soc. 2010, 132, 5820–5826.

(7)

Neklesa, T. K.; Winkler, J. D.; Crews, C. M. Targeted Protein Degradation by PROTACs. Pharmacol. Ther. 2017, 174, 138–144.

(8)

Collins, I.; Wang, H.; Caldwell, J. J.; Chopra, R. Chemical Approaches to Targeted Protein Degradation through Modulation of the Ubiquitin–proteasome Pathway. Biochem. J. 2017,

ACS Paragon Plus Environment

29

Journal of Chemical Information and Modeling 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 30 of 37

474, 1127–1147. (9)

Churcher, I. Protac-Induced Protein Degradation in Drug Discovery: Breaking the Rules or Just Making New Ones? J. Med. Chem. 2018, 61, 444–452.

(10)

Ottis, P.; Crews, C. M. Proteolysis-Targeting Chimeras: Induced Protein Degradation as a Therapeutic Strategy. ACS Chem. Biol. 2017, 12, 892–898.

(11)

Lai, A. C.; Crews, C. M. Induced Protein Degradation: An Emerging Drug Discovery Paradigm. Nat. Rev. Drug Discov. 2017, 16, 101–114.

(12)

Hughes, S. J.; Ciulli, A. Molecular Recognition of Ternary Complexes: A New Dimension in the Structure-Guided Design of Chemical Degraders. Essays Biochem. 2017, 61, 505516.

(13)

Tinworth, C. P.; Lithgow, H.; Churcher, I. Small Molecule-Mediated Protein Knockdown as a New Approach to Drug Discovery. Medchemcomm 2016, 7, 2206–2216.

(14)

Crews, C. M. Inducing Protein Degradation as a Therapeutic Strategy. J. Med. Chem. 2018, 61, 403–404.

(15)

Bondeson, D. P.; Smith, B. E.; Burslem, G. M.; Buhimschi, A. D.; Hines, J.; JaimeFigueroa, S.; Wang, J.; Hamman, B. D.; Ishchenko, A.; Crews, C. M. Lessons in PROTAC Design from Selective Degradation with a Promiscuous Warhead. Cell Chem. Biol. 2018, 25, 78–87.

(16)

Lai, A. C.; Toure, M.; Hellerschmied, D.; Salami, J.; Jaime-Figueroa, S.; Ko, E.; Hines, J.; Crews, C. M. Modular PROTAC Design for the Degradation of Oncogenic BCR-ABL. Angew. Chemie - Int. Ed. 2016, 55, 807–810.

ACS Paragon Plus Environment

30

Page 31 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

(17)

Burslem, G. M.; Smith, B. E.; Lai, A. C.; Jaime-Figueroa, S.; McQuaid, D. C.; Bondeson, D. P.; Toure, M.; Dong, H.; Qian, Y.; Wang, J.; Crew, A. P.; Hines, J.; Crews, C. M. The Advantages of Targeted Protein Degradation Over Inhibition: An RTK Case Study. Cell Chem. Biol. 2018, 25, 67–77.

(18)

Bondeson, D. P.; Mares, A.; Smith, I. E. D.; Ko, E.; Campos, S.; Miah, A. H.; Mulholland, K. E.; Routly, N.; Buckley, D. L.; Gustafson, J. L.; Zinn, N.; Grandi, P.; Shimamura, S.; Bergamini, G.; Faelth-Savitski, M.; Bantscheff, M.; Cox, C.; Gordon, D. A.; Willard, R. R.; Flanagan, J. J.; Casillas, L. N.; Votta, B. J.; den Besten, W.; Famm, K.; Kruidenier, L.; Carter, P. S.; Harling, J. D.; Churcher, I.; Crews, C. M. Catalytic in Vivo Protein Knockdown by Small-Molecule PROTACs. Nat. Chem. Biol. 2015, 11, 611–617.

(19)

Buhimschi, A. D.; Armstrong, H. A.; Toure, M.; Jaime-Figueroa, S.; Chen, T. L.; Lehman, A. M.; Woyach, J. A.; Johnson, A. J.; Byrd, J. C.; Crews, C. M. Targeting the C481S Ibrutinib-Resistance Mutation in Bruton’s Tyrosine Kinase Using PROTAC-Mediated Degradation. Biochemistry 2018, 57, 3564–3575.

(20)

Bassi, Z. I.; Fillmore, M. C.; Miah, A. H.; Chapman, T. D.; Maller, C.; Roberts, E. J.; Davis, L. C.; Lewis, D. E.; Galwey, N. W.; Waddington, K. E.; Parravicini, V.; Macmillan-Jones, A. L.; Gongora, C.; Humphreys, P. G.; Churcher, I.; Prinjha, R. K.; Tough, D. F. Modulating PCAF/GCN5 Immune Cell Function through a PROTAC Approach. ACS Chem. Biol. 2018, 13, 2862-2867.

(21)

Zengerle, M.; Chan, K. H.; Ciulli, A. Selective Small Molecule Induced Degradation of the BET Bromodomain Protein BRD4. ACS Chem. Biol. 2015, 10, 1770–1777.

(22)

DeGoey, D. A.; H.-J., C.; Cox, P. B.; Wendt, M. D. Beyond the Rule of 5: Lessons Learned

ACS Paragon Plus Environment

31

Journal of Chemical Information and Modeling 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 32 of 37

from AbbVie’s Drugs and Compound Collection. J. Med. Chem. 2018, 61, 2636–2651. (23)

Chan, K. H.; Zengerle, M.; Testa, A.; Ciulli, A. Impact of Target Warhead and Linkage Vector on Inducing Protein Degradation: Comparison of Bromodomain and Extra-Terminal (BET) Degraders Derived from Triazolodiazepine (JQ1) and Tetrahydroquinoline (IBET726) BET Inhibitor Scaffolds. J. Med. Chem. 2018, 61, 504–513.

(24)

Crew, A. P.; Raina, K.; Dong, H.; Qian, Y.; Wang, J.; Vigil, D.; Serebrenik, Y. V.; Hamman, B. D.; Morgan, A.; Ferraro, C.; Siu, K.; Neklesa, T. K.; Winkler, J. D.; Coleman, K. G.; Crews, C. M. Identification and Characterization of von Hippel-Lindau-Recruiting Proteolysis Targeting Chimeras (PROTACs) of TANK-Binding Kinase 1. J. Med. Chem. 2018, 61, 583–598.

(25)

Qin, C.; Hu, Y.; Zhou, B.; Fernandez-Salas, E.; Yang, C. Y.; Liu, L.; McEachern, D.; Przybranowski, S.; Wang, M.; Stuckey, J.; Meagher, J.; Bai, L.; Chen, Z.; Lin, M.; Yang, J.; Ziazadeh, D. N.; Xu, F.; Hu, J.; Xiang, W.; Huang, L.; Li, S.; Wen, B.; Sun, D.; Wang, S. Discovery of QCA570 as an Exceptionally Potent and Efficacious Proteolysis Targeting Chimera (PROTAC) Degrader of the Bromodomain and Extra-Terminal (BET) Proteins Capable of Inducing Complete and Durable Tumor Regression. J. Med. Chem. 2018, 61, 6685–6704.

(26)

Cyrus, K.; Wehenkel, M.; Choi, E. Y.; Lee, H.; Swanson, H.; Kim, K. B. Jostling for Position: Optimizing Linker Location in the Design of Estrogen Receptor-Targeting PROTACs. ChemMedChem 2010, 5, 979–985.

(27)

Maniaci, C.; Hughes, S. J.; Testa, A.; Chen, W.; Lamont, D. J.; Rocha, S.; Alessi, D. R.; Romeo, R.; Ciulli, A. Homo-PROTACs: Bivalent Small-Molecule Dimerizers of the VHL

ACS Paragon Plus Environment

32

Page 33 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

E3 Ubiquitin Ligase to Induce Self-Degradation. Nat. Commun. 2017, 8, 830. (28)

Ottis, P.; Toure, M.; Cromm, P. M.; Ko, E.; Gustafson, J. L.; Crews, C. M. Assessing Different E3 Ligases for Small Molecule Induced Protein Ubiquitination and Degradation. ACS Chem. Biol. 2017, 12, 2570–2578.

(29)

Gadd, M. S.; Testa, A.; Lucas, X.; Chan, K.-H.; Chen, W.; Lamont, D. J.; Zengerle, M.; Ciulli, A. Structural Basis of PROTAC Cooperative Recognition for Selective Protein Degradation. Nat. Chem. Biol. 2017, 13, 514–521.

(30)

Nowak, R. P.; Deangelo, S. L.; Buckley, D.; He, Z.; Donovan, K. A.; An, J.; Safaee, N.; Jedrychowski, M. P.; Ponthier, C. M.; Ishoey, M.; Zhang, T.; Mancias, J. D.; Gray, N. S.; Bradner, J. E.; Fischer, E. S. Plasticity in Binding Confers Selectivity in Ligand-Induced Protein Degradation. Nat. Chem. Biol. 2018, 14, 706–714.

(31)

Winter, G. E.; Buckley, D. L.; Paulk, J.; Roberts, J. M.; Souza, A.; Dhe-Paganon, S.; Bradner, J. E. Phthalimide Conjugation as a Strategy for in Vivo Target Protein Degradation. Science 2015, 348, 1376–1381.

(32)

Remillard, D.; Buckley, D. L.; Paulk, J.; Brien, G. L.; Sonnett, M.; Seo, H. S.; Dastjerdi, S.; Wühr, M.; Dhe-Paganon, S.; Armstrong, S. A.; Bradner, J. E. Degradation of the BAF Complex Factor BRD9 by Heterobifunctional Ligands. Angew. Chemie - Int. Ed. 2017, 56, 5738–5743.

(33)

Schiedel, M.; Herp, D.; Hammelmann, S.; Swyter, S.; Lehotzky, A.; Robaa, D.; Oláh, J.; Ovádi, J.; Sippl, W.; Jung, M. Chemically Induced Degradation of Sirtuin 2 (Sirt2) by a Proteolysis Targeting Chimera (PROTAC) Based on Sirtuin Rearranging Ligands

ACS Paragon Plus Environment

33

Journal of Chemical Information and Modeling 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 34 of 37

(SirReals). J. Med. Chem. 2018, 61, 482–491. (34)

Lebraud, H.; Wright, D. J.; Johnson, C. N.; Heightman, T. D. Protein Degradation by InCell Self-Assembly of Proteolysis Targeting Chimeras. ACS Cent. Sci. 2016, 2, 927–934.

(35)

Zorba, A.; Nguyen, C.; Xu, Y.; Starr, J.; Borzilleri, K.; Smith, J.; Zhu, H.; Farley, K. A.; Ding, W.; Schiemer, J.; Feng, X.; Chang, J. S.; Uccello, D. P.; Young, J. A.; Garcia-Irrizary, C. N.; Czabaniuk, L.; Schuff, B.; Oliver, R.; Montgomery, J.; Hayward, M. M.; Coe, J.; Chen, J.; Niosi, M.; Luthra, S.; Shah, J. C.; El-Kattan, A.; Qiu, X.; West, G. M.; Noe, M. C.; Shanmugasundaram, V.; Gilbert, A. M.; Brown, M. F.; Calabrese, M. F. Delineating the Role of Cooperativity in the Design of Potent PROTACs for BTK. Proc. Natl. Acad. Sci. 2018, 115, E7285-E7292.

(36)

Molecular Operating Environment (MOE) 2018.0101. Chemical Computing Group, ULC: Montreal, Quebec, H3A 2R7, Canada 2018.

(37)

Ito, T.; Ando, H.; Suzuki, T.; Ogura, T.; Hotta, K.; Imamura, Y.; Yamaguchi, Y.; Handa, H. Identification of a Primary Target of Thalidomide Teratogenicity. Science 2010, 327, 1345–1350.

(38)

Lopez-Girona, A.; Mendy, D.; Ito, T.; Miller, K.; Gandhi, A. K.; Kang, J.; Karasawa, S.; Carmel, G.; Jackson, P.; Abbasian, M.; Mahmoudi, A.; Cathers, B.; Rychak, E.; Gaidarova, S.; Chen, R.; Schafer, P. H.; Handa, H.; Daniel, T. O.; Evans, J. F.; Chopra, R. Cereblon Is a Direct Protein Target for Immunomodulatory and Antiproliferative Activities of Lenalidomide and Pomalidomide. Leukemia 2012, 26, 2326–2335.

(39)

Qi, J.; Picaud, S.; Shen, Y.; Smith, W. B.; Fedorov, O.; Morse, E. M.; Keates, T.; Hickman,

ACS Paragon Plus Environment

34

Page 35 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

T. T.; Felletar, I.; Philpott, M.; Munro, S.; McKeown, M. R.; Wang, Y.; Christie, A. L.; West, N.; Cameron, M. J.; Schwartz, B.; Heightman, T. D.; La Thangue, N.; French, C. A.; Wiest, O.; Kung, A. L.; Knapp, S.; Bradner, J. E. Selective Inhibition of BET Bromodomains. Nature 2010, 468, 1067–1073. (40)

Galdeano, C.; Gadd, M. S.; Soares, P.; Scaffidi, S.; Van Molle, I.; Birced, I.; Hewitt, S.; Dias, D. M.; Ciulli, A. Structure-Guided Design and Optimization of Small Molecules Targeting the Protein–Protein Interaction between the von Hippel–Lindau (VHL) E3 Ubiquitin Ligase and the Hypoxia Inducible Factor (HIF) Alpha Subunit with in Vitro Nanomolar Affinities. J. Med. Chem. 2014, 57, 8657–8663.

(41)

Huang, S. Y. Search Strategies and Evaluation in Protein-Protein Docking: Principles, Advances and Challenges. Drug Discov. Today 2014, 19, 1081–1096.

(42)

Douglass, E. F.; Miller, C. J.; Sparer, G.; Shapiro, H.; Spiegel, D. A. A Comprehensive Mathematical Model for Three-Body Binding Equilibria. J. Am. Chem. Soc. 2013, 135, 6092–6099.

(43)

Rodrigues, J. P. G. L. M.; Bonvin, A. M. J. J. Integrative Computational Modeling of Protein Interactions. FEBS J. 2014, 281, 1988–2003.

(44)

Huang, S. Y. Exploring the Potential of Global Protein-Protein Docking: An Overview and Critical Assessment of Current Programs for Automatic Ab Initio Docking. Drug Discov. Today 2015, 20, 969–977.

(45)

Pérez-Benito, L.; Henry, A.; Matsoukas, M.-T.; Lopez, L.; Pulido, D.; Royo, M.; Cordomí, A.; Tresadern, G.; Pardo, L. The Size Matters? A Computational Tool to Design Bivalent

ACS Paragon Plus Environment

35

Journal of Chemical Information and Modeling 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 36 of 37

Ligands. Bioinformatics 2018, 34, 3857-3863. (46)

Lucas, X.; Van Molle, I.; Ciulli, A. Surface Probing by Fragment-Based Screening and Computational Methods Identifies Ligandable Pockets on the von Hippel-Lindau (VHL) E3 Ubiquitin Ligase. J. Med. Chem. 2018, 61, 7387–7393.

(47)

Guharoy, M.; Bhowmick, P.; Sallam, M.; Tompa, P. Tripartite Degrons Confer Diversity and Specificity on Regulated Protein Degradation in the Ubiquitin-Proteasome System. Nat. Commun. 2016, 7, 1–13.

(48)

Riching, K. M.; Mahan, S.; Corona, C. R.; McDougall, M.; Vasta, J. D.; Robers, M. B.; Urh, M.; Daniels, D. L. Quantitative Live-Cell Kinetic Degradation and Mechanistic Profiling of PROTAC Mode of Action. ACS Chem. Biol. 2018, 13, 2758–2770.

ACS Paragon Plus Environment

36

Page 37 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Information and Modeling

For Table of Contents use only In silico Modeling of PROTAC-Mediated Ternary Complexes: Validation and Application Michael L. Drummond and Christopher I. Williams

ACS Paragon Plus Environment

37