1 Repurposing of Ring Framework through TR screening Surendra Kumar ǂ, Cheongyun Jang ǂ, Lalita Subedi, Sun Yeou Kim, Mi-hyun Kim * Gachon Institute of Pharmaceutical Science & Department of Pharmacy, College of Pharmacy, Gachon University, 191 Hambakmoeiro, Yeonsu-gu, Incheon, Republic of Korea ǂ The authors are co-first authors * Author for correspondence
2 Abstract In the drug repurposing approach, the chemically diverse and potentially safe molecules can be explored as therapeutic potential out of those originally targeted indications. However, the intellectual property rights, and competitive over-heating issue from sharing the same drug pipeline information are drawback of the approach. In current paradigm of drug discovery, despite approval of new chemical entity, the finding of a new ring system for drug discovery have stagnated. In addition, there are still rarely used ring systems from FDA approved drug that can be optimized and used for repurposing purpose. In the present work, these rarely used drug ring system were re-purposed with our bidirectional target and ring screening (TRdual screening) approach within our defined objectives. The TR-dual screening suggested RARbeta and cyproheptadine as a pair of target and ring system. The selected ring system was virtually grown with predefined criterion, ranking of molecular docking and top score poses made us select test compounds for synthesis and biological proof. Among the tested compounds, 6c proved RARbeta agonist activity. The concept of TR dual optimization can contribute to picking either a compound having acceptable ADMET profile or a target for a selected drug scaffold. Keywords: TR-dual screening, Scaffold (Ring Framework) Repurposing, ChOS (chemistry-oriented synthesis)
3 Introduction The paradigm of current drug R&D is selection & concentration and current massive investment in drug R&D has not been just massive but efficiently massive such as described in 5R-framework of Astrazeneca(1). In particular, Go/No Go decision in R&D stages means the willing to choose a more promising pipeline within available resource(2). However, the strategy can produce neither satisfied success rate nor a chance investigating black box information of the pipeline during R&D. Data generated at drug R&D process can t perfectly support assumed proof of concept (POC) to show the limitation in controlling hidden black box eg., exposure concentration of a drug near to a therapeutic target, POC, therapeutic window in clinical level (1). To overcome the productivity crisis, drug repurposing approaches have been developed. The chemically diverse and potentially safe molecules can be explored as therapeutic potential out of those originally targeted indications. However, the intellectual property rights, and competitive over-heating issue from sharing the same drug pipeline information are drawback of the approach. The described status of drug R&D made us highly motivated in thinking how to control the black box for overcoming current failure of drug R&D. Even though any researcher surely can t propose a breakthrough paradigm or strategy for controlling the black box, trials discriminative from current approaches are required and someday researchers can find any breakthrough or advance from the accumulated trials. As one discriminative drug discovery approach, our research group proposed chemistry-oriented synthesis (ChOS), the approach starting drug discovery from an unprecedented drug scaffold rather than from a specific target (3, 4). Through the ChOS approach, drug discovery driven by the novelty of a drug scaffold can compensate the general approach driven by a specific target to broaden the drug space of artificial drugs(4). In the view of selection & concentration, general approach focuses on the target and ChOS approach focuses on drug. After the study, now we really wonder transferring the point of view can make our discriminative insight and methods close to black box information of drug R&D. In the history of drug discovery, until 1980, synthetic chemists acquired novel drug scaffolds through serendipitous findings & intensive investigation during their active synthesis(5, 6). However, since the golden age, in the spite of annual approval of around 10 NMEs, it was reported that developing an unprecedented drug scaffold or a novel drug ring system, and finding a new ring for drugs have stagnated(7). In the analysis on frameworks/ring
4 systems/rings of FDA approved drugs, 237 ring systems were very rarely used (less than three times) among every ring system in FDA orange book. From the report, we thought (1) how much a ring system can make an effect on drug-like property of a whole drug structure, (2) how about repurposing a successful ring system from the indication of its approved drug to another indication, and (3) how to elicit a new promising drug scaffold from the successful ring system with desirable PK property. Even though a research corresponding to the described questions can t be an efficient strategy in drug R&D now, data generated from the research can be discriminating with the data driven from general approach. In this study, we chose one ring system through our bidirectional screening, target (T) & ring system (R) dual screening: TR- dual screening, generated virtual library of the chosen ring system, conducted chemo-centric target screening of the virtual library, and repurposed the chosen ring system from histamine receptor (anti-allergic) to retinoic acid receptor beta (anti-cancer/inflammation) with DMPK measurement. Figure 1: The general workflow for Target (T) Ring (R) screening.
5 Results Descriptor selection for ring system repurposing: For the scalable & bi-directional TR screening, descriptors were assigned to reported ring systems (349) and the descriptors were used to reduce the number of ring system. Properties of 349 ring systems extracted from every FDA-approved drug in Orange book (7) were calculated through CDK toolkit to analyze (1) fragment complexity, (2) hydrogen bond acceptors, (3) hydrogen bond donors, (4) molecular weight, and (5) VABC volume descriptor [ref]. Among all the properties calculated, VABC descriptor represents the approximate van der Waals volume (Å 3 /molecule) of molecules and can be calculated by the following formula: VvdW= all atom contributions 5.92NB 14.7RA 3.8RNR NB is the number of bonds, RA is the number of aromatic rings, and RNA is the number of nonaromatic rings. From 349 ring systems with less than three frequencies in clinically approved drug, 114 ring systems were chosen based on VABC Volume Descriptor > 140 (8) and 71 ring systems remained after the summation of HA+HB < 3. In this study, 3D structure of a core (a ring system) rather than side chains (substituents) was focused for maximizing the effect of a ring system on drug-target interaction. According to the criteria [Frequency in Orange book 2, MW 500, HA+HD 2, VABC 140] in the Figure 2, 71 ring systems were selected from total 349 ring systems (see supplementary file1; Table S1). Figure 2: Ring reduction workflow.
6 Target selection for ring system repurposing: In PDB bank, 38,529 targets were filtered under the conditions of (1) more than five available PDBs of a target, (2) existence of a ligand in the PDBs, (3) carbon construction of a ligand in the PDBs, and (4) desirable range of molecular weight (250 to 800) to produce 1,714 targets. The 292 targets were common in TTD database (3,261 targets), and ChEMBL version 21 (6,930 targets), and the refined PDB (1714 targets) so that every ChEMBL ligands of the 292 targets were collected and manipulated. Among 292 targets after the manipulation, the number of target was reduced into 97 targets (Figure 3). Figure 3: Target reduction workflow. Target (T)-Ring (R) screening Result: The 3D similarity calculation between 71 ring systems (query) and 3424 ligands of PDBs produced 246,512 paired data with several scores from 3 metrics and 2 color/shape (Figure 4). Expectedly, the complexity of the bi-direction screening as well as narrow distribution and deviation of the 3D similarity scores couldn t make us choose one target-ring system pair. In sequence, we chose the best target against each ring system (97 rows, see supplementary file1; Table S2) and the best ring system against each target (71 rows, see supplementary file1; Table S3) instead of choosing one uncertain best pair. PDBs annotated in the two best lists (best target, best ring system) were collected and after the annotation of accumulated
7 count of duplicated PDBs, the duplicated was removed to present 131 PDBs (97 target) for second screening (Docking based TR screening). Figure 4: 3D-similarity calculation between x-ray ligands and Rings. Docking based TR screening result: To select the best pair of target and rings, a second screening comprised of molecular docking was performed on 71 ring (R) system and 131 PDBs (97 target (T) ). Due to missing atom types for ring ( _2) system, it was excluded during docking. The Standard Precision (SP) and Xtra Precision (XP) mode is employed during the binding mode search. After docking, Uniprot accession id was assigned to PDBs and single entry among homologous protein was retained. A matrix of (97 row x 70 column) was prepared (see supplementary file 2) that consists of docking score values. The docking score values for each targets and each rings were transformed into ranking values and heatmap based on ranked targets and rings were created. The heatmap present the data in different colour intensity with high ranked as dark blue and low with yellow color. The heatmap for ranked target and ranked ring is shown in Figure 5 and 6.
8 Figure 5: Heatmap based on target (97) ranking according to docking score against Rings. Figure 6: Heatmap based on Ring (69) ranking according to docking score against target.
9 Overview on the selected target and ring system: The selected target P10826 belongs to the family of nuclear receptor denoted as RAR β (Retinoid Acid Receptor Beta). The natural modular of this receptor is retinoic acid, which specifically binds to retinoid acid receptor and then forms heterodimers with other receptors such as estrogen receptor alpha and AP-1 receptor and further activates downstream to regulate cell differentiation(9, 10). The retinoic acid receptor is divided into three classes of alpha, beta, and gamma, of which beta has been shown to have anti-cancer effects in epithelial cells. Multiple studies reported that RARβ play a role in breast cancer cell metastatic process and found to interact with ATRA (all trans retinoic acid)(11). In addition, the neuroprotective effect of retinoic acid as well as cancer cells is now being revealed(12). The selected ring system ( _1) has scaffold that abundantly found in ant-histamine drugs and known as cyproheptadine(13) but never was considered for RAR β in literature study. The selected ring system is rigid without rotatable bond. The nitrogen as the hetero atom is unique and it is easy to introduce the functional groups. Virtual library generation and molecular docking: Since the x-ray crystal structure in 4JYG has carboxyl functional group and it makes H-bond with Ser280 and Arg269 amino acid residues of retinoic acid receptor, thus the reagent library with carboxyl functional was selected to create the virtual libraries. During the virtual library generation, the cyproheptadine was set to core structure in demethylation state and libraries were grown to amine position. Additionally, the druggability of every compound in the virtual library was predicted to filter out undesirable compounds according to the pre-defined criteria and 4954 compounds were retained. These compounds were further docked into the selected target protein (PDB ID: 4JYG). During docking, the energy window for ring sampling was increased to enhance the conformer sampling and write out at most 2000 poses per ligand. From all the docking results, the poses were ranked according to docking scores and top scoring poses selected for further analysis. Comparing the binding mode of top ranked poses from virtual library and crystal bound ligand, a structural complementary in their binding mode was found (Figure 7, Table1). Among all the poses, compound 6c exhibit maximum binding affinity in terms of docking score of kcal/mol. Similarly, the compound 6b shows at par binding affinity of In addition, compound 8a, 6d, 6e, 6a, shows lower binding affinity and 6f shows lowest binding affinity among selected poses.
10 Figure 7: (A) The binding mode and superimposition of docked x-ray ligand and 6c (B) The structural superimposition of x-ray ligand and 6c Table 1: The selected N-arylated cyproheptadine derivatives. N R 1 R 5 R 2 R 4 R 3 Compound ID R 1 R 2 R 3 R 4 R 5 6a H COOH CH3 H H 6b H CH3 COOH H H 6c H H COOH H H 6d NO2 H COOH H H 6e H COOH H H H 6f COOH H H H H 8a NH2 H COOH H H 8b H COOH H NH2 H Discussion: The ring or ring systems plays an important role in medicinal chemistry either in the area of combinatorial chemistry, where ring systems are used as central scaffold or in scaffold hoping,
11 where central core is replaced with bioisosteric ring system. In bioactive molecules, heteroaromatic ring contributes approximately 97% and least by aliphatic ring (65%). Almost 73.3% of bioactive molecules have fused aromatic ring system ranging from one to three simple five or six membered rings(14). Despite having such diversity, the substituents grown on these cores have drug like properties. Any drug repurposing methodology can be challenged by intellectual property rights of drug molecules that limit the re-purposed drug molecule for its reuse in different target indication. The challenges can be overcome by the proposed bi-directional target and ring screening methodology. Rings and Target Selection for repurposing: A total of 349 ring systems were extracted from every FDA approved drug in orange book and based on predefined filter 71 rings were obtained. These rings were diverse in nature ranges from two to six fused ring. These fused ring system, open wide possibility for site modification, however considering the synthetic accessibility can further guide for optimization with the domain of application. Similarly, from 38,529 targets using predefined filtering criteria a total of 292 targets (3424 proteins) were retained and these targets cover wide range of biological activity. Bi-directional Target (T)-Ring(R) screening: Due to large number of protein structures, the initial bi-directional Target-Ring docking was not feasible, so a 3D similarity score based screening was performed. The initial screening result produces two list (best target and best ring system), which were accumulated and after removing the duplicates 97 targets (131 PDBs) were retained. The reduced targets and rings were further subjected for second target and ring screening based on docking. The bi-directional target and ring docking produces a matrix of docking score. Since our goal was to select best ring and target against each target and ring respectively, a heatmap was created. The heatmap has efficiently discriminated the poorly ranked targets for each rings and poorly ranked rings for each target, which were excluded from best pair selection criterion. After detailed analysis on two heatmaps (ranked targets and ranked rings), a matrix of best targets and best rings was created for selection of best pair of target and ring for repurposing. To select the best pair from matrix, the docking score is chosen as final selection criterion. After analysis of the docking matrix, it was found that, in most of the high docking
12 score T-R pairs; the target has similar x-ray crystal ligand as with docked ring. Thus these pairs in spite of having higher docking score were excluded to reduce the biasness in selection of best T-R pairs. Based on final selection criterion, the best T-R pair of P10826 (target) and _1 (ring) was chosen for repurpose. Molecular docking of virtual library of cyproheptadine: The binding pocket of 4JYG is hydrophobic in nature, as majority of amino acid belongs to aliphatic or hydrophobic amino acid residues. Despite having hydrophobic interaction of the drug molecules, the binding affinity is further strengthen by H-bond from hydrophilic amino acid residues. The list of binding amino acid residues for selected poses is shown in Table 2. In the binding pose, the compound 6b, 6c both makes the π-π interactions with Phe279 and hydrogen bond with Arg269 and Ser280 amino acid residues respectively. Comparing the binding affinity with, compound 6d and 8a, it is revealed that, despite having same functional group on the core, these compounds have lower binding affinity. The nature of binding pocket with in 3Å is Aromatic hydrophobic (Phe192, Trp218, Phe221, Phe279, Phe295), Aliphatic hydrophobic (Leu224, Ala225, Leu259, Leu262, Ile263, Ile266, Leu298, Ile380, Val288, Leu391, Ile403, Met407), Polar, uncharged (Cys228, Ser280) and positively charged (Arg269) amino acid residues. The hydrophobically enriched binding pocket is complementary with bound ligand and imparts the binding affinity, however, presence of substituents and its orientation contributes to the differences in binding affinity observed. Likewise, the carboxylic acid at R3 position imparts better binding affinity as compound containing such groups makes H-bond with Arg269 and Ser280 that stabilizes in the binding pocket. Changes in the position of carboxylic acid, greatly reduces the binding affinity (i.e. compounds 6a, 6e, and 6f). Moreover, the presence of polar groups like -NO2, or -NH2 at R1 position did not increase the binding affinity. Similarly, the low binding affinity has been shown by compound 6f, which contains carboxylic acid at the same position. Thus, it is quite evident from docking studies that, presence of hydrophobic groups increases the binding affinity, that further strengthen by the suitable positioning of hydrophilic groups.
14 Repurposed the chosen ring system from histamine receptor (anti-allergic) to retinoic acid receptor beta (anti-cancer/inflammation) with DMPK measurement: Among all the synthesized compounds, based on good docking scores, some compounds were selected for RAR β agonist activity (Figure 8). In all the tested compounds, most positions of substituents on phenyl ring at N-1 of cyproheptadine core were well tolerated, however, substitution of COOH group at R3 position has shown good RAR β agonist activity. Compound 6c has shown comparable agonist activity with reference compound RA (all trans-retinoic acid) and highest among all the tested compounds. Similarly, compound 8a also showed good activity as compared to rest of the tested compound. If we compare compound 6c (R1=H), 6b (R2=CH3), 6d (R1=NO2), and 8a (R1=NH2), all have docking scores of kcal/mol, kcal/mol, kcal/mol, and kcal/mol respectively, however, a very large difference in biological activity has been observed. Moreover, it seems, the common -COOH functionality at R3 position is complementary to binding amino acid residues (i.e. Ser280, Arg269) however the presence of other substituents, and their juxtaposition in the binding pocket might have caused such large difference in biological activity. Figure 8: (A) The human RAR β reporter bioassay for synthesized compound, (B) The concentration dependent RAR β bioassay for RA (all trans retinoic acid) and 6c.
15 Conclusion: Our ChOS approach made us propose the TR-dual screening to select a pair of ring system and target. The selected and identified ring system, cyproheptadine was abundantly found in ant-histamine drugs but was never considered for RAR β in literature study to present how to find a new scaffold from clinically used ring systems. This dual approach i.e. target to drug scaffold and drug scaffold to target can be combined for a novel approach in drug discovery. It is regretful that there are still majority of ring systems from FDA approved drug that are limited in use and can be exploited for drug-repurposing approach. Even though the pair of cyproheptadine and RAR β could be picked to validate our screening, the TR-dual screening can contribute to repurposing for current ring systems and future occurring ring systems (unprecedented drug scaffold). Experimental Section 1. Collection and preparation of dataset: The compound data containing ChEMBL Id and canonical smiles as 1,583,897 single entries was downloaded from the ChEMBL version 21 and imported into KNIME(15). The structural information carrying SMILES was converted and molecular properties were calculated by using CDK toolkit(16). All biological activity data were downloaded additionally and assigned to compound data. Similarly, the protein data containing protein number as 6,930 single entries was downloaded from the ChEMBL version 21 and imported into KNIME. From PDB bank, 219,393 ligands without hydrogen in 3D structures of deposited proteins were downloaded from the protein data bank downloader tool. As a criterion of protein filtration, protein accession id was assigned to PDB database downloaded from Uniprot database. Furthermore, to exclude protein containing marketed drug, protein data containing pharmaceutical development step information was downloaded from Therapeutic Target Database(17). The above following steps ensured in the final selection of research/clinical trial targets. 2. Preparation of ligands: Input to glide docking needs compounds in 3D type format. The conversion of 2D-to-3D was accomplished with the Schrodinger LigPrep program(18). The LigPrep generates multiple states tautomers, stereoisomers, desalt, and ionization at a ph range (7 ± 2) using Epik(19). Each ligand s specified chirality were retained and generated 32 conformers unit per ligand.
16 All 3D conformer structures were energy minimized with the OPLS_2005 force field(20). 3. Preparation of proteins: All protein structures were prepared using Protein Preparation Wizard(21). Water molecules were deleted, bond orders were assigned, hydrogens were added. Hydrogen bond was assigned using PROPKA (ph 7.0). Protein structures were refined through restrained minimization until converges within 0.30 Å RMSD for heavy atoms. 4. Prediction of ADMET properties Absorption, Distribution, Metabolism, Excretion and Toxicity properties of compounds were predicted by using QikProp module of Schrodinger(22). It predicts properties such as octanol/water partition, log BB, overall CNS activity, Caco-2 and MDCK cell permeability, logkhsa for human serum albumin binding and log IC50 for HERG K+ channel blockage. 5. Molecular Docking Simulations: The ligands were docked using SP (Standard Precision) and XP (Xtra Precision) protocol of Glide(23). The receptor grid had to be generated from protein structure before docking the ligands. An outer box was used for assigning ligands and an inner box was used to find a center site point of ligand in receptor grid generation. The inner boxes and outer boxes were set to the lengths of 10Å x 10Å x 10Å and original x-ray ligand respectively. The scaling factor of van der Waals radius was set to 1.0 and partial charge cutoff was set to The centroid of the outer box was selected by considering the original x-ray ligand and ligands were docked flexibly. In docking, the default setting conditions were applied. The post-docking minimization was performed with the threshold for rejecting minimized pose and number of pose per ligand to included was set to 0.50 kcal/mol and 10 respectively. 6. 3D similarity calculation 3D similarity calculation was performed using OpenEye s ROCS software(24, 25). The ROCS is a popular tool because of its efficient and effectiveness. Before running ROCS, ligand structures in x-ray protein set was filtered for druggability and it was accomplished by assigning TTD database information to ligand entries and excluding entries label with successful target. Additionally, inorganic compounds and compounds with molecular weight beyond the range of 250 to 800 were excluded. Similarly, the virtual libraries were filtered on the basis of predicted molecular properties (HERG score < -5, Oral > 80%, Rule of five < 2, Rule of three < 1, and protonated amide compound was deleted). The ROCS was designed to perform the Shape based similarity search and it was conducted by comparing the shape of
17 conformers of x-ray ligand as a reference set and ring structure as well as virtual libraries as query molecules(26,27). Shape similarity calculation was based on the Tanimoto combo score and during calculation a Tanimoto cut off score was set 0.1 to also include low score result. The maximum result number was set to 500 based on the high order of the Tanimoto combo score. For each similarity scores calculated for each x-ray ligands, the basic statistical values, such as the maximum, minimum and mean values were calculated. From shape similarity score matrix Ref Tversky Combo Max, Tanimoto Combo Mean, Ref Tversky Combo Mean, Tanimoto Combo Max, Ref Tversky Combo Min, Color Tanimoto Min, count scores were multiplied to each other, resulting in new scores (fused score) and sorted accordingly. Fused score = RefTverskyCombo_Max TanimotoCombo_Mean RefTverskyCombo_Mean TanimotoCombo_Max RefTverskyCombo_Min ColorTanimoto_Min count 7. Target (T) & ring system (R) dual screening In the adjusted scale of both targets and ring systems, a bi-directional TR dual screening (ring system based target screening, target system based ring screening) were conducted. As a preliminary screening, bi-directional docking was neither facile nor productive due to the size of the system (M: the number of PDBs on 97 targets = 3424 PDBs, N: the number of conformation of 71 ring systems = one conformer). Therefore, 3D similarity calculation of M x N matrix was expected for high speedy performance. At that time, chemo-centric assumption can permit the converted information from PDB structures of 97 targets into ligand conformers of the PDBs. The initial 3D similarity based screen helps to reduce the large pool of PDBs and feasible for bi-directional docking. 8. Virtual library design The virtual library was generated using combiglide (28). A reagent database (virtual library) limited to aryl or vinyl bromide and commercial available library in ZINC database was generated from Alfa-aesar database with multiple states tautomers, stereoisomers, ionization at a ph range (7 ± 2). The generated virtual libraries were untangled and minimized. 9. Synthetic chemistry General procedure of chemistry: All reactions were performed under nitrogen atmosphere
18 (anhydrous conditions). All of the reagents were purchased from Sigma Aldrich, TCI, Alfa Aeser with a purity of 95% or higher and used without further purification. All of the experimental instruments used in the reaction were thoroughly washed and then dried in a dry oven at 70 C. Anhydrous tetrahydrofuran (THF), toluene, diethyl ether (Et2O), tert-butyl methyl ether (TBME) and dimethylformamide (DMF) was purchased from Sigma Aldrich. Iridium catalysts and organolithium catalysts were stored in a refrigerator at -10 C under dark conditions. The reaction pattern was confirmed by UV irradiation, CAM stain (Cerium Ammonium Molybdate stain), PMA stain (Phospho molybdic acid stain) and Ninhydrin Stain on Silica TLC plate (60 F-254). BRUKER Ascend NMR instruments used for 1H-NMR (600 MHz) and 13C-NMR (150 MHz) spectra in CDCl3 solvent with regard to tetramethylsilane as the reference for chemical shift values. The chemical shifts (δ) values were expressed in parts per million (ppm) and coupling constants (J) in hertz (Hz). Synthesis of cyproheptadine derivatives and reagent library: The selected compounds (6a-f and 8a-b) (Table 1) were synthesized starting from 4-(5H-dibenzo[a,d]annulen-5-ylidene)piperidine. The 4-(5H-dibenzo[a,d]annulen- 5-ylidene)piperidine was obtained from commercially available cyproheptadine hydrochloride sesquihydrate (1), which was transformed into cyproheptadine (2) using 10% NaOH in ethanol. The demethylation was afforded by converting 2 into ethyl carboxylate derivative (3) followed by hydrolysis using KOH in ethanol to obtain 4-(5H-dibenzo[a,d]annulen-5-ylidene)piperidine (4) (Figure 9). Figure 9: General scheme for synthesis of 4-(5H-dibenzo[a,d]annulen-5-ylidene)piperidine (4). The obtained 4 is further N-arylated using various functional reagent library of methyl bromo benzoate derivative to obtain the target compound. Figure 10 outlined the reactions and condition for synthesis of target compound and detailed procedure is provided in
19 supplementary file 3. Moreover, due to unavailability of some functional reagent library, it was prepared by performing the nitration of bromo benzoic acid (see supplementary file 3). Synthesis of reagent library and target compound: As outlined in Scheme S1 (supplementary file 3), the reagent library was prepared starting from protecting the carboxylic group through methylation using H2SO4 and MeOH as reagent and solvent. Regardless of the ortho, meta, and para positions of the carboxyl group, the reaction proceeded at a high yield. The HNO3 and H2SO4 were used for nitrification. The Figure 9 outlined the synthesis of target compound. At first step, both Sphos and palladium catalysts were used 10 times more than the conventional reaction to maintain the anhydrous condition. The obtained substituted Nitro derivatives were reduced to amine using iron, and the reaction proceeded at a high yield regardless of the reagent. Finally, the final product was completed through demethylation using KOH. Figure 10: General scheme for synthesis of target compounds (6a-f and 8a-b). 10. Human Retenoic Acid Receptor Beta (RAR β) Reported Assay To measure the human RAR β agonist activity for our compounds, we have performed human RAR β assay (Cat No: #IB , PA, 16801, USA) in a 3*32 assays in a 96 well plate format. The assay was performed according to the kit protocol. Cells were seeded in a 96 well plate that has been provided with kit and treated with various concentrations of negative and positive control. Treated cell plate was incubated into a 37 C, humidified with 5% CO2
20 incubator for about hours. On next day the supernatant was removed and the Luciferase detection reagent was added. The light emission from each assay well has been quantified using luminescence system. The RAR β bioassay result for the synthesized compound is shown in Figure 8. Conflict of interest: The author(s) confirm that this article content has no conflicts of interest. Acknowledgment: This study was supported by the Basic Science Research Program of the National Research Foundation of Korea (NRF), which is funded by the Ministry of Education, Science and Technology (No.: 2017R1E1A1A ). Authors would like to thank OpenEye Scientific Software providing the academic free license. Author contributions: M.-h. K. and C. J. conceived and designed the study. Under M.-h. K. s plan, C. J. and S. K. carried out all modeling & data work and C. J. synthesized every tested compound. Under S.Y.K. guidance, L. S. performed in vitro experiment. M.-h. K., C. J. and S.K. analyzed the data and wrote the manuscript. M.-h. K. and S.Y.K. provided the molecular modeling lab and in vitro research work facility. All authors read and approved the final manuscript. References: 1. Morgan P, Brown DG, Lennard S, Anderton MJ, Barrett JC, Eriksson U, et al. Impact of a five-dimensional framework on R&D productivity at AstraZeneca. Nature reviews Drug discovery Mar;17(3): PubMed PMID: Epub 2018/01/20. eng. 2. Jekunen A. Decision-making in product portfolios of pharmaceutical research and development--managing streams of innovation in highly regulated markets. Drug design, development and therapy. 2014;8: PubMed PMID: Jang C, Yadav DK, Subedi L, Venkatesan R, Venkanna A, Afzal S, et al. Identification of novel acetylcholinesterase inhibitors designed by pharmacophore-based virtual screening, molecular docking and bioassay. Scientific reports /10/08;8(1): Venkanna A, Cho KH, Dhorma LP, Kumar DN, Hah JM, Park H-g, et al.
21 Chemistry-oriented synthesis (ChOS) and target deconvolution on neuroprotective effect of a novel scaffold, oxaza spiroquinone. European journal of medicinal chemistry. 2019;163: Ban TA. The role of serendipity in drug discovery. Dialogues in clinical neuroscience. 2006;8(3): PubMed PMID: Cha Y, Erez T, Reynolds I, Kumar D, Ross J, Koytiger G, et al. Drug repurposing from the perspective of pharmaceutical companies. British journal of pharmacology. 2018;175(2): Taylor RD, MacCoss M, Lawson AD. Rings in drugs. J Med Chem Jul 24;57(14): PubMed PMID: Epub 2014/01/30. eng. 8. Zhao YH, Abraham MH, Zissimos AM. Fast Calculation of van der Waals Volume as a Sum of Atomic and Bond Contributions and Its Application to Drug Compounds. The Journal of Organic Chemistry /09/01;68(19): Rochette-Egly C, Germain P. Dynamic and combinatorial control of gene expression by nuclear retinoic acid receptors (RARs). Nuclear receptor signaling. 2009;7(1):nrs Germain P, Iyer J, Zechel C, Gronemeyer H. Co-regulator recruitment and the mechanism of retinoic acid receptor synergy. Nature. 2002;415(6868): Nadendla E, Teyssier C, Delfosse V, Vivat V, Krishnasamy G, Gronemeyer H, et al. An unexpected mode of binding defines BMS948 as a full retinoic acid receptor β (RARβ, NR1B2) selective agonist. PloS one. 2015;10(5):e Kim JH, Yu KS, Jeong JH, Lee NS, Lee JH, Jeong YG, et al. All-trans-retinoic acid rescues neurons after global ischemia by attenuating neuroinflammatory reactions. Neurochemical research Dec;38(12): PubMed PMID: Epub 2013/10/26. eng. 13. Birknes B. The structure of an antihistamine: cyproheptadine hydrochloride sesquihydrate. Acta Crystallographica Section B. 1977;33(3): Ertl P, Jelfs S, Mühlbacher J, Schuffenhauer A, Selzer P. Quest for the rings. In silico exploration of ring universe to identify novel bioactive heteroaromatic scaffolds. Journal of medicinal chemistry. 2006;49(15): Berthold MR, Cebron N, Dill F, Gabriel TR, Kötter T, Meinl T, et al., editors. KNIME: The Konstanz Information Miner2008; Berlin, Heidelberg: Springer Berlin Heidelberg.
22 16. Steinbeck C, Han Y, Kuhn S, Horlacher O, Luttmann E, Willighagen E. The Chemistry Development Kit (CDK): An open-source Java library for chemo-and bioinformatics. Journal of chemical information and computer sciences. 2003;43(2): Li YH, Yu CY, Li XX, Zhang P, Tang J, Yang Q, et al. Therapeutic target database update 2018: enriched resource for facilitating bench-to-clinic research of targeted therapeutics. Nucleic acids research. 2017;46(D1):D1121-D Schrödinger Release , LigPrep, Schrödinger, LLC, New York, NY, Shelley JC, Cholleti A, Frye LL, Greenwood JR, Timlin MR, Uchimaya M. Epik: a software program for pk a prediction and protonation state generation for drug-like molecules. Journal of computer-aided molecular design. 2007;21(12): Jorgensen WL, Maxwell DS, Tirado-Rives J. Development and testing of the OPLS all-atom force field on conformational energetics and properties of organic liquids. Journal of the American Chemical Society. 1996;118(45): Sastry GM, Adzhigirey M, Day T, Annabhimoju R, Sherman W. Protein and ligand preparation: parameters, protocols, and influence on virtual screening enrichments. Journal of computer-aided molecular design. 2013;27(3): Schrödinger Release , QikProp, Schrödinger, LLC: New York, NY, USA, Friesner RA, Banks JL, Murphy RB, Halgren TA, Klicic JJ, Mainz DT, et al. Glide: a new approach for rapid, accurate docking and scoring. 1. Method and assessment of docking accuracy. J Med Chem Mar 25;47(7): PubMed PMID: Epub 2004/03/19. eng. 24. Hawkins PC, Skillman AG, Nicholls A. Comparison of shape-matching and docking as virtual screening tools. Journal of medicinal chemistry. 2007;50(1): Shape TK (2015) OpenEye scientific software (ver ), Santa Fe, NM. USA. 26. Gadhe CG, Lee E, Kim M-h. Finding new scaffolds of JAK3 inhibitors in public database: 3D-QSAR models & shape-based screening. Archives of pharmacal research. 2015;38(11): Kim H, Jang C, Yadav DK, Kim M-h. The comparison of automated clustering algorithms for resampling representative conformer ensembles with RMSD matrix. Journal of cheminformatics. 2017;9(1): Schrödinger , CombiGlide, Schrödinger, LLC, New York, NY, 2017.
Induced Fit Model For more than 100 years, the behaviour of enzymes had been explained by the "lock-and-key" mechanism developed by pioneering German chemist Emil Fischer. Fischer thought that the chemicals
The Schrödinger KNIME extensions Computational Chemistry and Cheminformatics in a workflow environment Jean-Christophe Mozziconacci Volker Eyrich Topics What are the Schrödinger extensions? Workflow application
Retrieving hits through in silico screening and expert assessment M.. Drwal a,b and R. Griffith a a: School of Medical Sciences/Pharmacology, USW, Sydney, Australia b: Charité Berlin, Germany Abstract:
PHARMA & LIFE SCIENCES WHITEPAPER In Silico Investigation of Off-Target Effects STREAMLINING IN SILICO PROFILING In silico techniques require exhaustive data and sophisticated, well-structured informatics
Virtual Libraries and Virtual Screening in Drug Discovery Processes using KNIME Iván Solt Solutions for Cheminformatics Drug Discovery Strategies for known targets High-Throughput Screening (HTS) Cells
Computational Chemistry in Drug Design Xavier Fradera Barcelona, 17/4/2007 verview Introduction and background Drug Design Cycle Computational methods Chemoinformatics Ligand Based Methods Structure Based
The Schrödinger KNIME extensions Computational Chemistry and Cheminformatics in a workflow environment Jean-Christophe Mozziconacci Volker Eyrich KNIME UGM, Berlin, February 2015 The Schrödinger Extensions
Using Phase for Pharmacophore Modelling 5th European Life Science Bootcamp March, 2017 Phase: Our Pharmacohore generation tool Significant improvements to Phase methods in 2016 New highly interactive interface
Introduction ntochem Providing drug discovery knowledge & small molecules... Supporting the task of medicinal chemistry Allows selecting best possible small molecule starting point From target to leads
Table of Contents Introduction... 2 CSD-CrossMiner Terminology... 2 Overview of CSD-CrossMiner... 3 Searching with a Pharmacophore... 4 Performing a Pharmacophore Search using CSD-CrossMiner Version 2.0
1 As you become familiar or continue to explore the Cresset technology and software applications, we encourage you to look through the user manual. This is accessible from the Help menu. However, don t
Dr. Sander B. Nabuurs Computational Drug Discovery group Center for Molecular and Biomolecular Informatics Radboud University Medical Centre The road to new drugs. How to find new hits? High Throughput
The PhilOEsophy There are only two fundamental molecular descriptors Where can we use shape? Virtual screening More effective than 2D Lead-hopping Shape analogues are not graph analogues Molecular alignment
Ligand Scout Tutorials Step : Creating a pharmacophore from a protein-ligand complex. Type ke6 in the upper right area of the screen and press the button Download *+. The protein will be downloaded and
Using AutoDock for Virtual Screening CUHK Croucher ASI Workshop 2011 Stefano Forli, PhD Prof. Arthur J. Olson, Ph.D Molecular Graphics Lab Screening and Virtual Screening The ultimate tool for identifying
FRAUNHOFER IME SCREENINGPORT Design of screening projects General remarks Introduction Screening is done to identify new chemical substances against molecular mechanisms of a disease It is a question of
Introduction to Chemoinformatics and Drug Discovery Irene Kouskoumvekaki Associate Professor February 15 th, 2013 The Chemical Space There are atoms and space. Everything else is opinion. Democritus (ca.
Similarity Search Uwe Koch Similarity Search The similar property principle: strurally similar molecules tend to have similar properties. However, structure property discontinuities occur frequently. Relevance
Table of Contents Creating a Pharmacophore Query from a Reference Molecule & Scaffold Hopping in CSD-CrossMiner Introduction... 2 CSD-CrossMiner Terminology... 2 Overview of CSD-CrossMiner... 3 Features
Identifying Interaction Hot Spots with SuperStar Version 1.0 November 2017 Table of Contents Identifying Interaction Hot Spots with SuperStar... 2 Case Study... 3 Introduction... 3 Generate SuperStar Maps
Supporting information Sensitive NMR Approach for Determining the Binding Mode of Tightly Binding Ligand Molecules to Protein Targets Wan-Na Chen, Christoph Nitsche, Kala Bharath Pilla, Bim Graham, Thomas
21 th /June/2018@CUGM Progress of Compound Library Design Using In-silico Approach for Collaborative Drug Discovery Kaz Ikeda, Ph.D. Keio University Self Introduction Keio University, Tokyo, Japan (Established
2 Viewing and Analyzing Proteins, Ligands and their Complexes 2 Overview Viewing the accessible surface Analyzing the properties of proteins containing thousands of atoms is best accomplished by representing
In silico pharmacology for drug discovery In silico drug design In silico methods can contribute to drug targets identification through application of bionformatics tools. Currently, the application of
Taller de docking y cribado virtual: Uso de herramientas computacionales en el diseño de fármacos Docking program GLIDE El programa de docking GLIDE Sonsoles Martín-Santamaría Shrödinger is a scientific
University of Silesia, Katowice, Poland 11 22 March 2013 Ligand-receptor interactions Dr. Pavel Polishchuk A.V. Bogatsky Physico-Chemical Institute of National Academy of Sciences of Ukraine Odessa, Ukraine
Structural Bioinformatics (C3210) Molecular Docking Molecular Recognition, Molecular Docking Molecular recognition is the ability of biomolecules to recognize other biomolecules and selectively interact
Applying Bioisosteric Transformations to Predict Novel, High Quality Compounds Dr James Chisholm,* Dr John Barnard, Dr Julian Hayward, Dr Matthew Segall*, Mr Edmund Champness*, Dr Chris Leeding,* Mr Hector
Supporting Material 2-Oxo-tetrahydro-1,8-naphthyridine-Based Protein Farnesyltransferase Inhibitors as Antimalarials Srinivas Olepu a, Praveen Kumar Suryadevara a, Kasey Rivas b, Christophe L. M. J. Verlinde
Hydrogen Bonding & Molecular Design Peter Kenny(firstname.lastname@example.org) Hydrogen Bonding in Drug Discovery & Development Interactions between drug and water molecules (Solubility, distribution, permeability,
Chapter 3 Molecular Modeling 3.1. Introduction In this study pharmacophore models will be created to screen a virtual library of FDA approved drugs for compounds that may inhibit MA-A and MA-B. The virtual
APPLICATION NOTE QSAR Modeling of ErbB1 Inhibitors Using Genetic Algorithm-Based Regression GAINING EFFICIENCY IN QUANTITATIVE STRUCTURE ACTIVITY RELATIONSHIPS ErbB1 kinase is the cell-surface receptor
DS QSAR and Library Design Webinar Practical QSAR and Library Design: Advanced tools for research teams Reservationless-Plus Dial-In Number (US): (866) 519-8942 Reservationless-Plus International Dial-In
SUPPLEMENTARY DATA, FIGURES AND TABLE BIOLOGICAL DATA Spheroids MARY-X size distribution, morphology and drug screening data Supplementary Figure S1: Spheroids MARY-X size distribution. Spheroid size was
EMPIRICAL VS. RATIONAL METHODS OF DISCOVERING NEW DRUGS PETER GUND Pharmacopeia Inc., CN 5350 Princeton, NJ 08543, USA email@example.com Empirical and theoretical approaches to drug discovery have often
Data Quality Issues That Can Impact Drug Discovery Sean Ekins 1, Joe Olechno 2 Antony J. Williams 3 1 Collaborations in Chemistry, Fuquay Varina, NC. 2 Labcyte Inc, Sunnyvale, CA. 3 Royal Society of Chemistry,
Supporting Information - I: Experimental Procedures and Characterization The Direct Reductive Amination of Electron-deficient Amines with Aldehydes: the Unique Reactivity of Re 2 O 7 Catalyst 1 Braja Gopal
Flexibility and Constraints in GOLD Version 2.1 August 2018 GOLD v5.6.3 Table of Contents Purpose of Docking... 3 GOLD s Evolutionary Algorithm... 4 GOLD and Hermes... 4 Handling Flexibility and Constraints
Electronic Supplementary Material (ESI) for Molecular BioSystems. This journal is The Royal Society of Chemistry 2015 Supporting Information Enhancing Specificity in the Janus Kinases: A Study on the Thienopyridine
est Drive K20 GPUs! Experience The Acceleration Run Computational Chemistry Codes on Tesla K20 GPU today Sign up for FREE GPU Test Drive on remotely hosted clusters www.nvidia.com/gputestd rive Shape Searching
CHEM 4170 Problem Set #1 0. Work problems 1-7 at the end of Chapter ne and problems 1, 3, 4, 5, 8, 10, 12, 17, 18, 19, 22, 24, and 25 at the end of Chapter Two and problem 1 at the end of Chapter Three
Data Mining in the Chemical Industry Glenn J. Myatt, Ph.D. Partner, Myatt & Johnson, Inc. firstname.lastname@example.org verview of presentation verview of the chemical industry Example of the pharmaceutical industry
Development of Pharmacophore Model for Indeno[1,2-b]indoles as Human Protein Kinase CK2 Inhibitors and Database Mining Samer Haidar 1, Zouhair Bouaziz 2, Christelle Marminon 2, Tiomo Laitinen 3, Anti Poso
User Guide for LeDock Hongtao Zhao, PhD Email: email@example.com Website: www.lephar.com Copyright 2017 Hongtao Zhao. All rights reserved. Introduction LeDock is flexible small-molecule docking software,
KIME-based scoring functions in Muse 3.0 KIME User Group Meeting 2013 Fabian Bös Certara Mission: End-to-End Model-Based Drug Development Certara was formed by acquiring and integrating Tripos, Pharsight,
Supporting Information Divergent Reactivity of gem-difluoro-enolates towards Nitrogen Electrophiles: Unorthodox Nitroso Aldol Reaction for Rapid Synthesis of -Ketoamides Mallu Kesava Reddy, Isai Ramakrishna,
Pharmacophore-based Drug design Ping-Chiang Lyu Institute of Bioinformatics and Structural Biology, Department of Life Science, National Tsing Hua University 96/08/07 Outline Part I: Analysis The analytical
Structure-Activity Modeling - QSAR Uwe Koch QSAR Assumption: QSAR attempts to quantify the relationship between activity and molecular strcucture by correlating descriptors with properties Biological activity
University of Wollongong Research Online Faculty of Science - Papers (Archive) Faculty of Science, Medicine and Health 2009 NMR study of complexes between low molecular mass inhibitors and the West Nile
Introduction Protein-Ligand Docking Evaluations Protein-ligand docking: Given a protein and a ligand, determine the pose(s) and conformation(s) minimizing the total energy of the protein-ligand complex
Enamine Golden Fragment Library 14 March 216 1794 compounds deliverable as entire set or as selected items. Fragment Based Drug Discovery (FBDD) [1,2] demonstrates remarkable results: more than 3 compounds
COMBINATORIAL CHEMISTRY: CURRENT APPROACH Dwivedi A. 1, Sitoke A. 2, Joshi V. 3, Akhtar A.K. 4* and Chaturvedi M. 1, NRI Institute of Pharmaceutical Sciences, Bhopal, M.P.-India 2, SRM College of Pharmacy,
CDK & Mass Spectrometry October 3, 2011 1/18 Stephan Beisken October 3, 2011 EBI is an outstation of the European Molecular Biology Laboratory. Chemistry Development Kit (CDK) An Open Source Java TM Library
LigandScout Automated Structure-Based Pharmacophore Model Generation Gerhard Wolber* and Thierry Langer * E-Mail: firstname.lastname@example.org Pharmacophores from LigandScout Pharmacophores & the Protein Data
Chapter 2 3-Bromoimidazo [1,2-a]pyridines: Bromination via in situ oxidation of HBr by nitric acid 2.1 INTRODUCTION Halogenated compounds play a very important role in getting higher value added pharmaceutical
COMPUTER AIDED DRUG DESIGN (CADD) AND DEVELOPMENT METHODS DRUG DEVELOPMENT Drug development is a challenging path Today, the causes of many diseases (rheumatoid arthritis, cancer, mental diseases, etc.)
Plan Day 1: What is Chemoinformatics and Drug Design? Methods and Algorithms used in Chemoinformatics including SVM. Cross validation and sequence encoding Example and exercise with herg potassium channel:
Docking Benzamidine Docking to Trypsin Relationship to Drug Design Ligand-based design QSAR Pharmacophore modeling Can be done without 3-D structure of protein Receptor/Structure-based design Molecular
Technical Note Characterization of Eleven 2,5-Dimethoxy-N-(2-methoxybenzyl)- phenethylamine (NBOMe) Derivatives and Differentiation from their 3- and 4- Methoxybenzyl Analogues - Part II Patrick A. Hays*,
Using Bayesian Statistics to Predict Water Affinity and Behavior in Protein Binding Sites Introduction J. Andrew Surface Hampden-Sydney College / Virginia Commonwealth University In the past several decades
Statistical concepts in QSAR. Computational chemistry represents molecular structures as a numerical models and simulates their behavior with the equations of quantum and classical physics. Available programs
MD Simulation in Pose Refinement and Scoring Using AMBER Workflows Yuan Hu (On behalf of Merck D3R Team) D3R Grand Challenge 2 Webinar Department of Chemistry, Modeling & Informatics Merck Research Laboratories,
Early Stages of Drug Discovery in the Pharmaceutical Industry Daniel Seeliger / Jan Kriegl, Discovery Research, Boehringer Ingelheim September 29, 2016 Historical Drug Discovery From Accidential Discovery
Dispensing Processes Profoundly Impact Biological, Computational and Statistical Analyses Sean Ekins 1, Joe Olechno 2 Antony J. Williams 3 1 Collaborations in Chemistry, Fuquay Varina, NC. 2 Labcyte Inc,
Figure 1: Schematic view of the ligand-binding domain from the vitamin D receptor (PDB file 1IE9). The crystallographic waters are shown as small spheres and the bound ligand is shown as a CPK model. HO
Pose and affinity prediction by ICM in D3R GC3 Max Totrov Molsoft Pose prediction method: ICM-dock ICM-dock: - pre-sampling of ligand conformers - multiple trajectory Monte-Carlo with gradient minimization
Structure-Based Drug Discovery An Overview Edited by Roderick E. Hubbard University of York, Heslington, York, UK and Vernalis (R&D) Ltd, Abington, Cambridge, UK RSC Publishing Contents Chapter 1 3D Structure
Introduction Quantitative Structure-Activity Relationships QSPR Quantitative Structure-Property Relationships What is? is a mathematical relationship between a biological activity of a molecular system
Protein structure prediction CS/CME/BioE/Biophys/BMI 279 Oct. 10 and 12, 2017 Ron Dror 1 Outline Why predict protein structure? Can we use (pure) physics-based methods? Knowledge-based methods Two major
Supplementary Note This section contains a detailed description of the chemical procedures and the characterization of products. The text is followed by a reaction scheme explaining the synthetic strategies
ombinatorial hemistry: molecular diversity "Synthesis and pplications of Small Molecule Libraries." Thompson, L..; llman, J.. hem. ev.,, -00. "esign, Synthesis, and valuation of Small-Molecule Libraries.
2013. 54, 6. 1115 1120 UDC 548.737:547.12 CHARACTERIZATION AND CRYSTAL STRUCTURES OF SOLVATED N -(4-HYDROXY-3-NITROBENZYLIDENE)-3-METHYLBENZOHYDRAZIDE AND N -(4-DIMETHYLAMINOBENZYLIDENE)-3-METHYLBENZOHYDRAZIDE
Introduction Chemical properties that affect binding of enzyme-inhibiting drugs to enzymes The production of new drugs requires time for development and testing, and can result in large prohibitive costs
ICM-Chemist-Pro How-To Guide Version 3.6-1h Last Updated 12/29/2009 ICM-Chemist-Pro ICM 3D LIGAND EDITOR: SETUP 1. Read in a ligand molecule or PDB file. How to setup the ligand in the ICM 3D Ligand Editor.
Chapter 4: Amino Acids All peptides and polypeptides are polymers of alpha-amino acids. lipid polysaccharide enzyme 1940s 1980s. Lipids membrane 1960s. Polysaccharide Are energy metabolites and many of
Navigation in Chemical Space Towards Biological Activity Peter Ertl Novartis Institutes for BioMedical Research Basel, Switzerland Data Explosion in Chemistry CAS 65 million molecules CCDC 600 000 structures
Synthetic chemistry ML5 and ML4 were identified as K P.(TREK-) activators using a combination of fluorescence-based thallium flux and automated patch-clamp assays. ML5, ML4, and ML5a were synthesized using
An Integrated Approach to in-silico Screening Joseph L. Durant Jr., Douglas. R. Henry, Maurizio Bronzetti, and David. A. Evans MDL Information Systems, Inc. 14600 Catalina St., San Leandro, CA 94577 Goals
AMRI COMPOUD LIBRARY COSORTIUM: A OVEL WAY TO FILL YOUR DRUG PIPELIE Muralikrishna Valluri, PhD & Douglas B. Kitchen, PhD Summary The creation of high-quality, innovative small molecule leads is a continual