Molecular Dynamics Simulation on the Suppression Mechanism of Phosphorylation to Ser222 by Allosteric Inhibitors Targeting MEK1/2 Kinase.

Allosteric inhibitors of mitogen-activated protein kinase 1 (MEK1) reveal distinct interactions with MEK1 activation loop residues. The structural analyses will determine whether, and how, distinct inhibitors suppress the phosphorylation of MEK1 and may guide future therapeutic development. In this study, we explored the suppression mechanism of the phosphorylation process in the presence of MEK allosteric inhibitors, such as selumetinib, trametinib, cobimetinib, and CH5126766, by employing molecular dynamics simulations accompanied by principal component analysis. The simulations of wildtype MEK1 show that Ser222 can come close to γ-phosphate but not Ser218. We have found the conformation where Ser222 is within 5 Šof distance, which makes Ser222 accessible for γ-phosphate. The conformation analysis from the simulations of MEK1 in the presence of allosteric inhibitors reveals that the inhibitor restricts the flexibility of Ser222 through strong interactions with the activation loop, Lys97, and water mediates interactions with amino acids in the vicinity. The results reveal that all the inhibitors act as screeners between the activation loop and Mg-ATP and restricting the flexibility of the activation loop through strong interaction causes the suppression of the phosphorylation process of MEK1. The results conclude that a strong interaction of allosteric inhibitors with the activation loop restricts the movement of Ser222 toward Mg-ATP, which could be the dominant factor for the suppression of phosphorylation in MEK1. This research will provide novel insights to design effective anticancer therapeutics for targeting MEK1 in the future.

Discovery of 2',6-Bis(4-hydroxybenzyl)-2-acetylcyclohexanone, a Novel FtsZ Inhibitor.

Multi-drug resistance is increasing in the pathogenic bacterium S. pneumoniae, which is mainly responsible for meningitis and community-acquired pneumonia (CAP), highlighting the need for new anti-pneumococcal agents. We have identified a potential anti-pneumococcal agent, enol 3, which acts by hindering the cell division process by perturbing Z-ring dynamics inside the cell. Enol 3 was also shown to inhibit FtsZ polymerization and induce its aggregation in vitro but does not affect the activity of tubulin and alkaline phosphatase. Docking studies show that 3 binds near the T7 loop, which is the catalytic site of FtsZ. Similar effects on Z-ring and FtsZ assembly were observed in B. subtilis, indicating that 3 could be a broad-spectrum anti-bacterial agent useful in targeting Gram-positive bacteria. In conclusion, compound 3 shows strong anti-pneumococcal activity, prompting further pre-clinical studies to explore its potential.

A Carbocyclic Curcumin Inhibits Proliferation of Gram-Positive Bacteria by Targeting FtsZ.

Inhibition of FtsZ assembly has been found to stall bacterial cell division. Here, we report the identification of a potent carbocyclic curcumin analogue (2d) that inhibits Bacillus subtilis 168 cell proliferation by targeting the assembly of FtsZ. 2d also showed potent inhibitory activity (minimum inhibitory concentrations of 2-4 mg/L) against several clinically important species of Gram-positive bacteria, including methicillin-resistant Staphylococcus aureus. In addition, 2d displayed a significantly reduced inhibitory effect on human cervical cancer cells in comparison to its effect on bacterial cells. Using live cell imaging of GFP-FtsZ by confocal microscopy, 2d was found to rapidly perturb the cytokinetic FtsZ rings in Bacillus subtilis cells. The immunofluorescence imaging of FtsZ also showed that 2d destroyed the Z-ring in bacteria within 5 min. Prolonged treatment with 2d produced filamentous bacteria, but 2d had no detectable effect either on the nucleoids or on the membrane potential of bacteria. 2d inhibited FtsZ assembly in vitro, whereas it had minimal effects on tubulin assembly. Interestingly, 2d strongly enhanced the GTPase activity of FtsZ and reduced the GTPase activity of tubulin. Furthermore, 2d bound to purified FtsZ with a dissociation constant of 4.0 ± 1.1 µM, and the binding of 2d altered the secondary structures of FtsZ. The results together suggested that the non-natural curcumin analogue 2d possesses powerful antibacterial activity against important pathogenic bacteria, and the evidence indicates that 2d inhibits bacterial proliferation by targeting FtsZ.

Potentiating the cellular targeting and anti-tumor activity of Dp44mT via binding to human serum albumin: two saturable mechanisms of Dp44mT uptake by cells.

Di-2-pyridylketone 4,4-dimethyl-3-thiosemicarbazone (Dp44mT) demonstrates potent anti-cancer activity. We previously demonstrated that 14C-Dp44mT enters and targets cells through a carrier/receptor-mediated uptake process. Despite structural similarity, 2-benzoylpyridine 4-ethyl-3-thiosemicarbazone (Bp4eT) and pyridoxal isonicotinoyl hydrazone (PIH) enter cells via passive diffusion. Considering albumin alters the uptake of many drugs, we examined the effect of human serum albumin (HSA) on the cellular uptake of Dp44mT, Bp4eT and PIH. Chelator-HSA binding studies demonstrated the following order of relative affinity: Bp4eT≈PIH>Dp44mT. Interestingly, HSA decreased Bp4eT and PIH uptake, potentially due to its high affinity for the ligands. In contrast, HSA markedly stimulated Dp44mT uptake by cells, with two saturable uptake mechanisms identified. The first mechanism saturated at 5-10 µM (B(max):1.20±0.04 × 107 molecules/cell; K(d):33±3 µM) and was consistent with a previously identified Dp44mT receptor/carrier. The second mechanism was of lower affinity, but higher capacity (B(max):2.90±0.12 × 107 molecules/cell; K(d):65±6 µM), becoming saturated at 100 µM and was only evident in the presence of HSA. This second saturable Dp44mT uptake process was inhibited by excess HSA and had characteristics suggesting it was mediated by a specific binding site. Significantly, the HSA-mediated increase in the targeting of Dp44mT to cancer cells potentiated apoptosis and could be important for enhancing efficacy.

Optimisation of LRRK2 inhibitors and assessment of functional efficacy in cell-based models of neuroinflammation.

LRRK2IN1 is a highly potent inhibitor of leucine-rich repeat kinase 2 (LRRK2, IC50 = 7.9 nM), an established target for treatment of Parkinson's disease. Two LRRK2IN1 analogues 1 and 2 were synthesised which retained LRRK2 inhibitory activity (1: IC50 = 72 nM; 2: IC50 = 51 nM), were predicted to have improved bioavailability and were efficacious in cell-based models of neuroinflammation. Analogue 1 inhibited IL-6 secretion from LPS-stimulated primary human microglia with EC50 = 4.26 µM. In order to further optimize the molecular properties of LRRK2IN1, a library of truncated analogues was designed based on docking studies. Despite lacking LRRK2 inhibitory activity, these compounds show anti-neuroinflammatory efficacy at micromolar concentration. The compounds developed were valuable tools in establishing a cell-based assay for assessing anti-neuroinflammatory efficacy of LRRK2 inhibitors. Herein, we present data that IL-1ß stimulated U87 glioma cell line is a reliable model for neuroinflammation, as data obtained in this model were consistent with results obtained using primary human microglia and astrocytes.

The development of CNS-active LRRK2 inhibitors using property-directed optimisation.

Mutations in PARK8/LRRK2 are the most common genetic cause of Parkinson's disease. Inhibition of LRRK2 kinase activity has neuroprotective benefits, and provides a means of addressing the underlying biochemical cause of Parkinson's disease for the first time. Initial attempts to develop LRRK2 inhibitors were largely unsuccessful and highlight shortcomings intrinsic to traditional, high throughput screening methods of lead discovery. Recently, amino-pyrimidine GNE-7915 was reported as a potent (IC50=9 nM) selective (1/187 kinases), brain-penetrant and non-toxic inhibitor of LRRK2. The use of in silico modelling, extensive in vitro assays and resource-efficient in vivo techniques to produce GNE-7915, reflects a trend towards the concerted optimisation of potency, selectivity and pharmacokinetic properties in early-stage drug development.

Design, Synthesis, and Pharmacological Evaluation of Fluorescent and Biotinylated Antagonists of ρ1 GABAC Receptors.

The ρ1 GABAC receptor is a ligand-gated chloride ion channel that shows promise as a therapeutic target for myopia, sleep disorders, memory and learning facilitation, and anxiety-related disorders. As such, there is a need for molecular probes to understand the role GABAC receptors play in physiological and pathological processes. To date, no labeled (either radioactive or fluorescent) GABAC selective ligand has been developed that can act as a marker for GABAC receptor visualization and localization studies. Herein, we report a series of fluorescent ligands containing different-sized linkers and fluorophores based around (S)-4-ACPBPA [(4-aminocyclopenten-1-yl)-butylphosphinic acid], a selective GABAC antagonist. One of these conjugates, (S)-4-ACPBPA-C5-BODIPY (13), displayed moderate potency (IC50 = 58.61 µM) and selectivity (>100 times) for ρ1 over α1ß2γ2L GABAA receptors. These conjugates are novel lead agents for the development of more potent and selective fluorescent probes for studying the localization and function of GABAC receptors in living cells.

Alanine scan of an immunosuppressive peptide (CP): analysis of structure-function relationships.

Core peptide is a hydrophobic peptide, the sequence of which is derived from the T-cell antigen receptor alpha-chain transmembrane region. Previous studies have shown that core peptide can inhibit T-cell-mediated immune responses both in vitro and in vivo. Here, we report the role each constituent amino acid plays within core peptide using an alanine scan and the amino acid effect on function using a biological antigen presentation assay. The biophysical behaviour of these analogues in model membranes was analysed using surface plasmon resonance studies and then binding correlated with T-cell function. Removal of any single hydrophobic amino acid between the two charged amino acids in core peptide (R, K) resulted in lower binding. Changing the overall net charge of core peptide, by removing either of the positively charged residues (R or K), had varying effects on peptide binding and IL-2 production. There was a direct correlation (ρ = 0.718) between peptide binding to model membranes and peptide ability to inhibit IL-2. Except for IL-2 inhibition, production of other T-cell cytokines such as GM-CSF, IFN-γ, IL-1α, IL-4, IL-5, IL-6, IL-10, IL-17 and T-cell antigen receptor alpha-chain was not detected using a fluorescent bead immunoassay. This study provides important structure-function relationships essential for further drug design.

The discovery of novel isoflavone pan peroxisome proliferator-activated receptor agonists.

Twenty three dual PPARα and γ molecules of natural product origin, previously reported by our group, were further investigated for pan PPAR transactivation against PPARδ. The in vitro cell toxicity profile, as well as, in silico study of the most active molecules within this new class of pan PPAR agonists are also described. 3',5' Dimethoxy-7 hydroxyisoflavone 6, Ψ-baptigenin 7, 4' fluoro-7 hydroxyisoflavone 8, and 3' methoxy-7 hydroxyisoflavone 9 were identified as the most potent molecules studied within the set compared to the commercially available pan PPAR agonist, bezafibrate 1. These novel active molecules may thus be useful as future leads in PPAR-related disorders, including type II diabetes mellitus and metabolic syndrome.

A σ(1) receptor pharmacophore derived from a series of N-substituted 4-azahexacyclo[5.4.1.0(2,6).0(3,10).0(5,9).0(8,11)]dodecan-3-ols (AHDs).

A library of N-substituted 4-azahexacyclo[5.4.1.0(2,6).0(3,10).0(5,9).0(8,11)]dodecan-3-ols (AHDs) was synthesized and subjected to competition binding assays at σ(1) and σ(2) receptors, as well as off-target screening of representative members at 44 other common central nervous system (CNS) receptors, transporters, and ion channels. Excluding 3 low affinity analogs, 31 ligands demonstrated nanomolar K(i) values for either σ receptor subtype. Several selective σ(1) and σ(2) ligands were discovered, with selectivities of up to 29.6 times for σ(1) and 52.4 times for σ(2), as well as several high affinity, subtype non-selective ligands. The diversity of structures and σ(1) affinities of the ligands allowed the generation of a σ(1) receptor pharmacophore that will enable the rational design of increasingly selective and potent σ(1) ligands for probing σ(1) receptor function.

Differentiating enantioselective actions of GABOB: a possible role for threonine 244 in the binding site of GABA(C) ρ(1) receptors.

Designing potent and subtype-selective ligands with therapeutic value requires knowledge about how endogenous ligands interact with their binding site. 4-Amino-3-hydroxybutanoic acid (GABOB) is an endogenous ligand found in the central nervous system in mammals. It is a metabolic product of GABA, the major inhibitory neurotransmitter. Homology modeling of the GABA(C) ρ(1) receptor revealed a potential H-bond interaction between the hydroxyl group of GABOB and threonine 244 (T244) located on loop C of the ligand binding site of the ρ(1) subunit. Using site-directed mutagenesis, we examined the effect of mutating T244 on the efficacy and pharmacology of GABOB and various ligands. It was found that mutating T244 to amino acids that lacked a hydroxyl group in their side chains produced GABA insensitive receptors. Only by mutating ρ(1)T244 to serine (ρ(1)T244S) produced a GABA responsive receptor, albeit 39-fold less sensitive to GABA than ρ(1)wild-type. We also observed changes in the activities of the GABA(C) receptor partial agonists, muscimol and imidazole-4-acetic acid (I4AA). At the concentrations we tested, the partial agonists antagonized GABA-induced currents at ρ(1)T244S mutant receptors (Muscimol: ρ(1)wild-type, EC(50) = 1.4 µM; ρ(1)T244S, IC(50) = 32.8 µM. I4AA: ρ(1)wild-type, EC(50) = 8.6 µM; ρ(1)T244S, IC(50) = 21.4 µM). This indicates that T244 is predominantly involved in channel gating. R-(-)-GABOB and S-(+)-GABOB are full agonists at ρ(1)wild-type receptors. In contrast, R-(-)-GABOB was a weak partial agonist at ρ(1)T244S (1 mM activates 26% of the current produced by GABA EC(50) versus ρ(1)wild-type, EC(50) = 19 µM; I(max) 100%), and S-(+)-GABOB was a competitive antagonist at ρ(1)T244S receptors (ρ(1)wild-type, EC(50) = 45 µM versus ρ(1)T244S, IC(50) = 417.4 µM, K(B) = 204 µM). This highlights that the interaction of GABOB with T244 is enantioselective. In contrast, the potencies of a range of antagonists tested, 3-aminopropyl(methyl)phosphinic acid (3-APMPA), 3-aminopropylphosphonic acid (3-APA), S- and R-(3-amino-2-hydroxypropyl)methylphosphinic acid (S-(-)-CGP44532 and R-(+)-CGP44533), were not altered. This suggests that T244 is not critical for antagonist binding. Receptor gating is dynamic, and this study highlights the role of loop C in agonist-evoked receptor activation, coupling agonist binding to channel gating.

Exploration of ring size in a series of cyclic vicinal diamines with σ1 receptor affinity.

Imidazolidine and 1,4-diazepane analogs of N-(2-benzofuranyl)methyl-N'-(4-alkoxybenzyl)piperazines were prepared to explore the effect of ring contraction and expansion on σ receptor affinity and subtype selectivity within a series of cyclic diamines. In vitro receptor binding assays revealed that all cyclic vicinal diamines possessed affinity and selectivity for σ(1) receptors. The imidazolidines possessed nanomolar σ(1) affinities (K(i)=6.45-53.5 nM), and relatively low levels of subtype selectivity (σ(2)/σ(1)=58-237). However, the piperazines and diazepanes achieved picomolar σ(1) interactions, with K(i) ranges of 0.05-10.28 and 0.10-0.194 nM, respectively. Moreover, the piperazines and diazepanes showed excellent discrimination over the σ(2) receptor, with σ(1) selectivities of 143-16140 and 220-11542, respectively.

Antiproliferative and antimigratory actions of synthetic long chain n-3 monounsaturated fatty acids in breast cancer cells that overexpress cyclooxygenase-2.

Cyclooxygenase-2 (COX-2) is overexpressed in many human cancers and converts the n-6 polyunsaturated fatty acid (PUFA) arachidonic acid to prostaglandin E(2) (PGE(2)), which drives tumorigenesis; in contrast, n-3 PUFA inhibit tumorigenesis. We tested the hypothesis that these antitumor actions of n-3 PUFA may involve the n-3 olefinic bond. n-3 Monounsaturated fatty acids (MUFAs) of chain length C16-C22 were synthesized and evaluated in MDA-MB-468 breast cancer cells that stably overexpressed COX-2 (MDA-COX-2 cells). Longer chain (C19-C22) n-3 MUFAs inhibited proliferation, activated apoptosis, decreased PGE(2) formation, and decreased cell invasion; C16-C18 analogues were less active. Molecular modeling showed that interactions of Arg120, Tyr355, and several hydrophobic amino acid residues in the COX-2 active site with C19-C22 MUFA analogues were favored. Thus, longer-chain n-3 MUFAs may be prototypes of novel anticancer agents that decrease the formation of PGE(2) in tumor cells that contain high levels of COX-2.

Role of human CYP3A4 in the biotransformation of sorafenib to its major oxidized metabolites.

The tyrosine kinase inhibitor drug sorafenib is used in the treatment of liver and renal cancers but adverse effects may necessitate dose interruption and under-dosage may lead to therapeutic failure. Sorafenib also undergoes cytochrome P450 (CYP)-dependent biotransformation to the N-oxide and other metabolites. However, although CYPs are major determinants of efficacy and toxicity the roles of these enzymes in the formation of multiple sorafenib metabolites are unclear. In the present study CYP-mediated pathways of sorafenib oxidation in human liver were evaluated. cDNA-expressed CYP3A4 was the major catalyst in the formation of the principal N-oxide and N-hydroxymethyl metabolites of sorafenib, as well as the minor N-desmethyl metabolite. In contrast, CYP3A5 exhibited only ~5% of the activity of CYP3A4 and eleven other CYPs and three flavin-containing monooxygenases were inactive. In human hepatic microsomes metabolite formation was correlated with CYP3A4-mediated midazolam 1'-hydroxylation, but not with other CYP-specific substrate oxidations. In accord with these findings the CYP3A4 inhibitor ketoconazole selectively inhibited microsomal sorafenib oxidation pathways. From computational modeling studies atoms in the structure of sorafenib that undergo biotransformation were within ~5.4 Å of the CYP3A4 heme. Important hydrogen bonding interactions between sorafenib and amino acids Ser-119 and Glu-374 in the active center of CYP3A4 were identified. These findings indicate that sorafenib is oxidized selectively by human CYP3A4. This information could be adapted in individualized approaches to optimize sorafenib safety and efficacy in cancer patients.

Consequences of linker length alteration of the α7 nicotinic acetylcholine receptor (nAChR) agonist, SEN12333.

A series of ligands based on SEN12333, containing either contracted or elongated alkyl chains, were synthesized and evaluated in molecular docking studies against a homology model of the α7 nicotinic acetylcholine receptor (nAChR) subtype. The predicted binding of all ligands was highly similar, with the exception of the analog containing a 5 methylene unit spacer. However, in vitro competition binding assays revealed that the ligands possessed dissimilar binding affinities, with a K(i) range of more than an order of magnitude (K(i)=0.50 to >10 µM), and only SEN12333 itself exhibited functional activity at the α7 nAChR.

The enantiomers of syn-2,3-difluoro-4-aminobutyric acid elicit opposite responses at the GABA(C) receptor.

The conformational behaviour and GABA receptor activity of the different stereoisomers of 2,3-difluoro-4-aminobutyric acid are described. Two enantiomeric GABA(C)-active ligands are identified, one of which is an agonist while the other is an antagonist. The results support an existing QSAR model of the bioactive geometry of GABA at GABA(C).

Synthesis and biological evaluation of 2,4-diaminoquinazoline derivatives as novel heat shock protein 90 inhibitors.

Novel 2,4-diaminoquinazoline derivatives originating from a virtual screening approach were designed, synthesized and their biological activities as heat shock protein 90 (Hsp90) inhibitors were evaluated. The prepared compounds exhibited significant anti-proliferative activities against DU-145, HT-29, HCT-116, A375P and MCF-7 cancer cell lines. The selected compounds were tested against Her2, a client protein of Hsp90, and showed significant reduction in Her2 protein expression. Compound 6b was found the most potent, reduced Her2 protein expression levels and induced Hsp70 protein expression levels significantly.

Structure based design of heat shock protein 90 inhibitors acting as anticancer agents.

Structure based drug design (SBDD) was used to discover heat shock protein 90 (HSP90) inhibitors useful in the treatment of cancer. By using the crystal structure of HSP90-ligand complex (1uyi), a docking model was prepared and was validated by external dataset containing known HSP90 inhibitors. This validated model was then used to virtually screen commercial databases, selected hits of which were bought and sent for real biological evaluation. Further as an alternative method, pharmacophores were generated using crystal structure conformations of ligands in HSP90 complexes (1uyi and 2bz5) and where used for virtual screening. Both cases yielded several hits containing novel scaffolds, particularly compound KHSP8 showed an IC(50) value of 0.902 µM in case of colon cancer (HT29), which is comparable to doxorubicin (0.828 µM). These compounds were being now used as leads for constructing small molecular libraries to get compounds with favourable pharmacokinetics and drug like properties.

Hybrid ortho/allosteric ligands for the adenosine A(1) receptor.

Many G protein-coupled receptors (GPCRs), including the adenosine A(1) receptor (A(1)AR), have been shown to be allosterically modulated by small molecule ligands. So far, in the absence of structural information, the exact location of the allosteric site on the A(1)AR is not known. We synthesized a series of bivalent ligands (4) with an increasing linker length between the orthosteric and allosteric pharmacophores and used these as tools to search for the allosteric site on the A(1)AR. The compounds were tested in both equilibrium radioligand displacement and functional assays in the absence and presence of a reference allosteric enhancer, (2-amino-4,5-dimethyl-3-thienyl)-[3-(trifluoromethyl)phenyl]methanone, PD81,723 (1). Bivalent ligand N(6)-[2-amino-3-(3,4-dichlorobenzoyl)-4,5,6,7-tetrahydrothieno[2,3-c]pyridin-6-yl-9-nonyloxy-4-phenyl]-adenosine 4h (LUF6258) with a 9 carbon atom spacer did not show significant changes in affinity or potency in the presence of 1, indicating that this ligand bridged both sites on the receptor. Furthermore, 4h displayed an increase in efficacy, but not potency, compared to the parent, monovalent agonist 2. From molecular modeling studies, we speculate that the allosteric site of the A(1)AR is located in the proximity of the orthosteric site, possibly within the boundaries of the second extracellular loop of the receptor.

Prospective validation of a comprehensive in silico hERG model and its applications to commercial compound and drug databases.

Ligand-based in silico hERG models were generated for 2 644 compounds using linear discriminant analysis (LDA) and support vector machines (SVM). As a result, the dataset used for the model generation is the largest publicly available (see Supporting Information). Extended connectivity fingerprints (ECFPs) and functional class fingerprints (FCFPs) were used to describe chemical space. All models showed area under curve (AUC) values ranging from 0.89 to 0.94 in a fivefold cross-validation, indicating high model consistency. Models correctly predicted 80 % of an additional, external test set; Y-scrambling was also performed to rule out chance correlation. Additionally models based on patch clamp data and radioligand binding data were generated separately to analyze their predictive ability when compared to combined models. To experimentally validate the models, 50 of the predicted hERG blockers from the Chembridge database and ten of the predicted non-hERG blockers from an in-house compound library were selected for biological evaluation. Out of those 50 predicted hERG blockers, tested at a concentration of 10 microM, 18 compounds showed more than 50 % displacement of [(3)H]astemizole binding to cell membranes expressing the hERG channel. K(i) values of four of the selected binders were determined to be in the micromolar and high nanomolar range (K(i) (VH01)=2.0 microM, K(i) (VH06)=0.15 microM, K(i) (VH19)=1.1 microM and K(i) (VH47)=18 microM). Of these four compounds, VH01 and VH47 showed also a second, even higher affinity binding site with K(i) values of 7.4 nM and 36 nM, respectively. In the case of non-hERG blockers, all ten compounds tested were found to be inactive, showing less than 50 % displacement of [(3)H]astemizole binding at 10 microM. These experimentally validated models were then used to virtually screen commercial compound databases to evaluate whether they contain hERG blockers. 109 784 (23 %) of Chembridge, 133 175 (38 %) of Chemdiv, 111 737 (31 %) of Asinex and 11 116 (18 %) of the Maybridge database were predicted to be hERG blockers by at least two of the models, a prediction which could, for example, be used as a pre-filtering tool for compounds with potential hERG liabilities.