==== Front Molecules Molecules molecules Molecules 1420-3049 MDPI 33276504 10.3390/molecules25235683 molecules-25-05683 Article Antiangiogenic Activity and in Silico Cereblon Binding Analysis of Novel Thalidomide Analogs https://orcid.org/0000-0002-4290-5268Peach Megan L. 1 Beedie Shaunna L. 23 Chau Cindy H. 2 https://orcid.org/0000-0002-9305-3560Collins Matthew K. 2 Markolovic Suzana 2 Luo Weiming 4 Tweedie David 4 https://orcid.org/0000-0001-5638-1955Steinebach Christian 5 Greig Nigel H. 4 https://orcid.org/0000-0002-9376-7897Gütschow Michael 5 https://orcid.org/0000-0001-8027-114XVargesson Neil 3 Nicklaus Marc C. 6 https://orcid.org/0000-0003-2428-5613Figg William D. 2* 1 Basic Science Program, Chemical Biology Laboratory, Frederick National Laboratory for Cancer Research, National Cancer Institute, Frederick, MD 21701, USA; megan.peach@nih.gov 2 Molecular Pharmacology Section, Genitourinary Malignancies Branch, Center for Cancer Research, National Cancer Institute, NIH, Bethesda, MD 20892, USA; shaunnabeedie90@gmail.com (S.L.B.); chauc@mail.nih.gov (C.H.C.); collinsmk@nih.gov (M.K.C.); suzana.markolovic@gmail.com (S.M.) 3 School of Medicine, Medical Sciences & Nutrition, Institute of Medical Sciences, University of Aberdeen, Aberdeen AB25 2ZD, UK; n.vargesson@abdn.ac.uk 4 Drug Design & Development Section, Translational Gerontology Branch, National Institute on Aging, NIH, Baltimore, MD 21224, USA; luowe@grc.nia.nih.gov (W.L.); tweedieda@grc.nia.nih.gov (D.T.); greign@grc.nia.nih.gov (N.H.G.) 5 Pharmaceutical Institute, University of Bonn, 53121 Bonn, Germany; c.steinebach@uni-bonn.de (C.S.); guetschow@uni-bonn.de (M.G.) 6 Chemical Biology Laboratory, Center for Cancer Research, National Cancer Institute, NIH, Frederick, MD 21701, USA; nicklaum@mail.nih.gov * Correspondence: figgw@mail.nih.gov; Tel.: +1-240-760-6179; Fax: +1-240-858-3020 02 12 2020 12 2020 25 23 568318 9 2020 27 11 2020 © 2020 by the authors.2020Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).Due to its antiangiogenic and anti-immunomodulatory activity, thalidomide continues to be of clinical interest despite its teratogenic actions, and efforts to synthesize safer, clinically active thalidomide analogs are continually underway. In this study, a cohort of 27 chemically diverse thalidomide analogs was evaluated for antiangiogenic activity in an ex vivo rat aorta ring assay. The protein cereblon has been identified as the target for thalidomide, and in silico pharmacophore analysis and molecular docking with a crystal structure of human cereblon were used to investigate the cereblon binding abilities of the thalidomide analogs. The results suggest that not all antiangiogenic thalidomide analogs can bind cereblon, and multiple targets and mechanisms of action may be involved. angiogenesiscereblondockingstructure–activity relationshipsthalidomide ==== Body 1. Introduction Despite its teratogenic toxicity, the antiangiogenic and anti-immunoinflammatory efficacy of thalidomide and its analogs have expanded its clinical use in the treatment of Hansen’s disease as well as multiple myeloma and other cancers [1,2,3]. As a result, there are intense efforts to improve the potency of thalidomide while decreasing its toxicity. Thalidomide was originally believed to be a multitarget drug due to its wide range of biological effects, including reduced TNF-α production, decreased or destabilized COX-2 expression, downregulation of VEGF and FGF, NF-κB inhibition (possibly by suppression of IκB kinase), inhibition of prostaglandin E2 secretion, and α1-acid glycoprotein binding [4,5]. In 2010, however, cereblon was identified as the single target of thalidomide binding, using an affinity purification assay with thalidomide-linked magnetic nanoparticles [6]. Cereblon is the substrate recognition component of a DDB1-CUL4-RBX1 E3 ubiquitin ligase complex. The current model for the bioactivity of thalidomide and its analogs is that binding to cereblon itself induces all downstream effects by triggering the ubiquitin-dependent proteasomal degradation of substrates for the E3 complex [7,8,9]. The chemical structure of thalidomide consists of two linked rings, a phthalimide and a glutarimide (Figure 1). Recent studies demonstrate that crystal structures of thalidomide, lenalidomide, or pomalidomide binding to cereblon [10,11,12] have shown that this interaction is mediated almost entirely by the glutarimide ring, whereas the phthalimide ring, along with a small surrounding region of the cereblon surface, is involved in binding to neo-substrates for ubiquitination in the ternary complex [13,14]. Thalidomide possesses one chiral center (the C3-carbon atom of the glutarimide ring) and consists of a racemic mixture of two optical isomers, (S)- and (R)-enantiomers that interconvert under physiological conditions [15]. Structural and biochemical studies established that the (S)-enantiomer exhibited a 10-fold stronger binding to cereblon and inhibition of self-ubiquitylation compared to the (R)-isomer with the teratogenic effects induced by the (S)-isomer [16,17]. Many previous efforts on designing thalidomide analogs to potentiate various biological effects have focused on modifications of the glutarimide ring [21,22,23]. However, since the interaction with cereblon is mediated by the glutarimide ring, it is possible that some of these analogs may not actually bind cereblon and may be exerting their effects by other mechanisms. Such a situation occurred with a series of thalidomide analogs developed by Celgene, initially based on the structure of thalidomide’s glutarimide ring-hydrolysis products [24]. These compounds, in which the glutarimide ring was replaced by a branched 3,4-dialkoxyphenyl-containing moiety, showed much more potent TNF-α inhibition relative to thalidomide. This was discovered to be due to a shift in the mechanism of action to binding and inhibiting phosphodiesterase 4 (PDE4), which elevates cAMP levels and in turn decreases TNF-α. Thalidomide itself does not bind PDE4 [25]. In this work, a cohort of 27 structurally diverse, recently identified thalidomide analogs [19,20] was tested in an ex vivo rat aorta ring (RAR) assay for angiogenesis inhibition. Notably, several of these analogs were found to have antiangiogenic, anti-inflammatory, and/or anticancer activity in in vivo models [18,26]. Given the chemical diversity of these thalidomide analogs, it is unknown which of them could bind to cereblon and which might have other targets. The recent structural studies of cereblon [10] enabled in silico protein docking experiments to analyze the chemical and structural requirements for cereblon binding and investigate their relationship to antiangiogenic activity. 2. Results and Discussion 2.1. Biological Testing of Thalidomide Analogs Recent in vitro and in vivo studies have identified the antiangiogenic and therapeutic potential of a novel series of heterocyclic and adamantyl ring-based [18,19] and tetrafluorinated [20,26] thalidomide analogs (Figure 1). To develop a model predictive of antiangiogenic activity, these analogs were screened in an ex vivo angiogenesis model, the rat aorta ring (RAR) assay. We previously demonstrated the RAR assay to be an accurate assessment of antiangiogenic activity, with cytotoxicity having no effect on microvessel outgrowth [27]. This assay, which is more representative than in vitro models and recapitulates the complexities of angiogenesis [28], was selected to assess antiangiogenic activity without in vivo metabolic activation. Using previously described conditions, rat aortic rings treated with TNP-470 (50 μM; positive control) showed little to no outgrowth (17.94% ± 5.42, n = 7) compared to vehicle-treated (0.5% DMSO; control) rings (98.05% ± 5.03, n = 11). Treatment with thalidomide (50 μM) itself had no effect on outgrowth (94.67% ± 8.25, n = 4). While pomalidomide treatment showed a similar response (98.14% ± 17.69, n = 4), lenalidomide treatment significantly decreased outgrowth (45.84% ± 5.65, n = 4). Compounds C4, C29, C46, C55, and C86 also significantly reduced microvessel outgrowth (C4 (18.54% ± 7.45, n = 4), C29 (37.43% ± 8.98, n = 3), C46 (20.72% ± 3.74, n = 7), C55 (22.95% ± 3.74, n = 3), and C86 (17.3% ± 9.34, n = 3)). The tetrafluorinated compounds were the most potent angiogenic inhibitors in this assay, where the inhibition of outgrowth was indistinguishable from the positive control-treated rings (Gu973 (5.91% ± 0.21, n = 3), Gu992 (11.37% ± 4.1, n = 3), Gu998 (6.17% ± 1.06, n = 3), and Gu1029 (9.02% ± 0.63, n = 4)). Figure 2 summarizes the RAR results of tested analogs. 2.2. In Silico Protein Docking Simulation Biochemical and structural studies have recently identified cereblon as a binding partner for thalidomide [6], and indicated that binding is mediated by the glutarimide ring in a relatively small and hydrophobic binding pocket composed of three tryptophan residues [10,11]. Since the series of thalidomide analogs under investigation here have been mainly designed with modifications to the glutarimide moiety, we investigated whether the cereblon binding pocket is able to accommodate these analogs and the extent to which computationally predicted cereblon binding affinities correlate with antiangiogenic activity in the RAR assay. The crystal structure of human cereblon in complex with lenalidomide (PDB entry 4TZ4) [10] was selected as the most appropriate for docking, even though it is of relatively low resolution (3.01Å), for two reasons: (1) it is the only available human ortholog structure, and (2) it is the full length cereblon structure rather than the truncated C-terminal cereblon-binding domain. Previous docking studies have shown that the truncated domain, which lacks part of two β-strands and a connecting hairpin loop, is not sufficient for successful docking [29]. In the full-length structure, these strands are stabilized by packing against the N-terminal domain, and the hairpin loop forms part of the binding site [10]. The ligand binding site in cereblon is quite small, and only the glutarimide ring of thalidomide or lenalidomide is buried in the pocket. The phthalimide ring hydrogen bonds to an asparagine residue on the pocket rim, but is otherwise fairly solvent exposed. The glutarimide imide group forms three hydrogen bonds on one side of the pocket, and the aliphatic carbons on the other side of the glutarimide ring are closely packed between three tryptophan residues [10,11]. As a consequence of this close packing, a preliminary docking run (results not shown) found that only compounds with an intact, unmodified glutarimide ring could dock correctly into the binding site. However, studies on the structural dynamics of the bacterial cereblon C-terminal domain have shown large flexibility and even partial unfolding of the domain in the absence of ligand, suggesting that cereblon might be able to adaptively bind native ligands of widely different sizes and chemistry compared to thalidomide [30]. We therefore developed a hybrid ligand-ensemble and induced-fit docking protocol. Ligand-ensemble docking is based on the assumption that ligands structurally similar to the cocrystal ligand should bind in a similar orientation, and even chemically divergent ligands should have some overlap in the types of interactions formed with residues in the binding site. Leveraging this assumption and taking advantage of the existing cocrystal cereblon structures can improve docking results for a series of congeneric analogs [31,32]. An initial docked pose for each thalidomide analog was generated by flexible alignment to the bound lenalidomide crystal ligand in its complex with human cereblon (4TZ4) [10]. Thus, each ligand was initially placed in the same binding mode, aligned via the glutarimide ring (when present) or the phathalimide ring otherwise. In most cases, this initial pose had serious steric clashes with the cereblon binding site, but each protein–ligand complex was then optimized with induced-fit docking. Induced-fit docking methods model flexibility and conformational changes in the protein binding site in various ways [33]. Here, protein sidechains in the neighborhood of the bound ligand pose were conformationally sampled, repacked, and energy minimized to optimize the binding site fit to the ligand [34,35]. With this protocol, the exact placement of the sidechains in a lower resolution protein structure is less of an issue because specific interactions for each individual ligand are independently optimized while the overall binding mode is maintained. In some cases, repacking of the cereblon binding site to accommodate a larger ligand required one of the tryptophan residues to flip away, out of the binding site, with a high protein conformational energy. We considered this to be evidence that it was not possible for that ligand to fit and bind cereblon. This scenario occurred with compounds having a large or very differently shaped glutarimide ring replacement: C53, C59, C64, and C65 as well as the methylene-bridged adamantyl compounds C70 and C74. The remaining compounds were able to be accommodated in the cereblon binding site without major sidechain conformational changes and proceeded to the final docking stage, where they were re-docked into their remodeled binding sites and a final induced-fit docking score was calculated as a function of the protein conformational energy and the nonbonded interactions between ligand and protein [35]. Induced fit docking scores are presented in Supplemental Table S1. As expected, thalidomide, pomalidomide, and lenalidomide themselves docked well and scored highly (Figure 3A). The thioxo analogs of thalidomide, C2–C19, generally fit easily into the cereblon binding site, though they appeared to bind more weakly than compounds with an unmodified glutarimide ring because sulfur is both a weaker hydrogen-bond acceptor and a bulkier atom than oxygen. When present at position Z3, the sulfur atom leads to a less preferred accommodation of the ligand than at position Z4 (Figure 1), and accordingly compounds C4 and C9 have the worst docking scores (Figure 3B). Compounds C44 and C46, with an unsaturated planar lactam ring, have a more pronounced tilt angle of the phthalimide skeleton relative to the 3,4-dihydropyridin-2(1H)-one moiety. Structures of cereblon complexed with the bound substrates CK1α [14] and GSPT1 [13] show that substrate binding also tilts the isoindolinone substructure in lenalidomide and the analog CC-885 in this same direction, so C44 and C46 may have bound conformations preoptimized to facilitate protein–substrate interaction (Figure 3C). Compound C55, in which the glutarimide ring has a different attachment geometry, docks well to cereblon. The (S)-stereoisomer appears to be preferred. Here, the isoindoline-1,3-diimine is also positioned differently relative to the substrate binding site, and it seems plausible that this may also be favorably affecting substrate binding (Figure 3D), although the exact mechanism for antiangiogenic activity remains to be determined. None of the adamantane-type compounds (C70–C93) fit well into the binding site or score well in docking. As mentioned above, the methylene-bridged compounds, C70 and C74, are completely unable to fit. The best-scoring one is the 2-adamantyl derivative C83 (Figure 3E), and the hydroxylated 1-adamantyl derivative C86 is barely able to fit into the binding site and shifts the isoindolinone out and away from its position in the other complexes (Figure 3F). Lastly, the tetrafluorinated compounds, where the glutarimide ring has been replaced with a substituted barbiturate ring, appear to fit well in the binding pocket and score as high as thalidomide itself. Here, however, the (R)-stereoisomers are predicted to bind, not the (S)-stereoisomers. The barbiturate ring has an additional carbonyl group which remains partially solvent accessible, and the carbonyl groups adjacent to R1 and R2 (Figure 1) make productive packing interactions with the binding site tryptophan residues (Figure 3G). Opening of the fused five-membered ring and replacing the CO unit, as in compounds Gu992 and Gu1029, slightly shifts the orientation of the tetrafluorobenzoyl ring (Figure 3H). Based on this docking pose, the fluorine atoms would not be involved with cereblon binding, but might affect the interaction with a substrate. In Figure 4, the induced-fit docking scores are plotted relative to % inhibition of angiogenesis in the RAR assay. Docking scores are generally at best only loosely correlated with target binding affinity, and it is known that the in vivo effects of thalidomide and its analogs do not necessarily correlate with cereblon binding affinity [36,37], perhaps due to the fact that they are exerting their effects by forming (or inhibiting the formation of) a ternary complex with cereblon, ligand, and substrate. Nevertheless, there was a general trend of increased cereblon binding leading to decreased angiogenic activity. The outliers are the most interesting compounds to consider here, and these are thalidomide, pomalidomide, and compound C17, which are predicted computationally to bind well and yet are inactive in the ex vivo assay. It has been argued that this is because thalidomide requires metabolic activation [38], but it may also be that thalidomide and pomalidomide, with two labile imide substructures, are quite vulnerable to spontaneous nonenzymatic hydrolysis [4,39,40], whereas compounds such as lenalidomide with an isoindolinone ring are hydrolyzed much more slowly and less extensively [41]. The compounds with a thioamide moiety within the five-membered ring (C29, C74) as well as those bearing a 6-thioxopiperidin-2-one (C2, C7, C14, C17, C19) or piperidine-2,6-dithione (C4, C9) are better stabilized and might be less prone to hydrolysis. Compound C17 can be compared to compound C14, which differs only in the position of the hydroxyl group on the phthalimide ring (R2 vs. R1, Figure 1). With the OH group as R2, outgrowth is not inhibited at all, whereas with OH as R1, outgrowth is reduced to 60%. This effect can also be seen in comparing compound C34 to lenalidomide, where positioning the amine group either as R2 or as R1 causes a fairly large shift in activity, from 71% outgrowth in the case of C34 to 46% outgrowth in the case of lenalidomide. The other outliers are compounds C4 and C86, which show excellent antiangiogenic activity, yet are not predicted to bind well to cereblon (C4) or even to fit into the binding site at all (C86). This suggests the possibility that another target or targets could be mediating these effects. We investigated this possibility with a ligand-based pharmacophore analysis. The set of aligned thalidomide analogs was used to develop a set of pharmacophore hypotheses, and multiple potential binding modes were searched for by clustering the active, antiangiogenic compounds according to which pharmacophore hypotheses they matched. The pharmacophore clustering revealed two potential binding modes. In this case, rather than different orientations in the same binding site, we suspect similar phthalimide binding sites on different targets. Compounds C4 and C86, along with other thioamide compounds C2, C7, C14, C19, and the lactam compounds C44 and C46, can be fit to a pharmacophore with an acceptor feature and an aromatic ring feature on the phthalimide or isoindolinone ring, and a hydrophobic feature along with a donor on the glutarimide or other ring. Thioamides are weak acceptors via their sulfur atoms, but stronger donors than amides, and the hydroxyl group in C86 can be aligned relatively closely to this donor feature as well (Figure 5A). The other potential binding mode has an aromatic ring feature on the phthalimide and the acceptor-donor-acceptor features on the imide group of the glutarimide ring, which matches lenalidomide, C29, C34, C55, Gu973, Gu998, Gu992, and Gu1029 (Figure 5B). The docking results suggest that the requirement for a match at the second acceptor site (Z4 in the glutarimide, Figure 1) is not necessary for favorable cereblon binding, and thus some compounds can match both pharmacophores. 3. Material and Methods 3.1. Thalidomide Analogs A broad series of novel thalidomide-based compounds were synthesized [19,20], dissolved in DMSO, and stored in stock concentration of 10 mM. The chemical structures of lead compounds of interest were confirmed by chemical characterization (purity > 99.5%). 3.2. Rat Aorta Ring (RAR) Assay of Angiogenesis The rat aortic ring assay was performed as previously described [23,26,38]. Briefly, 24-well tissue culture plates were covered with 250 µL of Matrigel (BD Biosciences) and allowed to set for 1 h at room temperature. Six- to eight-week old male Sprague Dawley rats were euthanized, and the descending aortas were dissected and cleaned in EBM media. The aorta was sliced to 1-mm cross-sections, placed on Matrigel-coated wells, and layered with additional Matrigel (250 µL). The rings were allowed to set, after which they were covered with endothelial cell growth media (EGM-II, Lonza, Walkersville, MD, USA) and incubated under 5% CO2 at 37 °C overnight. EGM-II consists of endothelial cell basal medium (EBM-II) and endothelial cell growth factors. The next day, media was replaced with EBM-II containing either the vehicle control (0.5% DMSO), 50 μM TNP-470 (a known angiogenesis inhibitor [42] as the positive control) or the test compounds at 50 μM. Rings were incubated for 4 days and then imaged on day 5 using an EVOS scope. The experiments were performed in triplicates using aortas from 3-4 different rats. The area of angiogenic sprouting, reported in square pixels, was quantified using Adobe Photoshop. Data was presented as percent growth based on the negative control (vehicle), which was normalized to 100% growth. 3.3. Molecular Modeling All molecular modeling was performed using Schrödinger software (Schrödinger, LLC; New York, NY, USA). version 2017-4. Conformations for the thalidomide analogs were generated in LigPrep using the MMFF94s forcefield [43]. The S chirality was retained for analogs with the same overall scaffold as thalidomide (C2–C53), and stereoisomers were enumerated for chiral centers in analogs where the scaffold varied (C55, C74, Gu973, Gu998, Gu992, and Gu1029). Two low energy ring conformations were generated for each analog because the glutarimide ring has two possible puckering conformations. Ligand conformers were aligned flexibly to the crystallized lenalidomide ligand in PDB structure 4TZ4 [10]. For compounds with an unmodified or minimally modified glutarimide ring (C2–C55), the atoms involved in hydrogen bonding interactions were matched: three carbonyl oxygen acceptors, and the NH donor in the glutarimide ring. For compounds lacking a glutarimide ring (C59–C93), the phthalimide ring was used for alignment. Both the R and S isomers of the barbiturate ring in compounds Gu973, Gu992, Gu998, and Gu1029 were aligned to hydrogen bonding atoms of the glutarimide. Upon alignment the isomers differ in whether the nitrogen aligned with the glutarimide donor NH is acyl-substituted. 3.4. Induced-Fit Docking As the ensemble of analog structures was now aligned to the lenalidomide-based pharmacophore in the frame of reference of the 4TZ4 crystal structure, these could now be used as an initial pose input for induced-fit docking. The 4TZ4 structure was prepared for docking [44] by building coordinates for missing residues, setting up zero-order bonds to the zinc atom, deleting all waters, protonating with simplified rules at neutral pH, and optimizing the hydrogen bond network. This was followed by a restrained minimization in the OPLS3 forcefield [45] to a heavy-atom RMSD convergence of 0.3 Å. This receptor structure, paired with individual ligands from the analog ensemble, underwent Prime protein–ligand complex refinement [34], using the OPLS3 forcefield and the VSGB solvation model [46]. Sidechains with atoms within 5 Å of any of the ensemble members (residues 58–60, 100, 156, 349–357, 374–389, 397, 400–402, 414, and 416 in 4TZ4) were refined with one pass of local optimization sampling, in the default environment with a dielectric of 80.0. Glide redocking was done with the set of compounds that were able to fit into the cereblon binding site without inducing large-scale unfolding or repacking of the binding side: C2–C46, C55, C72, C77, C83, C86, C93, Gu973, Gu998, Gu992, and Gu1029. Each ligand was redocked into its own induced-fit receptor structure. Glide grids were generated for each receptor structure with default parameters. Glide extra-precision (XP) docking runs sampled nitrogen inversions and ring conformations, penalizing nonpolar amides, with no core matching or constraints [47]. For each ligand, the pharmacophore-aligned conformation was docked as well as the set of original conformers generated by LigPrep. Docked poses were scored using the Prime induced-fit score of XPscore + 0.05 × PrimeEnergy [35]. 3.5. Pharmacophore Analysis A pharmacophore model was developed in Phase [48] based on the set of prealigned thalidomide analogs that were used for docking. Active compounds were defined as those with less than 75% outgrowth in the RAR assay. Pharmacophore hypotheses were required to match at least 50% of actives, and to have 3-5 features. Default feature types were used. Since as expected, no single hypothesis matched all the actives, the “Detect Binding Modes” tool was used to perform hierarchical clustering on both the active compounds and the hypotheses, each represented by bit strings. Clustering into 2, 3, or 4 possible binding modes was tested, and results were that active compounds could be consistently grouped into matches with two binding mode pharmacophore hypotheses. 4. Conclusions As a prodrug, thalidomide is appreciably more active in vivo than conveyed by in vitro data [23,38]. Given that thalidomide’s activity is in part due to its biologically active hydroxylated metabolite, it is of interest to identify the physicochemical features of potent thalidomide analogs that do not require hepatic activity. To identify such features, 27 chemically diverse thalidomide analogs were screened in the rat aorta ring (RAR) assay; 16 of these compounds exhibited inhibition of microvessel outgrowth, with seven of these being statistically significant (compounds C4, C46, C86, Gu973, Gu992, Gu998, and Gu1029; Figure 2). Here we analyzed the structural requirements for activity in this set of thalidomide analogs. Molecular modeling approaches have previously been used to examine the associations between the chemical structures of thalidomide analogs and their biological activities. A previous 3D-quantitive structure–activity relationship (3D-QSAR) study on thalidomide analogs used comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) procedures to identify structural components responsible for antiangiogenic activity [49]. Within this study, 29 structurally related thalidomide analogs were first examined for antiangiogenic activity in the RAR assay, and inhibition values were then used as the basis for CoMFA and CoMSIA analyses. These 3D-QSAR studies identified physicochemical features that affected antiangiogenic activity. In particular, perpendicular alignment of the glutarimide ring to the phthalimido ring, hydrophobicity, and steric bulk around the phthalimido ring, and hydrophilic regions within the side chain were found to promote activity; conversely, steric bulk and hydrogen bond donor groups around the glutarimide ring were identified to decrease activity [49]. These models offered insight into the structural requirements for antiangiogenic activity, speculative knowledge of the binding partner(s), and aid in the design of new thalidomide analogs [20,50,51,52]. Similar CoMFA and CoMSIA approaches were applied to identify trends in the properties of thalidomide analogs with anti-inflammatory activity [53]. 3D-QSAR approaches have the fundamental baseline assumption that all compounds are binding to the same target and making the same set of interactions, and further assumes that the correct alignment and superposition of the compounds can be correctly deduced [54,55]. However, trends in the antiangiogenic and anti-inflammatory potencies have led to speculation that more than one target may be mediating the effects of these thalidomide analogs [20,26]. Recently, pull-down assays and X-ray crystallography studies have shown that the target of thalidomide and its close analogs lenalidomide and pomalidomide is the E3 protein ligase complex component cereblon [6,10,11]. Protein docking of analogs can help elucidate if and how they interact with a known target of the parent compound. Using a crystal structure of human cereblon [10], we conducted in silico docking simulations to identify which compounds in the set of 27 thalidomide analogs are likely to bind to cereblon and to gain insight into their binding modes and predicted binding affinities. The small and closely packed nature of the cereblon binding pocket required the development of a hybrid ligand-ensemble and induced-fit docking protocol. This ensured that all compounds could be analyzed in the context of a similar binding mode, as expected for close analogs, and offered a feasible approach to the conformational flexibility of the binding site. Some structure–activity relationships could be observed directly from the RAR assay results. For example, all compounds with thalidomide-like glutarimide rings (C29, C34, C55) were able to reduce microvessel outgrowth (Figure 2). Of particular interest, the glutarimide of C55 is intact but is incorporated at position Z1 (Figure 1) of an isoindoline-1,3-diimine structure instead of substituting the phthalimide nitrogen (Figure 3D). When treated with C55, RAR microvessel outgrowth was reduced to 22.95% (±3.74), suggesting that the functional glutarimide moiety is sufficient for activity, and that modifications of the phthalimide are probably tolerated for cereblon binding. The glutarimide nitrogen of C53 is substituted with a space-filling residue (Figure 1). This compound did not show antiangiogenic activity in the RAR assay and was not predicted by docking to bind to cereblon. Previous studies have shown that thalidomide analogs without an intact glutarimide ring do not exhibit antiangiogenic activity in the RAR assay [49]. From the RAR assay data, it appears that other substructures can also be active. Docking results suggest that barbiturate moieties, as well as lactams and thiolactams, might fit into the cereblon binding site. However, compounds C4 and C9 bearing a piperidine-2,6-dithione, i.e., a doubly thionated glutarimide, along with compounds with large bicyclo or adamantyl rings (Figure 3E,F) were not predicted to be able to bind cereblon well. A recent paper that systematically characterized the cereblon binding abilities of probe compounds with a variety of 5- and 6-membered rings, using FRET pairing to the binding site tryptophans in a bacterial single-domain homolog of cereblon, found that barbiturate rings, as well as some but not all thiocarbonyl-containing 5-membered rings, did not bind cereblon [56]. This suggests that our induced-fit docking, which was designed to be as generous as possible in accommodating nonglutarimide rings into the binding site, may in fact be too lenient, and thus that even more of the thalidomide analogs examined here may be exerting antiangiogenic effects via a different mechanism. A pharmacophore analysis suggested that the strongly antiangiogenic, hydrophobic compounds C4 and C86 are likely to have another target. In fact, eight of the 13 nonfluorinated compounds with antiangiogenic activity in the RAR assay fit this potential target pharmacophore, although some of these compounds are also predicted to bind cereblon. C4 and C86 contain a piperidine-2,6-dithione and an adamantane group, respectively, in place of the glutarimide ring. Interestingly, both adamantane and thioamide derivatives have been reported as ion channel inhibitors [57,58,59]. Adamantane derivatives are largely used as antivirals [60] and in the treatment of Parkinson’s disease by indirect modulation of dopaminergic transmission [61]. All the tetrafluorinated thalidomide analogs (Gu973, Gu992, Gu998, and Gu1029) vastly reduced microvessel outgrowth and were predicted to bind cereblon as strongly as thalidomide via their alkylated barbiturate rings. Differences in the barbiturate ring puckering vs. glutarimide ring might shift the position of the phthalimide ring slightly toward the outside of the cereblon binding site (Figure 3G), and this difference might be magnified with the tetrafluorobenzamide derivatives Gu992 and Gu1029 (Figure 3H). In addition to this, the polyfluorination of the phenyl ring will nearly invert its quadrupole moment and strongly affect its intermolecular interactions [62]. This might suggest that cereblon complexes with tetrafluorinated analogs catalyze the ubiquitination of a new substrate or attract a new binding partner to inhibit angiogenesis by some other mechanism. A second possibility is that the tetrafluorinated analogs have a second target, specific to the fluorobenzene substructure. This is supported by studies with alkyl substituted tetrafluoro-phthalimides and tetrafluoro-benzamides that also show moderate antiangiogenic activity [63]. Overall, the in silico cereblon docking experiments suggest that cereblon binding alone does not account for the antiangiogenic activity of all the tested thalidomide analogs. Our recent findings demonstrate that loss of cereblon does not prevent thalidomide-induced antiangiogenesis, though downstream cereblon targets are upregulated [64]. Thalidomide itself appears to have one target [6], but analogs and metabolites of thalidomide could have multiple targets. Previous theories on thalidomide’s mechanism of action as a multitarget drug [22,65] have in a way become a self-fulfilling prophecy as modifications to the glutarimide moiety, which is directly involved in cereblon-binding of thalidomide, can lead to new target activities. The phthalimide moiety of thalidomide (and the similar isoindolinone moiety of lenalidomide) have been described as privileged scaffolds that can modulate many targets [66,67]. Though structure-based drug design has its limitations, computational modeling remains a powerful tool for preliminary drug screening and the design of new lead series. The herein described modeling studies must now be combined with the experimental determination of the affinity of this type of ligand to cereblon. Corresponding data will enable future work to identify thalidomide analogs that promote anti-inflammatory and antiangiogenic activity while reducing side effects and teratogenic activity, with the aim of making safer, more potent compounds available for clinical use. Acknowledgments The content of this publication does not necessarily reflect the views or policies of the Department of Health and Human Services, nor does mention of trade names, commercial products or organizations imply endorsement by the US Government. Sample Availability: Samples of the compounds are available from the authors. Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. Supplementary Materials The following are available online, Table S1: Induced-fit docking scores for compounds able to fit in the cereblon binding site. Click here for additional data file. Author Contributions Conceptualization, M.L.P., C.H.C., W.D.F. and M.C.N.; methodology, formal analysis, investigation, M.L.P., S.L.B., S.M., C.H.C., C.S. and M.K.C.; resources, N.H.G., M.G. and W.D.F.; writing—original draft preparation, S.M., C.H.C. and M.L.P.; writing—review and editing, M.L.P., S.L.B., C.H.C., M.K.C., S.M., W.L., D.T., C.S., N.H.G., M.G., N.V., M.C.N. and W.D.F.; supervision, W.D.F. and M.C.N.; funding acquisition, W.D.F., M.C.N. and N.V. All authors have read and agreed to the published version of the manuscript. Funding This research was supported in part by the Intramural Research Program of the Center for Cancer Research, National Cancer Institute (ZIA SC006538); in part with Federal funds from the Frederick National Laboratory for Cancer Research, National Institutes of Health, under contract HHSN261200800001E; the Intramural Research Program of the National Institute on Aging, National Institutes of Health; and a Wellcome Trust-NIH PhD Studentship to SB, WDF, and NV (Grant number 098252/Z/12/Z). Conflicts of Interest Authors W.D.F., N.H.G., W.L., S.B., N.V., D.T. have patents on compounds assessed in this study. Patent application No. US2016/054430. The authors declare no conflict of interest. Figure 1 Structures of thalidomide analogs tested and docked. Scheme 1: structures of thalidomide, lenalidomide, and pomalidomide. Scheme 2 comprises compounds where R1 is H, NO2, NH2, or OH; and R2 is H, NH2, OH, OAc, or NHCH(CH3)2; Z1 is O or CZ1 is CH2. Z3 is O or S; Z4 is O or S, or C-CZ4 is CH=CH. Scheme 3 comprises compounds where R1 is H or NH2; Z1 is O or CZ1 is CH2, or Z2 is iminoglutarimide; ring B is substituted glutarimide or a heterocyclic or polycarbocyclic moiety. Schemes 4 and 5 comprise compounds where R1 and R2 are alkyl groups. Compounds have been previously described: Scheme 2 [18], Scheme 3 [18,19], and Schemes 4 and 5 [20]. Outgrowth from rat aortic rings after a 5-day incubation is given as a percentage of control outgrowth. SEM is the standard error of the mean. A minimum of three rings was used per compound. Asterisk (*) denotes chiral carbon. Figure 2 Graphical representation of rat aorta ring (RAR) assay outgrowth. A minimum of three rings were used per treatment. Error bars are SEM. Statistical analysis was performed in Prism and is multiple comparisons (Dunnett’s test) one-way ANOVA (* p < 0.05; ** p < 0.01; *** p < 0.001). Figure 3 Thalidomide analog docking to cereblon. Hydrogen bonds are indicated by dashed orange lines, and aromatic ring stacking interactions are indicated by dashed light blue lines. (A) lenalidomide; (B) compound C4; (C) compound C44; (D) compound C55; (E) compound C83; (F) compound C86; (G) compound Gu998; (H) compound Gu992. The glutarimide ring makes a set of conserved hydrogen bonds to Trp 380 N and His 378 Nδ and O, and the phthalimide ring can hydrogen bond to Asn 351. Figure 4 Plot of induced-fit docking score against % outgrowth in the RAR assay. The induced-fit docking score is a function of conformational energy and nonbonded interactions in the receptor–ligand complex. Outliers from the general trend of increased cereblon binding (lower docking score) leading to decreased outgrowth are circled in red. Figure 5 Pharmacophores for two calculated target binding modes. Acceptor site points are colored pink, donors are colored light blue, hydrophobic regions are green, and aromatic rings are orange circles. (A) A potential non-cereblon binding pharmacophore for compounds with hydrophobic rings. (B) A cereblon-binding pharmacophore for glutarimide-containing compounds. ==== Refs References 1. Aragon-Ching J.B. Li H. Gardner E.R. Figg W.D. Thalidomide Analogues as Anticancer Drugs Recent Pat. Anticancer Drug Discov. 2007 2 167 174 10.2174/157489207780832478 17975653 2. Millrine D. Kishimoto T. A Brighter Side to Thalidomide: Its Potential Use in Immunological Disorders Trends. Mol. Med. 2017 23 348 361 10.1016/j.molmed.2017.02.006 28285807 3. Sherbet G.V. Therapeutic Potential of Thalidomide and Its Analogues in the Treatment of Cancer Anticancer Res. 2015 35 5767 5772 26503997 4. Lepper E.R. Smith N.F. Cox M.C. Scripture C.D. Figg W.D. Thalidomide Metabolism and Hydrolysis: Mechanisms and Implications Curr. Drug Metab. 2006 7 677 685 10.2174/138920006778017777 16918319 5. Paravar T. Lee D.J. Thalidomide: Mechanisms of Action Int. Rev. Immunol. 2008 27 111 135 10.1080/08830180801911339 18437602 6. 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 10.1126/science.1177319 20223979 7. Ito T. Handa H. Cereblon and Its Downstream Substrates as Molecular Targets of Immunomodulatory Drugs Int. J. Hematol. 2016 104 293 299 10.1007/s12185-016-2073-4 27460676 8. Kronke J. Fink E.C. Hollenbach P.W. MacBeth K.J. Hurst S.N. Udeshi N.D. Chamberlain P.P. Mani D.R. Man H.W. Gandhi A.K. Lenalidomide Induces Ubiquitination and Degradation of Ck1alpha in Del(5q) Mds Nature 2015 523 183 188 10.1038/nature14610 26131937 9. Matyskiela M.E. Couto S. Zheng X. Lu G. Hui J. Stamp K. Drew C. Ren Y. Wang M. Carpenter A. Sall4 Mediates Teratogenicity as a Thalidomide-Dependent Cereblon Substrate Nat. Chem. Biol. 2018 14 981 987 10.1038/s41589-018-0129-x 30190590 10. Chamberlain P.P. Lopez-Girona A. Miller K. Carmel G. Pagarigan B. Chie-Leon B. Rychak E. Corral L.G. Ren Y.J. Wang M. Structure of the Human Cereblon-Ddb1-Lenalidomide Complex Reveals Basis for Responsiveness to Thalidomide Analogs Nat. Struct. Mol. Biol. 2014 21 803 809 10.1038/nsmb.2874 25108355 11. Fischer E.S. Bohm K. Lydeard J.R. Yang H. Stadler M.B. Cavadini S. Nagel J. Serluca F. Acker V. Lingaraju G.M. Structure of the Ddb1-Crbn E3 Ubiquitin Ligase in Complex with Thalidomide Nature 2014 7512 49 53 10.1038/nature13527 12. Matyskiela M.E. Clayton T. Zheng X. Mayne C. Tran E. Carpenter A. Pagarigan B. McDonald J. Rolfe M. Hamann L.G. Crystal Structure of the Sall4-Pomalidomide-Cereblon-Ddb1 Complex Nat. Struct. Mol. Biol. 2020 27 319 322 10.1038/s41594-020-0405-9 32251415 13. Matyskiela M.E. Lu G. Ito T. Pagarigan B. Lu C.C. Miller K. Fang W. Wang N.Y. Nguyen D. Houston J. A Novel Cereblon Modulator Recruits Gspt1 to the Crl4(Crbn) Ubiquitin Ligase Nature 2016 535 252 257 10.1038/nature18611 27338790 14. Petzold G. Fischer E.S. Thoma N.H. Structural Basis of Lenalidomide-Induced Ck1alpha Degradation by the Crl4(Crbn) Ubiquitin Ligase Nature 2016 532 127 130 10.1038/nature16979 26909574 15. Eriksson T. Bjorkman S. Roth B. Fyge A. Hoglund P. Stereospecific Determination, Chiral Inversion in Vitro and Pharmacokinetics in Humans of the Enantiomers of Thalidomide Chirality 1995 7 44 52 10.1002/chir.530070109 7702998 16. Blaschke G. Kraft H.P. Fickentscher K. Kohler F. Chromatographic Separation of Racemic Thalidomide and Teratogenic Activity of Its Enantiomers (Author’s Transl) Arzneimittelforschung 1979 29 1640 1642 583234 17. Mori T. Ito T. Liu S. Ando H. Sakamoto S. Yamaguchi Y. Tokunaga E. Shibata N. Handa H. Hakoshima T. Structural Basis of Thalidomide Enantiomer Binding to Cereblon Sci. Rep. 2018 8 1294 10.1038/s41598-018-19202-7 29358579 18. Beedie S.L. Rore H.M. Barnett S. Chau C.H. Luo W. Greig N.H. Figg W.D. Vargesson N. In Vivo Screening and Discovery of Novel Candidate Thalidomide Analogs in the Zebrafish Embryo and Chicken Embryo Model Systems Oncotarget 2016 7 33237 33245 10.18632/oncotarget.8909 27120781 19. Luo W. Tweedie D. Beedie S.L. Vargesson N. Figg W.D. Greig N.H. Scerba M.T. Design, Synthesis and Biological Assessment of N-Adamantyl, Substituted Adamantyl and Noradamantyl Phthalimidines for Nitrite, Tnf-Alpha and Angiogenesis Inhibitory Activities Bioorg. Med. Chem. 2018 26 1547 1559 10.1016/j.bmc.2018.01.032 29472124 20. Ambrozak A. Steinebach C. Gardner E.R. Beedie S.L. Schnakenburg G. Figg W.D. Gutschow M. Synthesis and Antiangiogenic Properties of Tetrafluorophthalimido and Tetrafluorobenzamido Barbituric Acids Chem. Med. Chem. 2016 11 2621 2629 10.1002/cmdc.201600496 27805767 21. Gutschow M. Hecker T. Thiele A. Hauschildt S. Eger K. Aza Analogues of Thalidomide: Synthesis and Evaluation as Inhibitors of Tumor Necrosis Factor-Alpha Production in Vitro Bioorg. Med. Chem. 2001 9 1059 1065 10.1016/S0968-0896(00)00323-0 11354662 22. Hashimoto Y. Structural Development of Biological Response Modifiers Based on Retinoids and Thalidomide Mini Rev. Med. Chem. 2002 2 543 551 10.2174/1389557023405576 12370039 23. Ng S.S. Gutschow M. Weiss M. Hauschildt S. Teubert U. Hecker T.K. Luzzio F.A. Kruger E.A. Eger K. Figg W.D. Antiangiogenic Activity of N-Substituted and Tetrafluorinated Thalidomide Analogues Cancer Res. 2003 63 3189 3194 12810647 24. Muller G.W. Corral L.G. Shire M.G. Wang H. Moreira A. Kaplan G. Stirling D.I. Structural Modifications of Thalidomide Produce Analogs with Enhanced Tumor Necrosis Factor Inhibitory Activity J. Med. Chem. 1996 39 3238 3240 10.1021/jm9603328 8765505 25. Muller G.W. Shire M.G. Wong L.M. Corral L.G. Patterson R.T. Chen Y. Stirling D.I. Thalidomide Analogs and Pde4 Inhibition Bioorg Med. Chem. Lett. 1998 8 2669 2674 10.1016/S0960-894X(98)00475-2 9873600 26. Beedie S.L. Peer C.J. Pisle S. Gardner E.R. Mahony C. Barnett S. Ambrozak A. Gutschow M. Chau C.H. Vargesson N. Anticancer Properties of a Novel Class of Tetrafluorinated Thalidomide Analogues Mol. Cancer Ther. 2015 14 2228 2237 10.1158/1535-7163.MCT-15-0320 26269604 27. Bauer K.S. Cude K.J. Dixon S.C. Kruger E.A. Figg W.D. Carboxyamido-Triazole Inhibits Angiogenesis by Blocking the Calcium-Mediated Nitric-Oxide Synthase-Vascular Endothelial Growth Factor Pathway J. Pharmacol Exp. Ther. 2000 292 31 37 10604929 28. Iqbal F. Gratch Y.S. Szaraz P. Librach C.L. The Aortic Ring Co-Culture Assay: A Convenient Tool to Assess the Angiogenic Potential of Mesenchymal Stromal Cells in Vitro J. Vis. Exp. 2017 127 e56083 10.3791/56083 29. Murai T. Kawashita N. Tian Y.S. Takagi T. In Silico Analysis of Enantioselective Binding of Immunomodulatory Imide Drugs to Cereblon Springerplus 2016 5 1122 10.1186/s40064-016-2761-9 27478739 30. Hartmann M.D. Boichenko I. Coles M. Lupas A.N. Hernandez Alvarez B. Structural Dynamics of the Cereblon Ligand Binding Domain PLoS ONE 2015 10 e0128342 10.1371/journal.pone.0128342 26024445 31. Fu D.Y. Meiler J. Rosettaligandensemble: A Small-Molecule Ensemble-Driven Docking Approach Acs Omega 2018 3 3655 3664 10.1021/acsomega.7b02059 29732444 32. Huang S.Y. Li M. Wang J. Pan Y. Hybriddock: A Hybrid Protein-Ligand Docking Protocol Integrating Protein- and Ligand-Based Approaches J. Chem. Inf. Model. 2016 56 1078 1087 10.1021/acs.jcim.5b00275 26317502 33. Lexa K.W. Carlson H.A. Protein Flexibility in Docking and Surface Mapping Q Rev. Biophys. 2012 45 301 343 10.1017/S0033583512000066 22569329 34. Jacobson M.P. Friesner R.A. Xiang Z. Honig B. On the Role of the Crystal Environment in Determining Protein Side-Chain Conformations J. Mol. Biol. 2002 320 597 608 10.1016/S0022-2836(02)00470-9 12096912 35. Sherman W. Day T. Jacobson M.P. Friesner R.A. Farid R. Novel Procedure for Modeling Ligand/Receptor Induced Fit Effects J. Med. Chem. 2006 49 534 553 10.1021/jm050540c 16420040 36. Lopez-Girona A. Mendy D. Ito T. Miller K. Gandhi A.K. Kang J. Karasawa S. Carmel G. Jackson P. Abbasian M. Cereblon Is a Direct Protein Target for Immunomodulatory and Antiproliferative Activities of Lenalidomide and Pomalidomide Leukemia 2012 26 2326 2335 10.1038/leu.2012.119 22552008 37. Matyskiela M.E. Zhang W. Man H.W. Muller G. Khambatta G. Baculi F. Hickman M. LeBrun L. Pagarigan B. Carmel G. A Cereblon Modulator (Cc-220) with Improved Degradation of Ikaros and Aiolos J. Med. Chem. 2018 61 535 542 10.1021/acs.jmedchem.6b01921 28425720 38. Price D.K. Ando Y. Kruger E.A. Weiss M. Figg W.D. 5′-Oh-Thalidomide, a Metabolite of Thalidomide, Inhibits Angiogenesis Drug Monit. 2002 24 104 110 10.1097/00007691-200202000-00017 39. Hoffmann M. Kasserra C. Reyes J. Schafer P. Kosek J. Capone L. Parton A. Kim-Kang H. Surapaneni S. Kumar G. Absorption, Metabolism and Excretion of [14c]Pomalidomide in Humans Following Oral Administration Cancer Chemother. Pharm. 2013 71 489 501 10.1007/s00280-012-2040-6 40. Schumacher H. Smith R.L. Williams R.T. The Metabolism of Thalidomide: The Spontaneous Hydrolysis of Thalidomide in Solution Br. J. Pharm. Chemother. 1965 25 324 337 10.1111/j.1476-5381.1965.tb02053.x 41. Chen N. Zhou S. Palmisano M. Clinical Pharmacokinetics and Pharmacodynamics of Lenalidomide Clin. Pharm. 2017 56 139 152 10.1007/s40262-016-0432-1 42. Kruger E.A. Figg W.D. Tnp-470: An Angiogenesis Inhibitor in Clinical Development for Cancer Expert Opin. Investig. Drugs 2000 9 1383 1396 10.1517/13543784.9.6.1383 43. Halgren T.A. Mmff Vi. Mmff94s Option for Energy Minimization Studies J. Comput. Chem. 1999 20 720 729 10.1002/(SICI)1096-987X(199905)20:7<720::AID-JCC7>3.0.CO;2-X 44. Sastry G.M. Adzhigirey M. Day T. Annabhimoju R. Sherman W. Protein and Ligand Preparation: Parameters, Protocols, and Influence on Virtual Screening Enrichments J. Comput. Aided Mol. Des. 2013 27 221 234 10.1007/s10822-013-9644-8 23579614 45. Harder E. Damm W. Maple J. Wu C. Reboul M. Xiang J.Y. Wang L. Lupyan D. Dahlgren M.K. Knight J.L. Opls3: A Force Field Providing Broad Coverage of Drug-Like Small Molecules and Proteins J. Chem. Theory Comput. 2016 12 281 296 10.1021/acs.jctc.5b00864 26584231 46. Li J. Abel R. Zhu K. Cao Y. Zhao S. Friesner R.A. The Vsgb 2.0 Model: A Next Generation Energy Model for High Resolution Protein Structure Modeling Proteins 2011 79 2794 2812 10.1002/prot.23106 21905107 47. Friesner R.A. Murphy R.B. Repasky M.P. Frye L.L. Greenwood J.R. Halgren T.A. Sanschagrin P.C. Mainz D.T. Extra Precision Glide: Docking and Scoring Incorporating a Model of Hydrophobic Enclosure for Protein-Ligand Complexes J. Med. Chem. 2006 49 6177 6196 10.1021/jm051256o 17034125 48. Dixon S.L. Smondyrev A.M. Rao S.N. Phase: A Novel Approach to Pharmacophore Modeling and 3d Database Searching Chem. Biol. Drug Des. 2006 67 370 372 10.1111/j.1747-0285.2006.00384.x 16784462 49. Lepper E.R. Ng S.S. Gutschow M. Weiss M. Hauschildt S. Hecker T.K. Luzzio F.A. Eger K. Figg W.D. Comparative Molecular Field Analysis and Comparative Molecular Similarity Indices Analysis of Thalidomide Analogues as Angiogenesis Inhibitors J. Med. Chem. 2004 47 2219 2227 10.1021/jm0304820 15084120 50. Godin A.M. Araujo D.P. Menezes R.R. Brito A.M. Melo I.S. Coura G.M. Soares D.G. Bastos L.F. Amaral F.A. Ribeiro L.S. Activities of 2-Phthalimidethanol and 2-Phthalimidethyl Nitrate, Phthalimide Analogs Devoid of the Glutarimide Moiety, in Experimental Models of Inflammatory Pain and Edema Pharm. Biochem. Behav. 2014 122 291 298 10.1016/j.pbb.2014.04.008 51. Matijevic-Sosa J. Cvetnic Z. Antimicrobial Activity of N-Phthaloylamino Acid Hydroxamates Acta Pharm. 2005 55 387 399 16375828 52. Yamamoto T. Shibata N. Takashima M. Nakamura S. Toru T. Matsunaga N. Hara H. Enzymatic Resolution and Evaluation of Enantiomers of Cis-5’-Hydroxythalidomide Org. Biomol. Chem. 2008 6 1540 1543 10.1039/b802459f 18421383 53. Avila C.M. Romeiro N.C. da Silva G.M. Sant’Anna C.M. Barreiro E.J. Fraga C.A. Development of New Comfa and Comsia 3d-Qsar Models for Anti-Inflammatory Phthalimide-Containing Tnfalpha Modulators Bioorg. Med. Chem. 2006 14 6874 6885 10.1016/j.bmc.2006.06.042 16843662 54. Dean P.M. Molecular Similarity in Drug Design Blackie Academic & Professional New York, NY, USA 1995 55. Kubinyi H. Qsar and 3d Qsar in Drug Design.1. Methodology Drug Discov. 1997 2 457 467 10.1016/S1359-6446(97)01079-9 56. Boichenko I. Bar K. Deiss S. Heim C. Albrecht R. Lupas A.N. Hernandez Alvarez B. Hartmann M.D. Chemical Ligand Space of Cereblon ACS Omega. 2018 3 11163 11171 10.1021/acsomega.8b00959 31459225 57. Brunhofer G. Granig W.H. Studenik C.R. Erker T. A Journey from Benzanilides to Dithiobenzanilides: Synthesis of Selective Spasmolytic Compounds Bioorg. Med. Chem. 2011 19 994 1001 10.1016/j.bmc.2010.11.043 21185194 58. Fischer A. Schmidt C. Lachenicht S. Grittner D. Winkler M. Wrobel T. Rood A. Lemoine H. Frank W. Braun M. Synthesis of Benzofuran, Benzothiophene, and Benzothiazole-Based Thioamides and Their Evaluation as K(Atp) Channel Openers Chem. Med. Chem. 2010 5 1749 1759 10.1002/cmdc.201000297 20845362 59. Wanka L. Iqbal K. Schreiner P.R. The Lipophilic Bullet Hits the Targets: Medicinal Chemistry of Adamantane Derivatives Chem. Rev. 2013 113 3516 3604 10.1021/cr100264t 23432396 60. Cady S.D. Wang J. Wu Y. DeGrado W.F. Hong M. Specific Binding of Adamantane Drugs and Direction of Their Polar Amines in the Pore of the Influenza M2 Transmembrane Domain in Lipid Bilayers and Dodecylphosphocholine Micelles Determined by Nmr Spectroscopy J. Am. Chem. Soc. 2011 133 4274 4284 10.1021/ja102581n 21381693 61. Kaefer V. Semedo J.G. Silva Kahl V.F. Von Borowsky R.G. Gianesini J. Ledur Kist T.B. Pereira P. Picada J.N. DNA Damage in Brain Cells and Behavioral Deficits in Mice after Treatment with High Doses of Amantadine J. Appl. Toxicol. 2010 30 745 753 10.1002/jat.1550 20575081 62. Shimizu K. Costa Gomes M.F. Padua A.A. Rebelo L.P. Canongia Lopes J.N. On the Role of the Dipole and Quadrupole Moments of Aromatic Compounds in the Solvation by Ionic Liquids J. Phys. Chem. B. 2009 113 9894 9900 10.1021/jp903556q 19569634 63. Steinebach C. Ambrozak A. Dosa S. Beedie S.L. Strope J.D. Schnakenburg G. Figg W.D. Gutschow M. Synthesis, Structural Characterization, and Antiangiogenic Activity of Polyfluorinated Benzamides Chem. Med. Chem. 2018 13 2080 2089 10.1002/cmdc.201800263 30134015 64. Beedie S.L. Huang P.A. Harris E.M. Strope J.D. Mahony C. Chau C.H. Vargesson N. Figg W.D. Role of Cereblon in Angiogenesis and in Mediating the Antiangiogenic Activity of Immunomodulatory Drugs Faseb J. 2020 34 11395 11404 10.1096/fj.201903060RR 32677118 65. Hashimoto Y. Tanatani A. Nagasawa K. Miyachi H. Thalidomide as a Multitarget Drug and Its Application as a Template for Drug Design Drugs Future 2004 29 383 391 10.1358/dof.2004.029.04.792298 66. Sahn J.J. Su J.Y. Martin S.F. Facile and Unified Approach to Skeletally Diverse, Privileged Scaffolds Org. Lett. 2011 13 2590 2593 10.1021/ol200709h 21513290 67. Sharma U. Kumar P. Kumar N. Singh B. Recent Advances in the Chemistry of Phthalimide Analogues and Their Therapeutic Potential Mini. Rev. Med. Chem. 2010 10 678 704 10.2174/138955710791572442 20402635