Virtual screening of a natural compound library at orthosteric and allosteric binding sites of the neurotensin receptor

Abstract

Molecular dynamics (MD) simulation using the AMBER force field has been performed on the neurotensin (NT) receptor, a class A type G-protein-coupled receptor in its activated conformation co-crystallized with the non-peptide agonists. For structure-based hit molecule identification via natural chemical compound library, orthosteric sites on NT receptor have been mapped by docking using AutoDock4.0 and Vina with the known agonists and antagonists SR48692, SR142948, ML301 and ML314 of the receptor. Furthermore, clustering analysis on the MD trajectories by SIMULAID has been performed to filter receptor conformations for the allosteric binders from the Otava natural compound library. Comparative mappings of contrasting binding region patterns have been done between the crystal structure orthosteric sites as well as the binding regions in the SIMULAID-based cluster center conformations from MD trajectories with the FTmap server using the small organic molecule fragments as the probes. The distinct binding region in the cluster-based conformations in the extracellular region of the receptor has been identified for targeted docking by Otava natural chemical compound library using AutoDock4.0 and Vina docking suites to obtain putative allosteric binders. A group of compounds from the Otava library has been identified as showing high free energy in both AutoDock4.0 and Vina docking suites. Biophysical assessments on the natural compound computational hit molecules will be done to identify lead structures from the hit molecules.

Communicated by Ramaswamy H. Sarma     

New dipeptidyl peptidase 4 inhibitors among adamantane derivatives

Fifteen adamantane derivatives were synthesized. Preliminary evaluation of their potential as dipeptidyl peptidase 4 (DPP-4) inhibitors was performed in silico by the Microcosm informational technology, PASS system, and docking in AutoDock Vina. The DPP-4 inhibition was studied in vitro. The selectivity of action of the most active compounds was studied by the direct inhibition of human plasma DPP-4 and recombinant human DPP-8. The highest activity was found for the compounds containing a nitrogen atom in the β-position of the side chain, namely, derivatives of adamantane carboxylic acid and N-(3-adamantyl-allyl) thiourea. We demonstrated that the most active compound of the series, 3,5-dimethyladamantane 1-carboxamide, was a selective DPP-4 inhibitor with IC₅₀ 53.94 μM.     

Virtual Screening for HIV Protease Inhibitors: A Comparison of AutoDock 4 and Vina

Background

The AutoDock family of software has been widely used in protein-ligand docking research. This study compares AutoDock 4 and AutoDock Vina in the context of virtual screening by using these programs to select compounds active against HIV protease.

Methodology/Principal Findings

Both programs were used to rank the members of two chemical libraries, each containing experimentally verified binders to HIV protease. In the case of the NCI Diversity Set II, both AutoDock 4 and Vina were able to select active compounds significantly better than random (AUC = 0.69 and 0.68, respectively; p<0.001). The binding energy predictions were highly correlated in this case, with r = 0.63 and ι = 0.82. For a set of larger, more flexible compounds from the Directory of Universal Decoys, the binding energy predictions were not correlated, and only Vina was able to rank compounds significantly better than random.

Conclusions/Significance

In ranking smaller molecules with few rotatable bonds, AutoDock 4 and Vina were equally capable, though both exhibited a size-related bias in scoring. However, as Vina executes more quickly and is able to more accurately rank larger molecules, researchers should look to it first when undertaking a virtual screen.     

In silico investigation of ACE2 and the main protease of SARS-CoV-2 with phytochemicals from Myristica fragrans (Houtt.) for the discovery of a novel COVID-19 drug

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the virus that causes coronavirus disease 2019 (COVID-19), is a new coronavirus strain that was first reported in December 2019 in Wuhan, China. A specific treatment for COVID-19 has yet to be identified. Potential therapeutic targets include SARS-CoV-2 main protease (Mpro) and the SARS-CoV-2 spike-ACE2 interaction. Molecular docking, molecular dynamics (MD), solvent screening for the extraction of the specified compounds, and prediction of the drug properties of certain molecules were the methods used in this study to investigate compounds from the medicinal plant Myristica fragrans, which is one of twelve herbs in Prasachandaeng remedy (PSD). ArgusLab, AutoDock Vina, and AutoDock were used to perform docking tasks. The examined ligands were compared with panduratin A as a standard (Kanjanasirirat et al., 2020), which is a promising medicinal plant molecule for the treatment of COVID-19. Molecular docking revealed that malabaricones B and C and licarins A, B and C bound to SARS-CoV-2/ACE2 and SARS-CoV-2 Mpro with low binding energies compared to that of the standard ligand. Furthermore, appropriate solvent usage is important. Acetone was selected by COSMOquick software for compound extraction in this investigation because it can extract large amounts of all five of the abovementioned M. fragrans compounds. Furthermore, the drug-like properties of these compounds were studied utilizing the Lipinski, Veber, and Ghose criteria. The results revealed that these M. fragrans compounds have potential as effective medicines to combat the COVID-19 pandemic. However, to assess the therapeutic potential of these ligands, additional research is needed, which will use our findings as a foundation.     

Curcumin binds in silico to anti-cancer drug target enzyme MMP-3 (human stromelysin-1) with affinity comparable to two known inhibitors of the enzyme

In silico interaction of curcumin with the enzyme MMP-3 (human stromelysin-1) was studied by molecular docking using AutoDock 4.2 as the docking software application. AutoDock 4.2 software serves as a valid and acceptable docking application to study the interactions of small compounds with proteins. Interactions of curcumin with MMP-3 were compared to those of two known inhibitors of the enzyme, PBSA and MPPT. The calculated free energy of binding (ΔG binding) shows that curcumin binds with affinity comparable to or better than the two known inhibitors. Binding interactions of curcumin with active site residues of the enzyme are also predicted. Curcumin appears to bind in an extendended conformation making extensive VDW contacts in the active site of the enzyme. Hydrogen bonding and pi-pi interactions with key active site residues is also observed. Thus, curcumin can be considered as a good lead compound in the development of new inhibitors of MMP-3 which is a potential target of anticancer drugs. The results of these studies can serve as a starting point for further computational and experimental studies.     

An Integrated Study of Tyrosinase Inhibition by Rutin: Progress using a Computational Simulation

Abstract

Tyrosinase inhibition studies have recently gained the attention of researchers due to their potential application values. We simulated docking (binding energies for AutoDock Vina: ?9.1 kcal/mol) and performed a molecular dynamics simulation to verify docking results between tyrosinase and rutin. The docking results suggest that rutin mostly interacts with histidine residues located in the active site. A 10 ns molecular dynamics simulation showed that one copper ion at the tyrosinase active site was responsible for the interaction with rutin. Kinetic analyses showed that rutin-mediated inactivation followed a first-order reaction and mono- and biphasic rate constants occurred with rutin. The inhibition was a typical competitive type with K_i = 1.10 ± 0.25 mM. Measurements of intrinsic and ANS-binding fluorescences showed that rutin showed a relatively strong binding affinity for tyrosinase and one possible binding site that could be a copper was detected accompanying with a hydrophobic exposure of tyrosinase. Cell viability testing with rutin in HaCaT keratinocytes showed that no toxic effects were produced. Taken together, rutin has the potential to be a potent antipigment agent. The strategy of predicting tyrosinase inhibition based on hydroxyl group number and computational simulation may prove useful for the screening of potential tyrosinase inhibitors.     

Virtual screening to identify novel potential inhibitors for Glutamine synthetase of Mycobacterium tuberculosis

Abstract

Glutamine synthetase (GS) of Mycobacterium tuberculosis (Mtb) is an essential enzyme which is involved in nitrogen metabolism and cell wall synthesis. It is involved in the inhibition of phagosome-lysosome fusion by preventing acidification. Targeting GS can be helpful to control the infection of Mtb. In order to identify potential inhibitors, we screened chemical libraries (56,400 compounds of ChEMBL anti-mycobacterial, 1596 FDA approved drugs, 419 Natural product and 916 phytochemical) against this target using the virtual screening approach. Screening by molecular docking identified ten top-ranked compounds as GSMtb inhibitors and they were compared with known inhibitors (as control). Since GS enzyme (GSHs) is also present in human. We have compared the protein sequence of GS from Mtb and human using the P-BLAST in NCBI. We found ～27% identity in between these two sequences, so we also compared the binding affinity of inhibitor between Mtb and human. Finally, we identified top two compounds namely CHEMBL387509, CHEMBL226198 from ChEMBL anti-mycobacterial dataset, and Eriocitrin and Malvidin from phytochemical dataset which showed lees binding affinity towards GSHs whereas Pamidronate, and Phentermine from FDA approved drugs and (-)-Quinic Acid, Hexopyranuronic acid, Quebrachit, and Castanospermine from natural product showed protein-ligand interaction with Mtb protein while no interaction with GSHs. The top two docked complexes were subjected to molecular dynamic simulation to understand the stability of the molecule. Further, we calculated the binding free energy of the docked complex and analyzed hydrogen bond, salt bridge, pie stacking, and hydrophobic interaction in the docking region. These ligands exhibited very good binding affinity GSMtb enzymes. Therefore, these ligands are novel and drug-likeness compounds, and they may be potential inhibitors of M tuberculosis.

Communicated by Ramaswamy H. Sarma     

In silico approach towards identification of potential inhibitors of Helicobacter pylori DapE

Helicobacter pylori is a gastric mucosal pathogen and is associated with diseases like peptic ulcer and gastric cancer. To combat H. pylori infection, there is an urgent need for new class of antibiotics due to the emergence of drug-resistant strains. Enzymes involved in bacterial lysine biosynthetic pathways may be potential targets for antibacterial drug development, since lysine is an essential component of the bacterial peptidoglycan cell wall. No pathway exists for lysine biosynthesis in humans; hence, the inhibitors targeting bacterial enzymes may have selective toxicity. dapE-encoded N-succinyl-L,L-diaminopimelic acid desuccinylase (DapE) is a critical enzyme of this pathway and deletion of DapE gene is lethal to H. pylori, since the organism has no alternative pathway for lysine biosynthesis. In this study, we reported a 3D model structure of H. pylorie DapE, which consisted of a catalytic domain and a dimerization domain generated by MODELLER software. We also confirmed the stability of the modeled structure through 10 ns molecular dynamics simulation using GROMACS software. Next, to identify potential small molecule inhibitors of DapE, drug-like small molecule-screening library was generated. This was performed by Tanimoto-based similarity searching in the PubChem Database with DapE substrate L,L-SDAP as a query molecule, followed by fragment-based docking approach using GLIDE XP. This approach identified two potential substrate-competitive small molecule inhibitors of DapE. These new molecules may provide a starting point to search for novel therapeutics.     

Preventing the return of smallpox: molecular modeling studies on thymidylate kinase from Variola virus

Smallpox was one of the most devastating diseases in the human history and still represents a serious menace today due to its potential use by bioterrorists. Considering this threat and the non-existence of effective chemotherapy, we propose the enzyme thymidylate kinase from Variola virus (VarTMPK) as a potential target to the drug design against smallpox. We first built a homology model for VarTMPK and performed molecular docking studies on it in order to investigate the interactions with inhibitors of Vaccinia virus TMPK (VacTMPK). Subsequently, molecular dynamics (MD) simulations of these compounds inside VarTMPK and human TMPK (HssTMPK) were carried out in order to select the most promising and selective compounds as leads for the design of potential VarTMPK inhibitors. Results of the docking and MD simulations corroborated to each other, suggesting selectivity towards VarTMPK and, also, a good correlation with the experimental data.     

A computational approach to explore and identify potential herbal inhibitors for the p21-activated kinase 1 (PAK1)

Abstract

The oncogenic kinase PAK1 (p21-activated kinase 1) is involved in developing many diseases including cancers, neurofibromatosis, Alzheimer's disease, diabetes (type 2), and hypertension. Thus, it is thought to be a prominent therapeutic target, and its selective inhibitors have a huge market potential. Recently, herbal PAK1 inhibitors have gained immense interest over synthetic ones mainly due to their non-toxic effects. Till date, many herbal compounds have been suggested to inhibit PAK1, but their information on selectivity, bioavailability, ADMET (absorption, distribution, metabolism, excretion, and toxicity) properties, and molecular interactions with PAK1 has not been explored. Hence, this study was designed with computational approaches to explore and identify the best herbal PAK1-blockers showing good ADMET properties, druggable features and binding affinity with PAK1. Herbal inhibitors reported here were initially filtered with Lipinski’s rule of five (RO5). Then, molecular docking between these inhibitors and PAK1 catalytic sites was performed using AutoDock Vina and GOLD suite to determine the binding affinity and interactions. Finally, 200?ns molecular dynamics (MD) simulations on three top-ranked inhibitors including cucurbitacin I (C-I), nymphaeol A (NA), and staurosporine (SPN) were carried out. The binding free energies and interactions revealed that NA can strongly bind with the PAK1 catalytic cleft. PASS prediction and ADMET profiling supported that NA is appeared to be a more selective and safer inhibitor than C-I and SPN. These results conform to the previous experimental evidences, and therefore, NA from Okinawa propolis could be a promising inhibitor for treating PAK1-dependent illnesses.

Communicated by Ramaswamy H. Sarma     

rDock: A Fast,Versatile and Open Source Program for Docking Ligands to Proteins and Nucleic Acids

Identification of chemical compounds with specific biological activities is an important step in both chemical biology and drug discovery. When the structure of the intended target is available, one approach is to use molecular docking programs to assess the chemical complementarity of small molecules with the target; such calculations provide a qualitative measure of affinity that can be used in virtual screening (VS) to rank order a list of compounds according to their potential to be active. rDock is a molecular docking program developed at Vernalis for high-throughput VS (HTVS) applications. Evolved from RiboDock, the program can be used against proteins and nucleic acids, is designed to be computationally very efficient and allows the user to incorporate additional constraints and information as a bias to guide docking. This article provides an overview of the program structure and features and compares rDock to two reference programs, AutoDock Vina (open source) and Schrödinger''s Glide (commercial). In terms of computational speed for VS, rDock is faster than Vina and comparable to Glide. For binding mode prediction, rDock and Vina are superior to Glide. The VS performance of rDock is significantly better than Vina, but inferior to Glide for most systems unless pharmacophore constraints are used; in that case rDock and Glide are of equal performance. The program is released under the Lesser General Public License and is freely available for download, together with the manuals, example files and the complete test sets, at http://rdock.sourceforge.net/

This is a PLOS Computational Biology Software Article.

     

Design,synthesis and biological evaluation of novel N-phosphorylated and O-phosphorylated tacrine derivatives as potential drugs against Alzheimer’s disease

In this work, we designed, synthesised and biologically investigated a novel series of 14 N- and O-phosphorylated tacrine derivatives as potential anti-Alzheimer’s disease agents. In the reaction of 9-chlorotacrine and corresponding diamines/aminoalkylalcohol we obtained diamino and aminoalkylhydroxy tacrine derivatives. Next, the compounds were acid to give final products 6–13 and 16–21 that were characterised by ¹H, ¹³C, ³¹P NMR and MS. The results of the docking studies revealed that the designed phosphorus hybrids, in theory can bind to AChE and BChE. All compounds exhibited significantly lower AutoDock Vina scores compared to tacrine. The inhibitory potency evaluation was performed using the Ellman’s method. The most inhibitory activity against AChE exhibited compound 8 with an IC₅₀ value of 6.11 nM and against BChE 13 with an IC₅₀ value of 1.97 nM and they were 6- and 12-fold potent than tacrine. Compound 19 showed the lack of hepatocytotoxicity in MTT assay.     

Structure‐guided screening of chemical database to identify NS3‐NS2B inhibitors for effective therapeutic application in dengue infection

Dengue infection is the most common arthropod‐borne disease caused by dengue viruses, predominantly affecting millions of human beings annually. To find out promising chemical entities for therapeutic application in Dengue, in the current research, a multi‐step virtual screening effort was conceived to screen out the entire “screening library” of the Asinex database. Initially, through “Lipinski rule of five” filtration criterion almost 0.6 million compounds were collected and docked with NS3‐NS2B protein. Thereby, the chemical space was reduced to about 3500 compounds through the analysis of binding affinity obtained from molecular docking study in AutoDock Vina. Further, the “Virtual Screening Workflow” (VSW) utility of Schrödinger suite was used, which follows a stepwise multiple docking programs such as ‐ high‐throughput virtual screening (HTVS), standard precision (SP), and extra precision (XP) docking, and in postprocessing analysis the MM‐GBSA based free binding energy calculation. Finally, five potent molecules were proposed as potential inhibitors for the dengue NS3‐NS2B protein based on the investigation of molecular interactions map and protein‐ligand fingerprint analyses. Different pharmacokinetics and drug‐likeness parameters were also checked, which favour the potentiality of selected molecules for being drug‐like candidates. The molecular dynamics (MD) simulation analyses of protein‐ligand complexes were explained that NS3‐NS2B bound with proposed molecules quite stable in dynamic states as observed from the root means square deviation (RMSD) and root means square fluctuation (RMSF) parameters. The binding free energy was calculated using MM‐GBSA method from the MD simulation trajectories revealed that all proposed molecules possess such a strong binding affinity towards the dengue NS3‐NS2B protein. Therefore, proposed molecules may be potential chemical components for effective inhibition of dengue NS3‐NS2B protein subjected to experimental validation.     

Inhibition of Histo-blood Group Antigen Binding as a Novel Strategy to Block Norovirus Infections

Noroviruses (NoVs) are the most important viral pathogens that cause epidemic acute gastroenteritis. NoVs recognize human histo-blood group antigens (HBGAs) as receptors or attachment factors. The elucidation of crystal structures of the HBGA-binding interfaces of a number of human NoVs representing different HBGA binding patterns opens a new strategy for the development of antiviral compounds against NoVs through rational drug design and computer-aided virtual screening methods. In this study, docking simulations and virtual screening were used to identify hit compounds targeting the A and B antigens binding sites on the surface of the capsid P protein of a GII.4 NoV (VA387). Following validation by re-docking of the A and B ligands, these structural models and AutoDock suite of programs were used to screen a large drug-like compound library (derived from ZINC library) for inhibitors blocking GII.4 binding to HBGAs. After screening >2 million compounds using multistage protocol, 160 hit compounds with best predicted binding affinities and representing a number of distinct chemical classes have been selected for subsequent experimental validation. Twenty of the 160 compounds were found to be able to block the VA387 P dimers binding to the A and/or B HBGAs at an IC₅₀<40.0 µM, with top 5 compounds blocking the HBGA binding at an IC₅₀<10.0 µM in both oligosaccharide- and saliva-based blocking assays. Interestingly, 4 of the top-5 compounds shared the basic structure of cyclopenta [a] dimethyl phenanthren, indicating a promising structural template for further improvement by rational design.     

Identification of ALK5 inhibitor via structure-based virtual screening and ADMET prediction

Abstract

TGF-β plays a critical role in the initiation and progression of fibrosis in various organ systems such as kidney, heart, lung and liver. TGF-β and its receptors (ALK5 and TβR II) are able to control the cellular growth and promote several biological responses. To date, many pharmaceutical companies have employed virtual screening to identify potent inhibitors against ALK5. Nevertheless, none of these studies had involved the in silico ADMET evaluation and Raccoon filtering. In our experiment, all 57423 molecules were downloaded from TCM database and were filtered and converted to PDBQT formats by Raccoon software. Then 24?189 structures were run through AutoDock Vina in PyRx 0.8, 164 molecules were selected and further evaluated by ADMET Predictor 6.5, and 56 structures were selected and docked by Glide 6.2. Finally, the top 10 hits were identified as promising oral ALK5 inhibitors according to their Glide scores. The Glide scores of the best two compounds, 40686 and 33534, were ?10.75 and ?10.30?kcal/mol, respectively. This research provides a set of combined and detailed virtual screening protocol and is helpful for explaining the mechanism of receptor–ligand interactions.     

Evaluation of potential flavonoid inhibitors of glyoxalase-I based on virtual screening and in vitro studies

Glyoxalase-I (GLO-I) is a component of the ubiquitous detoxification system involved in the conversion of methylglyoxal (MG) to d-lactate in the glycolytic pathway. MG toxicity arises from its ability to form advanced glycation end products. GLO-I has been reported to be frequently overexpressed in various types of cancer cells. In this study, we performed structure-based virtual screening of focused flavonoids commercial library to identify potential and specific inhibitors of GLO-I. The compounds were ranked based on Glide extra precision docking score and five hits (curcumin, quercetin, morin, naringin and silibinin) were selected on the basis of their interaction with active site amino acid residues of GLO-I. Mixed mode QM/MM calculation was performed on the top-scoring hit to ascertain the role of zinc ion in ligand binding. In addition, the identified hits were subjected to MM/GBSA binding energy prediction, ADME prediction and similarity studies. The hits were tested in vitro for cell viability, and GLO-I inhibition. Naringin (ST072162) was found to be most potent inhibitor of GLO-I among the identified hits with highest glide XP dock score of ?14.906. These findings suggest that naringin could be a new scaffold for designing inhibitors against GLO-I with potential application as anticancer agents.     

Identification of 3′,4′,5′-trimethoxychalcone analogues as potent inhibitors of Helicobacter pylori-induced inflammation in human gastric epithelial cells

Efforts to identify potent small molecule inhibitors of Helicobacter pylori led to the evaluation of 23 3′,4′,5′-trimethoxychalcone analogues. Some of the compounds displayed potent antibacterial activity against H. pylori. Three most active and selective compounds 1, 7, and 13 also showed the bactericide activity against the reference as well as multidrug-resistant strains of H. pylori. Additionally, the aforementioned three compounds potentially inhibited the H. pylori adhesion and invasion to human gastric epithelial (AGS) cells. Furthermore, these selective compounds inhibited the H. pylori-induced gastric inflammation by reduced inflammatory mediator’s nuclear factor kappa B activation, and the secretion of interleukin-8.     

An integrated study of tyrosinase inhibition by rutin: progress using a computational simulation

Tyrosinase inhibition studies have recently gained the attention of researchers due to their potential application values. We simulated docking (binding energies for AutoDock Vina: -9.1 kcal/mol) and performed a molecular dynamics simulation to verify docking results between tyrosinase and rutin. The docking results suggest that rutin mostly interacts with histidine residues located in the active site. A 10 ns molecular dynamics simulation showed that one copper ion at the tyrosinase active site was responsible for the interaction with rutin. Kinetic analyses showed that rutin-mediated inactivation followed a first-order reaction and mono- and biphasic rate constants occurred with rutin. The inhibition was a typical competitive type with K(i) = 1.10±0.25 mM. Measurements of intrinsic and ANS-binding fluorescences showed that rutin showed a relatively strong binding affinity for tyrosinase and one possible binding site that could be a copper was detected accompanying with a hydrophobic exposure of tyrosinase. Cell viability testing with rutin in HaCaT keratinocytes showed that no toxic effects were produced. Taken together, rutin has the potential to be a potent anti-pigment agent. The strategy of predicting tyrosinase inhibition based on hydroxyl group number and computational simulation may prove useful for the screening of potential tyrosinase inhibitors.