Anti-tubercular activity and molecular docking studies of indolizine derivatives targeting mycobacterial InhA enzyme

A series of 1,2,3-trisubstituted indolizines (2a–2f, 3a–3d, and 4a–4c) were screened for in vitro whole-cell anti-tubercular activity against the susceptible H37Rv and multidrug-resistant (MDR) Mycobacterium tuberculosis (MTB) strains. Compounds 2b–2d, 3a–3d, and 4a–4c were active against the H37Rv-MTB strain with minimum inhibitory concentration (MIC) ranging from 4 to 32 µg/mL, whereas the indolizines 4a–4c, with ethyl ester group at the 4-position of the benzoyl ring also exhibited anti-MDR-MTB activity (MIC = 16–64 µg/mL). In silico docking study revealed the enoyl-acyl carrier protein reductase (InhA) and anthranilate phosphoribosyltransferase as potential molecular targets for the indolizines. The X-ray diffraction analysis of the compound 4b was also carried out. Further, a safety study (in silico and in vitro) demonstrated no toxicity for these compounds. Thus, the indolizines warrant further development and may represent a novel promising class of InhA inhibitors and multi-targeting agents to combat drug-sensitive and drug-resistant MTB strains.     

Larvicidal activity of novel anthraquinone analogues and their molecular docking studies

To investigate the larvicidal activities of novel anthraquinones (1a-1k) against Culex quinquefasciatus mosquito larvae. Novel anthraquinones (1a-1k) derivatives were synthesis via condensation method. The compounds were confirmed through FT-IR spectroscopy, ¹H & ¹³C NMR spectrum, and mass spectral studies. The larvicidal activity of compound 1c was highly active LD₅₀ 20.92 µg/mL against Culex quinquefasciatus compared standard permethrin with LD₅₀ 25.49 µg/mL. Molecular docking studies were carried out for compound 1c against Odorant-binding protein of Culex quinquefasciatus. The compound 1c (−9.8 Kcal/mol) was a potent larvicide with more binding energy than control permethrin (−9.7 Kcal/mol). Therefore, compound (1c) may be more significant inhibitors of mosquito larvicidal.     

In vitro anti-TB properties,in silico target validation,molecular docking and dynamics studies of substituted 1,2,4-oxadiazole analogues against Mycobacterium tuberculosis

The alarming increase in multi- and extensively drug-resistant (MDR and XDR) strains of Mycobacterium tuberculosis (MTB) has triggered the scientific community to search for novel, effective, and safer therapeutics. To this end, a series of 3,5-disubstituted-1,2,4-oxadiazole derivatives (3a–3i) were tested against H37Rv, MDR and XDR strains of MTB. Of which, compound 3a with para-trifluorophenyl substituted oxadiazole showed excellent activity against the susceptible H37Rv and MDR-MTB strain with a MIC values of 8 and 16 µg/ml, respectively.To understand the mechanism of action of these compounds (3a–3i) and identify their putative drug target, molecular docking and dynamics studies were employed against a panel of 20 mycobacterial enzymes reported to be essential for mycobacterial growth and survival. These computational studies revealed polyketide synthase (Pks13) enzyme as the putative target. Moreover, in silico ADMET predictions showed satisfactory properties for these compounds, collectively, making them, particularly compound 3a, promising leads worthy of further optimisation.     

Molecular modeling studies of atorvastatin analogues as HMGR inhibitors using 3D-QSAR,molecular docking and molecular dynamics simulations

3-Hydroxy-3-methylglutaryl coenzyme-A reductase (HMGR) is generally regarded as targets for the treatment of hypercholesterolemia. HMGR inhibitors (more commonly known as statins) are discovered as plasma cholesterol lowering molecules. In this work, 120 atorvastatin analogues were studied using a combination of molecular modeling techniques including three-dimensional quantitative structure–activity relationship (3D-QSAR), molecular docking and molecular dynamics (MD) simulation. The results show that the best CoMFA (comparative molecular field analysis) model has q² = 0.558 and r² = 0.977, and the best CoMSIA (comparative molecular similarity indices analysis) model has q² = 0.582 and r² = 0.919. Molecular docking and MD simulation explored the binding relationship of the ligand and the receptor protein. The calculation results indicated that the hydrophobic and electrostatic fields play key roles in QSAR model. After MD simulation, we found four vital residues (Lys735, Arg590, Asp690 and Asn686) and three hydrophobic regions in HMGR binding site. The calculation results show that atorvastatin analogues obtained by introduction of F atoms or gem-difluoro groups could obviously improve the inhibitory activity. The new HMGR inhibitor analogues design in this Letter had been submitted which is being currently synthesized by our laboratories.     

Pharmacophore modeling,molecular docking and molecular dynamics studies on natural products database to discover novel skeleton as non-purine xanthine oxidase inhibitors

Gout is a common inflammatory arthritis caused by the deposition of urate crystals within joints. It is increasingly in prevalence during the past few decades as shown by the epidemiological survey results. Xanthine oxidase (XO) is a key enzyme to transfer hypoxanthine and xanthine to uric acid, whose overproduction leads to gout. Therefore, inhibiting the activity of xanthine oxidase is an important way to reduce the production of urate. In the study, in order to identify the potential natural products targeting XO, pharmacophore modeling was employed to filter databases. Here, two methods, pharmacophore based on ligand and pharmacophore based on receptor-ligand, were constructed by Discovery Studio. Then GOLD was used to refine the potential compounds with higher fitness scores. Finally, molecular docking and dynamics simulations were employed to analyze the interactions between compounds and protein. The best hypothesis was set as a 3D query to screen database, returning 785 and 297 compounds respectively. A merged set of the above 1082 molecules was subjected to molecular docking, which returned 144 hits with high-fitness scores. These molecules were clustered in four main kinds depending on different backbones. What is more, molecular docking showed that the representative compounds established key interactions with the amino acid residues in the protein, and the RMSD and RMSF of molecular dynamics results showed that these compounds can stabilize the protein. The information represented in the study confirmed previous reports. And it may assist to discover and design new backbones as potential XO inhibitors based on natural products.     

Benzothiazole analogs as potential anti-TB agents: computational input and molecular dynamics

Biotin is very important for the survival of Mycobacterium tuberculosis. 7,8-Diamino pelargonic acid aminotransaminase (DAPA) is a transaminase enzyme involved in the biosynthesis of biotin. The benzothiazole title compounds were investigated for their in vitro anti-tubercular activity against two tubercular strains: H37Rv (ATCC 25,177) and MDR-MTB (multidrug-resistant M. tuberculosis, resistant to isoniazid, rifampicin, and ethambutol) by an agar incorporation method. The possible binding mode and predicted affinity were computed using a molecular docking study. Among the synthesized compounds in the series, the title compound {2-(benzo[d]thiazol-2-yl-methoxy)-5-fluorophenyl}-(4-chlorophenyl)-methanone was found to exhibit significant activity with minimum inhibitory concentrations of 1 μg/mL and 2 μg/mL against H37Rv and MDR-MTB, respectively; this compound showed the highest binding affinity (–24.75 kcal/mol) as well.     

Synthesis,anticancer activity and molecular docking studies on 1,2-diarylbenzimidazole analogues as anti-tubulin agents

Twenty-four 1,2-diarylbenzimidazole derivatives were designed, synthesized and biologically evaluated. It turned out that most of them were potential anticancer drugs. Among them, compound c24 showed the highest anti-tumor activity (GI₅₀ = 0.71–2.41 μM against HeLa, HepG2, A549 and MCF-7 cells), and low toxicity to normal cells (CC₅₀ > 100 μM against L02 cells). In the microtubule binding assay, c24 showed the most potent inhibition of microtubule polymerization (IC₅₀ = 8.47 μM). The binding ability of compound c24 to tubulin crystal was verified by molecular docking simulation experiment. Further studies on HepG2 and HeLa cells showed that compound c24 could cause mitotic arrest of tumor cells to G2/M phase then inducing apoptosis. To sum up, compound c24 is a promising microtubule assembly inhibitor.     

Pharmacophore generation,atom-based 3D-QSAR,molecular docking and molecular dynamics simulation studies on benzamide analogues as FtsZ inhibitors

FtsZ is an appealing target for the design of antimicrobial agent that can be used to defeat the multidrug-resistant bacterial pathogens. Pharmacophore modelling, molecular docking and molecular dynamics (MD) simulation studies were performed on a series of three-substituted benzamide derivatives. In the present study a five-featured pharmacophore model with one hydrogen bond acceptors, one hydrogen bond donors, one hydrophobic and two aromatic rings was developed using 97 molecules having MIC values ranging from .07 to 957 μM. A statistically significant 3D-QSAR model was obtained using this pharmacophore hypothesis with a good correlation coefficient (R² = .8319), cross validated coefficient (Q² = .6213) and a high Fisher ratio (F = 103.9) with three component PLS factor. A good correlation between experimental and predicted activity of the training (R² = .83) and test set (R² = .67) molecules were displayed by ADHRR.1682 model. The generated model was further validated by enrichment studies using the decoy test and MAE-based criteria to measure the efficiency of the model. The docking studies of all selected inhibitors in the active site of FtsZ protein showed crucial hydrogen bond interactions with Val 207, Asn 263, Leu 209, Gly 205 and Asn-299 residues. The binding free energies of these inhibitors were calculated by the molecular mechanics/generalized born surface area VSGB 2.0 method. Finally, a 15 ns MD simulation was done to confirm the stability of the 4DXD–ligand complex. On a wider scope, the prospect of present work provides insight in designing molecules with better selective FtsZ inhibitory potential.     

Synthesis,in vitro urease inhibitory activity,and molecular docking studies of thiourea and urea derivatives

The current study deals with the synthesis of urea and thiourea derivatives 1–37 which were characterized by various spectroscopic techniques including FAB-MS, ¹H-, and ¹³C NMR. The synthetic compounds were subjected to urease inhibitory activity and compounds exhibited good to moderate urease inhibitory activity having IC₅₀ values in range of 10.11–69.80 µM. Compound 1 (IC₅₀ = 10.11 ± 0.11 µM) was found to be most active and even better as compared to the standard acetohydroxamic acid (IC₅₀ = 27.0 ± 0.5 µM). A limited structure–activity relationship (SAR) was established and the compounds were also subjected to docking studies to confirm the binding interactions of ligands (compounds) with the active site of enzyme.     

Deciphering the GPER/GPR30-agonist and antagonists interactions using molecular modeling studies,molecular dynamics,and docking simulations

The G-protein coupled estrogen receptor 1 GPER/GPR30 is a transmembrane seven-helix (7TM) receptor involved in the growth and proliferation of breast cancer. Due to the absence of a crystal structure of GPER/GPR30, in this work, molecular modeling studies have been carried out to build a three-dimensional structure, which was subsequently refined by molecular dynamics (MD) simulations (up to 120 ns). Furthermore, we explored GPER/GPR30’s molecular recognition properties by using reported agonist ligands (G1, estradiol (E2), tamoxifen, and fulvestrant) and the antagonist ligands (G15 and G36) in subsequent docking studies. Our results identified the E2 binding site on GPER/GPR30, as well as other receptor cavities for accepting large volume ligands, through GPER/GPR30 π–π, hydrophobic, and hydrogen bond interactions. Snapshots of the MD trajectory at 14 and 70 ns showed almost identical binding motifs for G1 and G15. It was also observed that C107 interacts with the acetyl oxygen of G1 (at 14 ns) and that at 70 ns the residue E275 interacts with the acetyl group and with the oxygen from the other agonist whereas the isopropyl group of G36 is oriented toward Met141, suggesting that both C107 and E275 could be involved in the protein activation. This contribution suggest that GPER1 has great structural changes which explain its great capacity to accept diverse ligands, and also, the same ligand could be recognized in different binding pose according to GPER structural conformations.     

A combination of pharmacophore modeling,atom-based 3D-QSAR,molecular docking and molecular dynamics simulation studies on PDE4 enzyme inhibitors

Phosphodiesterases 4 enzyme is an attractive target for the design of anti-inflammatory and bronchodilator agents. In the present study, pharmacophore and atom-based 3D-QSAR studies were carried out for pyrazolopyridine and quinoline derivatives using Schrödinger suite 2014-3. A four-point pharmacophore model was developed using 74 molecules having pIC₅₀ ranging from 10.1 to 4.5. The best four feature model consists of one hydrogen bond acceptor, two aromatic rings, and one hydrophobic group. The pharmacophore hypothesis yielded a statistically significant 3D-QSAR model, with a high correlation coefficient (R^2?=?.9949), cross validation coefficient (Q^2?=?.7291), and Pearson-r (.9107) at six component partial least square factor. The external validation indicated that our QSAR model possessed high predictive power with R² value of .88. The generated model was further validated by enrichment studies using the decoy test. Molecular docking, free energy calculation, and molecular dynamics (MD) simulation studies have been performed to explore the putative binding modes of these ligands. A 10-ns MD simulation confirmed the docking results of both stability of the 1XMU–ligand complex and the presumed active conformation. Outcomes of the present study provide insight in designing novel molecules with better PDE4 inhibitory activity.     

Pharmacophore modeling,multiple docking,and molecular dynamics studies on Wee1 kinase inhibitors

Wee1-like protein kinase (Wee1) is a tyrosine kinase that regulates the G₂ checkpoint and prevents entry into mitosis in response to DNA damage. Based on a series of signaling pathways initiated by Wee1, Wee1 has been recognized as a potential target for cancer therapy. To discover potent Wee1 inhibitors with novel scaffolds, ligand-based pharmacophore model has been built based on 101 known Wee1 inhibitors. Then the best pharmacophore model, AADRRR.340, with good partial least square (PLS) statistics (R²?=?0.9212, Q²?=?0.7457), was selected and validated. The validated model was used as a three-dimensional (3D) search query for databases virtual screening. The filtered molecules were further analyzed and refined by Lipinski’s rule of 5, multiple docking procedures (high throughput virtual screening (HTVS), standard precision (SP), genetic optimization for ligand docking (GOLD), extra precision (XP), and unique quantum polarized ligand docking (QPLD)); absorption, distribution, metabolism, excretion, and toxicity (ADMET) screening; and the Prime/molecular mechanics generalized born surface area (MM-GBSA) method binding free energy calculations. Eight leads were identified as potential Wee1 inhibitors, and a 50?ns molecular dynamics (MD) simulation was carried out for top four inhibitors to predict the stability of ligand–protein complex. Molecular mechanics Poisson–Boltzmann surface area (MM-PBSA) based on MD simulation and the energy contribution per residue to the binding energy were calculated. In the end, three hits with good stabilization and affinity to protein were identified.

Communicated by Ramaswamy H. Sarma     

Synthesis of 1H-1,2,3-triazole derivatives as new α-glucosidase inhibitors and their molecular docking studies

Inhibition of α-glucosidase is an effective strategy for controlling the post-prandial hyperglycemia in diabetic patients. For the identification of new inhibitors of this enzyme, a series of new (R)-1-(2-(4-bromo-2-methoxyphenoxy) propyl)-4-(4-(trifluoromethyl) phenyl)-1H-1,2,3-triazole derivatives were synthesized (8a–d and 10a–e). The structures were confirmed by NMR, mass spectrometry and, in case of compound 8a, by single crystal X-ray crystallography. The α-glucosidase inhibitory activities were investigated in vitro. Most derivatives exhibited significant inhibitory activity against α-glucosidase enzyme. Their structure-activity relationship and molecular docking studies were performed to elucidate the active pharmacophore against this enzyme. Compound 10b was the most active analogue with IC₅₀ value of 14.2 µM, while compound 6 was found to be the least active having 218.1 µM. A preliminary structure-activity relationship suggested that the presence of 1H-1,2,3-triazole ring in 1H-1,2,3-triazole derivatives is responsible for this activity and can be used as anti-diabetic drugs. The molecular docking studies of all active compounds were performed, in order to understand the mode of binding interaction and the energy of this class of compounds.     

Synthesis,molecular docking,and antioxidant activity of new fluorescent tetrafluoroterphenyl analogues

Nucleophilic aromatic substitution (S_NAr) chemistry has been applied to develop many functionalized pentafluorobenzene derivatives. Those compounds are highly specific at the para position of the fluorinated ring. Therefore, they are typical adducts for the preparation of antioxidant molecular systems. In this context, we report the use of S_NAr chemistry as a suitable and simple approach for the synthesis of fluorescent antioxidant perfluorinated materials bearing ether bonds in various para-substituted alkoxy chains and with high purity and excellent yields. The fluoroterphenyl core was prepared via alkylation, Cu(I)-assisted decarboxylation, and cross-coupling using the potassium salt of fluorobenzoate, followed by the reaction with different alcohols. The structures of the synthesized fluoroterphenyl adducts were investigated using FT-IR, ¹H NMR, ¹³C NMR, and ¹⁹F NMR spectroscopy. The emission spectra and absorption spectra showed solvatochromism. The newly prepared tetrafluoroterphenyl analogues were investigated by antioxidant examination using the 2,2-diphenyl-1-picrylhydrazyl assay. Results were compared with ascorbic acid and butylated hydroxytoluene as references, and revealed that the tetrafluoroterphenyl analogues containing a decyl chain had the highest activity, with an IC₅₀ value of 22.36 ± 0.19 g/ml. The produced tetrafluoroterphenyl analogues were used in molecular docking strategies with a Protein Data Bank protein ID 5IKQ. The antioxidant investigations and docking results were convergent.     

Structural studies,homology modeling and molecular docking of novel non-competitive antagonists of GluK1/GluK2 receptors

Non-competitive ligands of kainate receptors have focused significant attention as medicinal compounds because they seem to be better tolerated than competitive antagonists and uncompetitive blocker of these receptors. Here we present structural studies (X-ray structure determination, NMR and MS spectra) of novel indole-derived non-competitive antagonists of GluK1/GluK2 receptors, homology models of GluK1 and GluK2 receptors based on novel AMPA receptor template as well as molecular docking of ligands to their molecular targets. We find that the allosteric site is in the receptor transduction domain, in one receptor subunit, not between the two subunits as it was indicated by our earlier studies.     

Identification of nonplanar small molecule for G-quadruplex grooves: molecular docking and molecular dynamic study

DNA G-quadruplex is an attractive drug target for anticancer therapy. Most G-quadruplex ligands have little selectivity, due to π-stacking interaction with common G-tetrads surface. Thanks to the varieties of G-quadruplex grooves, the groove-binding ligand is expected to create high selectivity. Therefore, developing novel molecular geometries that target G-quadruplex groove has been paid growing attention. In this work, steroid FG, a special nonplanar and nonaromatic small molecule, interacting with different conformations of G-quadruplexes has been studied by molecular docking and molecular dynamics simulations. The results showed the selectivity of the hydrophobic group of steroid FG for the wide groove of antiparallel G-quadruplex. The methyl groups on the tetracyclic ring of steroid represent the specific binding ability for the small hydrophobic cavity formed by reversed stacking of G-tetrads in antiparallel G-quadruplex groove. This work provides new insight for developing new classes of G-quadruplex groove-binding ligands.     

Copper(I)-catalysed synthesis of symmetrical perfluoroterphenyl analogues; fluorescence,antioxidant and molecular docking studies

Pentafluoroaryl analogues have been found to exhibit para specific nucleophilic aromatic substitution (S_NAr). Herein, we describe the use of S_NAr chemistry to create luminous perfluorinated symmetrical terphenyls. Both of S_NAr chemistry and copper(I)-catalysed decarboxylative cross-coupling were applied for the synthesis of the perfluorinated symmetrical terphenyls in high yields from the corresponding derivatives of aryl iodide and potassium salt of fluorobenzoate. A series of perfluorinated symmetrical terphenyls with different para alkoxy chains were synthesized. The synthesized perfluorinated terphenyl adducts were confirmed via elemental analysis, Fourier-transform infrared (FTIR), proton (¹H) carbon-13 (¹³C) and fluorine-19 (¹⁹F) nuclear magnetic resonance (NMR) spectra. The absorbance and fluorescence spectra showed solvatochromic activities. The new synthesized fluoroterphenyl hybrids were screened against antioxidant inspection over DPPH (2,2-diphenyl-1-picrylhydrazyl) performance, in assessment of vitamin C and butylated hydroxytoluene (BHT) as standard drugs exposed that fluoroterphenyl hybrid covering decyl hydrocarbons exhibited highest effectiveness through half maximal inhibitory concentration (IC₅₀) values of 21.74 μg/ml. Additionally, molecular docking procedures of the synthesized fluoroterphenyl hybrids were employed by using protein data bank (PDB ID: 5IKQ). The docking simulation displayed convenient and recognized findings with the antioxidant examination.     

Ligand-based pharmacophore modeling and docking studies on vitamin D receptor inhibitors

In recent years, pharmacophore modeling and molecular docking approaches have been extensively used to characterize the structural requirements and explore the conformational space of a ligand in the binding pocket of the selected target protein. Herein, we report a pharmacophore modeling and molecular docking of 45 compounds comprising of the indole scaffold as vitamin D receptor (VDR) inhibitors. Based on the selected best hypothesis (DRRRR.61), an atom-based three-dimensional quantitative structure-activity relationships model was developed to rationalize the structural requirement of biological activity modulating components. The developed model predicted the binding affinity for the training set and test set with R²_(training) = 0.8869 and R²_(test) = 0.8139, respectively. Furthermore, molecular docking and dynamics simulation were performed to understand the underpinning of binding interaction and stability of selected VDR inhibitors in the binding pocket. In conclusion, the results presented here, in the form of functional and structural data, agreed well with the proposed pharmacophores and provide further insights into the development of novel VDR inhibitors with better activity.     

Computational investigations of gram-negative bacteria phosphopantetheine adenylyltransferase inhibitors using 3D-QSAR,molecular docking and molecular dynamic simulations

Abstract

Phosphopantetheine adenylyltransferase (PPAT) has been recognized as a promising target to develop novel antimicrobial agents, which is a hexameric enzyme that catalyzes the penultimate step in coenzyme A biosynthesis. In this work, molecular modeling study was performed with a series of PPAT inhibitors using molecular docking, three-dimensional qualitative structure-activity relationship (3D-QSAR) and molecular dynamic (MD) simulations to reveal the structural determinants for their bioactivities. Molecular docking study was applied to understand the binding mode of PPAT with its inhibitors. Subsequently, 3D-QSAR model was constructed to find the features required for different substituents on the scaffolds. For the best comparative molecular field analysis (CoMFA) model, the Q² and R² values of which were calculated as 0.702 and 0.989, while they were calculated as 0.767 and 0.983 for the best comparative molecular similarity index analysis model. The statistical data verified the significance and accuracy of our 3D-QSAR models. Furthermore, MD simulations were carried out to evaluate the stability of the receptor–ligand contacts in physiological conditions, and the results were consistent with molecular docking studies and 3D-QSAR contour map analysis. Binding free energy was calculated with molecular mechanics generalized born surface area approach, the result of which coincided well with bioactivities and demonstrated that van der Waals accounted for the largest portion. Overall, our study provided a valuable insight for further research work on the recognition of potent PPAT inhibitors.

Communicated by Ramaswamy H. Sarma     

Molecular modeling studies of repandusinic acid as potent small molecule for hepatitis B virus through molecular docking and ADME analysis

Background: Hepatitis B virus (HBV) has affected over 300 million people worldwide which causes to induce mostly liver disease and liver cancer. It is a member of the family Hepadnaviridae which is a small DNA virus with unusual characters like retroviruses. Generally, hepatoprotective drugs provoke some side effects in human beings. For the reason, this study aims to identify alternative drug molecules from the natural source of medicinal plants with smaller quantity of side effects than those conventional drugs in treating HBV. Methods: We developed computational methods for calculating drug and target binding resemblance using the Maestro v10.2 of Schrodinger suite. The target and ligand molecules were obtained from recognized databases. Ligand molecules of 40 phytoconstituents were retrieved from variety of plants after we executed crucial analyses such as molecular docking and absorption, distribution, metabolism, and excretion (ADME) analysis.Results: In the docking analysis, the natural analogues repandusinic acid showed better docking scores of –14.768 with good binding contacts. The remaining bioactive molecules corilagin, furosin, nirurin, iso-quercetin and gallocatechin also showed better docking scores.Conclusion: This computational analysis reveals that repandusinic acid is a suitable drug candidate for HBV. Therefore, we recommend that this analogue is suitable in further exploration using in vitro studies.