Novel inhibitor design for hemagglutinin against H1N1 influenza virus by core hopping method

A deep learning and docking simulation-based virtual screening strategy enables the rapid identification of HIF-1α pathway activators from a marine natural product database

Artificial Intelligence is hailed as a cutting-edge technology for accelerating drug discovery efforts, and our goal was to validate its potential in predicting pharmacological inhibitors of EGLN1 using a deep learning-based architecture, one of its subsidiaries. Egl nine homolog 1 (EGLN1) inhibition prevents poly ubiquitination-mediated proteosomal destruction HIF-1α. The pharmacological interventions aimed at stabilizing HIF-1α have the potential to be a promising treatment option for a range of human diseases, including ischemic stroke. To unveil a novel EGLN1 inhibitor from marine natural products, a custom-based virtual screening was carried out using a Deep Convolutional Neural Network (DCNN) architecture, docking, and molecular dynamics simulation. The custom DCNN model was optimized and further employed to screen marine natural products from the CMNPD database. The docking was performed as a secondary strategy for screened hits. Molecular dynamics (MD) and molecular mechanics/generalized Born surface area (MM-GBSA) were used to analyze inhibitor binding and identify key interactions. The findings support the claim that deep learning-based virtual screening is a rapid, reliable and accurate method of identifying highly contributing drug candidates (EGLN1 inhibitors). This study demonstrates that deep learning architecture can significantly accelerate drug discovery and development, and provides a solid foundation for using (Z)-2-ethylhex-2-enedioic acid [(Z)-2-ethylhex-2-enedioic acid] as a potential EGLN1 inhibitor for treating various health complications.Communicated by Ramaswamy H. Sarma.

Use Chou's 5-Steps Rule With Different Word Embedding Types to Boost Performance of Electron Transport Protein Prediction Model

Living organisms receive necessary energy substances directly from cellular respiration. The completion of electron storage and transportation requires the process of cellular respiration with the aid of electron transport chains. Therefore, the work of deciphering electron transport proteins is inevitably needed. The identification of these proteins with high performance has a prompt dependence on the choice of methods for feature extraction and machine learning algorithm. In this study, protein sequences served as natural language sentences comprising words. The nominated word embedding-based feature sets, hinged on the word embedding modulation and protein motif frequencies, were useful for feature choosing. Five word embedding types and a variety of conjoint features were examined for such feature selection. The support vector machine algorithm consequentially was employed to perform classification. The performance statistics within the 5-fold cross-validation including average accuracy, specificity, sensitivity, as well as MCC rates surpass 0.95. Such metrics in the independent test are 96.82, 97.16, 95.76 percent, and 0.9, respectively. Compared to state-of-the-art predictors, the proposed method can generate more preferable performance above all metrics indicating the effectiveness of the proposed method in determining electron transport proteins. Furthermore, this study reveals insights about the applicability of various word embeddings for understanding surveyed sequences.

Molecular Docking and Dynamics Simulation Analysis of Thymoquinone and Thymol Compounds from Nigella sativa L. that Inhibit Cag A and Vac A Oncoprotein of Helicobacter pylori: Probable Treatment of H. pylori Infections

BACKGROUND

Helicobacter pylori infection is accountable for most of the peptic ulcer and intestinal cancers. Due to the uprising resistance towards H. pylori infection through the present and common proton pump inhibitors regimens, the investigation of novel candidates is the inevitable issue. Medicinal plants have always been a source of lead compounds for drug discovery. The research of the related effective enzymes linked with this gram-negative bacterium is critical for the discovery of novel drug targets.

OBJECTIVE

The aim of the study is to identify the best candidate to evaluate the inhibitory effect of thymoquinone and thymol against H. pylori oncoproteins, Cag A and Vac A in comparison to the standard drug, metronidazole by using a computational approach.

MATERIALS AND METHODS

The targeted oncoproteins, Cag A and Vac A were retrieved from RCSB PDB. Lipinski's rule and ADMET toxicity profiling were carried out on the phytoconstituents of the N. sativa. The two compounds of N. sativa were further analyzed by molecular docking and MD simulation studies. The reported phytoconstituents, thymoquinone and thymol present in N. sativa were docked with H. pylori Cag A and Vac A oncoproteins. Structures of ligands were prepared using ChemDraw Ultra 10 software and then changed into their 3D PDB structures using Molinspiration followed by energy minimization by using software Discovery Studio client 2.5.

RESULTS

The docking results revealed the promising inhibitory potential of thymoquinone against Cag A and Vac A with docking energy of -5.81 kcal/mole and -3.61kcal/mole, respectively. On the contrary, the inhibitory potential of thymol against Cag A and Vac A in terms of docking energy was -5.37 kcal/mole and -3.94kcal/mole as compared to the standard drug, metronidazole having docking energy of -4.87 kcal/mole and -3.20 kcal/mole, respectively. Further, molecular dynamic simulations were conducted for 5ns for optimization, flexibility prediction, and determination of folded Cag A and Vac A oncoproteins stability. The Cag A and Vac A oncoproteins-TQ complexes were found to be quite stable with the root mean square deviation value of 0.2nm.

CONCLUSION

The computational approaches suggested that thymoquinone and thymol may play an effective pharmacological role to treat H. pylori infection. Hence, it could be summarized that the ligands thymoquinone and thymol bound and interacted well with the proteins Cag A and Vac A as compared to the ligand MTZ. Our study showed that all lead compounds had good interaction with Cag A and Vac A proteins and suggested them to be a useful target to inhibit H. pylori infection.

Potent sialic acid inhibitors that target influenza A virus hemagglutinin

Eradicating influenza A virus (IAV) is difficult, due to its genetic drift and reassortment ability. As the infectious cycle is initiated by the influenza glycoprotein, hemagglutinin (HA), which mediates the binding of virions to terminal sialic acids moieties, HA is a tempting target of anti-influenza inhibitors. However, the complexity of the HA structure has prevented delineation of the structural characterization of the HA protein–ligand complex. Our computational strategy efficiently analyzed > 200,000 records of compounds held in the United States National Cancer Institute (NCI) database and identified potential HA inhibitors, by modeling the sialic acid (SA) receptor binding site (RBS) for the HA structure. Our modeling revealed that compound NSC85561 showed significant antiviral activity against the IAV H1N1 strain with EC₅₀ values ranging from 2.31 to 2.53 µM and negligible cytotoxicity (CC₅₀ > 700 µM). Using the NSC85561 compound as the template to generate 12 derivatives, robust bioassay results revealed the strongest antiviral efficacies with NSC47715 and NSC7223. Virtual screening clearly identified three SA receptor binding site inhibitors that were successfully validated in experimental data. Thus, our computational strategy has identified SA receptor binding site inhibitors against HA that show IAV-associated antiviral activity.

In Silico Design, Synthesis and Evaluation of Novel Series of Benzothiazole- Based Pyrazolidinediones as Potent Hypoglycemic Agents

Background:

The discovery of novel ligand binding domain (LBD) of peroxisome proliferator- activated receptor γ (PPARγ) has recently attracted attention to few research groups in order to develop more potent and safer antidiabetic agents.

Objective:

This study is focused on docking-based design and synthesis of novel compounds combining benzothiazole and pyrazolidinedione scaffold as potential antidiabetic agents.

Methods:

Several benzothiazole-pyrazolidinedione hybrids were synthesized and tested for their in vivo anti-hyperglycemic activity. Interactions profile of title compounds against PPARγ was examined through molecular modelling approach.

Results:

All tested compounds exhibited anti-hyperglycemic activity similar or superior to the reference drug Rosiglitazone. Introducing chlorine atom and alkyl group at position-6 and -5 respectively on benzothiazole core resulted in enhancing the anti-hyperglycemic effect. Docking study revealed that such groups demonstrated favorable hydrophobic interactions with novel LBD Ω- pocket of PPARγ protein.

Conclusion:

Among the tested compounds, N-(6-chloro-5-methylbenzo[d]thiazol-2-yl-4-(4((3,5- dioxopyrazolidin-4-ylidene)methyl)phenoxy)butanamide 5b was found to be the most potent compound and provided valuable insights to further develop novel hybrids as anti-hyperglycemic agents.

Marine compounds targeting the PI3K/Akt signaling pathway in cancer therapy

Cancer is a disease characterized by overproliferation, including that due to transformation, apoptosis disorders, proliferation, invasion, angiogenesis and metastasis, and is one of the deadliest diseases. Currently, conservative chemotherapy is used for cancer treatment due to a lack of effective drugs. The PI3K/Akt signaling pathway plays a very essential role in the pathogenesis of many cancers, and abnormal activation of this pathway leads to abnormal expression of a series of downstream proteins, which ultimately results in the excessive proliferation of cancer cells. Therefore, the PI3K/Akt signaling pathway is a critical target in cancer treatment. Marine drugs have attracted much attention in recent years, and studies have found that many extracts from oceanic animals, plants and microorganisms or their metabolites exert antitumor effects, including antiproliferative effects or the induction of cell cycle arrest, apoptosis or autophagy. However, most anticancer targets and the mechanisms of marine compounds remain unclear. The great potential of the development of marine drugs provides a new direction for cancer treatment. This review focuses on marine compounds that target the PI3K/Akt signaling pathway for the prevention and treatment of cancer and provides comprehensive information for those interested in research on marine drugs.

An Insightful 10-year Recollection Since the Emergence of the 5-steps Rule

OBJECTIVE

One of the most challenging and also the most difficult problems is how to formulate a biological sequence with a vector but considerably keep its sequence order information.

METHODS

To address such a problem, the approach of Pseudo Amino Acid Components or PseAAC has been developed.

RESULTS AND CONCLUSION

It has become increasingly clear via the 10-year recollection that the aforementioned proposal has been indeed very powerful.

Application of Machine Learning Approaches for the Design and Study of Anticancer Drugs

BACKGROUND

Globally the number of cancer patients and deaths are continuing to increase yearly, and cancer has, therefore, become one of the world's highest causes of morbidity and mortality. In recent years, the study of anticancer drugs has become one of the most popular medical topics.

OBJECTIVE

In this review, in order to study the application of machine learning in predicting anticancer drugs activity, some machine learning approaches such as Linear Discriminant Analysis (LDA), Principal components analysis (PCA), Support Vector Machine (SVM), Random forest (RF), k-Nearest Neighbor (kNN), and Naïve Bayes (NB) were selected, and the examples of their applications in anticancer drugs design are listed.

RESULTS

Machine learning contributes a lot to anticancer drugs design and helps researchers by saving time and is cost effective. However, it can only be an assisting tool for drug design.

CONCLUSION

This paper introduces the application of machine learning approaches in anticancer drug design. Many examples of success in identification and prediction in the area of anticancer drugs activity prediction are discussed, and the anticancer drugs research is still in active progress. Moreover, the merits of some web servers related to anticancer drugs are mentioned.

RAACBook: a web server of reduced amino acid alphabet for sequence-dependent inference by using Chou's five-step rule

By reducing amino acid alphabet, the protein complexity can be significantly simplified, which could improve computational efficiency, decrease information redundancy and reduce chance of overfitting. Although some reduced alphabets have been proposed, different classification rules could produce distinctive results for protein sequence analysis. Thus, it is urgent to construct a systematical frame for reduced alphabets. In this work, we constructed a comprehensive web server called RAACBook for protein sequence analysis and machine learning application by integrating reduction alphabets. The web server contains three parts: (i) 74 types of reduced amino acid alphabet were manually extracted to generate 673 reduced amino acid clusters (RAACs) for dealing with unique protein problems. It is easy for users to select desired RAACs from a multilayer browser tool. (ii) An online tool was developed to analyze primary sequence of protein. The tool could produce K-tuple reduced amino acid composition by defining three correlation parameters (K-tuple, g-gap, λ-correlation). The results are visualized as sequence alignment, mergence of RAA composition, feature distribution and logo of reduced sequence. (iii) The machine learning server is provided to train the model of protein classification based on K-tuple RAAC. The optimal model could be selected according to the evaluation indexes (ROC, AUC, MCC, etc.). In conclusion, RAACBook presents a powerful and user-friendly service in protein sequence analysis and computational proteomics. RAACBook can be freely available at http://bioinfor.imu.edu.cn/raacbook.

Database URL: http://bioinfor.imu.edu.cn/raacbook

The Cooperative Effect between Polybasic Region (PBR) and Polysialyltransferase Domain (PSTD) within Tumor-Target Polysialyltranseferase ST8Sia II

ST8Sia II (STX) is a highly homologous mammalian polysialyltransferase (polyST), which is a validated tumor-target in the treatment of cancer metastasis reliant on tumor cell polysialylation. PolyST catalyzes the synthesis of α2,8-polysialic acid (polySia) glycans by carrying out the activated CMP-Neu5Ac (Sia) to N- and O-linked oligosaccharide chains on acceptor glycoproteins. In this review article, we summarized the recent studies about intrinsic correlation of two polybasic domains, Polysialyltransferase domain (PSTD) and Polybasic region (PBR) within ST8Sia II molecule, and suggested that the critical amino acid residues within the PSTD and PBR motifs of ST8Sia II for polysialylation of Neural cell adhesion molecules (NCAM) are related to ST8Sia II activity. In addition, the conformational changes of the PSTD domain due to point mutations in the PBR or PSTD domain verified an intramolecular interaction between the PBR and the PSTD. These findings have been incorporated into Zhou's NCAM polysialylation/cell migration model, which will provide new perspectives on drug research and development related to the tumor-target ST8Sia II.

Progress in the Mechanism and Clinical Application of Cilostazol

Cilostazol is a unique platelet inhibitor that has been used clinically for more than 20 years. As a phosphodiesterase type III inhibitor, cilostazol is capable of reversible inhibition of platelet aggregation and vasodilation, has antiproliferative effects, and is widely used in the treatment of peripheral arterial disease, cerebrovascular disease, percutaneous coronary intervention, etc. This article briefly reviews the pharmacological mechanisms and clinical application of cilostazol.

Identifying Cancer Targets Based on Machine Learning Methods via Chou's 5-steps Rule and General Pseudo Components

In recent years, the successful implementation of human genome project has made people realize that genetic, environmental and lifestyle factors should be combined together to study cancer due to the complexity and various forms of the disease. The increasing availability and growth rate of 'big data' derived from various omics, opens a new window for study and therapy of cancer. In this paper, we will introduce the application of machine learning methods in handling cancer big data including the use of artificial neural networks, support vector machines, ensemble learning and naïve Bayes classifiers.

Investigation of Phospholipase Cγ1 Interaction with SLP76 Using Molecular Modeling Methods for Identifying Novel Inhibitors

The enzyme phospholipase C gamma 1 (PLCγ1) has been identified as a potential drug target of interest for various pathological conditions such as immune disorders, systemic lupus erythematosus, and cancers. Targeting its SH3 domain has been recognized as an efficient pharmacological approach for drug discovery against PLCγ1. Therefore, for the first time, a combination of various biophysical methods has been employed to shed light on the atomistic interactions between PLCγ1 and its known binding partners. Indeed, molecular modeling of PLCγ1 with SLP76 peptide and with previously reported inhibitors (ritonavir, anethole, daunorubicin, diflunisal, and rosiglitazone) facilitated the identification of the common critical residues (Gln805, Arg806, Asp808, Glu809, Asp825, Gly827, and Trp828) as well as the quantification of their interaction through binding energies calculations. These features are in agreement with previous experimental data. Such an in depth biophysical analysis of each complex provides an opportunity to identify new inhibitors through pharmacophore mapping, molecular docking and MD simulations. From such a systematic procedure, a total of seven compounds emerged as promising inhibitors, all characterized by a strong binding with PLCγ1 and a comparable or higher binding affinity to ritonavir (∆G_bind < -25 kcal/mol), one of the most potent inhibitor reported till now.

Advances in Electrochemistry for Monitoring Cellular Chemical Flux

The transport of organic and inorganic molecules, along with inorganic ions across the plasma membrane results in chemical fluxes that reflect the cellular function in healthy and diseased states. Measurement of these chemical fluxes enables the characterization of protein function and transporter stoichiometry, characterization of a single cell and embryo viability prior to implantation, and screening of pharmaceutical agents. Electrochemical sensors emerge as sensitive and non-invasive tools for measuring chemical fluxes immediately outside the cells in the boundary layer, that are capable of monitoring a diverse range of transported analytes including inorganic ions, gases, neurotransmitters, hormones, and pharmaceutical agents. Used on their own or in combination with other methods, these sensors continue to expand our understanding of the function of rare cells and small tissues. Advances in sensor construction and detection strategies continue to improve sensitivity under physiological conditions, diversify analyte detection, and increase throughput. These advances will be discussed in the context of addressing technical challenges to measuring chemical flux in the boundary layer of cells and measuring the resultant changes to the chemical concentration in the bulk media.

In silico minimalist approach to study 2D HP protein folding into an inhomogeneous space mimicking osmolyte effect: First trial in the search of foldameric backbones

We have employed our bioinformatics workbench, named Evolution, a Multi-Agent System based architecture with lattice-bead-models, evolutionary-algorithms, and correlated-networks as inhomogeneous spaces, with different correlation lengths, mimicking osmolyte effect (molecular crowding), to in silico survey protein folding. Resolution is with hydrophobic-polar (H-P) sequences in inhomogeneous 2D square lattices, since general biophysicochemical trends consider i) that the backbone is one of the major components responsible for protein folding and ii) osmolyte effect plays an important role to better folding kinetics and reach deeper optima. We have designed foldamers, as square n × n (n = 3, 4, 5, 6) arrays of hydrophobic cores stabilized by H⋯H contacts, attached through short PP (P₂) or long PPPP (P₄) loops, giving rise to 8 sequences (S₁ to S₈) with known optimal scores. Designed sequences were folded into different inhomogeneous spaces and indeed crowded media induced deeper optima, being crowding necessary to best fold, but the space should be enough constrained to induce folding without banning chain movement. The constrained space plays an important role to reach the optimal structure, depending on designed foldamer sequence size, for an optimal correlation length, implying that media affects the folding pathways as happens in real systems. Designed structures were found, moreover, they undergo to degenerated states, both folding states could survey considering i) backbone information and ii) osmolyte effect. In nature, the proteins fold in different structures aiming to reach a global minimum, but a local minimum could be enough to the protein to be functional or dysfunctional.

Metabolism of Oxalate in Humans: A Potential Role Kynurenine Aminotransferase/Glutamine Transaminase/Cysteine Conjugate Beta-lyase Plays in Hyperoxaluria

Hyperoxaluria, excessive urinary oxalate excretion, is a significant health problem worldwide. Disrupted oxalate metabolism has been implicated in hyperoxaluria and accordingly, an enzymatic disturbance in oxalate biosynthesis can result in the primary hyperoxaluria. Alanine glyoxylate aminotransferase-1 and glyoxylate reductase, the enzymes involving glyoxylate (precursor for oxalate) metabolism, have been related to primary hyperoxalurias. Some studies suggest that other enzymes such as glycolate oxidase and alanine glyoxylate aminotransferase-2 might be associated with primary hyperoxaluria as well, but evidence of a definitive link is not strong between the clinical cases and gene mutations. There are still some idiopathic hyperoxalurias, which require a further study for the etiologies. Some aminotransferases, particularly kynurenine aminotransferases, can convert glyoxylate to glycine. Based on biochemical and structural characteristics, expression level, subcellular localization of some aminotransferases, a number of them appear able to catalyze the transamination of glyoxylate to glycine more efficiently than alanine glyoxylate aminotransferase-1. The aim of this minireview is to explore other undermining causes of primary hyperoxaluria and stimulate research toward achieving a comprehensive understanding of underlying mechanisms leading to the disease. Herein, we reviewed all aminotransferases in the liver for their functions in glyoxylate metabolism. Particularly, kynurenine aminotransferase-I and III were carefully discussed regarding their biochemical and structural characteristics, cellular localization, and enzyme inhibition. Kynurenine aminotransferase-III is, so far, the most efficient putative mitochondrial enzyme to transaminate glyoxylate to glycine in mammalian livers, might be an interesting enzyme to look over in hyperoxaluria etiology of primary hyperoxaluria and should be carefully investigated for its involvement in oxalate metabolism.

Toward of Safer Phenylbutazone Derivatives by Exploration of Toxicity Mechanism

A drug design for safer phenylbutazone was been explored by reactivity and docking studies involving single electron transfer mechanism, as well as toxicological predictions. Several approaches about its structural properties were performed through quantum chemistry calculations at the B3LYP level of theory, together with the 6-31+G(d,p) basis sets. Molecular orbital and ionization potential were associated to electron donation capacity. The spin densities contribution showed a preferential hydroxylation at the para-positions of phenyl ring when compared to other positions. In addition, on electron abstractions the aromatic hydroxylation has more impact than alkyl hydroxylation. Docking studies indicate that six structures 1, 7, 8 and 13–15 have potential for inhibiting human as well as murine COX-2, due to regions showing similar intermolecular interactions to the observed for the control compounds (indomethacin and refecoxib). Toxicity can be related to aromatic hydroxylation. In accordance to our calculations, the derivatives here proposed are potentially more active as well safer than phenylbutazone and only structures 8 and 13–15 were the most promising. Such results can explain the biological properties of phenylbutazone and support the design of potentially safer candidates.

Novel antiviral drug discovery strategies to tackle drug-resistant mutants of influenza virus strains

INTRODUCTION

The emergence of drug-resistant influenza virus strains highlights the need for new antiviral therapeutics to combat future pandemic outbreaks as well as continuing seasonal cycles of influenza. Areas covered: This review summarizes the mechanisms of current FDA-approved anti-influenza drugs and patterns of resistance to those drugs. It also discusses potential novel targets for broad-spectrum antiviral drugs and recent progress in novel drug design to overcome drug resistance in influenza. Expert opinion: Using the available structural information about drug-binding pockets, research is currently underway to identify molecular interactions that can be exploited to generate new antiviral drugs. Despite continued efforts, antivirals targeting viral surface proteins like HA, NA, and M2, are all susceptible to developing resistance. Structural information on the internal viral polymerase complex (PB1, PB2, and PA) provides a new avenue for influenza drug discovery. Host factors, either at the initial step of viral infection or at the later step of nuclear trafficking of viral RNP complex, are being actively pursued to generate novel drugs with new modes of action, without resulting in drug resistance.

Hydrogen bond analysis of the EGFR-ErbB3 heterodimer related to non-small cell lung cancer and drug resistance

Lung cancer is the predominant cause of cancer deaths on a worldwide scale. A mutation in the epidermal growth factor receptor (EGFR) can cause non-small cell lung cancer (NSCLC). The L858R one-point mutation in exon 21 in EGFR is the most prevalent in NSCLC. For over 60% of EGFR-muted NSCLC, another mutation T790M can cause drug resistance. In this paper, we consider EGFR and ErbB3 heterodimers involving three structures of EGFR, wild-type, with L858R mutation, and with L858R and T790M mutations. We perform molecular dynamics (MD) simulations to analyze hydrogen bonds in all three instances. The hydrogen bonds contribute to the conformational stability of the protein and molecular recognition. Several other parameters are also investigated in the present study, which reveals significant changes in the dimer at different levels of mutation. The knowledge and results obtained from this study lead to useful insight into the mechanism of NSCLC drug resistance.

Small Molecule Inhibitor that Stabilizes the Autoinhibited Conformation of the Oncogenic Tyrosine Phosphatase SHP2

Genetic mutations in the phosphatase PTPN11 (SHP2) are associated with childhood leukemias. These mutations cause hyperactivation of SHP2 due to the disruption of the autoinhibitory conformation. By targeting the activation-associated protein conformational change, we have identified an SHP2 inhibitor ( E)-1-(1-(5-(3-(2,4-dichlorophenyl)acryloyl)-2-ethoxy-4-hydroxybenzyl)-1,2,5,6-tetrahydropyridin-3-yl)-1 H-benzo[ d]imidazol-2(3 H)-one (LY6, 1) using computer-aided drug design database screening combined with cell-based assays. This compound inhibited SHP2 with an IC₅₀ value of 9.8 μM, 7-fold more selective for SHP2 than the highly related SHP1. Fluorescence titration, thermal shift, and microscale thermophoresis quantitative binding assays confirmed its direct binding to SHP2. This compound was further verified to effectively inhibit SHP2-mediated cell signaling and proliferation. Furthermore, mouse and patient leukemia cells with PTPN11 activating mutations were more sensitive to this inhibitor than wild-type cells. This small molecule SHP2 inhibitor has a potential to serve as a lead compound for further optimization studies to develop novel anti-SHP2 therapeutic agents.

Exploration of synthetic multifunctional amides as new therapeutic agents for Alzheimer's disease through enzyme inhibition, chemoinformatic properties, molecular docking and dynamic simulation insights

A new series of multifunctional amides has been synthesized having moderate enzyme inhibitory potentials and mild cytotoxicity. 2-Furyl(1-piperazinyl)methanone (1) was coupled with 3,5-dichloro-2-hydroxybenzenesulfonyl chloride (2) to form {4-[(3,5-dichloro-2-hydroxyphenyl)sulfonyl]-1-piperazinyl}(2-furyl)methanone (3). Different elecrophiles were synthesized by the reaction of various un/substituted anilines (4a-o) with 2-bromoacetylbromide (5), 2‑bromo‑N-(un/substituted-phenyl)acetamides (6a-o). Further, equimolar ratios of 3 and 6a-o were allowed to react in the presence of K₂CO₃ in acetonitrile to form desired multifunctional amides (7a-o). The structural confirmation of all the synthesized compounds was carried out by their EI-MS, IR, ¹H NMR and ¹³C NMR spectral data. Enzyme inhibition activity was performed against acetyl and butyrylcholinestrase enzymes, whereby 7e showed very good activity having IC₅₀ value of 5.54 ± 0.03 and 9.15 ± 0.01 μM, respectively, relative to eserine, a reference standard. Hemolytic activity of the molecules was checked to asertain their cytotoxicity towards red blood cell membrance and it was observed that most of the compounds were not toxic up to certain range. Moreover, chemoinformatic protepties and docking simulation results also showed the significance of 7e as compared to other compounds. Based on in vitro and in silico analysis 7e could be used as a template for the development of new drugs against Alzheimer's disease.

A Comprehensive In Silico Method to Study the QSTR of the Aconitine Alkaloids for Designing Novel Drugs

A combined in silico method was developed to predict potential protein targets that are involved in cardiotoxicity induced by aconitine alkaloids and to study the quantitative structure⁻toxicity relationship (QSTR) of these compounds. For the prediction research, a Protein-Protein Interaction (PPI) network was built from the extraction of useful information about protein interactions connected with aconitine cardiotoxicity, based on nearly a decade of literature and the STRING database. The software Cytoscape and the PharmMapper server were utilized to screen for essential proteins in the constructed network. The Calcium-Calmodulin-Dependent Protein Kinase II alpha (CAMK2A) and gamma (CAMK2G) were identified as potential targets. To obtain a deeper insight on the relationship between the toxicity and the structure of aconitine alkaloids, the present study utilized QSAR models built in Sybyl software that possess internal robustness and external high predictions. The molecular dynamics simulation carried out here have demonstrated that aconitine alkaloids possess binding stability for the receptor CAMK2G. In conclusion, this comprehensive method will serve as a tool for following a structural modification of the aconitine alkaloids and lead to a better insight into the cardiotoxicity induced by the compounds that have similar structures to its derivatives.

Hyphenated 3D-QSAR statistical model-drug repurposing analysis for the identification of potent neuraminidase inhibitor

The Influenza A virus is one of the principle causes of respiratory illness in human. The surface glycoprotein of the influenza virus, neuraminidase (NA), has a vital role in the release of new viral particle and spreads infection in the respiratory tract. It has been long recognized as a valid drug target for influenza A virus infection. Oseltamivir is used as a standard drug of choice for the treatment of influenza. However, the emergence of mutants with novel mutations has increased the resistance to potent NA inhibitor. In the present investigation, we have employed computer-assisted combinatorial techniques in the screening of 8621 molecules from Drug Bank to find potent NA inhibitors. A three-dimensional pharmacophore model was generated from the previously reported 28 carbocylic influenza NA inhibitors along with oseltamivir using PHASE module of Schrödinger Suite. The model generated consists of one hydrogen bond acceptor (A), one hydrogen bond donors (D), one hydrophobic group (H), and one positively charged group (P), ADHP. The hypothesis was further validated for its integrity and significance using enrichment analysis. Subsequently, an atom-based 3D-QSAR model was built using the common pharmacophore hypothesis (CPH). The developed 3D-QSAR model was found to be statistically significant with R² value of 0.9866 and Q² value of 0.7629. Further screening was accomplished using three-stage docking process using the Glide algorithm. The resultant lead molecules were examined for its drug-like properties using the Qikprop algorithm. Finally, the calculated pIC₅₀ values of the lead compounds were validated by the AutoQSAR algorithm. Overall, the results from our analysis highlights that lisinopril (DB00722) is predicted to bind better with NA than currently approved drug. In addition, it has the best match in binding geometry conformations with the existing NA inhibitor. Note that the antiviral activity of lisinopril is reported in the literature. However, our paper is the first report on lisinopril activity against influenza A virus infection. These results are envisioned to help design the novel NA inhibitors with an increased antiviral efficacy.

Structural insights of dipeptidyl peptidase-IV inhibitors through molecular dynamics-guided receptor-dependent 4D-QSAR studies

Dipeptidyl peptidase-IV (DPP-IV) inhibitors are promising antidiabetic agents. Currently, several DPP-IV inhibitors have been approved for therapeutic use in diabetes mellitus. Receptor-dependent 4D-QSAR is comparatively a new approach which uses molecular dynamics simulations to generate conformational ensemble profiles of compounds representing a dynamic state of compounds at a target's binding site. This work describes a receptor-dependent 4D-QSAR study on triazolopiperazine derivatives. QSARINS multiple linear regression method was adopted to generate 4D-QSAR models. A model with 9 variables was found to have better predictive accuracy with [Formula: see text], [Formula: see text] (leave-one-out) = 0.592 and [Formula: see text] predicted = 0.597. The location of these 9 variables at the binding site of DPP-IV revealed the importance of the residues Val711, Tyr662, Tyr666, Val202, Asp200 and Thr199 in making critical interactions with DPP-IV inhibitors. The study of these critical interactions revealed the structural features required in DPP-IV inhibitors. Thus, in this study the importance of a halogen substituent on a phenyl ring, the extent of substitution on the triazolopiperazine ring, the presence of an ionizable amino group and the presence of a hydrophobic substituent that can bind deeper in binding pocket of DPP-IV were revealed.

Virtual screening, optimization, and identification of a novel specific PTP-MEG2 Inhibitor with potential therapy for T2DM

Megakaryocyte protein tyrosine phosphatase 2 (PTP-MEG2) is a tyrosine phosphatase expressed in megakaryocytic cells, and causes insulin sensitization when down regulated. Therefore, specific inhibitors of PTP-MEG2 are potential candidates for novel Type 2 Diabetes (T2DM)therapy. In this study, we discovered PTP-MEG2 inhibitors using high throughput and virtual screening (HTS/VS) and structural optimization in silicon. Eight compound-candidates were identified from the interactions with PTP-MEG2, protein tyrosine phosphatase 1B (PTP1B) and T cell protein tyrosine phosphatase (TCPTP). Results from enzymatic assays show compounds 4a and 4b inhibited PTP-MEG2 activity with an IC50 of 3.2 μM and 4.3 μM, respectively. Further, they showed a 7.5 and 5.5 fold change against PTP1B and TCPTP, respectively. We propose compounds 4a and 4b are PTP-MEG2 inhibitors with potential therapeutic use in T2DM treatment.

Small molecular floribundiquinone B derived from medicinal plants inhibits acetylcholinesterase activity

Being a neurodegenerative disorder, Alzheimer's disease (AD) is the one of the most terrible diseases. And acetylcholinesterase (AChE) is considered as an important target for treating AD. Acetylcholinesterase inhibitors (AChEI) are considered to be one of the effective drugs for the treatment of AD. The aim of this study is to find a novel potential AChEI as a drug for the treatment of AD. In this study, instead of using the synthetic compounds, we used those extracted from plants to investigate the interaction between floribundiquinone B (FB) and AChE by means of both the experimental approach such as fluorescence spectra, ultraviolet-visible (UV-vis) absorption spectrometry, circular dichroism (CD) and the theoretical approaches such as molecular docking. The findings reported here have provided many useful clues and hints for designing more effective and less toxic drugs against Alzheimer's disease.

Structure based design, synthesis and activity studies of small hybrid molecules as HDAC and G9a dual inhibitors

Aberrant enzymatic activities or expression profiles of epigenetic regulations are therapeutic targets for cancers. Among these, histone 3 lysine 9 methylation (H3K9Me2) and global de-acetylation on histone proteins are associated with multiple cancer phenotypes including leukemia, prostatic carcinoma, hepatocellular carcinoma and pulmonary carcinoma. Here, we report the discovery of the first small molecule capable of acting as a dual inhibitor targeting both G9a and HDAC. Our structure based design, synthesis, and screening for the dual activity of the small molecules led to the discovery of compound 14 which displays promising inhibition of both G9a and HDAC in low micro-molar range in cell based assays.

Molecular dynamics and protein interaction studies of lipopeptide (Iturin A) on α- amylase of Spodoptera litura

The small mottled willow moth (Spodoptera litura) is one of the best-known agricultural pest insects. To understand the insecticidal activity, we have selected iturin A compound produced by Bacillus amyloliquefaciens RHNK22 which showed the strongest and most common inhibitory effect on the Spodoptera litura protein. In this work we have identified the action of iturin A on α- amylase is a major digestive enzyme of Spodoptera litura using docking studies. A 3D model of α- amylase from Spodoptera litura was generated using 2HPH as a template with the help of Modeller7v7. With the aid of the molecular mechanics and molecular dynamics methods, the final model is obtained and is further checked by Procheck and Verify 3D graph programs, which showed that the final refined model is reliable. With this model, a adjustable docking study was performed with iturin A using GOLD software. The results indicated that ARG 18, THR15, LEU42 in α- amylase are important determinant residues in binding as they have strong hydrogen bonding interactions with iturin A. These hydrogen binding interactions play an important role for the stability of the complex.

The effect of glycosylation on the transferrin structure: A molecular dynamic simulation analysis

Transferrins have been defined by the highly cooperative binding of iron and a carbonate anion to form a Fe-CO3-Tf ternary complex. As such, the layout of the binding site residues affects transferrin function significantly; In contrast to N-lobe, C-lobe binding site of the transferrin structure has been less characterized and little research which surveyed the interaction of carbonate with transferrin in the C-lobe binding site has been found. In the present work, molecular dynamic simulation was employed to gain access into the molecular level understanding of carbonate binding site and their interactions in each lobe. Residues responsible for carbonate binding of transferrin structure were pointed out. In addition, native human transferrin is a glycoprotein that two N-linked complex glycan chains located in the C-lobe. Usually, in the molecular dynamic simulation for simplifying, glycan is removed from the protein structure. Here, we explore the effect of glycosylation on the transferrin structure. Glycosylation appears to have an effect on the layout of the binding site residue and transferrin structure. On the other hand, sometimes the entire transferrin formed by separated lobes that it allows the results to be interpreted in a straightforward manner rather than more parameters required for full length protein. But, it should be noted that there are differences between the separated lobe and full length transferrin, hence, a comparative analysis by the molecular dynamic simulation was performed to investigate such structural variations. Results revealed that separation in C-lobe caused a significant structural variation in comparison to N-lobe. Consequently, the separated lobes and the full length one are different, showing the importance of the interlobe communication and the impact of the lobes on each other in the transferrin structure.

Identification of Damaging nsSNVs in HumanERCC2 Gene

The hERCC2 gene is an important DNA repair molecule for initiating Cutaneous melanoma (CM). Therefore, it is advisable to study the possible functional SNVs in hERCC2. To achieve this goal, we collected total 2, 253 SNVs in hERCC2from the EMBL website, of which 303 are non-synonymous single nucleotide variants (nsSNVs). Then, SIFT and PolyPhen were used to predict the damaging nsSNVs, and four nsSNVs (rs368866996, rs377739017, rs370819591, and rs121913022) were suggested to be damaging mutations. Since I-Mutant2.0 showed a decrease in stability for the mutants containing each of the four nsSNVs, a 3D protein structure was modeled. Based on the comparison of the energy after minimization, RMSD and stabilizing residues between the native and mutant proteins' structure, rs121913022 was proposed to be the most damaging variant among the nsSNVs in hERCC2 gene by decreasing the stability of protein. The mutant G713R of hERCC2 protein caused by rs121913022 was found to have less expression level than native hERCC2 protein in melanoma cells. These results suggest that rs121913022 may have potentially important clinical and drug target implications.

Identification of CD4+ T-cell epitope and investigation of HLA distribution for the immunogenic proteins of Burkholderia pseudomallei using in silico approaches - A key vaccine development strategy for melioidosis

Melioidosis is a serious infectious diseases affecting multi-organ system in humans with high mortality rate. The disease is caused by the bacterium, Burkholderia pseudomallei and it is intrinsically resistant to many antibiotics. Thus, there is an urgent need for protective vaccine against B. pseudomallei; which may reduce morbidity and mortality in endemic areas. The identification of peptides that bind to major histocompatibility complex II class helps in understanding the nature of immune response and identifying T-cell epitopes for the design of new vaccines. Previous studies indicate that, ompA, bipB, fliC and groEL proteins of B. pseudomallei stimulate CD4+ T-cell immune response and act as protective immunogens. However, the data for CD4+ T-cell epitopes of these immunogenic proteins are very limited. Hence, in this present study we attempted to identify CD4+ T-cell epitopes in B. pseudomallei immunogenic proteins using in silico approaches. We did population coverage analysis for these identified epitopic core sequences to identify individuals in endemic areas expected to respond to a given set of these epitopes on the basis of HLA genotype frequencies. We observed that eight epitopic core sequences, two from each immunogenic protein, were associated with the maximum number of HLA-DR binding alleles. These eight peptides are found to be immunogenic in more than 90% of population in endemic areas considered. Thus, these eight peptides containing epitopic core sequences may act as probable vaccine candidates and they may be considered for the development of epitope-based vaccines for melioidosis.

In silico identification of putative bifunctional Plk1 inhibitors by integrative virtual screening and structural dynamics approach

Polo like kinase (Plk1) is a master regulator of cell cycle and considered as next generation antimitotic target in human. As Plk1 predominantly expresses in the dividing cells with a much higher expression in cancerous cells, it serves as a discriminative target for cancer therapeutics. Here we implied a novel and promising integrative strategy to identify "bifunctional" Plk1 inhibitors that compete simultaneously with ATP and substrate for their binding sites. We integrated structure-based virtual screening (SBVS) and molecular dynamics simulations with emphasis on unique structural properties of Plk1. Through screening of 20,000 compounds, nearly ~2000 hits were enriched and subjected to SBVS against ATP and substrate binding sites of Plk1. Subsequently, on the basis of their binding abilities to Plk1 kinase and polo box domains, filtration of candidate hits resulted in the isolation of 26 compounds. By exclusion of close analogs or isomers, 10 unique compounds were selected for detailed study. A representative compound was subjected to molecular dynamics simulation assay to have deep structural insights and to gauge critical structural crunch for inhibitor binding against kinase and polo box domains. Our integrative approach may complement high-throughput screening and identify bifunctional Plk1 inhibitors that may contribute in selective targeting of Plk1 to elicit desired biological process.