Identification of potential HIV-1 integrase strand transfer inhibitors: In silico virtual screening and QM/MM docking studies

HIV-1 integrase (IN) is a retroviral enzyme that catalyses integration of the reverse-transcribed viral DNA into the host genome, which is necessary for efficient viral replication. In this study, we have performed an in silico virtual screening for the identification of potential HIV-1 IN strand transfer (ST) inhibitors. Pharmacophore modelling and atom-based 3D-QSAR studies were carried out for a series of compounds belonging to 3-Hydroxypyrimidine-2,4-diones. Based on the ligand-based pharmacophore model, we obtained a five-point pharmacophore with two hydrogen bond acceptors (A), one hydrogen bond donor (D), one hydrophobic group (H) and one aromatic ring (R) as pharmacophoric features. The pharmacophore hypothesis AADHR was used as a 3D query in a sequential virtual screening study to filter small molecule databases Maybridge, ChemBridge and Asinex. Hits matching with pharmacophore hypothesis AADHR were retrieved and passed progressively through Lipinski’s rule of five filtering, molecular docking and hierarchical clustering. The five compounds with best hits with novel and diverse chemotypes were subjected to QM/MM docking, which showed improved docking accuracy. We further performed molecular dynamics simulation and found three compounds that form stable interactions with key residues. These compounds could be used as a leads for further drug development and rational design of HIV-1 IN inhibitors.     

Identification of potential influenza virus endonuclease inhibitors through virtual screening based on the 3D-QSAR model

Influenza endonucleases have appeared as an attractive target of antiviral therapy for influenza infection. With the purpose of designing a novel antiviral agent with enhanced biological activities against influenza endonuclease, a three-dimensional quantitative structure-activity relationships (3D-QSAR) model was generated based on 34 influenza endonuclease inhibitors. The comparative molecular similarity index analysis (CoMSIA) with a steric, electrostatic and hydrophobic (SEH) model showed the best correlative and predictive capability (q ² = 0.763, r ² = 0.969 and F = 174.785), which provided a pharmacophore composed of the electronegative moiety as well as the bulky hydrophobic group. The CoMSIA model was used as a pharmacophore query in the UNITY search of the ChemDiv compound library to give virtual active compounds. The 3D-QSAR model was then used to predict the activity of the selected compounds, which identified three compounds as the most likely inhibitor candidates.     

Discovery of novel indoleamine 2,3-dioxygenase 1 (IDO1) inhibitors by virtual screening

In this study, a combination of virtual screening methods were utilized to identify novel potential indoleamine 2,3-dioxygenase 1 (IDO1) inhibitors. A series of IDO1 potential inhibitors were identified by a combination of following steps: Lipinski's Rule of Five, Veber rules filter, molecular docking, HipHop pharmacophores, 3D-Quantitative structure activity relationship (3D-QSAR) studies and Pan-assay Interference Compounds (PAINS) filter. Three known categories of IDO1 inhibitors were used to constructed pharmacophores and 3D-QSAR models. Four point pharmacophores (RHDA) of IDO1 inhibitors were generated from the training set. The 3D-QSAR models were obtained using partial least squares (PLS) analyze based on the docking conformation alignment from the training set. The leave-one-out correlation (q²) and non-cross-validated correlation coefficient (r²_pred) of the best CoMFA model were 0.601 and 0.546, and the ones from the best CoMSIA model were 0.506 and 0.541, respectively. Six hits from Specs database were identified and analyzed to confirm their binding modes and key interactions to the amino acid residues in the protein. This work may provide novel backbones for new generation of inhibitors of IDO1.     

Structure-based virtual screening of natural products as potential stearoyl-coenzyme a desaturase 1 (SCD1) inhibitors

Stearyl coenzyme A desaturase enzyme 1 (SCD1) is a key enzyme that catalyzes the conversion of saturated fatty acids (SFA) into monounsaturated fatty acids (MUFA) and plays a vital role in lipid metabolism of tumor cells. SCD1 is overexpressed in a variety of malignant tumors, and its related inhibitors showed significant anti-tumor activity in vitro and in vivo experiments, which is a new target for tumor therapy. The focus of this study is to identify novel SCD1 inhibitors from natural products through computer simulations. First, 176,602 compounds from natural product databases were virtually screened. By molecular dynamics (MD) simulations, the ligand-protein interactions of 5 compounds with high docking manifestation were analyzed accurately. Then, MM-GBSA and MM-PBMA methods were used to verify the results. Finally, ADMET prediction was performed for the 5 compounds. As a result, two natural products with potential inhibition towards SCD1 were identified, which had the excellent docking manifestation, binding mode within SCD1 pocket and stability during molecular dynamics simulation. This study provides a meaningful model for the development and optimization of new inhibitors and anti-tumor drugs targeting SCD1.     

Identification of potential Gly/NMDA receptor antagonists by cheminformatics approach: a combination of pharmacophore modelling,virtual screening and molecular docking studies

The Gly/NMDA receptor has become known as potential target for the management of neurodegenerative diseases. Discovery of Gly/NMDA antagonists has thus attracted much attention in recent years. In the present research, a cheminformatics approach has been used to determine structural requirements for Gly/NMDA antagonism and to identify potential antagonists. Here, 37 quinoxaline derivatives were selected to develop a significant pharmacophore model with good certainty. The selected model was validated by leave-one-out cross-validation, an external test set, decoy set and Y-randomization test. Applicability domain was verified by the standardization approach. The validated 3D-QSAR model was used to screen virtual hits from the ZINC database by pharmacophore mapping. Molecular docking was used for assessment of receptor–ligand binding modes and binding affinities. The GlideScore and molecular interactions with critical amino acids were considered as crucial features to identify final hits. Furthermore, hits were analysed for in silico pharmacokinetic parameters and Lipinski’s rule of five, demonstrating their potential as drug-like candidates. The PubChem and SciFinder search tools were used to authenticate the novelty of leads retrieved. Finally, five different leads have been suggested as putative novel candidates for the exploration of potent Gly/NMDA receptor antagonists.     

Identification of potential inhibitors of human methionine aminopeptidase (type II) for cancer therapy: Structure-based virtual screening,ADMET prediction and molecular dynamics studies

Methionine Aminopeptidases MetAPs are divalent-cofactor dependent enzymes that are responsible for the cleavage of the initiator Methionine from the nascent polypeptides. MetAPs are classified into two isoforms: namely, MetAP1 and MetAP2. Several studies have revealed that MetAP2 is upregulated in various cancers, and its inhibition has shown to suppress abnormal or excessive blood vessel formation and tumor growth in model organisms. Clinical studies show that the natural product fumagillin, and its analogs are potential inhibitors of MetAP2. However, due to their poor pharmacokinetic properties and neurotoxicities in clinical studies, their further developments have received a great setback. Here, we apply structure-based virtual screening and molecular dynamics methods to identify a new class of potential inhibitors for MetAP2. We screened Otava’s Chemical Library, which consists of about 3 200 000 tangible-chemical compounds, and meticulously selected the top 10 of these compounds based on their inhibitory potentials against MetAP2. The top hit compounds subjected to ADMET predictor using 3 independent ADMET prediction programs, were found to be drug-like. To examine the stability of ligand binding mode, and efficacy, the unbound form of MetAP2, its complexes with fumagillin, spiroepoxytriazole, and the best promising compounds compound-3369841 and compound-3368818 were submitted to 100 ns molecular dynamics simulation. Like fumagillin, spiroepoxytriazole, and both compound-3369841 and compound-3368818 showed stable binding mode over time during the simulations. Taken together, these uninherited-fumagillin compounds may serve as new class of inhibitors or provide scaffolds for further optimization towards the design of more potent MetAP2 inhibitors -development of such inhibitors would be essential strategy against various cancer types.     

In silico design novel (5-imidazol-2-yl-4-phenylpyrimidin-2-yl)[2-(2-pyridylamino)ethyl]amine derivatives as inhibitors for glycogen synthase kinase 3 based on 3D-QSAR,molecular docking and molecular dynamics simulation

Glycogen Synthase Kinase 3 (GSK-3) is a member of cellular kinase with various functions, such as glucose regulation, cellular differentiation, neuronal function and cell apoptosis. It has been proved as an important therapeutic target in type 2 diabetes mellitus and Alzheimer's disease. To better understand their structure–activity relationships and mechanism of action, an integrated computational study, including three dimensional quantitative structure-activity relationship (3D-QSAR), molecular docking, and molecular dynamics (MD), was performed on 79 (5-Imidazol-2-yl-4-phenylpyrimidin-2-yl)[2-(2-pyridylamino)ethyl]amine GSK-3 inhibitors. In this paper, we constructed 3D-QSAR using comparative molecular field analysis (CoMFA) and comparative molecular similarity index analysis (CoMSIA) method. The results showed that the CoMFA model (q ² = 0.743, r² = 0.980) and the CoMSIA model (q² = 0.813, r² = 0.976) had stable and reliable predictive ability. The electrostatic and H-bond donor fields play important roles in the models. The contour maps of the model visually showed the relationship between the activity of compounds and their three-dimensional structure. Molecular docking was used to identify the key amino acid residues at the active site of GSK-3 and explore its binding mode with ligands. Based on 3D-QSAR models, contour maps and the binding feature between GSK-3 and inhibitor, we designed 10 novel compounds with good potential activity and ADME/T profile. Molecular dynamics simulation results validated that Ile62, Val70 and Lys85 located in the active site play a key role for GSK-3 complexed with inhibitors. These results might provide important information for designing GSK-3 inhibitors with high activity.     

High-throughput virtual screening and microsecond MD simulations to identify potential sugar mimic of the solute-binding protein BlAXBP of the ABC transporter from Bifidobacterium animalis subsp. Lactis

Xylotetraose is a prebiotic oligosaccharide can be utilized by the ABC transporter of the gut microbiota Bifidobacteria. BlAXBP is the solute binding protein of the ABC transporter, and its complex with xylotetraose has been solved by X-ray crystallography. Here, we have identified novel sugar mimic of BlAXBP by applying a high-throughput virtual screening of ZINC database containing a huge library with ∼22 M compounds. To begin with, we identified 18,571 ligands by a ligand-based virtual screening. Further, a total of 3968 compounds were selected for molecular docking due to their Tanimoto coefficient’s value were larger than a cutoff of 0.08. The molecular mechanics-generalized born surface area was used to evaluate the binding free energies, and the top 10 ligands with free energies below an energy threshold of -35.22 kcal/mol were selected. ZINC13783511 formed the most stable complex with BlAXBP and its recognition mechanism were further explored by microsecond MD simulations in explicit solvent. Free energy landscapes were used to evaluate conformational changes of BlAXBP in its ligand free and binding states. Collectively, this work identified potential novel sugar mimics to BlAXBP, providing novel atomic-level understanding of the binding mechanism.     

De novo design based identification of potential HIV-1 integrase inhibitors: A pharmacoinformatics study

In the present study, pharmacoinformatics paradigms include receptor-based de novo design, virtual screening through molecular docking and molecular dynamics (MD) simulation are implemented to identify novel and promising HIV-1 integrase inhibitors. The de novodrug/ligand/molecule design is a powerful and effective approach to design a large number of novel and structurally diverse compounds with the required pharmacological profiles. A crystal structure of HIV-1 integrase bound with standard inhibitor BI-224436 is used and a set of 80,000 compounds through the de novo approach in LigBuilder is designed. Initially, a number of criteria including molecular docking, in-silico toxicity and pharmacokinetics profile assessments are implied to reduce the chemical space. Finally, four de novo designed molecules are proposed as potential HIV-1 integrase inhibitors based on comparative analyses. Notably, strong binding interactions have been identified between a few newly identified catalytic amino acid residues and proposed HIV-1 integrase inhibitors. For evaluation of the dynamic stability of the protein-ligand complexes, a number of parameters are explored from the 100 ns MD simulation study. The MD simulation study suggested that proposed molecules efficiently retained their molecular interaction and structural integrity inside the HIV-1 integrase. The binding free energy is calculated through the Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) approach for all complexes and it also explains their thermodynamic stability. Hence, proposed molecules through de novo design might be critical to inhibiting the HIV-1 integrase.     

Multiple QSAR and molecular modelling for identification of potent human adenovirus inhibitors

This study has investigated docking-based 2D- and 3D-quantitative structure-activity relationships (QSARs) for a range of 53 hydroxybenzamide analogues as anti- Human adenoviruses (HAdVs). The best 3D-QSAR (Schrodinger, LLC, NY, 2020) and 2D-QSAR models were obtained for the training set and were found to be statistically significant, with cross-validated coefficients (q²) of 0.6775 and 0.7875, and coefficients of determination (r²) of 0.8106 and 0.8122, respectively. Our in-silico docking and virtual screening studies revealed significant higher binding affinity of dataset molecule 34 (-141.444 ​kcal/mol) and hit ZINC01088642 (-114.357 ​kcal/mol) with 4PIE protein than the standard drugs. In in-silico ADME/toxicity studies, molecule 34 and proposed hit ZINC01088642 were found safe with good intestinal absorption, aqueous solubility, medium blood–brain barrier (BBB), no eye corrosion, no skin irritancy, and non-mutagenic profiles. Molecular dynamics analysis showed good stability of complex, hit ZINC01088642 with protein, 4PIE over the simulation period of 20 ns. We believe that further experimental, as well as in-vitro investigation, will shed more lights on the identification of ZINC01088642 as a potential human adenovirus agent.     

Pharmacophore modelling,virtual screening and molecular docking studies on PLD1 inhibitors

Lipid metabolism plays a significant role in influenza virus replication and subsequent infection. The regulatory mechanism governing lipid metabolism and viral replication is not properly understood to date, but both Phospholipase D (PLD1 and PLD2) activities are stimulated in viral infection. In vitro studies indicate that chemical inhibition of PLD1 delays viral entry and reduction of viral loads. The current study reports a three-dimensional pharmacophore model based on 35 known PLD1 inhibitors. A sub-set of 25 compounds was selected as the training set and the remaining 10 compounds were kept in the test set. One hundred and twelve pharmacophore models were generated; a six-featured pharmacophore model (AADDHR.57) with survival score (2.69) produced a statistically significant three-dimensional quantitative structure–activity relationship model with r² = 0.97 (internal training set), r² = 0.71 (internal test set) and Q² = 0.64. The predictive power of the pharmacophore model was validated with an external test set (r² = 0.73) and a systematic virtual screening work-flow was employed showing an enrichment factor of 23.68 at the top 2% of the dataset (active and decoys). Finally, the model was used for screening of the filtered PubChem database to fetch molecules which can be proposed as potential PLD1 inhibitors for blocking influenza infection.     

Sub-structure-based screening and molecular docking studies of potential enteroviruses inhibitors

Rhinoviruses (RV), especially Human rhinovirus (HRVs) have been accepted as the most common cause for upper respiratory tract infections (URTIs). Pleconaril, a broad spectrum anti-rhinoviral compound, has been used as a drug of choice for URTIs for over a decade. Unfortunately, for various complications associated with this drug, it was rejected, and a replacement is highly desirable. In silico screening and prediction methods such as sub-structure search and molecular docking have been widely used to identify alternative compounds. In our study, we have utilised sub-structure search to narrow down our quest in finding relevant chemical compounds. Molecular docking studies were then used to study their binding interaction at the molecular level. Interestingly, we have identified 3 residues that is worth further investigation in upcoming molecular dynamics simulation systems of their contribution in stable interaction.     

Identification of potential Tpx inhibitors against pathogen-host interactions

Yersinia organisms cause many infectious diseases by invading human cells and delivering their virulence factors via the type three secretion system (T3SS). One alternative strategy in the fight against these pathogenic organisms is to interfere with their T3SS. Previous studies demonstrated that thiol peroxidase, Tpx is functional in the assembly of T3SS and its inhibition by salicylidene acylhydrazides prevents the secretion of pathogenic effectors. In this study, the aim was to identify potential inhibitors of Tpx using an integrated approach starting with high throughput virtual screening and ending with molecular dynamics simulations of selected ligands. Virtual screening of ZINC database of 500,000 compounds via ligand-based and structure-based pharmacophore models retrieved 10,000 hits. The structure-based pharmacophore model was validated using high-throughput virtual screening (HTVS). After multistep docking (SP and XP), common scaffolds were used to find common substructures and the ligand binding poses were optimized using induced fit docking. The stability of the protein–ligand complex was examined with molecular dynamics simulations and the binding free energy of the complex was calculated. As a final outcome eight compounds with different chemotypes were proposed as potential inhibitors for Tpx. The eight ligands identified by a detailed virtual screening protocol can serve as leads in future drug design efforts against the destructive actions of pathogenic bacteria.     

Computational insights into factor affecting the potency of diaryl sulfone analogs as Escherichia coli dihydropteroate synthase inhibitors

Dihydropteroate synthase (DHPS) is an alluring target for designing novel drug candidates to prevent infections caused by pathogenic Escherichia coli strains. Diaryl Sulfone (SO) compounds are found to inhibit DHPS competitively with respect to the substrate pABA (p-aminobenzoate). The extra aromatic ring of diaryl sulfone compounds found to stabilize them in highly flexible pABA binding loops. In this present study, a statistically significant 3D-QSAR model was developed using a data set of diaryl sulfone compounds. The favourable and unfavourable contributions of substitutions in sulfone compounds were illustrated by contour plot obtained from the developed 3D-QSAR model. Molecular docking calculations were performed to investigate the putative binding mode of diaryl sulfone compounds at the catalytic pocket. DFT calculations were carried out using SCF approach, B3LYP- 6-31 G (d) basis set to compute the HOMO, LUMO energies and their respective location at pABA binding pocket. Further, the developed model was validated by FEP (Free Energy Perturbation) calculations. The calculated relative free energy of binding between the highly potent and less potent sulfone compound was found to be −3.78 kcal/ mol which is comparable to the experimental value of −5.85 kcal/mol. A 10 ns molecular dynamics simulation of inhibitor and DHPS confirmed its stability at pABA catalytic site. Outcomes of the present work provide deeper insight in designing novel drug candidates for pathogenic Escherichia coli strains.     

Structure-based virtual screening of influenza virus RNA polymerase inhibitors from natural compounds: Molecular dynamics simulation and MM-GBSA calculation

The resistances of matrix protein 2 (M2) protein inhibitors and neuraminidase inhibitors for influenza virus have attracted much attention and there is an urgent need for new drug. The antiviral drugs that selectively act on RNA polymerase are less prone to resistance and possess fewer side effects on the patient. Therefore, there is increased interest in screening compounds that can inhibit influenza virus RNA polymerase. Three natural compounds were found by using molecular docking-based virtual screening, which could bind tightly within the polymerase acidic protein-polymerase basic protein 1 (PA-PB1) subunit of influenza virus polymerase. Firstly, their drug likeness properties were evaluated, which showed that the hepatotoxicity values of all the three compounds indicating they had less or no hepatotoxicity, and did not have the plasma protein biding (PPB) ability, the three compounds needed to be modified in some aspects, like bulky molecular size. The stability of the complexes of PA-hits was validated through molecular dynamics (MD) simulation, revealing compound 2 could form more stable complex with PA subunit. The torsional conformations of each rotatable bond of the ligands in PA subunit were also monitored, to investigate variation in the ligand properties during the simulation, compound 3 had fewer rotatable bonds, indicating that the molecule had stronger rigidity. The bar charts of protein–ligand contacts and contacts over the course of trajectory showed that four key residues (Glu623, Lys643, Asn703 and Trp706) of PA subunit that participated in hydrogen-bond, water bridge and hydrophobic interactions with the hit compounds. Finally, the binding free energy and contributed energies were calculated by using MM-GBSA method. Out of the three compounds, compound 1 showed the lowest total binding free energy. Among all the interactions, the contribution of the covalent binding and the van der Waals energy were more than other items, compound 1 formed more stable hydrogen bonds with the residues of PA subunit binding pocket. This study smoothed the path for the development of novel lead compounds with improved binding properties, high drug likeness, and low toxicity to humans for the treatment of influenza, which provided a good basis for further research on novel and effective influenza virus PA-PB1 interaction inhibitors.     

Evaluation of library ranking efficacy in virtual screening

We present the results of a comprehensive study in which we explored how the docking procedure affects the performance of a virtual screening approach. We used four docking engines and applied 10 scoring functions to the top-ranked docking solutions of seeded databases against six target proteins. The scores of the experimental poses were placed within the total set to assess whether the scoring function required an accurate pose to provide the appropriate rank for the seeded compounds. This method allows a direct comparison of library ranking efficacy. Our results indicate that the LigandFit/Ligscore1 and LigandFit/GOLD docking/scoring combinations, and to a lesser degree FlexX/FlexX, Glide/Ligscore1, DOCK/PMF (Tripos implementation), LigandFit1/Ligscore2 and LigandFit/PMF (Tripos implementation) were able to retrieve the highest number of actives at a 10% fraction of the database when all targets were looked upon collectively. We also show that the scoring functions rank the observed binding modes higher than the inaccurate poses provided that the experimental poses are available. This finding stresses the discriminatory ability of the scoring algorithms, when better poses are available, and suggests that the number of false positives can be lowered with conformers closer to bioactive ones.     

Discovery of potential Zika virus RNA polymerase inhibitors by docking-based virtual screening

Zika virus (ZIKV) infection has been associated with Guillain-Barre syndrome in adults and microcephaly in infants. The existence of insufficient structural data in most of the protein databases hinders the synthesis of anti-ZIKV pharmaceutics. In this work, we attempted to model the catalytic domain of the ZIKV RNA polymerase (RdRpC) along with a detailed assessment of conserved aspartates in ZIKV RdRpC palm domain as potential drug targets. The conserved and catalytically active aspartate residues present in the predicted RdRpC protein were virtually screened against a ZINC database for inhibitors, and the selected potential drug candidates were further filtered based on their ADMET profiles. One of the pharmacokinetically active compounds (Ligand 6) showed a remarkable docking profile against the strictly conserved aspartate residues of the RdRpC active site. We hypothesize that the Ligand 6 may form a potential drug candidate for RdRpC inhibition in the clinical treatment of ZIKV infection.     

基于药效团的BRD4靶向小分子抑制剂虚拟筛选

BRD4靶点和多种肿瘤密切相关,是具有良好成药性的热门靶点.本文选取活性较好且结构差异较大的BRD4小分子抑制剂作为训练集分子,基于配体小分子共同特征(HipHop)方法使用Discovery Studio 3.0分子模拟软件构建了药效团.药效团通过测试集验证、ROC曲线验证(SE(sensitivity)=0.937...     

Virtual screening of B-Raf kinase inhibitors: A combination of pharmacophore modelling,molecular docking, 3D-QSAR model and binding free energy calculation studies

B-Raf kinase has been identified as an important target in recent cancer treatment. In order to discover structurally diverse and novel B-Raf inhibitors (BRIs), a virtual screening of BRIs against ZINC database was performed by using a combination of pharmacophore modelling, molecular docking, 3D-QSAR model and binding free energy (ΔG_bind) calculation studies in this work. After the virtual screening, six promising hit compounds were obtained, which were then tested for inhibitory activities of A375 cell lines. In the result, five hit compounds show good biological activities (IC₅₀ < 50 μM). The present method of virtual screening can be applied to find structurally diverse inhibitors, and the obtained five structurally diverse compounds are expected to develop novel BRIs.     

Virtual screening and rational drug design method using structure generation system based on 3D-QSAR and docking

An efficient virtual and rational drug design method is presented. It combines virtual bioactive compound generation with 3D-QSAR model and docking. Using this method, it is possible to generate a lot of highly diverse molecules and find virtual active lead compounds. The method was validated by the study of a set of anti-tumor drugs. With the constraints of pharmacophore obtained by DISCO implemented in SYBYL 6.8, 97 virtual bioactive compounds were generated, and their anti-tumor activities were predicted by CoMFA. Eight structures with high activity were selected and screened by the 3D-QSAR model. The most active generated structure was further investigated by modifying its structure in order to increase the activity. A comparative docking study with telomeric receptor was carried out, and the results showed that the generated structures could form more stable complexes with receptor than the reference compound selected from experimental data. This investigation showed that the proposed method was a feasible way for rational drug design with high screening efficiency.