Molecular recognition of receptor sites using a genetic algorithm with a description of desolvation

Integrated ESI-MS/EPR/computational characterization of the binding of metal species to proteins: vanadium drug–myoglobin application

An integrated strategy based on ESI-MS spectrometry, EPR spectroscopy and docking/QM computational methods is applied to the systems formed by V^IVO²⁺ ions and four potential V^IVOL₂ drugs and myoglobin. This approach is generizable to other metals and proteins.

Methylation Products of 6β- N-Heterocyclic Substituted Naltrexamine Derivatives as Potential Peripheral Opioid Receptor Modulators

Two 6β- N-heterocyclic naltrexamine derivatives, NAP and NMP, have been identified as peripherally selective mu opioid receptor (MOR) antagonists. To further enhance the peripheral selectivity of both compounds, the 17-amino group and the nitrogen atom of the pyridine ring in both NAP and NMP were methylated to obtain dMNAP and dMNMP, respectively. Compared with NAP and NMP, the binding affinities of dMNAP and dMNMP shifted to MOR and KOR (kappa opioid receptor) dual selective and they acted as moderate efficacy partial agonists. The results from radioligand binding studies were further confirmed by molecular docking studies. In vivo studies demonstrated that dMNAP and dMNMP did not produce antinociception nor did they antagonize morphine's antinociceptive activity, indicating that these compounds did not act on the central nervous system. Meanwhile, both dMNAP and dMNMP significantly slowed down fecal excretion, which indicated that they were peripherally acting opioid receptor agonists. All together, these results suggested that dMNAP and dMNMP acted as peripheral mu/kappa opioid receptor modulators and may be applicable in the treatment of diarrhea in patients with bowel dysfunction.

Modeling and Simulation of hGAT1: A Mechanistic Investigation of the GABA Transport Process

Human γ-Aminobutyric acid transporter 1 (hGAT1) is a Na+/Cl- dependent co-transporter that plays a key role in the inhibitory neurotransmission of GABA in the brain. Due to the lack of structural data, the exact co-transport mechanism of GABA reuptake by hGAT1 remains unclear. To examine the roles of the co-transport ions and the nature of their interactions with GABA, homology modeling and molecular dynamics simulations of the hGAT1 in the open-to-out conformation were carried out. Our study focused on the sequential preloading of Na+ and Cl- ions, followed by GABA binding. Our simulations showed pre-loading of ions maintains the transport ready state of hGAT1 in the open-to-out conformation essential for GABA binding. Of the four putative preloaded states, GABA binding to the fully loaded state is most favored. Binding of Na+ ion to the Na1 site helps to maintain the open-to-out conformation for GABA binding as compared to the Na2 site. GABA binding to the mono-sodium or the di-sodium loaded states leads to destabilization of Na+ ions within their binding sites. The two most prominent interactions required for GABA binding include interaction between carboxylate group of GABA with the bound Na+ ion in Na1 binding site and the hydroxyl group of Y140. Overall our results support the fully loaded state as the predominate state for GABA binding. Our study further illustrates that Na+ ion within the Na1 site is crucial for GABA recognition. Therefore, a revised mechanism is proposed for the initial step of hGAT1 translocation cycle.

Molecular and computational approaches to understand resistance of New Delhi metallo β-lactamase variants (NDM-1, NDM-4, NDM-5, NDM-6, NDM-7)-producing strains against carbapenems

The discovery of NDM-1 and its variants has caused the emergence of antibiotic resistance in the community and hospital setting, causing major concern for health care across the globe. New Delhi Metallo-β-lactamase is known to hydrolyse almost all β-lactam antibiotics. Studies have shown the hydrolytic activates of NDM-1 and some of its variants, however a comparative study of these NDM variants has not been explored in detail. Hence, we proposed to check their catalytic activity by performing a comparative study between NDM-1 and its variants. The study was initiated to clone NDM variants (NDM-1, NDM-4, NDM-5, NDM-6 and NDM-7) followed by overexpression of the recombinant proteins to check their hydrolytic properties against β-lactam antibiotics. The minimum inhibitory concentration of carbapenems antibiotics for bla_NDM-5 clone was found four fold increased, whereas no change was observed in the clones having other variants. The hydrolytic activity of carbapenem with NDM-5 variant was found to be augmented as per the kinetics parameter where Km was decreased and kcat, kcat/Km values increased as compared to the NDM-1. Molecular docking studies were employed to identify the variations in the binding ability among all NDM variants with imipenem or meropenem. Simulation studies at 100 ns showed a good stability of NDM-5 with imipenem and meropenem as compared to NDM-1. CD spectroscopy data revealed significant changes in the secondary structure of NDM variants. We conclude that NDM-5 showed higher hydrolytic activity as compared to other variants. This study provides a comparative analysis of the severity of NDM producing strains.

Stereoselective Synthesis of the Isomers of Notoincisol A: Assigment of the Absolute Configuration of this Natural Product and Biological Evaluation

The total syntheses of all stereoisomers of notoincisol A, a recently isolated natural product with potential anti-inflammatory activity, are reported. The asymmetric synthesis was conducted employing a lipase-mediated kinetic resolution, which enables easy access to all required chiral building blocks with the aim of establishing the absolute configuration of the naturally occurring isomer. This was achieved by comparison of optical properties of the isolated compound with the synthetic derivatives obtained. Moreover, an assessment of the biological activity on PPARγ (peroxisome proliferator-activated receptor gamma) as a prominent receptor related to inflammation is reported. Only the natural isomer was found to activate the PPARγ receptor, and this phenomenon could be explained based on molecular docking studies. In addition, the pharmacological profiles of the isomers were determined using the GABA_A (gamma-aminobutyric acid A) ion channel receptor as a representative target for allosteric modulation related to diverse CNS activities. These compounds were found to be weak allosteric modulators of the α1β3 and α1β2γ2 receptor subtypes.

Validation and Applications of Protein-Ligand Docking Approaches Improved for Metalloligands with Multiple Vacant Sites

Decoding the interaction between coordination compounds and proteins is of fundamental importance in biology, pharmacy, and medicine. In this context, protein- ligand docking represents a particularly interesting asset to predict how small compounds could interact with biomolecules, but to date, very little information is available to adapt these methodologies to metal-containing ligands. Here, we assessed the predictive capability of a metal-compatible parameter set for the docking program GOLD for metallo ligands with multiple vacant sites and different geometries. The study first presents a benchmark of 25 well-characterized X-ray metallo ligand-protein adducts. In 100% of the cases, the docking solutions are superimposable to the X-ray determination, and in 92% the value of the root-mean-square deviation between the experimental and calculated structures is lower than 1.5 Å. After the validation step, we applied these methods to five case studies for the prediction of the binding of pharmacological active metal species to proteins: (i) the anticancer copper(II) complex [Cu^II(Br)(2-hydroxy-1-naphthaldehyde benzoyl hydrazine)(indazole)] to human serum albumin (HSA); (ii) one of the active species of antidiabetic and antitumor vanadium compounds, V^IVO²⁺ ion, to carboxypeptidase; (iii) the antiarthritic species [Au^I(PEt₃)]⁺ to HSA; (iv) the antitumor oxaliplatin to ubiquitin; (v) the antitumor ruthenium(II) compound RAPTA-PentaOH to cathepsin B. The calculations suggested that the binding modes are in good agreement with the partial information retrieved from spectroscopic and spectrometric analysis and allowed us, in certain cases, to propose additional hypotheses. This method is an important update in protein-metallo ligand docking, which could have a wide field of application, from biology and inorganic biochemistry to medicinal chemistry and pharmacology.

Nintedanib-cyclodextrin complex to improve bio-activity and intestinal permeability

Cyclodextrin complex of nintedanib was prepared aiming for increased bio-activity and improved transport across intestinal membrane with reduced p-glycoprotein (p-gp) efflux. Based on preliminary phase solubility studies and molecular modeling, sulfobutyl ether derivative of β-cyclodextrin (SBE-β-CD, Captisol^®) was selected to prepare inclusion complex. Complexation was confirmed using FTIR, 1H NMR, DSC, and XRD. Bioactivity of the formed complex was tested using lung fibroblast cells, WI-38 for anti-proliferative activity and effect on collagen deposition and cells migration. In-vitro permeability studies were performed using epiIntestinal tissue model to assess the effect of complexation on transport and p-gp efflux. Results of the study demonstrated that cyclodextrin complexation increased stability of nintedanib in PBS (pH 7.4) and simulated intestinal fluid (SIF). Further, bioactivity of nintedanib also improved. Interestingly, complexation has increased transport of nintedanib across intestinal membrane and reduced efflux ratio, suggesting the role of cyclodextrin complexation in modulating p-gp efflux.

Structure-guided design of pyridoclax derivatives based on Noxa / Mcl-1 interaction mode

Protein-protein interactions are attractive targets because they control numerous cellular processes. In oncology, apoptosis regulating Bcl-2 family proteins are of particular interest. Apoptotic cell death is controlled via PPIs between the anti-apoptotic proteins hydrophobic groove and the pro-apoptotic proteins BH3 domain. In ovarian carcinoma, it has been previously demonstrated that Bcl-x_L and Mcl-1 cooperate to protect tumor cells against apoptosis. Moreover, Mcl-1 is a key regulator of cancer cell survival and is a known resistance factor to Bcl-2/Bcl-x_L pharmacological inhibitors making it an attractive therapeutic target. Here, using a structure-guided design from the oligopyridine lead Pyridoclax based on Noxa/Mcl-1 interaction we identified a new derivative, active at lower concentration as compared to Pyridoclax. This new derivative selectively binds to the Mcl-1 hydrophobic groove and releases Bak and Bim from Mcl-1 to induce cell death and sensitize cancer cells to Bcl-2/Bcl-x_L targeting strategies.

Synthesis, molecular modeling and biological screening of some pyrazole derivatives as antileishmanial agents

Aim: Novel open chain and cyclized derivatives containing pyrazole scaffold were designed, synthesized and evaluated as antileishmanial compounds. Methodology & results: In silico reverse docking experiment suggested Leishmania major pteridine reductase (Lm-PTR1) as a putative target for the synthesized compounds. In vitro antileishmanial screening against L. major promastigotes and amastigotes using miltefosine and amphotericin B as references showed that the majority of the compounds displayed activity higher than miltefosine. Compounds 3i and 5 showed the highest antileishmanial activity with IC50 values of 1.45 ± 0.08 μM and 2.30 ± 0.09 μM, respectively, for the amastigote form. In silico drug-likeness and toxicity predictions showed acceptable profiles for most of the compounds, which were validated by experimental toxicity studies. Conclusion: This study offers promising entities for antileishmanial activity.

Anti-leprosy drug Clofazimine binds to human Raf1 kinase inhibitory protein and enhances ERK phosphorylation

Human Raf1 kinase inhibitory protein (hRKIP) is an important modulator of the Ras/Raf1/MEK/ERK signaling pathway. Here, we demonstrated that anti-leprosy drug Clofazimine can bind to hRKIP with a significantly stronger affinity than the endogenous substrate phosphatidylethanolamine (PE) by using Biolayer interference technology. Moreover, we identified that residues P74, S75, K80, P111, P112, V177, and P178 play crucial roles in the binding of hRKIP to Clofazimine by using a combination of Nuclear Magnetic Resonance spectroscopy and molecular docking approach. These residues are located at the conserved ligand-binding pocket of hRKIP. Furthermore, we found that 3.2 μM Clofazimine could significantly increase the ERK phosphorylation level by about 37%. Our results indicate that Clofazimine can enhance Ras/Raf1/MEK/ERK signaling transduction pathway via binding to hRKIP. This work provides valuable hints for exploiting Clofazimine as a potential lead compound to efficiently treat the diseases related to RKIP or the Ras/Raf/MEK/ERK pathway.

Dihydrochalcone Glucosides from the Subaerial Parts of Thonningia sanguinea and Their in Vitro PTP1B Inhibitory Activities

Six new and four known dihydrochalcone glucoside derivatives (1-10), the phenylpropanoid coniferin (11), and the lignans (+)-pinoresinol (12) and lariciresinol (13) were isolated from the subaerial plant parts of Thonningia sanguinea in the course of a screening campaign for new antidiabetic lead compounds. The structures of the new substances were elucidated by HRESIMS, NMR, GC-MS, and ECD data evaluation. 2'- O-(3-Galloyl-4,6- O- S_a-hexahydroxydiphenoyl-β-d-glucopyranosyl)-3-hydroxyphloretin (4), 2'- O-(4,6- O- S_a-hexahydroxydiphenoyl-β-d-glucopyranosyl)phloretin (5), 2'- O-(3- O-galloyl-4,6- O- S_a-hexahydroxydiphenoyl-β-d-glucopyranosyl)phloretin (6), and thonningianin B (9) showed moderate protein tyrosine phosphatase-1B inhibition in an enzyme assay (IC₅₀ values ranging from 19 to 25 μM), whereas thonningianin A (10) was identified as a more potent inhibitor (IC₅₀ = 4.4 μM). The observed activity differences could be explained by molecular docking experiments. The activity of 10 could further be confirmed in HEPG2 liver carcinoma cells, where the compound was able to increase the level of phosphorylated insulin receptors in a concentration-dependent manner.

An Evaluation of the Potential of NMR Spectroscopy and Computational Modelling Methods to Inform Biopharmaceutical Formulations

Protein-based therapeutics are considered to be one of the most important classes of pharmaceuticals on the market. The growing need to prolong stability of high protein concentrations in liquid form has proven to be challenging. Therefore, significant effort is being made to design formulations which can enable the storage of these highly concentrated protein therapies for up to 2 years. Currently, the excipient selection approach involves empirical high-throughput screening, but does not reveal details on aggregation mechanisms or the molecular-level effects of the formulations under storage conditions. Computational modelling approaches have the potential to elucidate such mechanisms, and rapidly screen in silico prior to experimental testing. Nuclear Magnetic Resonance (NMR) spectroscopy can also provide complementary insights into excipient–protein interactions. This review will highlight the underpinning principles of molecular modelling and NMR spectroscopy. It will also discuss the advancements in the applications of computational and NMR approaches in investigating excipient–protein interactions.

An update on the discovery and development of selective heat shock protein inhibitors as anti-cancer therapy

INTRODUCTION

Over the years, not a single HSP inhibitor has progressed into the post-market phase of drug development despite the success recorded in various pre-clinical and clinical studies. The inability of existing drugs to specifically target oncogenic HSPs has majorly accounted for these setbacks. Recent combinatorial strategies that incorporated computer-aided drug design (CADD) techniques are geared towards the development of highly specific HSP inhibitors with increased activities and minimal toxicities. Areas covered: In this review, strategic therapeutic approaches that have recently aided the development of selective HSP inhibitors were highlighted. Also, the significant contributions of CADD techniques over the years were discussed in detail. This article further describes promising computational paradigms and their applications towards the discovery of highly specific inhibitors of oncogenic HSPs. Expert opinion: The recent shift towards highly selective and specific HSP inhibition has shown great promise as evidenced by the development of paralog/isoform-selective HSP drugs. It could be further augmented with computer-aided drug design strategies, which incorporate reliable methods that would greatly enhance the design and optimization of novel inhibitors with improved activities and minimal toxicities.

Molecular Analysis and Modeling of Hepatitis E Virus Helicase and Identification of Novel Inhibitors by Virtual Screening

The hepatitis E virus- (HEV-) helicase as a novel drug-target was evaluated. While cell culture model was used for mutational characterization of helicase, in silico protein modeling and virtual screening were employed to identify helicase inhibitors. None of the saturation mutant replicons significantly affected RNA replication. Notably, mutants encompassing the Walker motifs replicated as wild-type, showing indispensability of nucleotides conservation in viability compared to known criticality of amino acids. A 3D modeling of HEV-helicase and screening of a compound dataset identified ten most promising inhibitors with drug likeness, notably, JFD02650, RDR03130, and HTS11136 that interacted with Walker A residues Gly975, Gly978, Ser979, and Gly980. Our model building and virtual identification of novel helicase inhibitors warrant further studies towards developing anti-HEV drugs.

Bridging Molecular Docking to Molecular Dynamics in Exploring Ligand-Protein Recognition Process: An Overview

Computational techniques have been applied in the drug discovery pipeline since the 1980s. Given the low computational resources of the time, the first molecular modeling strategies relied on a rigid view of the ligand-target binding process. During the years, the evolution of hardware technologies has gradually allowed simulating the dynamic nature of the binding event. In this work, we present an overview of the evolution of structure-based drug discovery techniques in the study of ligand-target recognition phenomenon, going from the static molecular docking toward enhanced molecular dynamics strategies.

Elastin is heterogeneously cross-linked

Elastin is an essential vertebrate protein responsible for the elasticity of force-bearing tissues such as those of the lungs, blood vessels, and skin. One of the key features required for the exceptional properties of this durable biopolymer is the extensive covalent cross-linking between domains of its monomer molecule tropoelastin. To date, elastin's exact molecular assembly and mechanical properties are poorly understood. Here, using bovine elastin, we investigated the different types of cross-links in mature elastin to gain insight into its structure. We purified and proteolytically cleaved elastin from a single tissue sample into soluble cross-linked and noncross-linked peptides that we studied by high-resolution MS. This analysis enabled the elucidation of cross-links and other elastin modifications. We found that the lysine residues within the tropoelastin sequence were simultaneously unmodified and involved in various types of cross-links with different other domains. The Lys-Pro domains were almost exclusively linked via lysinonorleucine, whereas Lys-Ala domains were found to be cross-linked via lysinonorleucine, allysine aldol, and desmosine. Unexpectedly, we identified a high number of intramolecular cross-links between lysine residues in close proximity. In summary, we show on the molecular level that elastin formation involves random cross-linking of tropoelastin monomers resulting in an unordered network, an unexpected finding compared with previous assumptions of an overall beaded structure.

Synthesis, anti-inflammatory screening, molecular docking, and COX-1,2/-5-LOX inhibition profile of some novel quinoline derivatives

New quinoline compounds comprising pyrazole scaffold through different amide linkages were synthesized. The synthesized compounds were evaluated for their anti-inflammatory activity. Eight compounds (5c, 11b,c, 12c, 14a,b, 20a and 21a) were found to exhibit promising anti-inflammatory profiles in acute and sub-acute inflammatory models. They were screened for their ulcerogenic activity and none of them showed significant ulcerogenic activity comparable to the reference drug celecoxib and are well tolerated by experimental animals with high safety margin (ALD50 > 0.3 g/kg). Compounds 5c, 11b,c, 12c, 14a,b, 20a and 21a showed significant in vitro LOX inhibitory activity higher than that of zileuton. In vitro COX-1/COX-2 inhibition study revealed that compounds 12c, 14a,b and 20a showed higher selectivity towards COX-2 than COX-1. Among the tested compounds, 12c, 14a and 14b showed the highest inhibitory activity against COX-2 with an IC50 values of 0.1, 0.11 and 0.11 μM respectively. The docking experiments attempted to postulate the binding mode for the most active compounds in the binding site of COX-2 enzymes and confirmed the high selectivity binding towards COX-2 enzyme over COX-1.

Importance of the second extracellular loop for melatonin MT1 receptor function and absence of melatonin binding in GPR50

BACKGROUND AND PURPOSE

Recent crystal structures of GPCRs have emphasized the previously unappreciated role of the second extracellular (E2) loop in ligand binding and gating and receptor activation. Here, we have assessed the role of the E2 loop in the activation of the melatonin MT1 receptor and in the inactivation of the closely related orphan receptor GPR50.

EXPERIMENTAL APPROACH

Chimeric MT1 -GPR50 receptors were generated and functionally analysed in terms of 2-[125 I]iodomelatonin binding, Gi /cAMP signalling and β-arrestin2 recruitment. We also used computational molecular dynamics (MD) simulations.

KEY RESULTS

MD simulations of 300 ns revealed (i) the tight hairpin structure of the E2 loop of the MT1 receptor (ii) the most suitable features for melatonin binding in MT1 receptors and (iii) major predicted rearrangements upon MT1 receptor activation, stabilizing interaction networks between Phe179 or Gln181 in the E2 loop and transmembrane helixes 5 and 6. Functional assays confirmed these predictions, because reciprocal replacement of MT1 and GPR50 residues/domains led to the predicted loss- and gain-of-melatonin action of MT1 receptors and GPR50 respectively.

CONCLUSIONS AND IMPLICATIONS

Our work demonstrated the crucial role of the E2 loop for MT1 receptor and GPR50 function by proposing a model in which the E2 loop is important in stabilizing active MT1 receptor conformations and by showing how evolutionary processes appear to have selected for modifications in the E2 loop in order to make GPR50 unresponsive to melatonin.

LINKED ARTICLES

This article is part of a themed section on Recent Developments in Research of Melatonin and its Potential Therapeutic Applications. To view the other articles in this section visit http://onlinelibrary.wiley.com/doi/10.1111/bph.v175.16/issuetoc.

Targeting Pin1 by inhibitor API-1 regulates microRNA biogenesis and suppresses hepatocellular carcinoma development

Hepatocellular carcinoma (HCC) is a leading cause of cancer death worldwide, but there are few effective treatments. Aberrant microRNA (miRNA) biogenesis is correlated with HCC development. We previously demonstrated that peptidyl-prolyl cis-trans isomerase NIMA-interacting 1 (Pin1) participates in miRNA biogenesis and is a potential HCC treatment target. However, how Pin1 modulates miRNA biogenesis remains obscure. Here, we present in vivo evidence that Pin1 overexpression is directly linked to the development of HCC. Administration with the Pin1 inhibitor (API-1), a specific small molecule targeting Pin1 peptidyl-prolyl isomerase domain and inhibiting Pin1 cis-trans isomerizing activity, suppresses in vitro cell proliferation and migration of HCC cells. But API-1-induced Pin1 inhibition is insensitive to HCC cells with low Pin1 expression and/or low exportin-5 (XPO5) phosphorylation. Mechanistically, Pin1 recognizes and isomerizes the phosphorylated serine-proline motif of phosphorylated XPO5 and passivates phosphorylated XPO5. Pin1 inhibition by API-1 maintains the active conformation of phosphorylated XPO5 and restores XPO5-driven precursor miRNA nuclear-to-cytoplasm export, activating anticancer miRNA biogenesis and leading to both in vitro HCC suppression and HCC suppression in xenograft mice.

CONCLUSION

Experimental evidence suggests that Pin1 inhibition by API-1 up-regulates miRNA biogenesis by retaining active XPO5 conformation and suppresses HCC development, revealing the mechanism of Pin1-mediated miRNA biogenesis and unequivocally supporting API-1 as a drug candidate for HCC therapy, especially for Pin1-overexpressing, extracellular signal-regulated kinase-activated HCC. (Hepatology 2018).

COX Inhibition Profiles and Molecular Docking Studies of the Lignan Hinokinin and Some Synthetic Derivatives

Encouraged by the anti-inflammatory activity of hinokinin in vivo, which is also observed for the analogues dinitrohinokinin and diidrocubebin, herein we used in vitro and in silico methods to assess their selectivity profiles and predict their binding modes with Cyclooxygenases (COX-1 and 2). The in vitro assays demonstrated dinitrohinokinin is about 13 times more selective for COX-2 than for COX-1, a similar profile observed for the drugs celecoxib (selective index ≈9) and meloxicam (selective index ≈11). Predictions of the binding modes suggested dinitrohinokinin interacts with COX-2 very similarly to rofecoxib, exploring residues at the hydrophilic pocket of the enzyme that accessible to ligands only in this isoform. This lignan also interacts with COX-1 in a similar mode to meloxicam, blocking the access of the substrate to the catalytic cleft. Therefore, dinitrohinokinin is a promising lead for the design of selective COX-2 inhibitors.

Binding Affinity via Docking: Fact and Fiction

In 1982, Kuntz et al. published an article with the title “A Geometric Approach to Macromolecule-Ligand Interactions”, where they described a method “to explore geometrically feasible alignment of ligands and receptors of known structure”. Since then, small molecule docking has been employed as a fast way to estimate the binding pose of a given compound within a specific target protein and also to predict binding affinity. Remarkably, the first docking method suggested by Kuntz and colleagues aimed to predict binding poses but very little was specified about binding affinity. This raises the question as to whether docking is the right tool to estimate binding affinity. The short answer is no, and this has been concluded in several comprehensive analyses. However, in this opinion paper we discuss several critical aspects that need to be reconsidered before a reliable binding affinity prediction through docking is realistic. These are not the only issues that need to be considered, but they are perhaps the most critical ones. We also consider that in spite of the huge efforts to enhance scoring functions, the accuracy of binding affinity predictions is perhaps only as good as it was 10–20 years ago. There are several underlying reasons for this poor performance and these are analyzed. In particular, we focus on the role of the solvent (water), the poor description of H-bonding and the lack of the systems’ true dynamics. We hope to provide readers with potential insights and tools to overcome the challenging issues related to binding affinity prediction via docking.

Machine intelligence decrypts β-lapachone as an allosteric 5-lipoxygenase inhibitor

Using machine learning, targets were identified for β-lapachone.

Using machine learning, targets were identified for β-lapachone. Resorting to biochemical assays, β-lapachone was validated as a potent, ligand efficient, allosteric and reversible modulator of 5-lipoxygenase (5-LO). Moreover, we provide a rationale for 5-LO modulation and show that inhibition of 5-LO is relevant for the anticancer activity of β-lapachone. This work demonstrates the power of machine intelligence to deconvolute complex phenotypes, as an alternative and/or complement to chemoproteomics and as a viable general approach for systems pharmacology studies.

Significance of catalase-peroxidase (KatG) mutations in mediating isoniazid resistance in clinical strains of Mycobacterium tuberculosis

OBJECTIVES

Isoniazid (INH) is still the most important first-line antitubercular drug. INH resistance is regarded as a major impediment to the tuberculosis (TB) control programme and contributes to the emergence of multidrug-resistant strains. Mutation at position 315 in the katG gene, encoding the catalase-peroxidase (KatG) enzyme, is the major cause of INH resistance in Mycobacterium tuberculosis. Therefore, investigation of the molecular mechanisms of INH resistance is the need of the hour.

METHODS

To understand the clinical importance of KatG mutants (MTs) leading to INH resistance, in this study five MTs (S315T, S315I, S315R, S315N and S315G) were modelled, docked and interacted with INH in dynamic state.

RESULTS

The binding affinity based on docking was found to be higher for MTs than for wild-type (WT) isolates, except for MT-S315R, indicating rigid binding of INH with MT proteins compared with the flexible binding seen in the WT. Analysis of molecular dynamics (MD) experiments suggested that fluctuations and deviations were higher at the INH binding residues for MTs than for the WT. Reduction in the hydrogen bond network after MD in all KatG enzymes implies an increase in the flexibility and stability of protein structures. Superimposition of MTs upon the WT structure showed a significant deviation that varies for the different MTs.

CONCLUSIONS

It can be inferred that the five KatG MTs affect enzyme activity in different ways, which could be attributed to conformational changes in MT KatG that result in altered binding affinity to INH and eventually to INH resistance.

Multifunctional cholinesterase inhibitors for Alzheimer's disease: Synthesis, biological evaluations, and docking studies of o/p-propoxyphenylsubstituted-1H-benzimidazole derivatives

This study indicates the synthesis, cholinesterase (ChE) inhibitory activity, and molecular modeling studies of 48 compounds as o- and p-(3-substitutedethoxyphenyl)-1H-benzimidazole derivatives. According to the ChE inhibitor activity results, generally, para series are more active against acetylcholinesterase (AChE) whereas ortho series are more active against butyrylcholinesterase (BuChE). The most active compounds against AChE and BuChE are compounds A12 and B14 with IC50 values of 0.14 and 0.22 μM, respectively. Additionally, the most active 16 compounds against AChE/BuChE were chosen to investigate the neuroprotective effects, and the results indicated that most of the compounds have free radical scavenging properties and show their effects by reducing free radical production; moreover, some of the compounds significantly increased the viability of SH-SY5Y cells exposed to H2 O2 . Overall, compounds A12 and B14 with potential AChE and BuChE inhibitory activities, high neuroprotection against H2 O2 -induced toxicity, free radical scavenging properties, and metal chelating abilities may be considered as lead molecules for the development of multi-target-directed ligands against Alzheimer's disease.

Docking of Covalent Ligands: Challenges and Approaches

Covalent ligands have recently regained considerable attention in drug discovery. The rational design of such ligands, however, is still faced with particular challenges, mostly related to the fact that covalent bond formation is a quantum mechanical phenomenon which cannot adequately be handled by the force fields or empirical approaches typically used for noncovalent protein-ligand interactions. Although the necessity for quantum chemical approaches is clear, they cannot yet routinely be applied on large data sets of ligands or for a broader exploration of binding modes in docking calculations. On the other hand, technical solutions for performing docking calculations with covalent ligands are available, but their scope is normally quite limited. Scoring functions typically neglect the contribution from covalent bond formation completely. In this situation, the question arises how to approach covalent ligands and which methods to choose for their docking and design.

Molecular principles for heparin oligosaccharide-based inhibition of neutrophil elastase in cystic fibrosis

Cystic fibrosis (CF) is a multifactorial disease in which dysfunction of protease-antiprotease balance plays a key role. The current CF therapy relies on dornase α, hypertonic saline, and antibiotics and does not address the high neutrophil elastase (NE) activity observed in the lung and sputum of CF patients. Our hypothesis is that variants of heparin, which potently inhibit NE but are not anticoagulant, would help restore the protease-antiprotease balance in CF. To realize this concept, we studied molecular principles governing the effectiveness of different heparins, especially 2-O,3-O-desulfated heparin (ODSH), in the presence of sputum components and therapeutic agents. Using sputa from CF patients and an NE activity assay, we found that heparins are ineffective if used in the absence of dornase. This is true even when mucolytics, such as DTT or N-acetylcysteine, were used. Computational modeling suggested that ODSH and DNA compete for binding to an overlapping allosteric site on NE, which reduces the anti-NE potential of ODSH. NE inhibition of both DNA and ODSH is chain length-dependent, but ODSH chains exhibit higher potency per unit residue length. Likewise, ODSH chains exhibit higher NE inhibition potential compared with DNA chains in the presence of saline. These studies suggest fundamental differences in DNA and ODSH recognition and inhibition of NE despite engaging overlapping sites and offer unique insights into molecular principles that could be used in developing antiprotease agents in the presence of current treatments, such as dornase and hypertonic saline.

A Diverse Benchmark Based on 3D Matched Molecular Pairs for Validating Scoring Functions

The prediction of protein-ligand interactions and their corresponding binding free energy is a challenging task in structure-based drug design and related applications. Docking and scoring is broadly used to propose the binding mode and underlying interactions as well as to provide a measure for ligand affinity or differentiate between active and inactive ligands. Various studies have revealed that most docking software packages reliably predict the binding mode, although scoring remains a challenge. Here, a diverse benchmark data set of 99 matched molecular pairs (3D-MMPs) with experimentally determined X-ray structures and corresponding binding affinities is introduced. This data set was used to study the predictive power of 13 commonly used scoring functions to demonstrate the applicability of the 3D-MMP data set as a valuable tool for benchmarking scoring functions.

Hierarchical Flexible Peptide Docking by Conformer Generation and Ensemble Docking of Peptides

Given the importance of peptide-mediated protein interactions in cellular processes, protein-peptide docking has received increasing attention. Here, we have developed a Hierarchical flexible Peptide Docking approach through fast generation and ensemble docking of peptide conformations, which is referred to as HPepDock. Tested on the LEADS-PEP benchmark data set of 53 diverse complexes with peptides of 3-12 residues, HPepDock performed significantly better than the 11 docking protocols of five small-molecule docking programs (DOCK, AutoDock, AutoDock Vina, Surflex, and GOLD) in predicting near-native binding conformations. HPepDock was also evaluated on the 19 bound/unbound and 10 unbound/unbound protein-peptide complexes of the Glide SP-PEP benchmark and showed an overall better performance than Glide SP-PEP+MM-GBSA and FlexPepDock in both bound and unbound docking. HPepDock is computationally efficient, and the average running time for docking a peptide is ∼15 min with the range from about 1 min for short peptides to around 40 min for long peptides.

Subtractive Genomics, Molecular Docking and Molecular Dynamics Simulation Revealed LpxC as a Potential Drug Target Against Multi-Drug Resistant Klebsiella pneumoniae

The emergence and dissemination of pan drug resistant clones of Klebsiella pneumoniae are great threat to public health. In this regard new therapeutic targets must be highlighted to pave the path for novel drug discovery and development. Subtractive proteomic pipeline brought forth UDP-3-O-[3-hydroxymyristoyl] N-acetylglucosamine deacetylase (LpxC), a Zn⁺² dependent cytoplasmic metalloprotein and catalyze the rate limiting deacetylation step of lipid A biosynthesis pathway. Primary sequence analysis followed by 3-dimensional (3-D) structure elucidation of the protein led to the detection of K. pneumoniae LpxC (KpLpxC) topology distinct from its orthologous counterparts in other bacterial species. Molecular docking study of the protein recognized receptor antagonist compound 106, a uridine-based LpxC inhibitory compound, as a ligand best able to fit the binding pocket with a Gold Score of 67.53. Molecular dynamics simulation of docked KpLpxC revealed an alternate binding pattern of ligand in the active site. The ligand tail exhibited preferred binding to the domain I residues as opposed to the substrate binding hydrophobic channel of subdomain II, usually targeted by inhibitory compounds. Comparison with the undocked KpLpxC system demonstrated ligand induced high conformational changes in the hydrophobic channel of subdomain II in KpLpxC. Hence, ligand exerted its inhibitory potential by rendering the channel unstable for substrate binding.

Symmetric dimeric adamantanes for exploring the structure of two viroporins: influenza virus M2 and hepatitis C virus p7

Background

Adamantane-based compounds have been identified to interfere with the ion-channel activity of viroporins and thereby inhibit viral infection. To better understand the difference in the inhibition mechanism of viroporins, we synthesized symmetric dimeric adamantane analogs of various alkyl-spacer lengths.

Methods

Symmetric dimeric adamantane derivatives were synthesized where two amantadine or rimantadine molecules were linked by various alkyl-spacers. The inhibitory activity of the compounds was studied on two viroporins: the influenza virus M2 protein, expressed in Xenopus oocytes, using the two-electrode voltage-clamp technique, and the hepatitis C virus (HCV) p7 channels for five different genotypes (1a, 1b, 2a, 3a, and 4a) expressed in HEK293 cells using whole-cell patch-clamp recording techniques.

Results

Upon testing on M2 protein, dimeric compounds showed significantly lower inhibitory activity relative to the monomeric amantadine. The lack of channel blockage of the dimeric amantadine and rimantadine analogs against M2 wild type and M2-S31N mutant was consistent with previously proposed drug-binding mechanisms and further confirmed that the pore-binding model is the pharmacologically relevant drug-binding model. On the other hand, these dimers showed similar potency to their respective monomeric analogs when tested on p7 protein in HCV genotypes 1a, 1b, and 4a while being 700-fold and 150-fold more potent than amantadine in genotypes 2a and 3a, respectively. An amino group appears to be important for inhibiting the ion-channel activity of p7 protein in genotype 2a, while its importance was minimal in all other genotypes.

Conclusion

Symmetric dimeric adamantanes can be considered a prospective class of p7 inhibitors that are able to address the differences in adamantane sensitivity among the various genotypes of HCV.

Discovery of novel dual acetylcholinesterase inhibitors with antifibrillogenic activity related to Alzheimer's disease

AIM

Alzheimer's disease is a progressive and neurodegenerative disorder of the CNS, affecting elderly people. The current pharmacological approach is based on the improvement of cholinergic neurotransmission by inhibiting acetylcholinesterase (AChE) with AChE inhibitors. The disease is also characterized by the accelerated accumulation of β-amyloid plaques around neurons. Furthermore, in vitro studies revealed that AChE can induce β-amyloid peptide (Aβ) aggregation.

METHODOLOGY

Computer-aided molecular design by virtual screening was here employed to discover novel potential AChE inhibitors, with antifibrillogenic properties, in other words, inhibiting Aβ aggregation.

RESULTS

Compounds 1, 4 and 6 showed interesting AChE inhibition. In addition, they particularly inhibit Aβ aggregation in vitro, indicating to be promising novel anti-Alzheimer agents.

Computational Fragment-Based Drug Design: Current Trends, Strategies, and Applications

Fragment-based drug design (FBDD) has become an effective methodology for drug development for decades. Successful applications of this strategy brought both opportunities and challenges to the field of Pharmaceutical Science. Recent progress in the computational fragment-based drug design provide an additional approach for future research in a time- and labor-efficient manner. Combining multiple in silico methodologies, computational FBDD possesses flexibilities on fragment library selection, protein model generation, and fragments/compounds docking mode prediction. These characteristics provide computational FBDD superiority in designing novel and potential compounds for a certain target. The purpose of this review is to discuss the latest advances, ranging from commonly used strategies to novel concepts and technologies in computational fragment-based drug design. Particularly, in this review, specifications and advantages are compared between experimental and computational FBDD, and additionally, limitations and future prospective are discussed and emphasized.

Methylindoles and Methoxyindoles are Agonists and Antagonists of Human Aryl Hydrocarbon Receptor

Novel methylindoles were identified as endobiotic and xenobiotic ligands of the human aryl hydrocarbon receptor (AhR). We examined the effects of 22 methylated and methoxylated indoles on the transcriptional activity of AhRs. Employing reporter gene assays in AZ-AHR transgenic cells, we determined full agonist, partial agonist, or antagonist activities of tested compounds, having substantially variable EC₅₀, IC₅₀, and relative efficacies. The most effective agonists (E_MAX relative to 5 nM dioxin) of the AhR were 4-Me-indole (134%), 6-Me-indole (91%), and 7-MeO-indole (80%), respectively. The most effective antagonists of the AhR included 3-Me-indole (IC₅₀; 19 μM), 2,3-diMe-indole (IC₅₀; 11 μM), and 2,3,7-triMe-indole (IC₅₀; 12 μM). Reverse transcription polymerase chain reaction analyses of CYP1A1 mRNA in LS180 cells confirmed the data from gene reporter assays. The compound leads, 4-Me-indole and 7-MeO-indole, induced substantial nuclear translocation of the AhR and enriched binding of the AhR to the CYP1A1 promoter, as observed using fluorescent immunohistochemistry and chromatin immunoprecipitation assays, respectively. Molecular modeling and docking studies suggest the agonists and antagonists likely share the same binding pocket but have unique binding modes that code for their affinity. Binding pocket analysis further revealed that 4-methylindole and 7-methoxyindole can simultaneously bind to the pocket and produce synergistic interactions. Together, these data show a dependence on subtle and specific chemical indole structures as AhR modulators and furthermore underscore the importance of complete evaluation of indole compounds as nuclear receptor ligands.

Identification of Novel Allosteric Modulators of Glutamate Transporter EAAT2

Dysfunction of excitatory amino acid transporters (EAATs) has been implicated in the pathogenesis of various neurological disorders, such as stroke, brain trauma, epilepsy, and neurodegenerative diseases, among others. EAAT2 is the main subtype responsible for glutamate clearance in the brain, having a key role in regulating transmission and preventing excitotoxicity. Therefore, compounds that increase the expression or activity of EAAT2 have therapeutic potential for neuroprotection. Previous studies identified molecular determinants for EAAT2 transport stimulation in a structural domain that lies at the interface of the rigid trimerization domain and the central substrate binding transport domain. In this work, a hybrid structure based approach was applied, based on this molecular domain, to create a high-resolution pharmacophore. Subsequently, virtual screening of a library of small molecules was performed, identifying 10 hit molecules that interact at the proposed domain. Among these, three compounds were determined to be activators, four were inhibitors, and three had no effect on EAAT2-mediated transport in vitro. Further characterization of the two best ranking EAAT2 activators for efficacy, potency, and selectivity for glutamate over monoamine transporters subtypes and NMDA receptors and for efficacy in cultured astrocytes is demonstrated. Mutagenesis studies suggest that the EAAT2 activators interact with residues forming the interface between the trimerization and transport domains. These compounds enhance the glutamate translocation rate, with no effect on substrate interaction, suggesting an allosteric mechanism. The identification of these novel positive allosteric modulators of EAAT2 offers an innovative approach for the development of therapies based on glutamate transport enhancement.

Design, synthesis, biological evaluation and in vivo testing of dual phosphodiesterase 5 (PDE5) and histone deacetylase 6 (HDAC6)-selective inhibitors for the treatment of Alzheimer's disease

We have identified chemical probes that act as dual phosphodiesterase 5 (PDE5) and histone deacetylase 6 (HDAC6)-selective inhibitors (>1 log unit difference versus class I HDACs) to decipher the contribution of HDAC isoforms to the positive impact of dual-acting PDE5 and HDAC inhibitors on mouse models of Alzheimer's disease (AD) and fine-tune this systems therapeutics approach. Structure- and knowledge-based approaches led to the design of first-in-class molecules with the desired target compound profile: dual PDE5 and HDAC6-selective inhibitors. Compound 44b, which fulfilled the biochemical, functional and ADME-Tox profiling requirements and exhibited adequate pharmacokinetic properties, was selected as pharmacological tool compound and tested in a mouse model of AD (Tg2576) in vivo.

Diverse signalling of the platelet P2Y1 receptor leads to a dichotomy in platelet function

Platelet P2Y₁ receptor signalling via RhoGTPases is necessary for platelet-dependent leukocyte recruitment, where no platelet aggregation is observed. We investigated signalling cascades involved in distinct P2Y₁-dependent platelet activities in vitro, using specific inhibitors for phospholipase C (PLC) (U73122, to inhibit the canonical pathway), and RhoGTPases: Rac1 (NSC23766) and RhoA (ROCK inhibitor GSK429286). Human platelet rich plasma (for platelet aggregation) or isolated washed platelets (for chemotaxis assays) was treated with U73122, GSK429286 or NSC23766 prior to stimulation with adenosine diphosphate (ADP) or the P2Y₁ specific agonist MRS2365. Aggregation, chemotaxis (towards f-MLP), or platelet-induced human neutrophil chemotaxis (PINC) towards macrophage derived chemokine (MDC) was assessed. Molecular docking of ADP and MRS2365 to P2Y₁ was analysed using AutoDock Smina followed by GOLD molecular docking in the Accelrys Discovery Studio software. Inhibition of PLC, but not Rac1 or RhoA, suppressed platelet aggregation induced by ADP and MRS2365. In contrast, platelet chemotaxis and PINC, were significantly attenuated by inhibition of platelet Rac1 or RhoA, but not PLC. MRS2365, compared to ADP had a less pronounced effect on P2Y₁-induced aggregation, but a similar efficacy to stimulate platelet chemotaxis and PINC, which might be explained by differences in molecular interaction of ADP compared to MRS2365 with the P2Y₁ receptor. Platelet P2Y₁ receptor activation during inflammation signals through alternate pathways involving Rho GTPases in contrast to canonical P2Y₁ receptor induced PLC signalling. This might be explained by selective molecular interactions of ligands within the orthosteric site of the P2Y₁ receptor.

Molecular recognition by synthetic receptors: Application in field-effect transistor based chemosensing

Molecular recognition, i.e., ability of one molecule to recognize another through weak bonding interactions, is one of the bases of life. It is often implemented to sensing systems of high merits. Preferential recognition of the analyte (guest) by the receptor (host) induces changes in physicochemical properties of the sensing system. These changes are measured by using suitable signal transducers. Because of possibility of miniaturization, fast response, and high sensitivity, field-effect transistors (FETs) are more frequently being used for that purpose. A FET combined with a biological material offers the potential to overcome many challenges approached in sensing. However, low stability of biological materials under measurement conditions is a serious problem. To circumvent this problem, synthetic receptors were integrated with the gate surface of FETs to provide robust performance. In the present critical review, the approach utilized to devise chemosensors integrating synthetic receptors and FET transduction is discussed in detail. The progress in this field was summarized and important outcome was provided.