Current Computer - Aided Drug Design - Volume 16, Issue 5, 2020

Background: The computational studies on 2-phenazinamines with their protein targets have been carried out to design compounds with potential anticancer activity and selectivity over specific BCR-ABL Tyrosine kinase. Methods: This has been achieved through G-QSAR and molecular docking studies. Computational chemistry was done by using VLife MDS 4.3 and Autodock 4.2. 2D and structures of ligands were drawn by using Chemdraw 2D Ultra 8.0 and were converted into 3D. These were optimized by using semi-empirical method called MOPAC. The protein structure was downloaded as PDB file from RCSC protein data bank. PYMOL was used for studying the binding interactions. The G-QSAR models generated were found to possess training (r2=0.8074), cross-validation (q2=0.6521), and external validation (pred_r2=0.5892) which proved their statistical significance. Accordingly, the newly designed series of 2-phenazinamines viz., 3-chloro-4-aryl-1-(phenazin-7-yl) azetidin-2-ones (4a-4e) were subjected to wet lab synthesis. Alternatively, docking studies were also conducted which showed binding interactions of some derivatives with > 30% higher binding energy values than the standard anticancer drug imatinib. The lower energy values obtained for these derivatives indicate energetically favorable interaction with protein binding site as compared to standard imatinib. Results: G-QSAR and molecular docking studies predicted better anticancer activity for the synthesized azitidine derivatives of 2-phenazinamines (4a-4e) as compared to standard drug. Conclusion: It is therefore surmised that the molecular manipulations at appropriate sites of these derivatives suggested by structure activity relationship data will prove to be beneficial in raising anticancer potential.

Background: A series of new six thiazolyl-2-amine-based Schiff base derivatives (4a-4f) were synthesized by a sequential multistep reaction starting with Salicylaldehyde. Methods: All the Schiff base derivatives were screened in-vitro for their antibacterial activity against Mycobacterium tuberculosis (H37RV strain) ATCC No-27294. The synthesized compounds were characterized by FTIR, 1H-NMR, 13C-NMR and Mass spectrometry. Results: Among the compounds tested, 4c and 4f derivatives exhibited potent antitubercular activity against M. tuberculosis at MIC 6.25 μg/mL. Conclusion: We extended our study to explore the inhibition mechanism by conducting molecular docking analysis by using Schrodinger’s molecular modeling software. All the newly synthesized compounds were found to be in-silico AMES test non-toxic and non-carcinogens. The good Qikprop’s Absorption, Distribution, Metabolism and Excretion (ADMET) would definitely help the researchers in order to make more potent Anti-TB agents.

Background: Shikimate pathway is essential for tubercular bacillus but it is absent in mammals. Therefore, Shikimate kinase and other enzymes in the pathway are potential targets for the development of novel anti-tuberculosis drugs. Objective: In the present study, Shikimate kinase is selected as the target for in silico screening of phytochemicals with an aim to discover a novel herbal drug against Mycobacterium tuberculosis (Mtb). Methods: A structure-based drug discovery approach is undertaken for the execution of the objective. Virtual screening of phytochemical database NPACT against the target, Shikimate kinase (PDB ID 3BAF), is carried out followed by toxicity and drug-likeness filtration. Finally, a lead, narirutin was selected for in vitro anti-tubercular study. Results: Narirutin, present in citrus fruits, emerges as the lead. It is considered to be non-toxic with predicted high LD50 value, 12000 mg/kg body weight. The phytochemical is tested for its antitubercular activity in vitro. It has MIC99 62.5 μg/mL against the MtbH37Rv strain. Conclusion: This is the first-ever report to show anti-tuberculosis potency of narirutin.

Background: Benzothiazoles are reported to have bioorganic and pharmaceutical chemistry applications. Introduction: A series of substituted N-bromoamido-2-aminobenzothiazoles was synthesized from substituted anilines via 2-aminobenzothiazoles and it was further evaluated for its antimicrobial activity. Methods: All the newly synthesized compounds were characterized by FT-IR, NMR and mass spectra and purity profiles were studied by HPLC analysis. The antimicrobial testing (MIC determination) was newly performed with agar micro-broth dilution method for these analogs. Results: Among the synthesized compound 3b showed the highest activity with MIC value of 3.12 μg/mL against Bacillus, E. coli, S. aureus and Klebsiella and 6.25 μg/mL against C. albicans. The ADME properties as calculated by using Qikprop were found within acceptable range. Derivatives shows a good-moderate binding affinity towards target Cytochrome P450 14 alpha-sterol demethylase (CYP51) (PDB ID: 1EA1). Conclusion: Our in-silico and in-vitro studies on a series of substituted aminobenzothiazoles may be helpful for further designing of more potent antimicrobials in future.

Introduction: The enzyme Glycogen Synthase Kinase 3-β (GSK-3β) is related to neuronal cell degeneration, representing a promising target to treat Alzheimer’s Disease (AD). Methods: In this work, we performed a molecular modeling study of existing GSK-3β inhibitors by means of evaluation of their IC50 values, derivation of a pharmacophore model, molecular docking simulations, ADME/Tox properties predictions, molecular modifications and prediction of synthetic viability. Results: In this manner, inhibitor 15 (CID 57399952) was elected a template molecule, since it demonstrated to bear relevant structural groups able to interact with GSK-3β, and also presented favorable ADME/Tox predicted properties, except for mutagenicity. Based on this inhibitor chemical structure we proposed six analogues that presented the absence of alerts for mutagenic and carcinogenic activity, both for rats and mouse; likewise they all presented low risk alerts for inhibition of hERG and medium prediction of synthetic viability. Conclusion: It is concluded that the analogues of GSK-3β inhibitors were optimized in relation to the toxicity endpoint of the template molecule, being, therefore, presented as novel and promising drug candidates for AD treatment.

Background: H+/K+ ATPase a protein present in the gastric parietal cells is a better target for the prevention and treatment of gastric ulcer. Plant flavonoids have been reported to elicit anti-ulcer activity by inhibiting the proton pump as well as by antioxidant defense mechanism. Methods: Chloroform fraction of hydro-alcoholic extract of passion fruit was screened for proton pump inhibitory assay using goat parietal cell. In-silico computational docking studies were carried out using Glide program in order to validate the inhibitory action of selected constituents. Results: The flavonoid rich fruit possess a promising radical scavenging activity against DPPH. 10.41μg/mL is sufficient to inhibit 50% of ATPase enzyme activity. A synergistic activity was also achieved by the fruit with sub-effective doses of lansoprazole. Fenton's oxidation induced by H2O2 was also blunted by the fruit extract. Conclusion: The in-vitro and in-silico findings indicated that, passion fruit can be a good dietary supplement for the prevention and management of ulcer.

Background: Pyrazol-5-amine derivatives are an important class of heterocyclic compounds. However, there are less 4-alkyl substituted pyrazoles reported. Objective: Here reported are the design, synthesis and biological evaluation of 3-aryl-4- alkylpyrazol-5-amines derivatives. Methods: A serials of 3-aryl-4-alkylpyrazol-5-amines were designed and the biological action targets were screened by target fishing function of Discovery Studio software. The synthesis route involved 3-oxo-3-arylpropanenitrile formation, alkylation, pyrazole formation, and amides formation. The antitumor activities of these compounds were carried out by thiazolyl blue tetrazolium bromide (MTT) method using U-2 OS (osteosarcoma) and A549 (lung cancer) tumor cells. Results: Eight 3-aryl-4-alkylpyrazol-5-amines were synthesized, and their structures were verified by 1H NMR, 13C NMR, and HRMS. Thirteen pharmacophores were mapped out by target fishing. Compound 5h showed anti-proliferation activities against U-2 OS and A549 tumor cell with IC50 value of 0.9 μM and 1.2 μM, respectively. Conclusion: Compound 5h might represent a promising scaffold for the further development of novel antitumor drugs.

Aim and Objective: Cycloalkanes have been largely used in the field of medicine, components of food, pharmaceutical drugs, and they are mainly used to produce fuel. In present study the relationship between molecular descriptors and thermodynamic properties such as the standard enthalpies of formation (ΔH°f), the standard enthalpies of fusion (ΔH°fus), and the standard Gibbs free energy of formation (ΔG°f)of the cycloalkanes is represented. Materials and Methods: The Genetic Algorithm (GA) and multiple linear regressions (MLR) were successfully used to predict the thermodynamic properties of cycloalkanes. A large number of molecular descriptors were obtained with the Dragon program. The Genetic algorithm and backward method were used to reduce and select suitable descriptors. Results: QSPR models were used to delineate the important descriptors responsible for the properties of the studied cycloalkanes. The multicollinearity and autocorrelation properties of the descriptors contributed in the models were tested by calculating the Variance Inflation Factor (VIF), Pearson Correlation Coefficient (PCC) and the Durbin–Watson (DW) statistics. The predictive powers of the MLR models were discussed using Leave-One-Out Cross-Validation (LOOCV) and test set validation methods. The statistical parameters of the training, and test sets for GA–MLR models were calculated. Conclusion: The results of the present study indicate that the predictive ability of the models was satisfactory and molecular descriptors such as: the Functional group counts, Topological indices, GETAWAY descriptors, Constitutional indices, and molecular properties provide a promising route for developing highly correlated QSPR models for prediction the studied properties.

Background: Human African trypanosomiasis is a fatal disease prevalent in approximately 36 sub-Saharan countries. Emerging reports of drug resistance in Trypanosoma brucei are a serious cause of concern as only limited drugs are available for the treatment of the disease. Pteridine reductase is an enzyme of Trypanosoma brucei. Methods: It plays a critical role in the pterin metabolic pathway that is absolutely essential for its survival in the human host. The success of finding a potent inhibitor in structure-based drug design lies within the ability of computational tools to efficiently and accurately dock a ligand into the binding cavity of the target protein. Here we report the computational characterization of Trypanosoma brucei pteridine reductase (Tb-PR) active-site using twenty-four high-resolution co-crystal structures with various drugs. Structurally, the Tb-PR active site can be grouped in two clusters; one with high Root Mean Square Deviation (RMSD) of atomic positions and another with low RMSD of atomic positions. These clusters provide fresh insight for rational drug design against Tb-PR. Henceforth, the effect of several factors on docking accuracy, including ligand and protein flexibility were analyzed using Fred. Results: The online server was used to analyze the side chain flexibility and four proteins were selected on the basis of results. The proteins were subjected to small-scale virtual screening using 85 compounds, and statistics were calculated using Bedroc and roc curves. The enrichment factor was also calculated for the proteins and scoring functions. The best scoring function was used to understand the ligand protein interactions with top common compounds of four proteins. In addition, we made a 3D structural comparison between the active site of Tb-PR and Leishmania major pteridine reductase (Lm- PR). We described key structural differences between Tb-PR and Lm-PR that can be exploited for rational drug design against these two human parasites. Conclusion: The results indicated that relying just on re-docking and cross-docking experiments for virtual screening of libraries isn’t enough and results might be misleading. Hence it has been suggested that small scale virtual screening should be performed prior to large scale screening.

Background: Small interfering RNAs (siRNAs) are known as commonly used targeting mRNAs tools for suppressing gene expression. Since Signal Transducer and Activator of Transcription 4 (STAT4) is considered as a significant transcription factor for generation and differentiation of Th1 cells during vascular dysfunction and atherosclerosis, suppressing STAT4 could represent novel immunomodulatory therapies against atherosclerosis. Objective: Therefore, the current study was conducted to design efficient siRNAs specific for STAT4 and to evaluate different criteria affecting their functionality. Methods: In the present study, all related sequences of STAT4 gene were retrieved from Gen Bank database. Multiple sequence alignment was carried out to recognize Open Reading Frame (ORF) and conserved region. Then, siDirect 2.0 server was applied for the development of candidate siRNA molecules and confirmation of predicted molecules was performed using Dharma siRNA technology and GeneScript siRNA targetfinder. In addition, BLAST tool was used against whole Genebank databases to identify potential off-target genes. DNA/RNA GC content calculator and mfold server were used to calculate GC content and secondary structure prediction of designed siRNA, respectively. Finally, IntaRNA program was used to study the thermodynamics of interaction between predicted siRNA and target gene. Results: Based on the obtained results, three efficient siRNA molecules were designed and validated for STAT4 gene silencing using computational methods, which may result in suppressing STAT4 gene expression. Conclusion: According to our results, this study shows that siRNA targeting STAT4 can be considered as a therapeutic agent in many Th1-mediated pathologic conditions specially atherosclerosis.

Background: Peptidyl-prolyl cis-trans isomerase NIMA-interacting 1 (Pin1) is an enzyme that isomerizes phosphorylated serine or threonine motifs adjacent to proline residues. Pin1 has important roles in several cellular signaling pathways, consequently impacting the development of multiple types of cancers. Methods: Based on the previously reported inhibitory activity of pentacyclic triterpenoids isolated from the gum resin of Boswellia genus against Pin1, we designed a computational experiment using molecular docking, pharmacophore filtering, and structural clustering allied to molecular dynamics (MD) simulations and binding free energy calculations to explore the inhibitory activity of new triterpenoids against Pin1 structure. Results: Here, we report different computational evidence that triterpenoids from neem (Azadirachta indica A. Juss), such as 6-deacetylnimbinene, 6-Oacetylnimbandiol, and nimbolide, replicate the binding mode of the Pin1 substrate peptide, interacting with high affinity with the binding site and thus destabilizing the Pin1 structure. Conclusions: Our results are supported by experimental data, and provide interesting structural insights into their molecular mechanism of action, indicating that their structural scaffolds could be used as a start point to develop new inhibitors against Pin1.

Background: The treatment of multiple-drug-resistant tuberculosis (MDR-TB) with currently available marketed drugs remains a global health concern. The cases of resistant tuberculosis patients are increasing day by day. Objective: The objective of this study is to highlight the need of developing shorter, simpler and tolerable drug regimens. Methods: In the present study, we synthesized various halo-substituted 2-aryloxyacetohydrazones via a series of reactions from halo-substituted phenols. All the compounds were characterized by using various spectroscopic methods, such as NMR, FT-IR, UV spectroscopy, etc. Results: All the synthesized hydrazones showed theoretically good interactions with enzyme enoyl reductase (pdb id: 4tzk). All the synthesized compounds (5a-5o) showed moderate to good activity (3.125-100 μg/mL) against Mycobacteria tuberculosis, H37RV strain. Conclusion: Our results would pave a new way for the development of more effective Anti-TB agents in the future.

Background: Cholestatic liver disease, a serious chronic condition that develops progressive hepatic degeneration through free radicals. Objective: The present study was designed to extract and identify two flavonoids in Phlomoides hyoscyamoides plant, native to Iran and evaluate the role of quercetin identified on the liver injury among bile ductligated rats. Methods: This study was conducted on 25 male Wistar rats within three groups of sham control, mere bile duct-ligated, and bile duct-ligated with quercetin. The bile duct-ligated animals received quercetin at a dose of 50 mg/kg/day for 10 days, followed by biochemical tests, oxidative stress markers, activity of antioxidant enzymes and hematoxylin and eosin staining. Molecular docking was used to explore the interactive behavior of quercetin with glutathione peroxidase. Results: According to analyses of the obtained extract, two main active ingredients of P. hyoscyamoides were rutin and quercetin. Bile duct-ligated group showed a significant liver necrosis, a clear increase in plasma and tissue oxidative stress parameters, and a decrease in glutathione peroxidase activity as compared to sham control group. Quercetin injection in bile duct-ligated rats resulted in significant decrease in hydroxyproline, protein carbonyl and histopathologic indexes and significant increase in glutathione peroxidase activity (P-value≤0.05). Based on the molecular docking, the quercetin was able to regulate the glutathione peroxidase activity. Conclusion: The quercetin acts as an enzyme inducer by renewing the glutathione peroxidase activity and inhibiting the oxidation of proteins and hence decreases the oxidative stress. These results could be a sign of confirming the positive role of quercetin in attenuating the liver damage and degeneration.

Background: Multicentric Castleman Disease (MCD) is a confrontational lymphoproliferative disorder described by symptoms such as lymph node proliferation, unwarranted secretion of inflammatory cytokines, hyperactive immune system, and in severe cases, multiple organ dysfunction. Interleukin-6 (IL-6) is a pleiotropic cytokine which is involved in a large range of physiological processes in our body such as pro-inflammation, anti-inflammation, differentiation of T-cells and is reported to be a key pathological factor in MCD. In the case of MCD, it was observed that IL-6 is overproduced from T-cells and macrophages which disturb Hepcidin, a vital regulator of iron trafficking in macrophage. The present study endeavour to expound the inhibitor which binds to IL-6 protein receptor with high affinity. Methods: MolegroVirtual Docker software was employed to find the best-established drug from the list of selected inhibitors of IL-6. This compound was subjected to virtual screening against PubChem database to get inhibitors with a very similar structure. These inhibitors were docked to obtain a compound binding with high affinity to the target protein. The established compound and the virtual screened compound were subjected to relative analysis of interactivity energy variables and ADMET profile studies. Results: Among all the selected inhibitors, the virtual screened compound PubChem CID: 101119084 is seen to possess the highest affinity with the target protein. Comparative studies and ADMET analysis further implicate this compound as a better inhibitor of the IL-6 protein. Conclusion: Hence, this compound recognized in the study possesses high potential as an IL-6 inhibitor which might assist in the treatment of Multicentric Castleman Disease and should be examined for its efficiency by in vivo studies.

Background: HIV-1 Integrase (IN) is an important target for the development of the new anti-AIDS drugs. HIV-1 LEDGF/p75 inhibitors, which block the integrase and LEDGF/p75 interaction, have been validated for reduction in HIV-1 viral replicative capacity. Methods: In this work, computational Quantitative Structure-Activity Relationship (QSAR) models were developed for predicting the bioactivity of HIV-1 integrase LEDGF/p75 inhibitors. We collected 190 inhibitors and their bioactivities in this study and divided the inhibitors into nine scaffolds by the method of T-distributed Stochastic Neighbor Embedding (TSNE). These 190 inhibitors were split into a training set and a test set according to the result of a Kohonen’s self-organizing map (SOM) or randomly. Multiple Linear Regression (MLR) models, support vector machine (SVM) models and two consensus models were built based on the training sets by 20 selected CORINA Symphony descriptors. Results: All the models showed a good prediction of pIC50. The correlation coefficients of all the models were more than 0.7 on the test set. For the training set of consensus Model C1, which performed better than other models, the correlation coefficient(r) achieved 0.909 on the training set, and 0.804 on the test set. Conclusion: The selected molecular descriptors show that hydrogen bond acceptor, atom charges and electronegativities (especially π atom) were important in predicting the activity of HIV-1 integrase LEDGF/p75-IN inhibitors.