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

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

Discovery of novel mGluR1 antagonists: a multistep virtual screening approach based on an SVM model and a pharmacophore hypothesis significantly increases the hit rate and enrichment factor

Development of glutamate non-competitive antagonists of mGluR1 (Metabotropic glutamate receptor subtype 1) has increasingly attracted much attention in recent years due to their potential therapeutic application for various nervous disorders. Since there is no crystal structure reported for mGluR1, ligand-based virtual screening (VS) methods, typically pharmacophore-based VS (PB-VS), are often used for the discovery of mGluR1 antagonists. Nevertheless, PB-VS usually suffers a lower hit rate and enrichment factor. In this investigation, we established a multistep ligand-based VS approach that is based on a support vector machine (SVM) classification model and a pharmacophore model. Performance evaluation of these methods in virtual screening against a large independent test set, M-MDDR, show that the multistep VS approach significantly increases the hit rate and enrichment factor compared with the individual SB-VS and PB-VS methods. The multistep VS approach was then used to screen several large chemical libraries including PubChem, Specs, and Enamine. Finally a total of 20 compounds were selected from the top ranking compounds, and shifted to the subsequent in vitro and in vivo studies, which results will be reported in the near future.     

Docking- and pharmacophore-based virtual screening for the identification of novel Mycobacterium tuberculosis protein tyrosine phosphatase B (MptpB) inhibitor with a thiobarbiturate scaffold

Mycobacterium tuberculosis (Mtb) protein tyrosine phosphatase B (MptpB) is an important virulence factor for Mtb that contributes to survival of the bacteria in macrophages. The absence of a human ortholog makes MptpB an attractive target for new therapeutics to treat tuberculosis. MptpB inhibitors could be an effective treatment to overcome emerging TB drug resistance. Adopting a structure-based virtual screening strategy, we successfully identified thiobarbiturate-based drug-like MptpB inhibitor 15 with an IC₅₀ of 22.4 μM, and as a non-competitive inhibitor with a K_i of 24.7 μM. Importantly, not only did it exhibit moderate cell membrane permeability, compound 15 also displayed potent inhibition of intracellular TB growth in the macrophage, making it an excellent lead compound for anti-TB drug discovery. To the best of our knowledge, this novel thiobarbiturate is the first class of MptpB inhibitor reported so far that leveraged docking- and pharmacophore-based virtual screening approaches. The results of preliminary structure-activity relationship demonstrated that compound 15 identified herein was not a singleton and may inspire the design of novel selective and drug-like MptpB inhibitors.     

Virtual screening of small molecules databases for discovery of novel PARP-1 inhibitors: combination of in silico and in vitro studies

Poly(ADP-ribose) polymerase-1 (PARP-1) enzyme has critical roles in DNA replication repair and recombination. Thus, PARP-1 inhibitors play an important role in the cancer therapy. In the current study, we have performed combination of in silico and in vitro studies in order to discover novel inhibitors against PARP-1 target. Structure-based virtual screening was carried out for an available small molecules database. A total of 257,951 ligands from Otava database were screened at the binding pocket of PARP-1 using high-throughput virtual screening techniques. Filtered structures based on predicted binding energy results were then used in more sophisticated molecular docking simulations (i.e. Glide/standard precision, Glide/XP, induced fit docking – IFD, and quantum mechanics polarized ligand docking – QPLD). Potential high binding affinity compounds that are predicted by molecular simulations were then tested by in vitro methods. Computationally proposed compounds as PARP-1 inhibitors (Otava Compound Codes: 7111620047 and 7119980926) were confirmed by in vitro studies. In vitro results showed that compounds 7111620047 and 7119980926 have IC₅₀ values of 0.56 and 63 μM against PARP-1 target, respectively. The molecular mechanism analysis, free energy perturbation calculations using long multiple molecular dynamics simulations for the discovered compounds which showed high binding affinity against PARP-1 enzyme, as well as structure-based pharmacophore development (E-pharmacophore) studies were also studied.     

Evaluation in zebrafish model of the toxicity of rhodamine B‐conjugated crotamine,a peptide potentially useful for diagnostics and therapeutics

Crotamine is defensin‐like cationic peptide from rattlesnake venom that possesses anticancer, antimicrobial, and antifungal properties. Despite these promising biological activities, toxicity is a major concern associated with the development of venom‐derived peptides as therapeutic agents. In the present study, we used zebrafish as a system model to evaluate the toxicity of rhodamine B‐conjugated (RhoB) crotamine derivative. The lethal toxic concentration of RhoB‐crotamine was as low as 4 μM, which effectively kill zebrafish larvae in less than 10 min. With non‐lethal concentrations (<1 μM), crotamine caused malformation in zebrafish embryos, delayed or completely halted hatching, adversely affected embryonic developmental programming, decreased the cardiac functions, and attenuated the swimming distance of zebrafish. The RhoB‐crotamine translocated across vitelline membrane and accumulated in zebrafish yolk sac. These results demonstrate the sensitive responsivity of zebrafish to trial crotamine analogues for the development of novel therapeutic peptides with improved safety, bioavailability, and efficacy profiles.     

A convenient screening test system and a model for thermal germination responses of wild plant seeds: behaviour of model and real seeds in the system

Abstract A convenient test system for screening the thermal germination behaviour of seeds was developed for both basic and applied research in seed germination ecophysiology. Only two temperature-controlled facilities, a test period of about a month and a relatively small number of sample seeds arc needed to obtain information on the thermal-germination parameters of individual seed populations, such as lower or higher limit temperatures, and thermal times required for germination. In the test system, the germination performances of sample seed populations were compared under two temperature regimes: a gradually increasing temperature regime and a gradually decreasing temperature regime, in which the seeds were subjected to gradually changing temperatures in the range of 4 36°C. In order to assess the effects of various values for thermal-germination parameters on the patterns of germination performance in the system, the behaviour of model seeds characterized by a definite set of thermal germination parameters were investigated. Referring to the results of the simulations, the actual germination patterns of some wild-seed populations in the test system were interpreted in terms of thermal-germination parameters.     

A structural model for the inhibition of calpain by calpastatin: crystal structures of the native domain VI of calpain and its complexes with calpastatin peptide and a small molecule inhibitor

The Ca(2+)-dependent cysteine protease calpain along with its endogenous inhibitor calpastatin is widely distributed. The interactions between calpain and calpastatin have been studied to better understand the nature of calpain inhibition by calpastatin, which can aid the design of small molecule inhibitors to calpain. Here we present the crystal structure of a complex between a calpastatin peptide and the calcium-binding domain VI of calpain. DIC19 is a 19 residue peptide, which corresponds to one of the three interacting domains of calpastatin, which is known to interact with domain VI of calpain. We present two crystal structures of DIC19 bound to domain VI of calpain, determined by molecular replacement methods to 2.5A and 2.2A resolution. In the process of crystallizing the inhibitor complex, a new native crystal form was identified which had the homodimer 2-fold axis along a crystallographic axis as opposed to the previously observed dimer in the asymmetric unit. The crystal structures of the native domain VI and its inhibitor PD150606 (3-(4-iodophenyl)-2-mercapto-(Z)-2-propenoic acid) complex were determined with the help of molecular replacement methods to 2.0A and 2.3A resolution, respectively. In addition, we built a homology model for the complex between domain IV and DIA19 peptide of calpastatin. Finally, we present a model for the calpastatin-inhibited calpain.     

The RecA protein as a model molecule for molecular systematic studies of bacteria: Comparison of trees of RecAs and 16S rRNAs from the same species

The evolution of the RecA protein was analyzed using molecular phylogenetic techniques. Phylogenetic trees of all currently available complete RecA proteins were inferred using multiple maximum parsimony and distance matrix methods. Comparison and analysis of the trees reveal that the inferred relationships among these proteins are highly robust. The RecA trees show consistent subdivisions corresponding to many of the major bacterial groups found in trees of other molecules including the , , , , proteobacteria, cyanobacteria, high-GC gram-positives, and the Deinococcus-Thermus group. However, there are interesting differences between the RecA trees and these other trees. For example, in all the RecA trees the proteins from gram-positive species are not monophyletic. In addition, the RecAs of the cyanobacteria consistently group with those of the high-GC gram-positives. To evaluate possible causes and implications of these and other differences phylogenetic trees were generated for small-subunit rRNA sequences from the same (or closely related) species as represented in the RecA analysis. The trees of the two molecules using these equivalent species-sets are highly congruent and have similar resolving power for close, medium, and deep branches in the history of bacteria. The implications of the particular similarities and differences between the trees are discussed. Some of the features that make RecA useful for molecular systematics and for studies of protein evolution are also discussed.     

Molecular modeling,dynamics studies and density functional theory approaches to identify potential inhibitors of SIRT4 protein from Homo sapiens : a novel target for the treatment of type 2 diabetes

Type 2 diabetes is one of the biggest health challenges in the world and WHO projects it to be the 7th leading cause of death in 2030. It is a chronic condition affecting the way our body metabolizes sugar. Insulin resistance is high risk factor marked by expression of Lipoprotein Lipases and Peroxisome Proliferator-Activated Receptor that predisposes to type 2 diabetes. AMP-dependent protein kinase in AMPK signaling pathway is a central sensor of energy status. Deregulation of AMPK signaling leads to inflammation, oxidative stress, and deactivation of autophagy which are implicated in pathogenesis of insulin resistance. SIRT4 protein deactivates AMPK as well as directly inhibits insulin secretion. SIRT4 overexpression leads to dyslipidimeia, decreased fatty acid oxidation, and lipogenesis which are the characteristic features of insulin resistance promoting type 2 diabetes. This makes SIRT4 a novel therapeutic target to control type 2 diabetes. Virtual screening and molecular docking studies were performed to obtain potential ligands. To further optimize the geometry of protein–ligand complexes Quantum Polarized Ligand Docking was performed. Binding Free Energy was calculated for the top three ligand molecules. In view of exploring the stereoelectronic features of the ligand, density functional theory approach was implemented at B3LYP/6-31G* level. 30 ns MD simulation studies of the protein–ligand complexes were done. The present research work proposes ZINC12421989 as potential inhibitor of SIRT4 with docking score (?7.54 kcal/mol), docking energy (?51.34 kcal/mol), binding free energy (?70.21 kcal/mol), and comparatively low energy gap (?0.1786 eV) for HOMO and LUMO indicating reactivity of the lead molecule.