Structure-based screening of chemical libraries to identify small molecules that are likely to bind with the SET and RING-associated (SRA) domain of Ubiquitin-like, PHD and Ring Finger-containing 1 (UHRF1)
BMC Research Notes volume 13, Article number: 254 (2020)
UHRF1 is a multi-domain protein that recognizes both histone and DNA modification marks on chromatin. UHRF1 is involved in various cellular processes that lead to tumorigenesis and thus attracted considerable attention as a potential anti-cancer drug target. The SRA domain is a unique to the UHRF family. SRA domain recognizes 5-methylcytosine in hemimethylated DNA and necessary for maintenance DNA methylation mediated by DNMT1. Small molecules capable of interacting with the SRA domain may reduce aberrant methylation levels by preventing the interaction of 5-methylcytosine with the SRA domain and thereby blocking substrate access to the catalytic center of DNMT1. The data were collected to identify and predict an initial set of small molecules that are expected to bind to the SRA domain.
Nearly 2.4 million molecules from various chemical libraries were screened with the SRA domain of UHRF1 using Schrodinger’s Small Molecule Drug Discovery Suite. The data is available in the form of a methodology presentation, MS Excel files listing the top hits, and Maestro pose viewer files that provide visualization of how the identified ligands interact with the SRA domain.
UHRF1 functions as an epigenomic controller and is involved in various cellular mechanisms that lead to tumorigenesis . UHRF1 has been proven to increase the activity and specificity of DNMT1 . The SRA domain of UHRF1 is a DNA-binding domain and recognizes 5-methylcytosine (5mC) in hemimethylated CpG dinucleotides [3,4,5,6,7]. Due to the 5mC binding epitope architecture, the SRA domain is a highly promising site for small molecules targeting . The SRA domain of UHRF1 interacts directly with DNMT1 and thereby provides improved substrate (hemimethylated DNA) access to the catalytic center of DNMT1, leading to an increase of DNA methylation activity . In vitro studies have shown that UHRF1 can cause a fivefold increase in DNMT1 activity, and the SRA domain on its own can lead to a 1.9-fold increase in the activity of DNMT1. The interaction between UHRF1 and DNMT1 causes a nearly two-fold increase in the preferential targeting of hemimethylated DNA by DNMT1 . Significantly, the expression levels of UHRF1 were described to be 5- to 70-folds lesser than those observed for HDAC1 and DNMT1 in healthy tissues. Thus, any potential adverse effects that may result due to the inhibition of UHRF1 expression or function are expected to be reasonably manageable when compared with consequences that are caused by the direct inhibition of DNMT1 . Thus, preventing the interaction between the SRA domain and hemimethylated DNA via small molecules is a viable strategy to prevent aberrant DNA methylation . Additional information about targeting the SRA domain for anti-cancer drug development was published earlier .
The identification of small molecules that are predicted to bind to the SRA domain of UHRF1 was performed via virtual screening using Schrodinger’s Small Molecule Drug Discovery Suite. The crystal structures of UHRF1 is available in the public domain. The structure of the SRA domain and its interaction with hemimethylated DNA has been published [3, 5, 6]. The small molecule libraries were screened using the SRA domain (PDB Id: 3DWH) . The downloaded PDB structure was prepared using the protein preparation wizard, which confirmed structural correctness at the start of the screening work. The Asp469 residue, which forms a hydrogen bond with the methylcytosine , was chosen as the active site, and a primary grid was prepared 10 A0 from the Asp469 residue . The other residues that were selected to define the grid were Tyr466 and Tyr478 that sandwich 5-methylcytosine, and also Thr479, that is known to play a crucial role in the preferential recognition of cytosine .
A personal computer with the i7-4700MQ quad-core processor and 32 GB memory was used for this work. The small molecule libraries in the SDF format were prepared with LigPrep, to generate precise 3D molecular models for virtual screening. Epik was utilized for the consistent estimation of pKa values and to return chemically functional structures. The compounds were subjected to a filter to eliminate reactive compounds and analyzed via QIKPROP for the reliable projection of the ADME properties of the small molecules. The structure-based screening was performed using Schrodinger’s virtual screening workflow, which involves sequentially running Glide HTVS, Glide SP, and Glide XP on the prepared compound libraries. The virtual screening workflow removed 90% of the compounds at each phase, thus permitting only the top 10% of the small molecules on to the next step .
Nearly 2.4 million small molecules were screened using the SDF files of compound libraries from ChemDiv (San Diego, CA) and Timtec (Newark, DE). The numbers mentioned in parenthesis is the number of small molecules of the library. TIMTEC’s libraries include the Actimol collection (127,937), HTS part I, and HTS part II (400,000 & 491,349). ChemDiv libraries that were screened were Discovery Chemistry 1, 2 and 3 (350,000, 350,000 and 277,772) and New Chemistry 1 and 2 (250,000 and 206,249). The focused libraries from ChemDiv that were screened include bromodomain (6114), cancer stem cells (19,956), 3D mimetics (9461), soluble diversity (9624), targeted diversity (46,817), and methyltransferase (11,647) libraries. The specific libraries were chosen to facilitate the identification of diverse drug-like molecules that are likely to interact with an anti-cancer drug target with a crucial role in epigenomic regulation.
The data is available in the form of Maestro pose viewer files that is output by Glide. Glide is a sophisticated numerical algorithm optimized for docking accuracy and database enrichment. The pose viewer file contains a set of selected entries in Maestro in which the first entry is the protein (SRA domain), and all the other entries are poses of the docked ligand. After entering into the Pose Viewing Mode, the ligand poses can be navigated. The output files thus provide information about the identified molecules and visualize the predicted interactions with the SRA domain of UHRF1 (Table 1).
The present investigation is limited to the selected small molecule libraries from ChemDiv and Timtec.
The structure-based virtual screening was carried out using most of the default parameters of the Schrodinger’s Small Molecule Drug Discovery Suite.
The small molecule hits that were identified in the present study only narrow down the number of compounds that needs to be evaluated initially in an in vitro assay.
The small molecules identified in this study have not been evaluated in a biochemical or biophysical assay. Some of the identified small molecules may not show a binding response to the SRA domain of UHRF1 in a biochemical or biophysical assay. If a successful binding interaction is detected in an in vitro assay, the molecules need to be validated further in a series of biochemical, biophysical, and cell-based assays.
Ubiquitin-like containing PHD Ring Finger 1
SET and RING-associated domain
Su(var)3-9, enhancer-of-zeste and Trithorax
Really interesting new gene
DNA methyltransferase 1
Histone deacetylase 1
Protein Data Bank
Patnaik D, Esteve PO, Pradhan S. Targeting the SET and RING-associated (SRA) domain of ubiquitin-like, PHD and ring finger-containing 1 (UHRF1) for anti-cancer drug development. Oncotarget. 2018;9(40):26243–58.
Bashtrykov P, Jankevicius G, Jurkowska RZ, Ragozin S, Jeltsch A. The UHRF1 protein stimulates the activity and specificity of the maintenance DNA methyltransferase DNMT1 by an allosteric mechanism. J Biol Chem. 2014;289(7):4106–15.
Hashimoto H, Horton JR, Zhang X, Bostick M, Jacobsen SE, Cheng X. The SRA domain of UHRF1 flips 5-methylcytosine out of the DNA helix. Nature. 2008;455(7214):826–9.
Frauer C, Hoffmann T, Bultmann S, Casa V, Cardoso MC, Antes I, Leonhardt H. Recognition of 5-hydroxymethylcytosine by the Uhrf1 SRA domain. PLoS ONE. 2011;6(6):e21306.
Arita K, Ariyoshi M, Tochio H, Nakamura Y, Shirakawa M. Recognition of hemi-methylated DNA by the SRA protein UHRF1 by a base-flipping mechanism. Nature. 2008;455(7214):818–21.
Avvakumov GV, Walker JR, Xue S, Li Y, Duan S, Bronner C, Arrowsmith CH, Dhe-Paganon S. Structural basis for recognition of hemi-methylated DNA by the SRA domain of human UHRF1. Nature. 2008;455(7214):822–5.
Qian C, Li S, Jakoncic J, Zeng L, Walsh MJ, Zhou MM. Structure and hemimethylated CpG binding of the SRA domain from human UHRF1. J Biol Chem. 2008;283(50):34490–4.
Myrianthopoulos V, Cartron PF, Liutkeviciute Z, Klimasauskas S, Matulis D, Bronner C, Martinet N, Mikros E. Tandem virtual screening targeting the SRA domain of UHRF1 identifies a novel chemical tool modulating DNA methylation. Eur J Med Chem. 2016;114:390–6.
Berkyurek AC, Suetake I, Arita K, Takeshita K, Nakagawa A, Shirakawa M, Tajima S. The DNA methyltransferase Dnmt1 directly interacts with the SET and RING finger-associated (SRA) domain of the multifunctional protein Uhrf1 to facilitate accession of the catalytic center to hemi-methylated DNA. J Biol Chem. 2014;289(1):379–86.
Unoki M. Current and potential anti-cancer drugs targeting members of the UHRF1 complex including epigenetic modifiers. Recent Pat Anti-Cancer Drug Discov. 2011;6(1):116–30.
Patnaik D. Structure-based screening of chemical libraries to identify small molecules that are likely to bind with the SET and RING-associated (SRA) domain of Ubiquitin-like, PHD and Ring Finger–containing 1 (UHRF1). figshare; 2020. https://doi.org/10.6084/m9.figshare.12086727.
DP wishes to thank Dr. Sriharsa Pradhan, New England Biolabs Inc., for scientific advice. DP also wishes to thank Salma Rafi, Aliz Axmann, and Schrodinger Inc. for the free trial of the Schrodinger’s Small Molecule Drug Discovery Suite; and letter of support for the grant application 1 R43 CA196070-01.
No external funding was received.
Ethics approval and consent to participate
Consent for publication
DP provides consent to publish the data.
The author declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Patnaik, D. Structure-based screening of chemical libraries to identify small molecules that are likely to bind with the SET and RING-associated (SRA) domain of Ubiquitin-like, PHD and Ring Finger-containing 1 (UHRF1). BMC Res Notes 13, 254 (2020). https://doi.org/10.1186/s13104-020-05103-4