Identification of a novel selective small-molecule inhibitor of protein arginine methyltransferase 5 (PRMT5) by virtual screening,resynthesis and biological evaluations

As one of the most promising anticancer target in protein arginine methyltransferase (PRMT) family, PRMT5 has been drawing more and more attentions, and many efforts have been devoted to develop its inhibitors. In this study, three PRMT5 inhibitors (9, 16, and 23) with novel scaffolds were identified by performing pharmacophore- and docking-based virtual screening combined with in vitro radiometric-based scintillation proximity assay (SPA). Substructure search based on the scaffold of the most active 9 afforded 26 additional analogues, and SPA results indicated that two analogues (9–1 and 9–2) showed increased PRMT5 inhibitory activity compared with the parental compound. Resynthesis of 9, 9–1, and 9–2 confirmed their PRMT5 enzymatic inhibition activity. In addition, compound 9–1 displayed selectivity against PRMT5 over other key homological members (PRMT1 and CARM1 (PRMT4)). While the structure–activity relationship (SAR) of this series of compounds was discussed to provide clues for further structure optimization, the probable binding modes of active compounds were also probed by molecular docking and molecular dynamics simulations. Finally, the antiproliferative effect of 9–1 on MV4-11 leukemia cell line was confirmed and its impact on regulating the target gene of PRMT5 was also validated. The hit compounds identified in this work have provided more novel scaffolds for future hit-to-lead optimization of small-molecule PRMT5 inhibitors.     

Understanding protein arginine methyltransferase 1 (PRMT1) product specificity from molecular dynamics

Protein arginine methyltransferases (PRMTs) catalyze the post-translational methylation of specific arginyl groups within targeted proteins to regulate fundamental biological responses in eukaryotic cells. The major Type I PRMT enzyme, PRMT1, strictly generates monomethyl arginine (MMA) and asymmetric dimethylarginine (ADMA), but not symmetric dimethylarginine (SDMA). Multiple diseases can arise from the dysregulation of PRMT1, including heart disease and cancer, which underscores the need to elucidate the origin of product specificity. Molecular dynamics (MD) simulations were carried out for WT PRMT1 and its M48F, H293A, H293S, and H293S-M48F mutants bound with S-adenosylmethionine (AdoMet) and the arginine substrate in an unmethylated or methylated form. Experimental site-directed mutagenesis and analysis of the resultant products were also performed. Two specific PRMT1 active site residues, Met48 and His293, have been determined to play a key role in dictating product specificity, as: (1) the single mutation of Met48 to Phe enabled PRMT1 to generate MMA, ADMA, and a limited amount of SDMA; (2) the single mutation of His293 to Ser formed the expected MMA and ADMA products only; whereas (3) the double mutant H293S-M48F-PRMT1 produced SMDA as the major product with limited amounts of MMA and ADMA. Calculating the formation of near-attack conformers resembling S_N2 transition states leading to either the ADMA or SDMA products finds that Met48 and His293 may enable WT PRMT1 to yield ADMA exclusively by precluding MMA from binding in an orientation more conducive to SDMA formation, i.e., the methyl group bound at the arginine N_η2 position.     

RioK1, a new interactor of protein arginine methyltransferase 5 (PRMT5), competes with pICln for binding and modulates PRMT5 complex composition and substrate specificity

Protein arginine methylation plays a critical role in differential gene expression through modulating protein-protein and protein-DNA/RNA interactions. Although numerous proteins undergo arginine methylation, only limited information is available on how protein arginine methyltransferases (PRMTs) identify their substrates. The human PRMT5 complex consists of PRMT5, WD45/MEP50 (WD repeat domain 45/methylosome protein 50), and pICln and catalyzes the symmetrical arginine dimethylation of its substrate proteins. pICln recruits the spliceosomal Sm proteins to the PRMT5 complex for methylation, which allows their subsequent loading onto snRNA to form small nuclear ribonucleoproteins. To understand how the PRMT5 complex is regulated, we investigated its biochemical composition and identified RioK1 as a novel, stoichiometric component of the PRMT5 complex. We show that RioK1 and pICln bind to PRMT5 in a mutually exclusive fashion. This results in a PRMT5-WD45/MEP50 core structure that either associates with pICln or RioK1 in distinct complexes. Furthermore, we show that RioK1 functions in analogy to pICln as an adapter protein by recruiting the RNA-binding protein nucleolin to the PRMT5 complex for its symmetrical methylation. The exclusive interaction of PRMT5 with either pICln or RioK1 thus provides the first mechanistic insight into how a methyltransferase can distinguish between its substrate proteins.     

结合分子相似性、药效团和分子对接筛选新的HIV-1蛋白酶抑制剂

结合分子相似性、药效团和分子对接建立兼顾计算效率和预测准确度的HIV-1蛋白酶抑制剂筛选方法。首先通过对现有HIV-1蛋白酶抑制剂分子进行相似性分析,选取代表性的HIV-1蛋白酶抑制剂作为模板分子,构建和优化药效团模型,并从1万个化合物中优先筛选出500个化合物。而后采用分子对接方法进一步考察化合物与HIV-1蛋白酶结合情况,得到4个新的活性候选化合物,并进行其结合自由能计算和抗突变性分析。结果表明新候选化合物ST025723和HIV-1蛋白酶表现出较好的相互作用和抗突变性,具有深入研究的价值,同时也证明分子相似性、药效团和分子对接相结合能够快速有效地发现新颖活性候选化合物。     

Discovery of novel HCV polymerase inhibitors using pharmacophore-based virtual screening

We report the use of pharmacophore-based virtual screening as an efficient tool for the discovery of novel HCV polymerase inhibitors. A three-dimensional pharmacophore model for the HCV-796 binding site, NNI site IV inhibitor, to the enzyme was built by means of the structure-based focusing module in Cerius2 program. Using these models as a query for virtual screening, we produced a successful example of using pharmacophore-based virtual screening to identify novel compounds with HCV replicon assay through inhibition of HCV polymerization. Among the hit compounds, compounds 1 and 2 showed 56% and 48% inhibition of NS5B polymerization activity at 20 μM, respectively. In addition, compound 1 also exhibited replicon activity with EC₅₀ value of 2.16 μM. Following up the initial hit, we obtained derivatives of compound 1 and evaluated polymerization inhibition activity and HCV replicon assay. These results provide information necessary for the development of more potent NS5B inhibitors.     

Discovery of new potent inhibitors for carbonic anhydrase IX by structure-based virtual screening

Through structure-based virtual screening, some dozen of benzene sulfonamides with novel scaffolds are identified as potent inhibitors against carbonic anhydrase (CA) IX with IC₅₀ values ranging from 2.86 to 588.34 nM. Among them, compounds 1 and 9 show high selectivity against tumor-target CA IX over CA II (the selectivity ratios are 21.3 and 136.6, respectively). The possible binding poses of hit compounds are also explored and the selectivity is elucidated by molecular docking simulations. The hit compounds discovered in this work would provide novel scaffolds for further hit-to-lead optimization.     

Kinetic mechanism of protein arginine methyltransferase 6 (PRMT6)

The protein arginine methyltransferases (PRMTs) are a family of enzymes that catalyze the mono- and dimethylation of arginine residues in a variety of proteins. Although these enzymes play important roles in a variety of cellular processes, aberrant PRMT activity is associated with several disease states, including heart disease and cancer. In an effort to guide the development of inhibitors targeting individual PRMTs, we initiated studies to characterize the molecular mechanisms of PRMT catalysis. Herein, we report studies on the kinetic mechanism of PRMT6. Initial velocity, product inhibition, and dead-end analog inhibition studies with the AcH4-21 and R1 peptides, as well as their monomethylated versions, indicate, in contrast to a previous report, that PRMT6 utilizes a rapid equilibrium random mechanism with dead-end EAP and EBQ complexes.     

Pyrazole inhibitors of coactivator associated arginine methyltransferase 1 (CARM1)

This study reports the identification and Hits to Leads optimization of inhibitors of coactivator associated arginine methyltransferase (CARM1). Compound 7b is a potent, selective inhibitor of CARM1.     

Isoxazol-5(2H)-one: a new scaffold for potent human neutrophil elastase (HNE) inhibitors

Human neutrophil elastase (HNE) is an important target for the development of novel and selective inhibitors to treat inflammatory diseases, especially pulmonary pathologies. Here, we report the synthesis, structure–activity relationship analysis, and biological evaluation of a new series of HNE inhibitors with an isoxazol-5(2H)-one scaffold. The most potent compound (2o) had a good balance between HNE inhibitory activity (IC₅₀ value =20?nM) and chemical stability in aqueous buffer (t_1/2=8.9?h). Analysis of reaction kinetics revealed that the most potent isoxazolone derivatives were reversible competitive inhibitors of HNE. Furthermore, since compounds 2o and 2s contain two carbonyl groups (2-N-CO and 5-CO) as possible points of attack for Ser195, the amino acid of the active site responsible for the nucleophilic attack, docking studies allowed us to clarify the different roles played by these groups.     

Regulation of protein arginine methyltransferase 8 (PRMT8) activity by its N-terminal domain

Human protein arginine methyltransferase PRMT8 has been recently described as a type I enzyme in brain that is localized to the plasma membrane by N-terminal myristoylation. The amino acid sequence of human PRMT8 is almost 80% identical to human PRMT1, the major protein arginine methyltransferase activity in mammalian cells. However, the activity of a recombinant PRMT8 GST fusion protein toward methyl-accepting substrates is much lower than that of a GST fusion of PRMT1. We show here that both His-tagged and GST fusion species lacking the initial 60 amino acid residues of PRMT8 have enhanced enzymatic activity, suggesting that the N-terminal domain may regulate PRMT8 activity. This conclusion is supported by limited proteolysis experiments showing an increase in the activity of the digested full-length protein, consistent with the loss of the N-terminal domain. In contrast, the activity of the N-terminal truncated protein was slightly diminished by limited proteolysis. Significantly, we detect automethylation at two sites in the N-terminal domain, as well as binding sites for SH3 domain-containing proteins. We suggest that the N-terminal domain may function as an autoregulator that may be displaced by interaction with one or more physiological inducers.     

Discovery of potential sclerostin inhibitors from plants with loop2 region of sclerostin inhibition by interacting with residues outside Pro-Asn-Ala-Ile-Gly motif

Abstract

Sclerostin, an antagonist of the Wnt/β-catenin signaling pathway, was discovered as a potential therapeutic target for stimulating bone formation in osteoporosis. In this study, molecular docking was employed to predict the binding of 29 herbal compounds, which were reported as bone formation stimulators, to the loop2 region of sclerostin. Then, the 50 ns molecular dynamics (MD) simulation of the complexes between sclerostin and the top 10 hits obtained from molecular docking were carried out. Root mean square deviations (RMSDs) analysis of MD trajectories pointed out that all ligands-complexes remain stable throughout the duration of MD simulations. In addition, the molecular mechanics/generalized born surface area (MM/GBSA) binding free energy and energy decomposition analyses were determined. The results here suggested that baicalin is the most promising inhibitor of sclerostin. Interestingly, baicalin binds to sclerostin via the hydrophobic interaction with the amino acid residues on loop2 region but outside the Pro-Asn-Ala-Ile-Gly (PNAIG) motif, particularly the Arg-Gly-Lys-Trp-Trp-Arg (RGKWWR) motif. This finding could be a novel strategy for developing new sclerostin inhibitors in the future.

Communicated by Ramaswamy H. Sarma     

Expansion of the scaffold diversity for the development of highly selective butyrylcholinesterase (BChE) inhibitors: Discovery of new hits through the pharmacophore model generation,virtual screening and molecular dynamics simulation

Butyrylcholinesterase (BChE) is recently considered as a new target for the treatment of Alzheimer’s disease (AD). There is an increasing interest in the development of BChE inhibitors. In the present study, a set of pharmacophore models for BChE was developed and validated. Based on the models, virtual screening was performed on five compound collections, from which seventeen potential hits were retained for biological investigation. In total, eight of these seventeen potential hits showed selective BChE inhibitory activity. Moreover, four compounds displayed IC₅₀ values in sub-micromolar range on eqBChE and three displayed IC₅₀ values < 2 μM on huBChE. The diverse scaffolds of the active compounds provided good starting point further development of selective BChE inhibitors. As far as we concerned, here we disclose the first selective pharmacophore model targeting BChE. The high rate of the model in the identification of active hits indicates it is a valuable tool for the development of selective BChE inhibitors, which may benefit the treatment of AD.     

Discovery of novel 2-(aminoheteroaryl)-thiazole-5-carboxamides as potent and orally active Src-family kinase p56(Lck) inhibitors

A series of substituted 2-(aminoheteroaryl)-thiazole-5-carboxamide analogs have been synthesized as novel, potent inhibitors of the Src-family kinase p56^Lck. Among them, compound 2 displayed superior in vitro potency and excellent in vivo efficacy.     

Discovery of ZAP70 inhibitors by high-throughput docking into a conformation of its kinase domain generated by molecular dynamics

Very few selective inhibitors of the zeta-chain associated protein kinase 70 kDa (ZAP70) have been reported despite its importance in autoimmune diseases. Here, to induce a fit of the so-called gatekeeper residue (Met414) and hydrophobic pocket next to it, a potent Janus kinase 2 (JAK2) inhibitor was first docked into the ATP binding site of ZAP70 by structural alignment of the kinase domains. The resulting model of the complex between ZAP70 and the JAK2 inhibitor was then relaxed by an explicit solvent molecular dynamics simulation with restraints on the backbone atoms. High-throughput docking into the induced-fit conformation of ZAP70 generated by molecular dynamics has revealed 10 low-micromolar inhibitors which correspond to six distinct chemotypes. One of these ZAP70 inhibitors has an IC₅₀ of 110 nM for JAK2.     

Structural features of the reprolysin atrolysin C and tissue inhibitors of metalloproteinases (TIMPs) interaction

Atrolysin C is a P-I snake venom metalloproteinase (SVMP) from Crotalus atrox venom, which efficiently degrades capillary basement membranes, extracellular matrix, and cell surface proteins to produce hemorrhage. The tissue inhibitors of metalloproteinases (TIMPs) are effective inhibitors of matrix metalloproteinases which share some structural similarity with the SVMPs. In this work, we evaluated the inhibitory profile of TIMP-1, TIMP-2, and the N-terminal domain of TIMP-3 (N-TIMP-3) on the proteolytic activity of atrolysin C and analyzed the structural requirements and molecular basis of inhibitor-enzyme interaction using molecular modeling. While TIMP-1 and TIMP-2 had no inhibitory activity upon atrolysin C, the N-terminal domain of TIMP-3 (N-TIMP-3) was a potent inhibitor with a K(i) value of approximately 150nM. The predicted docking structures of atrolysin C and TIMPs were submitted to molecular dynamics simulations and the complex atrolysin C/N-TIMP-3 was the only one that maintained the inhibitory conformation. This study is the first to shed light on the structural determinants required for the interaction between a SVMP and a TIMP, and suggests a structural basis for TIMP-3 inhibitory action and related proteins such as the ADAMs.     

Discovery of novel potent HCV NS5B polymerase non-nucleoside inhibitors bearing a fused benzofuran scaffold

This letter describes the discovery of a fused benzofuran scaffold viable for preparing a series of novel potent HCV NS5B polymerase non-nucleoside inhibitors. Designed on the basis of the functionalized benzofuran derivative nesbuvir (HCV-796), these compounds presumably bind similarly to the allosteric binding site in the “palm” domain of HCV NS5B protein. SAR of each potential hydrogen-bonding interaction site of this novel scaffold is discussed along with some preliminary genotypic profile and PK data of several advanced compounds.