Structural basis for the activity of drugs that inhibit phosphodiesterases

Phosphodiesterases (PDEs) comprise a large family of enzymes that catalyze the hydrolysis of cAMP or cGMP and are implicated in various diseases. We describe the high-resolution crystal structures of the catalytic domains of PDE4B, PDE4D, and PDE5A with ten different inhibitors, including the drug candidates cilomilast and roflumilast, for respiratory diseases. These cocrystal structures reveal a common scheme of inhibitor binding to the PDEs: (i) a hydrophobic clamp formed by highly conserved hydrophobic residues that sandwich the inhibitor in the active site; (ii) hydrogen bonding to an invariant glutamine that controls the orientation of inhibitor binding. A scaffold can be readily identified for any given inhibitor based on the formation of these two types of conserved interactions. These structural insights will enable the design of isoform-selective inhibitors with improved binding affinity and should facilitate the discovery of more potent and selective PDE inhibitors for the treatment of a variety of diseases.     

Pharmacophore Modeling and Virtual Screening for the Discovery of New type 4 cAMP Phosphodiesterase (PDE4) Inhibitors

Type 4 cAMP phosphodiesterase (PDE4) inhibitors show a broad spectrum of anti-inflammatory effects in almost all kinds of inflamed cells, by an increase in cAMP levels which is a pivotal second messenger responsible for various biological processes. These inhibitors are now considered as the potential drugs for treatment of chronic inflammatory diseases. However, some recently marketed inhibitors e.g., roflumilast, have shown adverse effects such as nausea and emesis, thus restricting its use. In order to identify novel PDE4 inhibitors with improved therapeutic indexes, a highly correlating (r = 0.963930) pharmacophore model (Hypo1) was established on the basis of known PDE4 inhibitors. Validated Hypo1 was used in database screening to identify chemical with required pharmacophoric features. These compounds are further screened by using the rule of five, ADMET and molecular docking. Finally, twelve hits which showed good results with respect to following properties such as estimated activity, calculated drug-like properties and scores were proposed as potential leads to inhibit the PDE4 activity. Therefore, this study will not only assist in the development of new potent hits for PDE4 inhibitors, but also give a better understanding of their interaction with PDE4. On a wider scope, this will be helpful for the rational design of novel potent enzyme inhibitors.     

In-silico analysis of Sirt2 from Schistosoma mansoni: structures,conformations, and interactions with inhibitors

Sirtuins are NAD+-dependent lysine deacetylases member of the class III HDAC family. These are demonstrated to be therapeutic targets in parasitic diseases like schistosomiasis. Observations suggested that sirtuin enzyme is necessary for the functionality of fe/male reproductive system, due to which SmSirt2 is treated as a potential therapeutic target. There are no structural and molecular features of SmSirt2 have been reported yet. In this study, homology modeling has been used to determine the three-dimensional features of the SmSITRT2. Further, structure validation has been performed by energy minimization and Ramachandran plot. Validated structures are further subjected to molecular docking and virtual screening to find the best lead molecules for downstream analysis. Ten lead molecules were selected while comparing virtual screening of hSirt2 and SmSirt2 both. These leads are further compared with AKG2 which is known inhibitor of hSirt2 (?8.8 kcal/mol). Out of selected 10 leads, four of them (ZINC23995485 (?9.5 kcal/mol), ZINC53298162 (?9.4 kcal/mol), ZINC70927268 (?10.0 kcal/mol), ZINC89878705 (?11.2 kcal/mol)) have shown better interaction with SmSirt2, in which ZINC89878705 (?11.2 kcal/mol) shows a more compact packing as compared to AKG2 and rest of ligands. These molecules could be further subject to in vitro study and model of SmSirt2 has been proposed for further structure-based drug design projects concerning sirtuins from Schistosoma mansoni.     

Identification and evaluation of bioactive natural products as potential inhibitors of human microtubule affinity-regulating kinase 4 (MARK4)

Microtubule affinity-regulating kinase 4 (MARK4) has recently been identified as a potential drug target for several complex diseases including cancer, diabetes and neurodegenerative disorders. Inhibition of MARK4 activity is an appealing therapeutic option to treat such diseases. Here, we have performed structure-based virtual high-throughput screening of 100,000 naturally occurring compounds from ZINC database against MARK4 to find its potential inhibitors. The resulted hits were selected, based on the binding affinities, docking scores and selectivity. Further, binding energy calculation, Lipinski filtration and ADMET prediction were carried out to find safe and better hits against MARK4. Best 10 compounds bearing high specificity and binding efficiency were selected, and their binding pattern to MARK4 was analyzed in detail. Finally, 100 ns molecular dynamics simulation was performed to evaluate; the dynamics stability of MARK4-compound complex. In conclusion, these selected natural compounds from ZINC database might be potential leads against MARK4, and can further be exploited in drug design and development for associated diseases.     

High throughput screening,docking, and molecular dynamics studies to identify potential inhibitors of human calcium/calmodulin-dependent protein kinase IV

Calcium/calmodulin-dependent protein kinase IV (CAMKIV) is associated with many diseases including cancer and neurodegenerative disorders and thus being considered as a potential drug target. Here, we have employed the knowledge of three-dimensional structure of CAMKIV to identify new inhibitors for possible therapeutic intervention. We have employed virtual high throughput screening of 12,500 natural compounds of Zinc database to screen the best possible inhibitors of CAMKIV. Subsequently, 40 compounds which showed significant docking scores (?11.6 to ?10.0?kcal/mol) were selected and further filtered through Lipinski rule and drug likeness parameter to get best inhibitors of CAMKIV. Docking results are indicating that ligands are binding to the hydrophobic cavity of the kinase domain of CAMKIV and forming a significant number of non-covalent interactions. Four compounds, ZINC02098378, ZINC12866674, ZINC04293413, and ZINC13403020, showing excellent binding affinity and drug likeness were subjected to molecular dynamics simulation to evaluate their mechanism of interaction and stability of protein-ligand complex. Our observations clearly suggesting that these selected ligands may be further employed for therapeutic intervention to address CAMKIV associated diseases.

Communicated by Ramaswamy H. Sarma     

Comparative docking of dual conformations in human fatty acid synthase thioesterase domain reveals potential binding cavity for virtual screening of ligands

Human fatty acid synthase (hFASN), a homo dimeric lipogenic enzyme with seven catalytic domains, is an important clinical target in cancer, metabolic syndrome and infections. Here, molecular modelling and docking methods were implemented to examine the inter-molecular interactions of thioesterase (TE) domain in hFASN with its physiological substrate, and to identify potential chemical inhibitors. TE catalyses the hydrolysis of thioester bond between palmitate and the 4’ phosphopantetheine of acyl carrier protein, releasing 16-carbon palmitate. The crystal structure of hFASN TE in two inhibitory conformations (A and B) were geometry-optimized and used for molecular docking with palmitate, orlistat (a known FASN inhibitor) and virtual screening against compounds from National Cancer Institute (NCI) database. Relatively, low binding affinity was observed during the complex formation of palmitate with A (?.164 kcal/mol) and B (?.332 kcal/mol) forms of TE, when compared with orlistat-docked TE (A form: ?5.872 kcal/mol and B form: ?5.484 kcal/mol), clearly indicating that the native inhibited conformation (crystal structure) was unfavourable for substrate binding. We used these orlistat dual binding modes as positive controls for prioritizing the ligands during virtual screening. From 2, 31,617 molecules in the NCI database, 916 high-scoring compounds (hit ligands) were obtained for A-form and 4582 for B-form of the TE-domain, which were then ranked according to glide docking score, XP H bond score, absorption, distribution, metabolism and excretion and binding free energy (Prime/MM-GBSA). Consequently, two top scoring ligands (NSC: 319661 and NSC: 153166) emerged as promising drug candidates that may be tested in FASN-over-expressing diseases.     

Identification of novel natural inhibitors of Opisthorchis felineus cytochrome P450 using structure-based screening and molecular dynamic simulation

Opisthorchis felineus is the etiological agent of opisthorchiasis in humans. O. felineus cytochrome P450 (OfCYP450) is an important enzyme in the parasite xenobiotic metabolism. To identify the potential anti-opisthorchid compound, we conducted a structure-based virtual screening of natural compounds from the ZINC database (n = 1,65,869) against the OfCYP450. The ligands were screened against OfCYP450 in four sequential docking modes that resulted in 361 ligands having better docking score. These compounds were evaluated for Lipinski and ADMET prediction, and 10 compounds were found to fit well with re-docking studies. After refinement by docking and drug-likeness analyses, four potential inhibitors (ZINC2358298, ZINC8790946, ZINC70707116, and ZINC85878789) were identified. These ligands with reference compounds (itraconazole and fluconazole) were further subjected to molecular dynamics simulation (MDS) and binding energy analyses to compare the dynamic structure of protein after ligand binding and the stability of the OfCYP450 and bound complexes. The binding energy analyses were also calculated. The results suggested that the compounds had a negative binding energy with ?259.41, ?110.09, ?188.25, ?163.30, ?202.10, and ?158.79 kJ mol^?1 for itraconazole, fluconazole, and compounds with IDs ZINC2358298, ZINC8790946, ZINC70707116, and ZINC85878789, respectively. These lead compounds displayed significant pharmacological and structural properties to be drug candidates. On the basis of MDS results and binding energy analyses, we concluded that ZINC8790946, ZINC70707116, and ZINC85878789 have excellent potential to inhibit OfCYP450.     

Molecular docking and virtual screening for novel protein tyrosine phosphatase 1B (PTP1B) inhibitors

Protein tyrosine phosphatase 1B (PTP1B) functions as major negative regulator of insulin and leptin signaling pathways. In view of this, PTP1B is an significant target for drug development against cancer, diabetes and obesity. The aim of the current study is to identify PTP1B inhibitors by means of virtual screening with docking. 523,366 molecules from ZINC database have been screened and based on DOCK grid scores and hydrogen bonding interactions five new potential inhibitors were identified. ZINC12502589, ZINC13213457, ZINC25721858, ZINC31392733 and ZINC04096400 were identified as potential lead molecules for inhibition of PTP1B. The identified molecules were subjected to Lipinski''s rule of five parameters and found that they did not violate any rule. More specific analysis of pharmacological parameters may be scrutinized through a complete ADME/Tox evaluation. Pharma algorithm was used to Calculate ADME–Tox profiles for such molecules. In general, all the molecules presented advantages and as well as disadvantages when compared to each other. No marked difference in health effects and toxicity profiles were observed among these molecules.     

The Unrecognized Effects of Phosphodiesterase 4 on Epithelial Cells in Pulmonary Inflammation

Acute pulmonary inflammation is characterized by migration of polymorphonuclear neutrophils (PMNs) into the different compartments of the lung, passing an endothelial and epithelial barrier. Recent studies showed evidence that phosphodiesterase (PDE)4-inhibitors stabilized endothelial cells. PDE4B and PDE4D subtypes play a pivotal role in inflammation, whereas blocking PDE4D is suspected to cause gastrointestinal side effects. We thought to investigate the particular role of the PDE4-inhibitors roflumilast and rolipram on lung epithelium. Acute pulmonary inflammation was induced by inhalation of LPS. PDE4-inhibitors were administered i.p. or nebulized after inflammation. The impact of PDE4-inhibitors on PMN migration was evaluated in vivo and in vitro. Microvascular permeability, cytokine levels, and PDE4B and PDE4D expression were analyzed. In vivo, both PDE4-inhibitors decreased transendothelial and transepithelial migration even when administered after inflammation, whereas roflumilast showed a superior effect compared to rolipram on the epithelium. Both inhibitors decreased TNFα, IL6, and CXCL2/3. CXCL1, the strong PMN chemoattractant secreted by the epithelium, was significantly more reduced by roflumilast. In vitro assays with human epithelium also emphasized the pivotal role of roflumilast on the epithelium. Additionally, LPS-induced stress fibers, an essential requirement for a direct migration of PMNs into the alveolar space, were predominantly reduced by roflumilast. Expression of PDE4B and PDE4D were both increased in the lungs by LPS, PDE4-inhibitors decreased mainly PDE4B. The topical administration of PDE4-inhibitors was also effective in curbing down PMN migration, further highlighting the clinical potential of these compounds. In pulmonary epithelial cells, both subtypes were found coexistent around the nucleus and the cytoplasm. In these epithelial cells, LPS increased PDE4B and, to a lesser extend, PDE4D, whereas the effect of the inhibitors was prominent on the PDE4B subtype. In conclusion, we determined the pivotal role of the PDE4-inhibitor roflumilast on lung epithelium and emphasized its main effect on PDE4B in hyperinflammation.     

PDE4 inhibitor, roflumilast protects cardiomyocytes against NO-induced apoptosis via activation of PKA and Epac dual pathways

Myocyte apoptosis plays an important role in myocardial infarction and cAMP is crucial in the regulation of myocyte apoptosis. Phosphodiesterase-4 (PDE4) inhibitor blocks the hydrolysis of cAMP via inhibition of PDE4 and is attractive candidate for novel anti-inflammatory drugs. However, its function in cardiovascular diseases and cardiomyocyte apoptosis is unclear. Therefore, we investigated whether roflumilast, a PDE4 inhibitor, exerts protective effect against NO-induced apoptosis in both of H9c2 cells and neonatal rat cardiomyocytes (NRCMs), focusing on cAMP downstream molecules such as protein kinase A (PKA) and exchange protein directly activated by cAMP (Epac). According to our data, intracellular cAMP was increased by roflumilast treatment in H9c2 cells and NRCMs. Roflumilast inhibited SNP-induced apoptosis and this effect was reversed by PKA specific inhibitor H-89 and KT-5720. In addition, PKA specific activator N(6)-benzoyladenosine 3',5-cyclic monophosphate (N(6)Bz-cAMP) mimicked the effects of roflumilast. CREB phosphorylation by roflumilast was also inhibited by H-89, indicating that roflumilast protects SNP-induced apoptosis via PKA-dependent pathway. Roflumilast increased Epac1/GTP-Rap1 and the protective effect was abolished by Epac1 siRNA transfection, demonstrating that Epac signaling was also involved in this protective response. In support, Epac specific activator 8-(4-chlrorophenylthio)-2'-O-methyladenosine-3',5'-cyclic monophosphate (8CPT-2Me-cAMP) protected SNP-induced apoptosis. PI3K/Akt inhibitor LY294002 blocked roflumilast-induced Akt phosphorylation and protective effect. Furthermore, inhibition of Epac1 with siRNA had no effect on roflumilast-induced CREB phosphorylation, whereas inhibited Akt phosphorylation, implicating that Akt phosphorylation was regulated by Epac pathway. In addition, it was also observed that rolipram and cilomilast exert similar effects as roflumilast. In summary, our data indicate that roflumilast protects NO-induced apoptosis via both cAMP-PKA/CREB and Epac/Akt-dependent pathway. Our study suggests a possibility of PDE4 inhibitor roflumilast as a potential therapeutic agent against myocardial ischemia/reperfusion (I/R) injury.     

Computational screening of dipeptidyl peptidase IV inhibitors from micoroalgal metabolites by pharmacophore modeling and molecular docking

Dipeptidyl peptidase IV (DPP‐IV) catalyzes conversion of GLP‐1 (glucagon like peptide 1) to inert structure, which results in insufficient secretion of insulin and increase in postprandial blood glucose level. The present study attempts to identify novel inhibitors from bioactive metabolites present in microalgae against DPP‐IV through virtual screening, molecular docking, and pharmacophore modeling for the active target. Possible binding modes of all 60 ligands against DPP‐IV receptor were constructed using MTiOpenScreen virtual screening server. Pharmacophore model was built based on identified 38 DPP‐IV test ligands by using the web‐based PharmaGist program which encompasses hydrogen‐bond acceptors, hydrophobic groups, spatial features, and aromatic rings. The pharmacophore model having highest scores was selected to screen active DPP‐IV ligands. Highest scoring model was used as a query in ZincPharmer screening. All identified ligands were filtered, based on the Lipinski's rule‐of‐five and were subjected to docking studies. In the process of docking analyses, we considered different bonding modes of one ligand with multiple active cavities of DPP‐IV with the help of AutoDock 4.0. The docking analyses indicate that the bioactive constituents, namely, β‐stigmasterol, barbamide, docosahexaenoic acid, arachidonic acid, and harman showed the best binding energies on DPP‐IV receptor and hydrogen bonding with ASP545, GLY741, TYR754, TYR666, ARG125, TYR547, SER630, and LYS554 residues. This study concludes that docosahexaenoic acid, arachidonic acid, β‐stigmasterol, barbamide, harman, ZINC58564986, ZINC56907325, ZINC69432950, ZINC69431828, ZINC73533041, ZINC84287073, ZINC69849395, and ZINC10508406 act as possible DPP‐IV inhibitors.     

Combined anti-inflammatory effects of β2-adrenergic agonists and PDE4 inhibitors on astrocytes by upregulation of intracellular cAMP

Inflammation is an important hallmark of all neurodegenerative diseases and activation of different glial populations may be involved in the progression of some of these disorders. Especially, the activation of astroglia can lead to long-term detrimental morphological changes, such as scar formation. Therefore, improved strategies to modulate inflammation in these cells are currently being investigated. We investigated the interaction of phosphodiesterase (PDE) 4 inhibitors, such as rolipram, with other agents raising cellular cAMP levels. When used alone, none of the PDE4 inhibitors increased cAMP levels. The adenylate cyclase activator forskolin, the β₂-adrenergic agonist clenbuterol and the mixed β₁/β₂-adrenergic agonist isoproterenol increased intracellular cAMP levels of cortical murine astrocytes. This increase was synergistically elevated by rolipram or the PDE4 inhibitor RO-201724, but not by inhibition of PDE3. Inflammatory stimulation of the cells with the cytokines TNF-α, IL-1β and IFN-γ strongly induced PDE4B and augmented overall PDE4 activity, while PDE3 activity was low. Clenbuterol and forskolin caused downregulation of cytokines and chemokines such as IL-6 and MCP-1. This effect was further enhanced by rolipram, but not by the PDE3 inhibitor milrinone. The cAMP-raising drug combinations attenuated the upregulation of TNF-α and IL-6 mRNA and the secretion of IL-6, but did not affect initial NF-κB signalling triggered by the stimulating cytokines. These results indicate that PDE4 may be a valuable anti-inflammatory target in brain diseases, especially under conditions associated with stimulation of cAMP-augmenting astrocyte receptors as is observed by clenbuterol treatment.     

Investigation of PDE5/PDE6 and PDE5/PDE11 selective potent tadalafil-like PDE5 inhibitors using combination of molecular modeling approaches,molecular fingerprint-based virtual screening protocols and structure-based pharmacophore development

The essential biological function of phosphodiesterase (PDE) type enzymes is to regulate the cytoplasmic levels of intracellular second messengers, 3′,5′-cyclic guanosine monophosphate (cGMP) and/or 3′,5′-cyclic adenosine monophosphate (cAMP). PDE targets have 11 isoenzymes. Of these enzymes, PDE5 has attracted a special attention over the years after its recognition as being the target enzyme in treating erectile dysfunction. Due to the amino acid sequence and the secondary structural similarity of PDE6 and PDE11 with the catalytic domain of PDE5, first-generation PDE5 inhibitors (i.e. sildenafil and vardenafil) are also competitive inhibitors of PDE6 and PDE11. Since the major challenge of designing novel PDE5 inhibitors is to decrease their cross-reactivity with PDE6 and PDE11, in this study, we attempt to identify potent tadalafil-like PDE5 inhibitors that have PDE5/PDE6 and PDE5/PDE11 selectivity. For this aim, the similarity-based virtual screening protocol is applied for the “clean drug-like subset of ZINC database” that contains more than 20 million small compounds. Moreover, molecular dynamics (MD) simulations of selected hits complexed with PDE5 and off-targets were performed in order to get insights for structural and dynamical behaviors of the selected molecules as selective PDE5 inhibitors. Since tadalafil blocks hERG1 K channels in concentration dependent manner, the cardiotoxicity prediction of the hit molecules was also tested. Results of this study can be useful for designing of novel, safe and selective PDE5 inhibitors.     

Virtual screening for novel COX-2 inhibitors using the ZINC database

Cyclooxygenase-2 (COX-2) enzyme binds to arachidonic acid and releases metabolites that are used to induce pain and inflammation. COX-2 selective inhibitors such as celecoxib, rofecoxib and valdecoxib are currently used to reduce inflammatory response. However, they lack anti-thrombotic activity and hence lead to cardiovascular and renal liabilities apart from gastrointestinal irritation. Therefore, there is still a need to develop more potent COX-2 inhibitors. In this paper, we report the screening of various compounds from the ZINC database (contains 4.6 million small molecule compounds) using the eHiTS (electronic High Throughput Screening) software tool against the COX-2 protein. The strategy employed can be conveniently split into two categories, viz. screening and docking, respectively. Screening was performed using molecular constraints tool to filter compounds with physico-chemical properties similar to the 6COX bound ligand SC-558. The analysis resulted in 1042 Lipinski compliant hits which are docked and scored to identify structurally novel ligands that make similar interactions to those of known ligands or may have different interactions with other parts of the binding site. Our screening approach identified two molecules ZINC00663976 (eHITS score of -7.135 kcal/mol) and ZINC02062094 (eHITS score of -7.242 kcal/mol) from the ZINC database. Their energy scores are better than the 6COX bound co-crystallized ligand SC-558 with an eHiTS score of -6.559 kcal/mol. Both the ligands were docked within the binding pocket forming interactions with Leu352, Phe518, Met522, Val523, Ala527 and Ser353. Visual inspection suggested similar orientation and binding mode for ZINC02062094 with SC-558 ligand. The NH group of the ligand formed hydrogen bond interactions with the backbone NH of Ala527.     

Discovery of triazines as selective PDE4B versus PDE4D inhibitors

In this study we report a series of triazine derivatives that are potent inhibitors of PDE4B. We also provide a series of structure activity relationships that demonstrate the triazine core can be used to generate subtype selective inhibitors of PDE4B versus PDE4D. A high resolution co-crystal structure shows that the inhibitors interact with a C-terminal regulatory helix (CR3) locking the enzyme in an inactive ‘closed’ conformation. The results show that the compounds interact with both catalytic domain and CR3 residues. This provides the first structure-based approach to engineer PDE4B-selective inhibitors.     

Virtual screening to identify <Emphasis Type="Italic">Leishmania braziliensis N</Emphasis>-myristoyltransferase inhibitors: pharmacophore models,docking, and molecular dynamics

Leishmaniasis is caused by several protozoa species belonging to genus Leishmania that are hosted by humans and other mammals. Millions of new cases are recorded every year and the drugs available on the market do not show satisfactory efficacy and safety. A hierarchical virtual screening approach based on the pharmacophore model, molecular docking, and molecular dynamics was conducted to identify possible Leishmania braziliensis N-misristoyltransferase (LbNMT) inhibitors. The adopted pharmacophore model had three main features: four hydrophobic centers, four hydrogen-bond acceptor atoms, and one positive nitrogen center. The molecules (n=15,000) were submitted to alignment with the pharmacophore model and only 27 molecules aligned to model. Six molecules were submitted to molecular docking, using receptor PDB ID 5A27. After docking, the ZINC35426134 was a top-ranked molecule (? 64.61 kcal/mol). The molecule ZINC35426134 shows hydrophobic interactions with Phe82, Tyr209, Val370, and Leu391 and hydrogen bonds with Asn159, Tyr318, and Val370. Molecular dynamics simulations were performed with the protein in its APO and HOLO forms for 37 ns in order to assess the stability of the protein–ligand complex. Results showed that the HOLO form was more stable than the APO one, and it suggests that the ZINC35426134 binding stabilizes the enzyme. Therefore, the selected molecule has the potential to meet the herein proposed target.     

In silico docking and molecular dynamics simulation of 3-dehydroquinate synthase (DHQS) from <Emphasis Type="Italic">Mycobacterium tuberculosis</Emphasis>

The shikimate pathway is as an attractive target because it is present in bacteria, algae, fungi, and plants but does not occur in mammals. In Mycobacterium tuberculosis (MTB), the shikimate pathway is integral to the biosynthesis of naphthoquinones, menaquinones, and mycobactin. In these study, novel inhibitors of 3-dehydroquinate synthase (DHQS), an enzyme that catalyzes the second step of the shikimate pathway in MTB, were determined. 12,165 compounds were selected from two public databases through virtual screening and molecular docking analysis using PyRx 8.0 and Autodock 4.2, respectively. A total of 18 compounds with the best binding energies (?13.23 to ?8.22 kcal/mol) were then selected and screened for absorption, distribution, metabolism, excretion, and toxicity (ADMET) analysis, and nine of those compounds were found to satisfy all of the ADME and toxicity criteria. Among those nine, the three compounds—ZINC633887?(binding energy =??10.29 kcal/mol), ZINC08983432?(?9.34 kcal/mol), and PubChem73393?(?8.61 kcal/mol)—with the best binding energies were further selected for molecular dynamics (MD) simulation analysis. The results of the 50-ns MD simulations showed that the two compounds ZINC633887 and PubChem73393 formed stable complexes with DHQS and that the structures of those two ligands remained largely unchanged at the ligand-binding site during the simulations. These two compounds identified through docking and MD simulation are potential candidates for the treatment of TB, and should undergo validation in vivo and in vitro.     

Identification potential inhibitors against the Streptococcus quorum-sensing signal pathway

Abstract

Streptococcal infections are common in human and antibiotics are frequently prescribed in clinical practice. However, infections caused by drug-resistant strains are particularly difficult to treat using common antibiotics. Hence, there is an urgent need for new antibiotics. Quorum sensing is a regulatory mechanism involving cell communication that is thought to play an important role in various bacterial infections, including those caused by Streptococcus. The ATP-binding cassette transporter ComA of Streptococcus is essential for quorum-sensing signal production. The inhibition of the ComA peptidase domain (ComA PEP) suppresses the quorum-sensing pathway and resulting changes in phenotype and/or behavior. Using virtual screening and molecular dynamics simulations, two promising candidate compounds, ZINC32918029 and ZINC6751571, were found. These compounds had similar binding modes and interactions to the experimentally determined reference inhibitor 6CH. However, a significantly stronger negative binding energy was achieved (?113.501?±?15.312?KJ/mol and ?103.153?±?11.912?KJ/mol for ZINC32918029 and ZINC6751571, respectively). Molecular dynamics simulations also revealed that ZINC32918029 and ZINC6751571 had a strong affinity for ComA PEP. These results indicate that ZINC32918029 and ZINC6751571 are promising candidate inhibitors of the Streptococcus quorum-sensing pathway.

Communicated by Ramaswamy H. Sarma