Discovering novel carbonic anhydrase type IX (CA IX) inhibitors from seven million compounds using virtual screening and in vitro analysis

Carbonic anhydrase type IX (CA IX) enzyme is mostly over expressed in different cancer cell lines and tumor tissues. Potent CA IX inhibitors can be effective for adjusting the pH imbalance in tumor cells. In the present work, we represented the successful application of high throughput virtual screening (HTVS) of large dataset from ZINC database included of ～7 million compounds to discover novel inhibitors of CA IX. HTVS and molecular docking were performed using consequence Glide/standard precision (SP), extra precision (XP) and induced fit docking (IFD) molecular docking protocols. For each compound, docking code calculates a set of low-energy poses and then exhaustively scans the binding pocket of the target with small compounds. Novel CA IX inhibitor candidates were suggested based on molecular modeling studies and a few of them were tested using in vitro analysis. These compounds were determined as good inhibitors against human CA IX target with K_i in the range of 0.85–1.58?μM. In order to predict the pharmaceutical properties of the selected compounds, ADME (absorption, distribution, metabolism and excretion) analysis was also carried out.     

基于处方挖掘与分子动力学模拟筛选严重急性呼吸综合征冠状病毒2潜在抑制剂分子的研究

目的 从中药筛选具有潜在抑制严重急性呼吸综合征冠状病毒2 （SARS-CoV-2） 活性的成分，进一步从原子水平揭 示其抑制SARS-CoV-2 表面刺突蛋白（S 蛋白） 受体结合域（RBD） 与血管紧张素转化酶2 （ACE2） 结合的内在机制。 方法 检索新型冠状病毒（简称“新冠肺炎”） 治疗中药处方，构建“新冠肺炎中药候选活性成分数据库”。用具有ACE2 抑制活性的小分子化合物构建HipHop药效团模型，并对“新冠肺炎中药候选活性成分数据库”中活性成分筛选。采用分子 对接和分子动力学模拟方法研究候选活性成分与ACE2 的结合方式及其对SARS-CoV-2 S 蛋白与ACE2 识别的影响。 结果 本文通过中药处方挖掘和分子动力学模拟，从143 个新冠肺炎治疗中药处方中筛选出10 种可与SARS-CoV-2 S 蛋白/ 人源ACE2 识别位点结合的中药成分。其中，枇杷叶主要活性成分23-trans-p-coumaryhormentic acid 与ACE2 具有最高的亲和 力，且23-trans-p-coumaryhormentic acid 的结合可有效阻断SARS-CoV-2 S蛋白与宿主细胞ACE2 的结合。结论 本文通过虚 拟筛选发现了SARS-CoV-2 潜在抑制剂分子23-trans-p-coumaryhormentic acid，同时从原子水平预测了其抑制SARS-CoV-2 S 蛋白与ACE2 结合的内在机制，这将为SARS-CoV-2 特异性抗病毒药物的研发提供理论依据。     

Kinetic and in silico studies of hydroxy-based inhibitors of carbonic anhydrase isoforms I and II

A series of hydroxy and phenolic compounds have been assayed for the inhibition of two physiologically relevant carbonic anhydrase (CA, EC 4.2.1.1) isozymes, the cytosolic human isozymes I and II. The investigated molecules showed inhibition constants in the range of 1.07–4003 and 0.09–31.5?μM at the hCA I and hCA II enzymes, respectively. In order to investigate the binding mechanisms of these inhibitors, in silico studies were also applied. Molecular docking scores of the studied compounds are compared using three different scoring algorithms, namely Glide/SP, Glide/XP and Glide/IFD. In addition, different ADME (absorption, distribution, metabolism and excretion) analysis was performed. All the examined compounds were found within the acceptable range of pharmacokinetic profiles.     

Shape-based virtual screening,docking, and molecular dynamics simulations to identify Mtb-ASADH inhibitors

Aspartate β-semialdehyde dehydrogenase (ASADH) is a key enzyme for the biosynthesis of essential amino acids and several important metabolites in microbes. Inhibition of ASADH enzyme is a promising drug target strategy against Mycobacterium tuberculosis (Mtb). In this work, in silico approach was used to identify potent inhibitors of Mtb-ASADH. Aspartyl β-difluorophosphonate (β-AFP), a known lead compound, was used to understand the molecular recognition interactions (using molecular docking and molecular dynamics analysis). This analysis helped in validating the computational protocol and established the participation of Arg99, Glu224, Cys130, Arg249, and His256 amino acids as the key amino acids in stabilizing ligand–enzyme interactions for effective binding, an essential feature is H-bonding interactions with the two arginyl residues at the two ends of the ligand. Best binding conformation of β-AFP was selected as a template for shape-based virtual screening (ZINC and NCI databases) to identify compounds that competitively inhibit the Mtb-ASADH. The top rank hits were further subjected to ADME and toxicity filters. Final filter was based on molecular docking analysis. Each screened molecule carries the characteristics of the highly electronegative groups on both sides separated by an average distance of 6?Å. Finally, the best predicted 20 compounds exhibited minimum three H-bonding interactions with Arg99 and Arg249. These identified hits can be further used for designing the more potent inhibitors against ASADH family. MD simulations were also performed on two selected compounds (NSC4862 and ZINC02534243) for further validation. During the MD simulations, both compounds showed same H-bonding interactions and remained bound to key active residues of Mtb-ASADH.     

Ligand‐based pharmacophore modelling and screening of DNA minor groove binders targeting Staphylococcus aureus

The recognition of DNA by small molecules is of special importance in the design of new drugs. Many natural and synthetic compounds have the ability to interact with the minor groove of DNA. In the present study, identification of minor groove binding compounds was attained by the combined approach of pharmacophore modelling, virtual screening and molecular dynamics approach. Experimentally reported 32 minor groove binding compounds were used to develop the pharmacophore model. Based on the fitness score, best three pharmacophore hypotheses were selected and used as template for screening the compounds from drug bank database. This pharmacophore‐based screening provides many compounds with the same pharmacological properties. All these compounds were subjected to four phases of docking protocols with combined Glide‐quantum‐polarized ligand docking approach. Molecular dynamics results indicated that selected compounds are more active and showed good interaction in the binding site of DNA. Based on the scoring parameters and energy values, the best compounds were selected, and antibacterial activity of these compounds was identified using in vitro antimicrobial techniques. Copyright © 2014 John Wiley & Sons, Ltd.     

Prediction of estrogen receptor β ligands potency and selectivity by docking and MM-GBSA scoring methods using three different scaffolds

This study aimed to identify the docking and molecular mechanics-generalized born surface area (MM-GBSA) re-scoring parameters which can correlate the binding affinity and selectivity of the ligands towards oestrogen receptor β (ERβ). Three different series of ERβ ligands were used as dataset and the compounds were docked against ERβ (protein data bank (PDB) ID: 1QKM) using Glide and ArgusLab. Glide docking showed superior results when compared with ArgusLab. Docked poses were then rescored using Prime-MM-GBSA to calculate free energy binding. Correlations were made between observed activities of ERβ ligands with computationally predicted values from docking, binding energy parameters. ERβ ligands experimental binding affinity/selectivity did not correlate well with Glide and ArgusLab score. Whereas calculated Glide energy (coulomb-van der Waal interaction energy) correlated significantly with binding affinity of ERβ ligands (r²?=?0.66). MM-GBSA re-scoring showed correlation of r²?=?0.74 with selectivity of ERβ ligands. These results will aid the discovery of novel ERβ ligands with isoform selectivity.     

Discovery of novel inhibitors of Mycobacterium tuberculosis MurG: homology modelling,structure based pharmacophore,molecular docking,and molecular dynamics simulations

MurG (Rv2153c) is a key player in the biosynthesis of the peptidoglycan layer in Mycobacterium tuberculosis (Mtb). This work is an attempt to highlight the structural and functional relationship of Mtb MurG, the three-dimensional (3D) structure of protein was constructed by homology modelling using Discovery Studio 3.5 software. The quality and consistency of generated model was assessed by PROCHECK, ProSA and ERRAT. Later, the model was optimized by molecular dynamics (MD) simulations and the optimized model complex with substrate Uridine-diphosphate-N-acetylglucosamine (UD1) facilitated us to employ structure-based virtual screening approach to obtain new hits from Asinex database using energy-optimized pharmacophore modelling (e-pharmacophore). The pharmacophore model was validated using enrichment calculations, and finally, validated model was employed for high-throughput virtual screening and molecular docking to identify novel Mtb MurG inhibitors. This study led to the identification of 10 potential compounds with good fitness, docking score, which make important interactions with the protein active site. The 25 ns MD simulations of three potential lead compounds with protein confirmed that the structure was stable and make several non-bonding interactions with amino acids, such as Leu290, Met310 and Asn167. Hence, we concluded that the identified compounds may act as new leads for the design of Mtb MurG inhibitors.     

Design,synthesis and molecular modelling studies of some pyrazole derivatives as carbonic anhydrase inhibitors

Abstract

In this study, newly synthesised compounds 6, 8, 10 and other compounds (1–5, 7 and 9) and their inhibitory properties against the human isoforms hCA I and hCA II were reported for the first time. Compounds 1–10 showed effective inhibition profiles with K _I values in the range of 5.13–16.9?nM for hCA I and of 11.77–67.39?nM against hCA II, respectively. Molecular docking studies were also performed with Glide XP to get insight into the inhibitory activity and to evaluate the binding modes of the synthesised compounds to hCA I and II. More rigorous binding energy calculations using MM-GBSA protocol which agreed well with observed activities were then performed to improve the docking scores. Results of in silico calculations showed that all compounds obey drug likeness properties. The new compounds reported here might be promising lead compounds for the development of new potent inhibitors as alternatives to classical hCA inhibitors.     

Binding mode exploration of LuxR-thiazolidinedione analogues,e-pharmacophore-based virtual screening in the designing of LuxR inhibitors and its biological evaluation

Master quorum sensing (QS) regulator LuxR of Vibrio harveyi is a unique member of the TetR protein superfamily. Recent studies have demonstrated the contribution of thiazolidinedione analogues in blocking QS by decreasing the DNA-binding ability of LuxR. However, the precise mechanism of thiazolidinedione analogues binding to LuxR is still unclear. In the present study, molecular docking combined with molecular dynamics (MD) simulations was performed to understand the mechanism of ligand binding to the protein. The binding pattern of thiazolidinedione analogues showed strong hydrogen bonding interactions with the amine group (NH) of polar amino acid residue Asn133 and carbonyl (C=O) interaction with negatively charged amino acid residue Gln137 in the binding site of LuxR. The stability of the protein–ligand complexes was confirmed by running 50 ns of MD simulations. Further, the four-featured pharmacophore hypothesis (AHHD) consists of one acceptor (A), two hydrophobic regions (HH) and one donor (D) group was used to screen compounds from ChemBridge database. The identified hit molecules were shown to have excellent pharmacokinetic properties under the acceptable range. Based on the computational studies, ChemBridge_5343641 was selected for in vitro assays. The 1-(4-chlorophenoxy)-3-[(4,6-dimethyl-2-pyrimidinyl)thio]-2-propanol (ChemBridge_5343641) showed significant reduction in bioluminescence in a dose-dependent manner. In addition, ChemBridge_5343641 inhibits biofilm formation and motility in V. harveyi. The result from the study suggests that ChemBridge_5343641 could serve as an anti-QS molecule.     

Identification of potential inhibitors of Fasciola gigantica thioredoxin1: computational screening,molecular dynamics simulation,and binding free energy studies

Fasciola gigantica is the causative organism of fascioliasis and is responsible for major economic losses in livestock production globally. F. gigantica thioredoxin1 (FgTrx1) is an important redox-active enzyme involved in maintaining the redox homeostasis in the cell. To identify a potential anti-fasciolid compound, we conducted a structure-based virtual screening of natural compounds from the ZINC database (n = 1,67,740) against the FgTrx1 structure. The ligands were docked against FgTrx1 and 309 ligands were found to have better docking score. These compounds were evaluated for Lipinski and ADMET prediction, and 30 compounds were found to fit well for re-docking studies. After refinement by molecular docking and drug-likeness analysis, three potential inhibitors (ZINC15970091, ZINC9312362, and ZINC9312661) were identified. These three ligands were further subjected to molecular dynamics simulation (MDS) to compare the dynamics and stability of the protein structure after binding of the ligands. The binding free energy analyses were calculated to determine the intermolecular interactions. The results suggested that the two compounds had a binding free energy of –82.237, and –109.52 kJ.mol^?1 for compounds with IDs ZINC9312362 and ZINC9312661, respectively. These predicted compounds displayed considerable pharmacological and structural properties to be drug candidates. We concluded that these two compounds could be potential drug candidates to fight against F. gigantica parasites.     

Design of inhibitors of thymidylate kinase from Variola virus as new selective drugs against smallpox: part II

Abstract

Acknowledging the importance of studies toward the development of measures against terrorism and bioterrorism, this study aims to contribute to the design of new prototypes of potential drugs against smallpox. Based on a former study, nine synthetic feasible prototypes of selective inhibitors for thymidylate kinase from Variola virus (VarTMPK) were designed and submitted to molecular docking, molecular dynamics simulations and binding energy calculations. The compounds are simplifications of two more complex scaffolds, with a guanine connected to an amide or alcohol through a spacer containing ether and/or amide groups, formerly suggested as promising for the design of selective inhibitors of VarTMPK. Our study showed that, despite the structural simplifications, the compounds presented effective energy values in interactions with VarTMPK and HssTMPK and that the guanine could be replaced by a simpler imidazole ring linked to a –NH₂ group, without compromising the affinity for VarTMPK. It was also observed that a positive charge in the imidazole ring is important for the selectivity toward VarTMPK and that an amide group in the spacer does not contribute to selectivity. Finally, prototype 3 was pointed as the most promising to be synthesized and experimentally evaluated.

Communicated by Ramaswamy H. Sarma     

Identification of peptidomimetic compounds as potential inhibitors against MurA enzyme of Mycobacterium tuberculosis

Abstract

Increasing prevalence of resistance to anti-tubercular drugs has become the foremost challenge now. According to WHO, over half a million of multidrug resistance cases (rifampicin, isoniazid, etc.) were reported in 2017, mostly emerging from countries such as China, India, and Russia. Therefore, developing new drugs or repurposing existing ones is need of the hour. The Mycobacterium cell wall biogenesis pathway offers many attractive targets for drug discovery against Tuberculosis (TB). MurA, a transferase enzyme that catalyzes the initial step of peptidoglycan (PG) biosynthesis, is one among them. A peptidoglycan layer resides over the plasma membrane and is an integral component of the bacterial cell wall. Therefore, disruption of their formation through inhibition of MurA enzyme should lead to deficiency in Mycobacterium cell synthesis. Based on this strategy, we have designed this study where two libraries of peptidomimetic compounds (Asinex & ChemDiv) were first screened against our modeled MurA structure and then validated through molecular dynamic simulations. From our virtual screening, top four compounds (ChemDiv: D675-0102, D675-0217; Asinex: BDE25373574, BDE 26717803) were selected based on their docking scores, binding energies, and interactions with catalytic site residues, for further evaluation. Results revealed stable ligand-MurA interactions throughout 50?ns of MD simulation and also druggability acceptable pharmacokinetic profile for all four compounds. Thus, based on our findings, these compounds could be considered as potential inhibitors of Mycobacterium MurA enzyme and hence be further tested for in vitro experimental validation as TB therapeutic drug candidate.

Communicated by Ramaswamy H. Sarma     

Structure-based screening and molecular dynamics simulations offer novel natural compounds as potential inhibitors of Mycobacterium tuberculosis isocitrate lyase

Mycobacterium tuberculosis is the etiological agent of tuberculosis in humans and is responsible for more than two million deaths annually. M. tuberculosis isocitrate lyase (MtbICL) catalyzes the first step in the glyoxylate cycle, plays a pivotal role in the persistence of M. tuberculosis, which acts as a potential target for an anti-tubercular drug. To identify the potential anti-tuberculosis compound, we conducted a structure-based virtual screening of natural compounds from the ZINC database (n = 1,67,748) against the MtbICL structure. The ligands were docked against MtbICL in three sequential docking modes that resulted in 340 ligands having better docking score. These compounds were evaluated for Lipinski and ADMET prediction, and 27 compounds were found to fit well with re-docking studies. After refinement by molecular docking and drug-likeness analyses, three potential inhibitors (ZINC1306071, ZINC2111081, and ZINC2134917) were identified. These three ligands and the reference compounds were further subjected to molecular dynamics simulation and binding energy analyses to compare the dynamic structure of protein after ligand binding and the stability of the MtbICL and bound complexes. The binding free energy analyses were calculated to validate and capture the intermolecular interactions. The results suggested that the three compounds had a negative binding energy with ?96.462, ?143.549, and ?122.526 kJ mol^?1 for compounds with IDs ZINC1306071, ZINC2111081, and ZINC2134917, respectively. These lead compounds displayed substantial pharmacological and structural properties to be drug candidates. We concluded that ZINC2111081 has a great potential to inhibit MtbICL and would add to the drug discovery process against tuberculosis.     

Ugi Bis-Amide Derivatives as Tyrosinase Inhibitor; Synthesis,Biology Assessment,and in Silico Analysis

Herein, a straightforward synthetic strategy mediated by Ugi reaction was developed to synthesize novel series of compounds as tyrosinase inhibitors. The structures of all compounds were confirmed by FT-IR, ¹H-NMR, ¹³C-NMR, and CHNOS techniques. The tyrosinase inhibitory activities of all synthesized derivatives 5a – m were determined against mushroom tyrosinase and it was found that derivative 5c possesses the best inhibition with an IC₅₀ value of 69.53±0.042 μM compared to the rest of the synthesized derivatives. Structure–activity relationships (SARs) showed that the presence of 4-MeO or 4-NO₂ at the R² position plays a key role in tyrosinase inhibitory activities. The enzyme kinetics studies showed that compound 5c is an noncompetitive inhibitor. For in silico study, the allosteric site detection was first applied to find the appropriate binding site and then molecular docking and molecular dynamic studies were performed to reveal the position and interactions of 5c as the most potent inhibitor within the tyrosinase active site. The results showed that 5c bind well with the proposed binding site and formed a stable complex with the target protein.     

Screening of hepatitis C NS5B polymerase inhibitors containing benzothiadiazine core: a steered molecular dynamics

Hepatic C virus (HCV) is a global health problem, resulting in liver cirrhosis and inflammation that can develop to hepatocellular carcinoma and fatality. The NS5B polymerase of HCV plays an important role in viral RNA replication process, making it an attractive therapeutic target for design and development of anti-HCV drugs. To search new potent compounds against the HCV NS5B polymerase, the molecular docking and the steered molecular dynamics (SMD) simulation techniques were performed. The potential potent inhibitors of the NS5B polymerase were screened out from the ZINC database using structural similarity search and molecular docking technique. Five top-hit compounds (the ZINC compounds 49888724, 49054741, 49777239, 49793673, and 49780355) were then studied by the SMD simulations based on the hypothesis that a high rupture force relates to a high binding efficiency. The results demonstrated that the ZINC compound 49888724 had a greater maximum rupture force, reflecting a good binding strength and inhibitory potency than known inhibitors and the rest four ZINC compounds. Therefore, our finding indicated that the ZINC compound 49888724 is a potential candidate to be a novel NS5B inhibitor for further design. Besides, the van der Waals interaction could be considered as the main contribution for stabilizing the NS5B-ligand complex.     

In silico investigation of PARP-1 catalytic domains in holo and apo states for the design of high-affinity PARP-1 inhibitors

The rational design of high-affinity inhibitors of poly-ADP-ribose polymerase-1 (PARP-1) is at the heart of modern anti-cancer drug design. While relevance of enzyme to DNA repair processes in cellular environment is firmly established, the structural and functional understanding of the main determinants for high-affinity ligands controlling PARP-1 activity is still lacking. The conserved active site of PARP-1 represents an ideal target for inhibitors and may offer a novel target at the treatment of breast cancer. To fill the gap in the structural knowledge, we report on the combination of molecular dynamics (MD) simulations, principal component analysis (PCA), and conformational analysis that analyzes in great details novel binding mode for a number of inhibitors at the PARP-1. While optimization of the binding affinity for original target is an important goal in the drug design, many of the promising molecules for treatment of the breast cancer are plagued by significant cardiotoxicity. One of the most common side-effects reported for a number of polymerase inhibitors is its off-target interactions with cardiac ion channels and hERG1 channel, in particular. Thus, selected candidate PARP-1 inhibitors were also screened in silico at the central cavities of hERG1 potassium ion channel.     

Campesterol Semi-Synthetic Derivatives as Potential Antibacterial: in vitro and in silico Evaluation

In this study, twelve campesterol derivatives ( 2 – 13 ) were prepared by esterification reaction at the hydroxy group in C-3 and catalytic hydrogenation at the carbon-carbon double bond in C-5(6). All obtained compounds were characterized by IR, ¹H-NMR, ¹³C-NMR, and MS spectra. Campesterol ( 1 ) and its derivatives ( 2 – 13 ) were evaluated in vitro against Staphylococcus aureus (ATCC 6538), Streptococcus mutans (ATCC 0046), Escherichia coli (ATCC 10536), Pseudomonas aeruginosa (ATCC 15442), and Klebsiella pneumoniae (ATCC 10031) using the microdilution method. Among tested compounds, 4 , 6 , 9 , 11 , 12 , and 13 displayed the best antibacterial activity. Moreover, to support the antibacterial activity experiments, the investigation of molecular interactions of more active compounds, and also compound 1 and neomycin, used as starting material and positive control, respectively, at the binding site of the target proteins was performed using molecular docking simulations. Four compounds ( 7 , 9 , 10 and 11 ) are herein described for the first time.