Identification of novel EED-EZH2 PPI inhibitors using an in silico fragment mapping method.

Enhancer of zeste homolog 2 (EZH2) is a histone lysine methyltransferase that is overexpressed in many cancers. Numerous EZH2 inhibitors have been developed as anticancer agents, but recent studies have also focused on protein-protein interaction (PPI) between embryonic ectoderm development (EED) and EZH2 as a novel drug discovery target. Because EED indirectly enhances EZH2 enzymatic activity, EED-EZH2 PPI inhibitors suppress the methyltransferase activity and inhibit cancer growth. By contrast to the numerous promising EZH2 inhibitors, there are a paucity of EED-EZH2 PPI inhibitors reported in the literature. Here, we aimed to discover novel EED-EZH2 PPI inhibitors by first identifying possible binders of EED using an in-house knowledge-based in silico fragment mapping method. Next, 3D pharmacophore models were constructed from the arrangement pattern of the potential binders mapped onto the EED surface. In all, 16 compounds were selected by 3D pharmacophore-based virtual screening followed by docking-based virtual screening. In vitro evaluation revealed that five of these compounds exhibited inhibitory activities. This study has provided structural insights into the discovery and the molecular design of novel EED-EZH2 PPI inhibitors using an in silico fragment mapping method.

[Identifying Inhibitors of USP7-HDM2 Protein-Protein Interaction (PPI) by the in Silico Fragment-mapping Method].

Proteolysis mediated by the ubiquitin-proteome system plays an important role in cancer. Recently, a deubiquitinating enzyme, ubiquitin-specific protease 7 (USP7) has attracted attention as a key regulator of the p53-human double minute 2 (HDM2) pathway in cancer cells. Although some USP7 enzyme inhibitors have been identified, issues related to activity and selectivity prevent their therapeutic application. In this study, we aimed to search for novel USP7-HDM2 protein-protein interaction (PPI) inhibitors that do not affect the USP7 enzyme activity. Using the fragment-mapping program Fsubsite and the canonical subsite-fragment database (CSFDB) developed in our laboratory, we mapped a variety of fragments onto USP7 protein and constructed 3D-pharmacophore models based on the arrangement patterns of the mapped fragments. Finally, we performed 3D pharmacophore-based virtual screening of a commercial compound database and successfully selected promising USP7-HDM2 PPI inhibitor candidates.

Validation of molecular force field parameters for peptides including isomerized amino acids.

Recently, stereoinversions and isomerizations of amino acid residues in the proteins of living beings have been observed. Because isomerized amino acids cause structural changes and denaturation of proteins, isomerizations of amino acid residues are suspected to cause age-related diseases. In this study, AMBER molecular force field parameters were tested by using computationally generated nonapeptides and tripeptides including stereoinverted and/or isomerized amino acid residues. Energy calculations by using density functional theory were also performed for comparison. Although the force field parameters were developed by parameter fitting for l-α-amino acids, the accuracy of the computational results for d-amino acids and ß-amino acids was comparable to those for l-α-amino acids. The conformational energies for tripeptides calculated by using density functional theory were reproduced more accurately than those for nonapeptides calculated by using the molecular mechanical force field. The evaluations were performed for the ff99SB, ff03, ff12SB, and the latest ff14SB force field parameters.

Multi-step virtual screening to develop selective DYRK1A inhibitors.

Developing selective inhibitors for a particular kinase remains a major challenge in kinase-targeted drug discovery. Here we performed a multi-step virtual screening for dual-specificity tyrosine-phosphorylation-regulated kinase 1A (DYRK1A) inhibitors by focusing on the selectivity for DYRK1A over cyclin-dependent kinase 5 (CDK5). To examine the key factors contributing to the selectivity, we constructed logistic regression models to discriminate between actives and inactives for DYRK1A and CDK5, respectively, using residue-based binding free energies. The residue-based parameters were calculated by molecular mechanics-generalized Born surface area (MM-GBSA) decomposition methods for kinase-ligand complexes modeled by computer ligand docking. Based on the findings from the logistic regression models, we built a three-dimensional (3D) pharmacophore model and chose filter criteria for the multi-step virtual screening. The virtual hit compounds obtained from the screening were assessed for their inhibitory activities against DYRK1A and CDK5 by in vitro assay. Our screening identified two novel selective DYRK1A inhibitors with IC50 values of several µM for DYRK1A and >100µM for CDK5, which can be further optimized to develop more potent selective DYRK1A inhibitors.

Prediction of three-dimensional structures and structural flexibilities of wild-type and mutant cytochrome P450 1A2 using molecular dynamics simulations.

In this study, we investigated the effect of genetic polymorphism on the three-dimensional (3D) conformation of cytochrome P450 1A2 (CYP1A2) using molecular dynamics (MD) simulations. CYP1A2, a major drug-metabolizing enzyme among cytochrome P450 enzymes (CYPs), is known to have many variant alleles. The genetic polymorphism of CYP1A2 may cause individual differences in the pharmacokinetics of medicines. By performing 100ns or longer MD simulations, we investigated the influence of amino acid mutation on the 3D structures and the dynamic properties of proteins. The results show that the static structures were changed by the mutations of amino acid residues, not only near the mutated residues but also in distant portions of the proteins. Moreover, the mutation of only one amino acid was shown to change the structural flexibility of proteins, which may influence the substrate recognition and enzymatic activity. Our results clearly suggest that it is necessary to investigate the dynamic property as well as the static 3D structure for understanding the change of the enzymatic activity of mutant CYP1A2.

AMF-26, a novel inhibitor of the Golgi system, targeting ADP-ribosylation factor 1 (Arf1) with potential for cancer therapy.

ADP-ribosylation factor 1 (Arf1) plays a major role in mediating vesicular transport. Brefeldin A (BFA), a known inhibitor of the Arf1-guanine nucleotide exchange factor (GEF) interaction, is highly cytotoxic. Therefore, interaction of Arf1 with ArfGEF is an attractive target for cancer treatment. However, BFA and its derivatives have not progressed beyond the pre-clinical stage of drug development because of their poor bioavailability. Here, we aimed to identify novel inhibitors of the Arf1-ArfGEF interaction that display potent antitumor activity in vivo but with a chemical structure distinct from that of BFA. We exploited a panel of 39 cell lines (termed JFCR39) coupled with a drug sensitivity data base and COMPARE algorithm, resulting in the identification of a possible novel Arf1-ArfGEF inhibitor AMF-26, which differed structurally from BFA. By using a pulldown assay with GGA3-conjugated beads, we demonstrated that AMF-26 inhibited Arf1 activation. Subsequently, AMF-26 induced Golgi disruption, apoptosis, and cell growth inhibition. Computer modeling/molecular dynamics (MD) simulation suggested that AMF-26 bound to the contact surface of the Arf1-Sec7 domain where BFA bound. AMF-26 affected membrane traffic, including the cis-Golgi and trans-Golgi networks, and the endosomal systems. Furthermore, using AMF-26 and its derivatives, we demonstrated that there was a significant correlation between cell growth inhibition and Golgi disruption. In addition, orally administrated AMF-26 (83 mg/kg of body weight; 5 days) induced complete regression of human breast cancer BSY-1 xenografts in vivo, suggesting that AMF-26 is a novel anticancer drug candidate that inhibits the Golgi system, targeting Arf1 activation.

Brazilian propolis-derived components inhibit TNF-α-mediated downregulation of adiponectin expression via different mechanisms in 3T3-L1 adipocytes.

BACKGROUND: Previous reports suggest that Brazilian propolis has multiple biological functions and may help to restore adiponectin expression and insulin sensitivity. However, little is known about the molecular mechanisms by which these compounds inhibit the downregulation of adiponectin. METHODS: The effect of various Brazilian propolis-derived components on inhibition of tumor necrosis factor-α (TNF-α)-mediated downregulation of adiponectin expression in 3T3-L1 adipocytes and molecular mechanism was investigated. RESULTS AND CONCLUSIONS: Pretreatment with either artepillin C (C3) or its derivative (C4) significantly inhibited TNF-α-mediated downregulation of adiponectin expression in 3T3-L1 adipocytes. Interestingly, C3 strongly activated peroxisome proliferator-activated receptor γ (PPARγ) transcriptional activity. Treatment of adipocytes with C3 resulted in the upregulation of adiponectin and fatty acid-binding protein 4 expression, but C4 did not significantly induce PPARγ transactivation. C4 did, however, inhibit the TNF-α-induced c-Jun-NH(2)-terminal kinase (JNK) signaling that is involved in adiponectin expression. Molecular docking studies based on hPPARγ with C3 and JNK1 with C4 clearly supported our experimental results. These data demonstrate that 1) both C3 and C4 significantly inhibit the TNF-α-mediated downregulation of adiponectin in adipocytes, 2) C3 functions as a PPARγ agonist, and its inhibition of the effect of TNF-α is due to this PPARγ transactivation, and 3) C4 is an effective inhibitor of JNK activation, thus inhibiting the TNF-α-mediated downregulation of adiponectin. GENERAL SIGNIFICANCE: Brazilian propolis-derived components (C3 and C4) can significantly inhibit TNF-α-mediated downregulation of adiponectin in adipocytes, although they do so via different mechanisms.

New AMBER force field parameters of heme iron for cytochrome P450s determined by quantum chemical calculations of simplified models.

The heme protein, cytochrome P450, is an oxidoreductase that plays an important role in drug metabolism. To model P450s using molecular mechanics methods and classical molecular dynamics simulations, force field parameters and atomic charges are required. Because these parameters are generally obtained by quantum chemical methods, an appropriate simplified model for the iron-porphyrin system was needed. In this study, two models with a five-coordinated Fe(III) mimicking the sextet spin state of P450s are proposed, which are optimized by semiempirical and ab initio unrestricted Hartree-Fock methods. The results produced using the simpler of the two models were similar to those of the more complex model; therefore, the more simplified model of P450 can be used without a loss of accuracy. Furthermore, several quantum chemical calculations were carried out on the simpler model to investigate which method was most suitable for iron-porphyrin systems. The results calculated by hybrid density functional theory (DFT), with the MIDI basis set for iron, reproduced the three-dimensional structures determined by X-ray diffraction and extended X-ray absorption fine-structure experiments. From these results, atomic charges and force-field parameters for molecular mechanics and molecular dynamics calculations were obtained.

Use of photoaffinity labeling and site-directed mutagenesis for identification of the key residue responsible for extraordinarily high affinity binding of UCN-01 in human alpha1-acid glycoprotein.

7-Hydroxystaurosporine (UCN-01) is a protein kinase inhibitor anticancer drug currently undergoing a phase II clinical trial. The low distribution volumes and systemic clearance of UCN-01 in human patients have been found to be caused in part by its extraordinarily high affinity binding to human alpha1-acid glycoprotein (hAGP). In the present study, we photolabeled hAGP with [3H]UCN-01 without further chemical modification. The photolabeling specificity of [3H]UCN-01 was confirmed by findings in which other hAGP binding ligands inhibited formation of covalent bonds between hAGP and [3H]UCN-01. The amino acid sequence of the photolabeled peptide was concluded to be SDVVYTDXK, corresponding to residues Ser-153 to Lys-161 of hAGP. No PTH derivatives were detected at the 8th cycle, which corresponded to the 160th Trp residue. This strongly implies that Trp-160 was photolabeled by [3H]UCN-01. Three recombinant hAGP mutants (W25A, W122A, and W160A) and wild-type recombinant hAGP were photolabeled by [3H]UCN-01. Only mutant W160A showed a marked decrease in the extent of photoincorporation. These results strongly suggest that Trp-160 plays a prominent role in the high affinity binding of [3H]UCN-01 to hAGP. A docking model of UCN-01 and hAGP around Trp-160 provided further details of the binding site topology.