In silico discovery of novel Retinoic Acid Receptor agonist structures
© Schapira et al, licensee BioMed Central Ltd. 2001
Received: 6 April 2001
Accepted: 4 June 2001
Published: 4 June 2001
Several Retinoic Acid Receptors (RAR) agonists have therapeutic activity against a variety of cancer types; however, unacceptable toxicity profiles have hindered the development of drugs. RAR agonists presenting novel structural and chemical features could therefore open new avenues for the discovery of leads against breast, lung and prostate cancer or leukemia.
We have analysed the induced fit of the active site residues upon binding of a known ligand. The derived binding site models were used to dock over 150,000 molecules in silico (or virtually) to the structure of the receptor with the Internal Coordinates Mechanics (ICM) program. Thirty ligand candidates were tested in vitro.
Two novel agonists resulting from the predicted receptor model were active at 50 nM. One of them displays novel structural features which may translate into the development of new ligands for cancer therapy.
The retinoic acid receptors (RAR-α, -β, and -γ) are transcription factors regulating a variety of endocrine metabolic pathways. Unlike anti-estrogens, such as tamoxifen or raloxifene, ligands targeted against the RAR isoforms can present anticancer activity against both estrogen receptor positive and negative breast tumor cells . As a result, such molecules could constitute a novel generation of drugs against breast cancer. For reasons not yet clear, both agonists and antagonists of RAR can present anti-tumor activity against breast, prostate, lung cancer or leukemia [1,2,3,4,5,6,7]. The development of both types of ligands could therefore have important biomedical implications. We have recently demonstrated that antagonists could be discovered rationally, based on a model of the antagonist-bound conformation of the receptor . Our goal here is to discover innovative molecular structures with RAR agonist activity.
We applied a flexible virtual screening algorithm (Molsoft ICM, virtual library screening module ) which rapidly docks hundreds of thousands of flexible compound structures into the ligand binding pocket of RAR, and discovered two novel RAR-β selective agonists. One of these ligands displays original structural and chemical characteristics, which could be used in the development of novel compounds for cancer prevention and therapy.
Results and discussion
We first built a model of the RAR-α agonist binding pocket from the crystal structure of the RAR-γ ligand binding domain (RAR-γ LBD) / all-trans RA complex . All but three amino acids in the vicinity of the ligand are conserved between the two isoforms. These three non-identical residues -A234, M272, and A397- were changed to the RAR-α isoform -S234, I272 and V397- and the energy of the system was minimized (see "Materials and Methods").
A high throughput virtual screening was carried out on the Available Chemicals Directory (MDL Information Systems, San Leandro, CA), a compound structure database of over 150,000 molecules. Each compound was automatically docked into a grid representation of RAR-α, as previously described for Am580, and assigned a score according to the quality of the fit [8, 17, 20]. The 5364 ligand candidates which scored better (i.e. lower) than -32 kcal/mol were preselected for a more refined energy minimization procedure, with flexible receptor side chains, and the energy of complexation was predicted as previously described  (see "Materials and Methods"). After careful visual examination of the 300 candidates displaying the lowest predicted binding energy, 30 molecules were selected and purchased to be experimentally tested in vitro.
The structure of agonist 2 is almost similar to that of Ch55 (Fig. 1e), another RAR-β selective agonist  with 0.44, 0.04 and 1.7 nM Kd for RAR-α, RAR-β and RAR-γ respectively . However, agonist 2 has an additional hydroxyl group on the aromatic ring. The observation that Ch55 displays some potency against the growth of human non-small cell lung carcinoma cells  as well as anti-angiogenic activity  suggests that agonist 2 is a good lead candidate for cancer therapy. Since this compound was discovered through a single virtual library screening, this is a very encouraging result which also constitutes a validation of our method.
The structure of agonist 1 is entirely novel; it is the first RAR ligand described so far with either a thiazole ring or a trifluoro group. It presents a very limited level of flexibility and fits tightly into the receptor's binding pocket. This compound illustrates the benefits of the in silico screening procedure, a rational approach which is based solely on the structure of the receptor, and is not biased towards existing ligands. Indeed, structure-based drug design usually derives novel compound structures from ligands previously described. Our approach enabled the discovery ab initio of RAR agonists displaying some activity at 50 nM with a 7% success rate.
The specificity profile of the two RAR ligands described here are similar to that of Ch55, a compound with ten times higher affinity for RAR-β than RAR-α, and which can block the growth of certain cancer types . It will be interesting to test these molecules for anti-proliferative activity, both in vitro and in vivo. The observation that agonist 1 does not belong to any of the series of RAR ligands described so far makes this compound particularly interesting for further development.
While a model of RAR-α was used to conduct the virtual screening, agonists 1 and 2 have an EC-50 of 200 nM for RAR-β, and an EC-50 of 4 and 2 μM respectively for RAR-α. This observation emphasizes the primary goal of structure-based screening procedures, which is to significantly accelerate the initial steps of lead identification projects by automatically discriminating between binders and non-binders, rather than ranking ligands according to their affinity for the receptor. The discovery of 2 RAR-α agonists out of 30 molecules tested is a good illustration of this approach.
This report details the rapid discovery of RAR agonists with novel structural features, thanks to a powerful virtual ligand screening approach, and a research strategy where considerations on existing ligands are avoided. One of the molecules presented here constitute a good framework for the development of a novel series of RAR ligands very different from all structures described so far. Such ligands could present more favorable specificity and toxicity profiles, and have important applications in cancer therapy.
Materials and Methods
Modeling of RAR-alpha ligand binding pocket
The crystal structure of RAR-γ-LBD complexed to all-trans RA was used as a template  and the three residues in the vicinity of the ligand which are not conserved between the two isoforms were modified accordingly: A234, M272, and A397 were changed to S234, I272 and V397 respectively. The rotation variables of the side chains within 3.5 Å of the modified residues were unfixed and the energy of the system was minimized in the internal coordinate space, according to the ICM method [26, 29].
Docking of AM580 into RAR-alpha
The flexible ligand was docked into a combination of five potential map representations of the RAR-α ligand binding pocket, which account for two Van der Waals boundaries, hydrophobicity, electrostatics and hydrogen bonding profiles [20,27]. This rapid docking procedure was followed by a more refined energy minimization of the complex, with a full atom representation of the receptor, and flexible receptor side chains, according to the ICM stochastic global optimization algorithm [26, 27] as implemented in the Molsoft ICM 2.7 program .
Virtual screening of the compound structure database
The ICM program  was used to perform both receptor modeling and virtual ligand screening. The procedure followed was similar than previously described : each flexible ligand of the Available Chemicals Directory (MDL Information Systems, San Leandro) was docked automatically into the combination of potential maps described above, and assigned a score according to its fit with the receptor. The scoring function included continuum as well as discreet electrostatics, hydrophobicity and entropy parameters . The screening of the database of over 150,000 ligands took less than a month on 10 "194 MHZ IP25" processors. Using the same computing power, the 5364 compounds which scored better (i.e. lower) than -32 were preselected for a second automatic round of selection: they were all docked in two days into a full atom representation of the receptor, with flexible receptor side chains, according to a global energy optimization in internal coordinates [17, 26,27,28]. The structure generated by the initial library screening was used as a starting point for this second, more refined, docking procedure. The binding energy of the compounds was then evaluated with a boundary elements implementation of solvation electrostatics  and the 300 compounds showing the lowest binding energy were selected for further examination. After careful visual inspection for shape complementarity, hydrogen bonding network, compound flexibility, and potential Van der Waals clashes, 30 molecules out of 300 were selected and purchased to be experimentally tested in vitro. We should stress out here that the Van der Waals term is too noisy to be added to the binding energy function, as discussed in details elsewhere . As a consequence, some compounds with a low predicted binding energy were not retained, because they seemed likely to clash sterically with the receptor.
Biological activity of the ligand candidates
HeLa cells were transfected by calcium phosphate precipitation using 1 μg of plasmid expressing the full length receptor isoform and 1 μg of a ΔMTV-IR-CAT reporter gene, as previously described [22, 23]. Ligands were dissolved in DMSO at 20 mM final concentration. Cell cultures were supplemented with indicated ligands immediately after addition of the calcium phosphate/DNA precipitate. Media and ligands were replaced after 24 h and cells were harvested and assayed for CAT activity 24 h later: acetylated and unreacted [14C] chloramphenicol was excised from a thin layer chromatography plate and quantitated in a liquid scintillation counter. The amount of cellular protein extract after two days of incubation was measured to determine the compounds toxicity. Agonist 1 was purchased from Bionet Research (catalog number 1G-433S). Agonist 2 was purchased from Sigma-Aldrich (Sigma Aldrich library of rare chemicals. Catalog number S08503-1).
We would like to thank Maxim Totrov for his help and support throughout the project, Claudio Cavasotto for helping with the manuscript, and the Department of Defense (Grants DAMD179818133, DAMD179919318 and DAMO179919323) for supporting the project.
This work was also supported by NIH grant DK16636 and NYS Empire Award C015710 to HHS.
- Fitzgerald P, Teng M, Chandraratna RA, Heyman RA, Allegretto EA: Retinoic acid receptor alpha expression correlates with retinoid-induced growth inhibition of human breast cancer cells regardless of estrogen receptor status. Cancer Res. 1997, 57: 2642–2645.PubMedGoogle Scholar
- Fanjul AN, Piedrafita FJ, Al-Shamma H, Pfahl M: Apoptosis induction and potent antiestrogen receptor-negative breast cancer activity in vivo by a retinoid antagonist. Cancer Res. 1998, 58: 4607–4610.PubMedGoogle Scholar
- Lu XP, Fanjul A, Picard N, Shroot B, Pfahl MA: Selective retinoid with high activity against an androgen-resistant prostate cancer cell type. Int. J. Cancer 1999, 80: 272–278. 10.1002/(SICI)1097-0215(19990118)80:2<272::AID-IJC17>3.3.CO;2-OView ArticlePubMedGoogle Scholar
- Sun SY, Yue P, Dawson MI, Shroot B, Michel S, Lamph WW, Heyman RA, Teng M, Chandraratna RA, Shudo K, Hong WK, Lotan R: Differential effects of synthetic nuclear retinoid receptor-selective retinoids on the growth of human non-small cell lung carcinoma cells. Cancer Res. 1997, 57: 4931–4939.PubMedGoogle Scholar
- Sun SY, Yue P, Wu GS, El-Deiry WS, Shroot B, Hong WK, Lotan R: Implication of p53 in growth arrest and apoptosis induced by the synthetic retinoid CD437 in human lung cancer cells. Cancer Res. 1999, 59: 2829–2833.PubMedGoogle Scholar
- Mologni L, Ponzanelli I, Bresciani F, Sardiello G, Bergamaschi D, Gianni M, Reichert U, Rambaldi A, Terao M, Garattini E: The novel synthetic retinoid 6-[3-adamantyl-4-hydroxyphenyl]-2-naphthalene carboxylic acid (CD437) causes apoptosis in acute promyelocytic leukemia cells through rapid activation of caspases. Blood 1999, 93: 1045–1061.PubMedGoogle Scholar
- Chandaratna RA: Future trends: a new generation of retinoids. J. Am. Acad. Dermatol. 1998, 39(4 Pt 2): S149-S152.View ArticleGoogle Scholar
- Schapira M, Raaka B, Samuels HH, Abagyan R: Rational discovery of nuclear receptor antagonists. Proc. Natl. Acad. Sci. U S A 2000, 97: 1008–1013. 10.1073/pnas.97.3.1008PubMed CentralView ArticlePubMedGoogle Scholar
- Lippman SM, Benner SE, Hong WK: Retinoid chemoprevention studies in upper aerodigestive tract and lung carcinogenesis. Cancer Res. 1994, 54: 2025s-2028s.PubMedGoogle Scholar
- Costa A, De Palo G, Decensi A, Formelli F, Chiesa F, Nava M, Camerini T, Marubini E, Veronesi U: Retinoids in cancer chemoprevention clinical trials with the synthetic analogue fenretinide. Ann. N. Y. Acad. Sci. 1995, 768: 148–162.View ArticlePubMedGoogle Scholar
- Liu Y, Lee MO, Wang HG, Li Y, Hashimoto Y, Klaus M, Reed JC, Zhang X: Retinoic acid receptor beta mediates the growth-inhibitory effect of retinoic acid by promoting apoptosis in human breast cancer cells. Mol. Cell. Biol. 1996, 16: 1138–1149.PubMed CentralPubMedGoogle Scholar
- Widschwendter M, Berger J, Daxenbichler G, Muller-Holzner E, Widschwendter A, Mayr A, Marth C, Zeimet AG: Loss of retinoic acid receptor beta expression in breast cancer and morphologically normal adjacent tissue but not in the normal breast tissue distant from the cancer. Cancer Res. 1997, 57: 4158–4161.PubMedGoogle Scholar
- Li Y, Dawson MI, Agadir A, Lee MO, Long L, Hobbs PD, Zhang XK: Regulation of RAR-β expression by RAR- and RXR-selective retinoids in human lung cancer cell lines: effect on growth inhibition and apoptosis induction. Int. J. Cancer 1998, 75: 88–95. 10.1002/(SICI)1097-0215(19980105)75:1<88::AID-IJC14>3.0.CO;2-9View ArticlePubMedGoogle Scholar
- Xu XC, Liu X, Tahara E, Lippman SM, Lotan R: Expression and up-regulation of retinoic acid receptor-beta is associated with retinoid sensitivity and colony formation in esophageal cancer cell lines. Cancer Res. 1999, 59: 2477–2483.PubMedGoogle Scholar
- Lin B, Chen G, Xiao D, Kolluri SK, Cao X, Su H, Zhang XK: Orphan receptor COUP-TF is required for induction of retinoic acid receptor beta growth inhibition and apoptosis by retinoic acid in cancer cells. Mol. Cell. Biol. 2000, 20: 957–970. 10.1128/MCB.20.3.957-970.2000PubMed CentralView ArticlePubMedGoogle Scholar
- Sun SY, Wan H, Yue P, Hong WK, Lotan R: Evidence that retinoic acid receptor beta induction by retinoids is important for tumor cell growth inhibition. J. Biol. Chem. 2000, 275: 17149–17153. 10.1074/jbc.M000527200View ArticlePubMedGoogle Scholar
- Molsoft LLC: ICM 2.7 Program Manual. 1998.Google Scholar
- Renaud JP, Rochel N, Ruff M, Vivat V, Chambon P, Gronemeyer H, Moras D: Crystal structure of the RAR-γ ligand-binding domain bound to all-trans retinoic acid. Nature 1995, 378: 681–689. 10.1038/378681a0View ArticlePubMedGoogle Scholar
- Delescluse C, Cavey MT, Martin B, Bernard BA, Reichert U, Maignan J, Darmon M, Shroot B: Selective high affinity retinoic acid receptor alpha or beta-gamma ligands. Mol. Pharmacol. 1991, 40: 556–562.PubMedGoogle Scholar
- Totrov M, Abagyan R: Derivation of sensitive discrimination potential for virtual ligand screening. Proceedings of the 3rd annual international conference on computational molecular biology (RECOMB 99) Lyon France, ACM Press 1999, 312–317.View ArticleGoogle Scholar
- Schapira M, Totrov M, Abagyan R: Prediction of the binding energy for small molecules peptides and proteins. J. Mol. Recognit. 1999, 12: 177–190. 10.1002/(SICI)1099-1352(199905/06)12:3<177::AID-JMR451>3.3.CO;2-QView ArticlePubMedGoogle Scholar
- Qi JS, Desai-Yajnik V, Greene ME, Raaka BM, Samuels HH: The ligand-binding domains of the thyroid hormone/retinoid receptor gene subfamily function in vivo to mediate heterodimerization gene silencing and transactivation. Mol. Cell. Biol. 1995, 15: 1817–1825.PubMed CentralPubMedGoogle Scholar
- Li D, Desai-Yajnik V, Lo E, Schapira M, Abagyan R, Samuels HH: NRIF3 is a novel coactivator mediating functional specificity of nuclear hormone receptors. Mol. Cell. Biol. 1999, 19: 7191–7202.PubMed CentralPubMedGoogle Scholar
- Ogiso Y, Kitagawa K, Nishino H, Iwashima A, Shudo K: Suppression of c-mos expression in teratocarcinoma cells with a new type of inducer of differentiation 35-di-tert-butylchalcone 4'-carboxylic acid. Exp. Cell. Res. 1987, 173: 262–266.View ArticlePubMedGoogle Scholar
- Oikawa T, Hirotani K, Nakamura O, Shudo K, Hiragun A, Iwaguchi T: A highly potent antiangiogenic activity of retinoids. Cancer Lett. 1989, 48: 157–162.View ArticlePubMedGoogle Scholar
- Abagyan RA: Protein structure prediction by global energy optimization Computer Simulation of Biomolecular Systems. In: Theoretical and Experimental Applications Volume 3 (Edited by van Gunsteren WF et al.) Leiden, ESCOM Science Publishers BV 1997.Google Scholar
- Totrov M, Abagyan R: Flexible protein-ligand docking by global energy optimization in internal coordinates. Proteins 1997, Suppl 1: 215–220. Publisher Full Text 10.1002/(SICI)1097-0134(1997)1+<215::AID-PROT29>3.3.CO;2-IGoogle Scholar
- Abagyan R, Totrov M: Biased probability Monte Carlo conformational searches and electrostatic calculations for peptides and proteins. J. Mol. Biol. 1994, 235: 983–1002. 10.1006/jmbi.1994.1052View ArticlePubMedGoogle Scholar
- Abagyan RA, Totrov MM, Kuznetsov DN: ICM - a new method for protein modelling and design applications to docking and structure prediction from the distorted native conformation. J. Comp. Chem. 1994, 15: 488–506.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article: verbatim copying and redistribution of this article are permitted in all media for any purpose, provided this notice is preserved along with the article's original URL.