Combined structure-based pharmacophore and 3D-QSAR studies on phenylalanine series compounds as TPH1 inhibitors

AI-based 3D-QSAR model of FDA-approved repurposed drugs for inhibiting sclerostin

BACKGROUND

Wnt activation promotes bone formation and prevents bone loss. The Wnt pathway antagonist sclerostin and additional anti-sclerostin antibodies were discovered as a result of the development of the monoclonal antibody romosozumab. These monoclonal antibodies greatly increase the risk of cardiac arrest. Three-dimensional quantitative structure-activity relationships (3D-QSAR) predicts biological activities of ligands based on their three-dimensional features by employing powerful chemometric investigations such as artificial neural networks (ANNs) and partial least squares (PLS).

OBJECTIVE

In this study, ligand-receptor interactions were investigated using 3D-QSAR Comparative molecular field analysis (CoMFA). Estimates of steric and electrostatic characteristics in CoMFA are made using Lennard-Jones and Coulomb potentials.

METHODS

To identify the conditions necessary for the activity of these molecules, fifty Food and Drug Administration (FDA)-approved medications were chosen for 3D-QSAR investigations and done by CoMFA. For QSAR analysis, there are numerous tools available. This study employed Open 3D-QSAR for analysis due to its simplicity of use and capacity to produce trustworthy results. Four tools were used for the analysis on this platform: Py-MolEdit, Py-ConfSearch, and Py-CoMFA.

RESULTS

Maps that were generated were used to determine the screen's r2 (Coefficient of Multiple Determinations) value and q2 (correlation coefficient). These numbers must be fewer than 1, suggesting a good, trustworthy model. Cross-validated (q2) 0.532 and conventional (r2) correlation values of 0.969 made the CoMFA model statistically significant. The model showed that hydroxamic acid inhibitors are significantly more sensitive to the steric field than the electrostatic field (70%) (30%). This hypothesis states that steric (43.1%), electrostatic (26.4%), and hydrophobic (20.3%) qualities were important in the design of sclerostin inhibitors.

CONCLUSION

With 3D-QSAR and CoMFA, statistically meaningful models were constructed to predict ligand inhibitory effects. The test set demonstrated the model's robustness. This research may aid in the development of more effective sclerostin inhibitors that are synthesised using FDA-approved medications.

Inhibition of Tryptophan Hydroxylases and Monoamine Oxidase-A by the Proton Pump Inhibitor, Omeprazole-In Vitro and In Vivo Investigations

Serotonin (5-HT) is a hormone and neurotransmitter that modulates neural activity as well as a wide range of other physiological processes including cardiovascular function, bowel motility, and platelet aggregation. 5-HT synthesis is catalyzed by tryptophan hydroxylase (TPH) which exists as two distinct isoforms; TPH1 and TPH2, which are responsible for peripheral and central 5-HT, respectively. Due to the implication of 5-HT in a number of pathologies, including depression, anxiety, autism, sexual dysfunction, irritable bowel syndrome, inflammatory bowel disease, and carcinoid syndrome, there has been a growing interest in finding modulators of these enzymes in recent years. We thus performed high-throughput screening (HTS) using a fluorescence-based thermal shift assay (DSF) to search the Prestwick Chemical Library containing 1,280 compounds, mostly FDA-approved drugs, for TPH1 binders. We here report the identification of omeprazole, a proton pump inhibitor, as an inhibitor of TPH1 and TPH2 with low micromolar potency and high selectivity over the other aromatic amino acid hydroxylases. The S-enantiomer of omeprazole, esomeprazole, has recently also been described as an inhibitor of monoamine oxidase-A (MAO-A), the main enzyme responsible for 5-HT degradation, albeit with lower potency compared to the effect on TPH1 and TPH2. In order to investigate the net effect of simultaneous inhibition of TPH and MAO-A in vivo, we administered high-dose (100 mg/kg) omeprazole to CD-1 mice for 4 days, after which the animals were subjected to the tail suspension test. Finally, central (whole brain) and peripheral (serum) 5-HT content was measured using liquid chromatography-mass spectrometry (LC-MS). Omeprazole treatment significantly increased 5-HT concentrations, both in brain and in serum, and reduced the time spent immobile in the tail suspension test relative to vehicle control. Thus, the MAO-A inhibition afforded by high-dose omeprazole appears to overcome the opposing effect on 5-HT produced by inhibition of TPH1 and TPH2. Further modification of proton pump inhibitor scaffolds may yield more selective modulators of 5-HT metabolism.

Discovery and biological characterization of a novel scaffold for potent inhibitors of peripheral serotonin synthesis

Aim: Tryptophan hydroxylase 1 (TPH1) catalyzes serotonin synthesis in peripheral tissues. Selective TPH1 inhibitors may be useful for treating disorders related to serotonin dysregulation. Results & methodology: Screening using a thermal shift assay for TPH1 binders yielded Compound 1 (2-(4-methylphenyl)-1,2-benzisothiazol-3(2 H)-one), which showed high potency (50% inhibition at 98 ± 30 nM) and selectivity for inhibiting TPH over related aromatic amino acid hydroxylases in enzyme activity assays. Structure–activity relationships studies revealed several analogs of 1 showing comparable potency. Kinetic studies suggested a noncompetitive mode of action of 1, with regards to tryptophan and tetrahydrobiopterin. Computational docking studies and live cell assays were also performed. Conclusion: This TPH1 inhibitor scaffold may be useful for developing new therapeutics for treating elevated peripheral serotonin.

Design, synthesis, and biological evaluation of polo-like kinase 1/eukaryotic elongation factor 2 kinase (PLK1/EEF2K) dual inhibitors for regulating breast cancer cells apoptosis and autophagy

Both PLK1 and EEF2K are serine⁄threonine kinases that play important roles in the proliferation and programmed cell death of various types of cancer. They are highly expressed in breast cancer tissues. Based on the multiple-complexes generated pharmacophore models of PLK1 and homology models of EEF2K, the integrated virtual screening is performed to discover novel PLK1/EEF2K dual inhibitors. The top ten hit compounds are selected and tested in vitro, and five of them display PLK1 and EEF2K inhibition in vitro. Based on the docking modes of the most potent hit compound, a series of derivatives are synthesized, characterized and biological assayed on the PLK1, EEF2K as well as breast cancer cell proliferation models. Compound 18i with satisfied inhibitory potency are shifted to molecular mechanism studies contained molecular dynamics simulations, cell cycles, apoptosis and autophagy assays. Our results suggested that these novel PLK1/EEF2K dual inhibitors can be used as lead compounds for further development breast cancer chemotherapy.

3D-QSAR, molecular docking, and molecular dynamic simulations for prediction of new Hsp90 inhibitors based on isoxazole scaffold

Heat shock protein 90(Hsp90), as a molecular chaperone, play a crucial role in folding and proper function of many proteins. Hsp90 inhibitors containing isoxazole scaffold are currently being used in the treatment of cancer as tumor suppressers. Here in the present studies, new compounds based on isoxazole scaffold were predicted using a combination of molecular modeling techniques including three-dimensional quantitative structure-activity relationship (3D-QSAR), molecular docking and molecular dynamic (MD) simulations. Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) were also done. The steric and electrostatic contour map of CoMFA and CoMSIA were created. Hydrophobic, hydrogen bond donor and acceptor of CoMSIA model also were generated, and new compounds were predicted by CoMFA and CoMSIA contour maps. To investigate the binding modes of the predicted compounds in the active site of Hsp90, a molecular docking simulation was carried out. MD simulations were also conducted to evaluate the obtained results on the best predicted compound and the best reported Hsp90 inhibitors in the 3D-QSAR model. Findings indicate that the predicted ligands were stable in the active site of Hsp90.

Crystal structure-based discovery of a novel synthesized PARP1 inhibitor (OL-1) with apoptosis-inducing mechanisms in triple-negative breast cancer

Poly (ADP-ribose) polymerase-1 (PARP1) is a highly conserved enzyme focused on the self-repair of cellular DNA damage. Until now, numbers of PARP inhibitors have been reported and used for breast cancer therapy in recent years, especially in TNBC. However, developing a new type PARP inhibitor with distinctive skeleton is alternatively promising strategy for TNBC therapy. In this study, based on co-crystallization studies and pharmacophore-docking-based virtual screening, we discovered a series of dihydrodibenzo[b,e]-oxepin compounds as PARP1 inhibitors. Lead optimization result in the identification of compound OL-1 (2-(11-(3-(dimethylamino)propylidene)-6,11- dihydrodibenzo[b,e]oxepin )-2-yl)acetohydrazide), which has a novel chemical scaffold and unique binding interaction with PARP1 protein. OL-1 demonstrated excellent potency (inhibiting PARP1 enzyme activity with IC₅₀ = 0.079 μM), as well as inhibiting PARP-modulated PARylation and cell proliferation in MDA-MB-436 cells (BRAC1 mutation). In addition, OL-1 also inhibited cell migration that closely related to cancer metastasis and displayed remarkable anti-tumor efficacy in MDA-MB-436 xenograft model without apparent toxicities. These findings highlight a new small-molecule PAPR1 inhibitor (OL-1) that has the potential to impact future TNBC therapy.

Combining structure-based pharmacophore modeling, virtual screening, and in silico ADMET analysis to discover novel tetrahydro-quinoline based pyruvate kinase isozyme M2 activators with antitumor activity

Compared with normal differentiated cells, cancer cells upregulate the expression of pyruvate kinase isozyme M2 (PKM2) to support glycolytic intermediates for anabolic processes, including the synthesis of nucleic acids, amino acids, and lipids. In this study, a combination of the structure-based pharmacophore modeling and a hybrid protocol of virtual screening methods comprised of pharmacophore model-based virtual screening, docking-based virtual screening, and in silico ADMET (absorption, distribution, metabolism, excretion and toxicity) analysis were used to retrieve novel PKM2 activators from commercially available chemical databases. Tetrahydroquinoline derivatives were identified as potential scaffolds of PKM2 activators. Thus, the hybrid virtual screening approach was applied to screen the focused tetrahydroquinoline derivatives embedded in the ZINC database. Six hit compounds were selected from the final hits and experimental studies were then performed. Compound 8 displayed a potent inhibitory effect on human lung cancer cells. Following treatment with Compound 8, cell viability, apoptosis, and reactive oxygen species (ROS) production were examined in A549 cells. Finally, we evaluated the effects of Compound 8 on mice xenograft tumor models in vivo. These results may provide important information for further research on novel PKM2 activators as antitumor agents.

Tryptophan hydroxylase 1 (Tph-1)-targeted bone anabolic agents for osteoporosis

Tryptophan hydroxylase 1 (Tph-1), the principal enzyme for peripheral serotonin biosynthesis, provides a novel target to design anabolic agents for osteoporosis. Here, we present a design, synthesis of a novel series of ursolic acid derivatives under the guidance of docking technique, and bioevaluation of the derivatives using RBL2H3 cells and ovariectomized (OVX) rats. Of the compounds, 9a showed a potent inhibitory activity on serotonin biosynthesis. Further investigations revealed that 9a, as an efficient Tph-1 binder identified by SPR (estimated KD: 6.82 μM), suppressed the protein and mRNA expressions of Tph-1 and lowered serotonin contents in serum and gut without influence on brain serotonin. Moreover, oral administration of 9a elevated serum level of N-terminal propeptide of procollagen type 1 (P1NP), a bone formation marker, and improved bone microarchitecture without estrogenic side effects in ovariectomized rats. Collectively, 9a may serve as a new candidate for bone anabolic drug discovery.

Design, synthesis, antifungal activities and 3D-QSAR of new N,N'-diacylhydrazines containing 2,4-dichlorophenoxy moiety

A series of new N,N'-diacylhydrazine derivatives were designed and synthesized. Their structures were verified by 1H-NMR, mass spectra (MS) and elemental analysis. The antifungal activities of these N,N'-diacylhydrazines were evaluated. The bioassay results showed that most of these N,N'-diacylhydrazines showed excellent antifungal activities against Cladosporium cucumerinum, Corynespora cassiicola, Sclerotinia sclerotiorum, Erysiphe cichoracearum, and Colletotrichum orbiculare in vivo. The half maximal effective concentration (EC50) of one of the compounds was also determined, and found to be comparable with a commercial drug. To further investigate the structure-activity relationship, comparative molecular field analysis (CoMFA) was performed on the basis of antifungal activity data. Both the steric and electronic field distributions of CoMFA are in good agreement in this study.

Novel hybrid virtual screening protocol based on molecular docking and structure-based pharmacophore for discovery of methionyl-tRNA synthetase inhibitors as antibacterial agents

Methione tRNA synthetase (MetRS) is an essential enzyme involved in protein biosynthesis in all living organisms and is a potential antibacterial target. In the current study, the structure-based pharmacophore (SBP)-guided method has been suggested to generate a comprehensive pharmacophore of MetRS based on fourteen crystal structures of MetRS-inhibitor complexes. In this investigation, a hybrid protocol of a virtual screening method, comprised of pharmacophore model-based virtual screening (PBVS), rigid and flexible docking-based virtual screenings (DBVS), is used for retrieving new MetRS inhibitors from commercially available chemical databases. This hybrid virtual screening approach was then applied to screen the Specs (202,408 compounds) database, a structurally diverse chemical database. Fifteen hit compounds were selected from the final hits and shifted to experimental studies. These results may provide important information for further research of novel MetRS inhibitors as antibacterial agents.

Molecular dynamics simulation of tryptophan hydroxylase-1: binding modes and free energy analysis to phenylalanine derivative inhibitors

Serotonin is a neurotransmitter that modulates many central and peripheral functions. Tryptophan hydroxylase-1 (TPH1) is a key enzyme of serotonin synthesis. In the current study, the interaction mechanism of phenylalanine derivative TPH1 inhibitors was investigated using molecular dynamics (MD) simulations, free energy calculations, free energy decomposition analysis and computational alanine scanning. The predicted binding free energies of these complexes are consistent with the experimental data. The analysis of the individual energy terms indicates that although the van der Waals and electrostatics interaction contributions are important in distinguishing the binding affinities of these inhibitors, the electrostatic contribution plays a more crucial role in that. Moreover, it is observed that different configurations of the naphthalene substituent could form different binding patterns with protein, yet lead to similar inhibitory potency. The combination of different molecular modeling techniques is an efficient way to interpret the interaction mechanism of inhibitors and our work could provide valuable information for the TPH1 inhibitor design in the future.

Discovery of novel focal adhesion kinase inhibitors using a hybrid protocol of virtual screening approach based on multicomplex-based pharmacophore and molecular docking

Focal adhesion kinase (FAK) is a tyrosine kinase that functions as a key orchestrator of signals leading to invasion and metastasis. In the current study, the multicomplex-based pharmacophore (MCBP)-guided method has been suggested to generate a comprehensive pharmacophore of FAK kinase based on seven crystal structures of FAK-inhibitor complexes. In this investigation, a hybrid protocol of virtual screening methods, comprising of pharmacophore model-based virtual screening (PB-VS) and docking-based virtual screening (DB-VS), is used for retrieving new FAK inhibitors from commercially available chemical databases. This hybrid virtual screening approach was then applied to screen several chemical databases, including the Specs (202,408 compounds) database. Thirty-five compounds were selected from the final hits and should be shifted to experimental studies. These results may provide important information for further research of novel FAK inhibitors.