Characterization of a putative orexin receptor in Ciona intestinalis sheds light on the evolution of the orexin/hypocretin system in chordates.

Tunicates are evolutionary model organisms bridging the gap between vertebrates and invertebrates. A genomic sequence in Ciona intestinalis (CiOX) shows high similarity to vertebrate orexin receptors and protostome allatotropin receptors (ATR). Here, molecular phylogeny suggested that CiOX is divergent from ATRs and human orexin receptors (hOX1/2). However, CiOX appears closer to hOX1/2 than to ATR both in terms of sequence percent identity and in its modelled binding cavity, as suggested by molecular modelling. CiOX was heterologously expressed in a recombinant HEK293 cell system. Human orexins weakly but concentration-dependently activated its Gq signalling (Ca2+ elevation), and the responses were inhibited by the non-selective orexin receptor antagonists TCS 1102 and almorexant, but only weakly by the OX1-selective antagonist SB-334867. Furthermore, the 5-/6-carboxytetramethylrhodamine (TAMRA)-labelled human orexin-A was able to bind to CiOX. Database mining was used to predict a potential endogenous C. intestinalis orexin peptide (Ci-orexin-A). Ci-orexin-A was able to displace TAMRA-orexin-A, but not to induce any calcium response at the CiOX. Consequently, we suggested that the orexin signalling system is conserved in Ciona intestinalis, although the relevant peptide-receptor interaction was not fully elucidated.

KNeMAP: a network mapping approach for knowledge-driven comparison of transcriptomic profiles.

MOTIVATION: Transcriptomic data can be used to describe the mechanism of action (MOA) of a chemical compound. However, omics data tend to be complex and prone to noise, making the comparison of different datasets challenging. Often, transcriptomic profiles are compared at the level of individual gene expression values, or sets of differentially expressed genes. Such approaches can suffer from underlying technical and biological variance, such as the biological system exposed on or the machine/method used to measure gene expression data, technical errors and further neglect the relationships between the genes. We propose a network mapping approach for knowledge-driven comparison of transcriptomic profiles (KNeMAP), which combines genes into similarity groups based on multiple levels of prior information, hence adding a higher-level view onto the individual gene view. When comparing KNeMAP with fold change (expression) based and deregulated gene set-based methods, KNeMAP was able to group compounds with higher accuracy with respect to prior information as well as is less prone to noise corrupted data. RESULT: We applied KNeMAP to analyze the Connectivity Map dataset, where the gene expression changes of three cell lines were analyzed after treatment with 676 drugs as well as the Fortino et al. dataset where two cell lines with 31 nanomaterials were analyzed. Although the expression profiles across the biological systems are highly different, KNeMAP was able to identify sets of compounds that induce similar molecular responses when exposed on the same biological system. AVAILABILITY AND IMPLEMENTATION: Relevant data and the KNeMAP function is available at: https://github.com/fhaive/KNeMAP and 10.5281/zenodo.7334711.

Azulene as a biphenyl mimetic in orexin/hypocretin receptor agonists.

Azulene is a rare ring structure in drugs, and we investigated whether it could be used as a biphenyl mimetic in known orexin receptor agonist Nag 26, which is binding to both orexin receptors OX1 and OX2 with preference towards OX2. The most potent azulene-based compound was identified as an OX1 orexin receptor agonist (pEC50 = 5.79 ± 0.07, maximum response = 81 ± 8% (s.e.m. of five independent experiments) of the maximum response to orexin-A in Ca2+ elevation assay). However, the azulene ring and the biphenyl scaffold are not identical in their spatial shape and electron distribution, and their derivatives may adopt different binding modes in the binding site.

Expanding the Paradigm of Structure-Based Drug Design: Molecular Dynamics Simulations Support the Development of New Pyridine-Based Protein Kinase C-Targeted Agonists.

Protein kinase C (PKC) modulators hold therapeutic potential for various diseases, including cancer, heart failure, and Alzheimer's disease. Targeting the C1 domain of PKC represents a promising strategy; the available protein structures warrant the design of PKC-targeted ligands via a structure-based approach. However, the PKC C1 domain penetrates the lipid membrane during binding, complicating the design of drug candidates. The standard docking-scoring approach for PKC lacks information regarding the dynamics and the membrane environment. Molecular dynamics (MD) simulations with PKC, ligands, and membranes have been used to address these shortcomings. Previously, we observed that less computationally intensive simulations of just ligand-membrane interactions may help elucidate C1 domain-binding prospects. Here, we present the design, synthesis, and biological evaluation of new pyridine-based PKC agonists implementing an enhanced workflow with ligand-membrane MD simulations. This workflow holds promise to expand the approach in drug design for ligands targeted to weakly membrane-associated proteins.

Exploring cooperative molecular contacts using a PostgreSQL database system.

Cooperative molecular contacts play an important role in protein structure and ligand binding. Here, we constructed a PostgreSQL database that stores structural information in the form of atomic environments and allows flexible mining of molecular contacts. Taking the Ser-His-Asp/Glu catalytic triad as a first test case, we demonstrate that the presence of a carboxylate oxygen atom in the vicinity of a His is associated with shorter Ser-OH..N-His bond in the PDB30 subset. We prospectively mine catalytic triads in unannotated proteins, suggesting catalytic functions for unannotated proteins. As a second test case, we demonstrate that this database system can include ligand atoms, represented by Sybyl atom types, by evaluating the proportion of counter-ions for ligand carboxylate oxygens.

Analysis of Biologics Molecular Descriptors towards Predictive Modelling for Protein Drug Development Using Time-Gated Raman Spectroscopy.

Pharmaceutical proteins, compared to small molecular weight drugs, are relatively fragile molecules, thus necessitating monitoring protein unfolding and aggregation during production and post-marketing. Currently, many analytical techniques take offline measurements, which cannot directly assess protein folding during production and unfolding during processing and storage. In addition, several orthogonal techniques are needed during production and market surveillance. In this study, we introduce the use of time-gated Raman spectroscopy to identify molecular descriptors of protein unfolding. Raman spectroscopy can measure the unfolding of proteins in-line and in real-time without labels. Using K-means clustering and PCA analysis, we could correlate local unfolding events with traditional analytical methods. This is the first step toward predictive modeling of unfolding events of proteins during production and storage.

Macrodomain Binding Compound MRS 2578 Inhibits Alphavirus Replication.

Alphaviruses are positive-strand RNA viruses causing febrile disease. Macrodomain-containing proteins, involved in ADP-ribose-mediated signaling, are encoded by both host cells and several virus groups, including alphaviruses. In this study, compound MRS 2578 that targets the human ADP-ribose glycohydrolase MacroD1 inhibited Semliki Forest virus production as well as viral RNA replication and replicase protein expression. The inhibitor was similarly active in alphavirus trans-replication systems, indicating that it targets the viral RNA replication stage.

Exploration of Pyrazolo[1,5-a]pyrimidines as Membrane-Bound Pyrophosphatase Inhibitors.

Inhibition of membrane-bound pyrophosphatase (mPPase) with small molecules offer a new approach in the fight against pathogenic protozoan parasites. mPPases are absent in humans, but essential for many protists as they couple pyrophosphate hydrolysis to the active transport of protons or sodium ions across acidocalcisomal membranes. So far, only few nonphosphorus inhibitors have been reported. Here, we explore the chemical space around previous hits using a combination of screening and synthetic medicinal chemistry, identifying compounds with low micromolar inhibitory activities in the Thermotoga maritima mPPase test system. We furthermore provide early structure-activity relationships around a new scaffold having a pyrazolo[1,5-a]pyrimidine core. The most promising pyrazolo[1,5-a]pyrimidine congener was further investigated and found to inhibit Plasmodium falciparum mPPase in membranes as well as the growth of P. falciparum in an ex vivo survival assay.

Multisubstituted pyrimidines effectively inhibit bacterial growth and biofilm formation of Staphylococcus aureus.

Biofilms are multicellular communities of microorganisms that generally attach to surfaces in a self-produced matrix. Unlike planktonic cells, biofilms can withstand conventional antibiotics, causing significant challenges in the healthcare system. Currently, new chemical entities are urgently needed to develop novel anti-biofilm agents. In this study, we designed and synthesized a set of 2,4,5,6-tetrasubstituted pyrimidines and assessed their antibacterial activity against planktonic cells and biofilms formed by Staphylococcus aureus. Compounds 9e, 10d, and 10e displayed potent activity for inhibiting the onset of biofilm formation as well as for killing pre-formed biofilms of S. aureus ATCC 25923 and Newman strains, with half-maximal inhibitory concentration (IC50) values ranging from 11.6 to 62.0 µM. These pyrimidines, at 100 µM, not only decreased the number of viable bacteria within the pre-formed biofilm by 2-3 log10 but also reduced the amount of total biomass by 30-50%. Furthermore, these compounds were effective against planktonic cells with minimum inhibitory concentration (MIC) values lower than 60 µM for both staphylococcal strains. Compound 10d inhibited the growth of S. aureus ATCC 25923 in a concentration-dependent manner and displayed a bactericidal anti-staphylococcal activity. Taken together, our study highlights the value of multisubstituted pyrimidines to develop novel anti-biofilm agents.

Mitoxantrone, pixantrone and mitoxantrone (2-hydroxyethyl)piperazine are toll-like receptor 4 antagonists, inhibit NF-κB activation, and decrease TNF-alpha secretion in primary microglia.

Toll-like receptor 4 (TLR4) recognizes various endogenous and microbial ligands and is an essential part in the innate immune system. TLR4 signaling initiates transcription factor NF-κB and production of proinflammatory cytokines. TLR4 contributes to the development or progression of various diseases including stroke, neuropathic pain, multiple sclerosis, rheumatoid arthritis and cancer, and better therapeutics are currently sought for these conditions. In this study, a library of 140 000 compounds was virtually screened and a resulting hit-list of 1000 compounds was tested using a cellular reporter system. The topoisomerase II inhibitor mitoxantrone and its analogues pixantrone and mitoxantrone (2-hydroxyethyl)piperazine were identified as inhibitors of TLR4 and NF-κB activation. Mitoxantrone was shown to bind directly to the TLR4, and pixantrone and mitoxantrone (2-hydroxyethyl)piperazine were shown to inhibit the production of proinflammatory cytokines such as tumor necrosis factor alpha (TNFα) in primary microglia. The inhibitory effect on NF-κB activation or on TNFα production was not mediated through cytotoxity at ≤ 1 µM concentration for pixantrone and mitoxantrone (2-hydroxyethyl)piperazine treated cells, as assessed by ATP counts. This study thus identifies a new mechanism of action for mitoxantrone, pixantrone, and mitoxantrone (2-hydroxyethyl)piperazine through the TLR4.

Binding Site Interactions of Modulators of Breast Cancer Resistance Protein, Multidrug Resistance-Associated Protein 2, and P-Glycoprotein Activity.

ATP-binding cassette (ABC)-transporters protect tissues by pumping their substrates out of the cells in many physiological barriers, such as the blood-brain barrier, intestine, liver, and kidney. These substrates include various endogenous metabolites, but, in addition, ABC transporters recognize a wide range of compounds, therefore affecting the disposition and elimination of clinically used drugs and their metabolites. Although numerous ABC-transporter inhibitors are known, the underlying mechanism of inhibition is not well characterized. The aim of this study is to deepen our understanding of transporter inhibition by studying the molecular basis of ligand recognition. In the current work, we compared the effect of 44 compounds on the active transport mediated by three ABC transporters: breast cancer resistance protein (BCRP and ABCG2), multidrug-resistance associated protein (MRP2 and ABCC2), and P-glycoprotein (P-gp and ABCB1). Eight compounds were strong inhibitors of all three transporters, while the activity of 36 compounds was transporter-specific. Of the tested compounds, 39, 25, and 11 were considered as strong inhibitors, while 1, 4, and 11 compounds were inactive against BCRP, MRP2, and P-gp, respectively. In addition, six transport-enhancing stimulators were observed for P-gp. In order to understand the observed selectivity, we compared the surface properties of binding cavities in the transporters and performed structure-activity analysis and computational docking of the compounds to known binding sites in the transmembrane domains and nucleotide-binding domains. Based on the results, the studied compounds are more likely to interact with the transmembrane domain than the nucleotide-binding domain. Additionally, the surface properties of the substrate binding site in the transmembrane domains of the three transporters were in line with the observed selectivity. Because of the high activity toward BCRP, we lacked the dynamic range needed to draw conclusions on favorable interactions; however, we identified amino acids in both P-gp and MRP2 that appear to be important for ligand recognition.

Discovery of Membrane-Bound Pyrophosphatase Inhibitors Derived from an Isoxazole Fragment.

Membrane-bound pyrophosphatases (mPPases) regulate energy homeostasis in pathogenic protozoan parasites and lack human homologues, which makes them promising targets in e.g. malaria. Yet only few nonphosphorus inhibitors have been reported so far. Here, we explore an isoxazole fragment hit, leading to the discovery of small mPPase inhibitors with 6-10 µM IC50 values in the Thermotoga maritima test system. Promisingly, the compounds retained activity against Plasmodium falciparum mPPase in membranes and inhibited parasite growth.

Deciphering Imidazoline Off-targets by Fishing in the Class A of GPCR field.

Based on the finding that a central antihypertensive agent with high affinity for I1-type imidazoline receptors - rilmenidine, shows cytotoxic effects on cultured cancer cell lines, it has been suggested that imidazoline receptors agonists might have a therapeutic potential in the cancer therapy. Nevertheless, potential rilmenidine side effects caused by activation of α-adrenoceptors, or other associated receptors and enzymes, might hinder its therapeutic benefits. Considering that human α-adrenoceptors belong to the rhodopsin-like class A of G-protein-coupled receptors (GPCRs) it is reasonable to assume that imidazolines might have the affinity for other receptors from the same class. Therefore, to investigate possible off-target effects of imidazoline ligands we have prepared a reverse docking protocol on class A GPCRs, using imidazoline ligands and their decoys. To verify our in silico results, three ligands with high scores and three ligands with low scores were tested for antagonistic activity on α2 - adrenoceptors.

Modeling membrane proteins: The importance of cysteine amino-acids.

Computational modeling of membrane proteins is critical to understand biochemical systems and to support chemical biology. In this work, we use a dataset of 448 non-redundant membrane protein chains to expose a "rule" that governs membrane protein structure: free cysteine thiols are not found accessible to oxidative compartments such as the extracellular space, but are rather involved in disulphide bridges. Taking as examples the 1018 three-dimensional models produced during the GPCR Dock 2008, 2010 and 2013 competitions and 390 models for a GPCR target in CASP13, we show that this rule was not accounted for by the modeling community. We thus highlight a new direction for model development that should lead to more accurate membrane protein models, especially in the loop domains.

Screening for Thermotoga maritima Membrane-Bound Pyrophosphatase Inhibitors.

Membrane-bound pyrophosphatases (mPPases) are dimeric enzymes that occur in bacteria, archaea, plants, and protist parasites. These proteins cleave pyrophosphate into two orthophosphate molecules, which is coupled with proton and/or sodium ion pumping across the membrane. Since no homologous proteins occur in animals and humans, mPPases are good candidates in the design of potential drug targets. Here we present a detailed protocol to screen for mPPase inhibitors utilizing the molybdenum blue reaction in a 96 well plate system. We use mPPase from the thermophilic bacterium Thermotoga maritima (TmPPase) as a model enzyme. This protocol is simple and inexpensive, producing a consistent and robust result. It takes only about one hour to complete the activity assay protocol from the start of the assay until the absorbance measurement. Since the blue color produced in this assay is stable for a long period of time, subsequent assay(s) can be performed immediately after the previous batch, and the absorbance can be measured later for all batches at once. The drawback of this protocol is that it is done manually and thus can be exhausting as well as require good skills of pipetting and time keeping. Furthermore, the arsenite-citrate solution used in this assay contains sodium arsenite, which is toxic and should be handled with necessary precautions.

Synthesis, Identification, and Structure-Activity Relationship Analysis of GATA4 and NKX2-5 Protein-Protein Interaction Modulators.

Transcription factors GATA4 and NKX2-5 directly interact and synergistically activate several cardiac genes and stretch-induced cardiomyocyte hypertrophy. Previously, we identified phenylisoxazole carboxamide 1 as a hit compound, which inhibited the GATA4-NKX2-5 transcriptional synergy. Here, the chemical space around the molecular structure of 1 was explored by synthesizing and characterizing 220 derivatives and structurally related compounds. In addition to the synergistic transcriptional activation, selected compounds were evaluated for their effects on transcriptional activities of GATA4 and NKX2-5 individually as well as potential cytotoxicity. The structure-activity relationship (SAR) analysis revealed that the aromatic isoxazole substituent in the southern part regulates the inhibition of GATA4-NKX2-5 transcriptional synergy. Moreover, inhibition of GATA4 transcriptional activity correlated with the reduced cell viability. In summary, comprehensive SAR analysis accompanied by data analysis successfully identified potent and selective inhibitors of GATA4-NKX2-5 transcriptional synergy and revealed structural features important for it.

Asymmetry in catalysis by Thermotoga maritima membrane-bound pyrophosphatase demonstrated by a nonphosphorus allosteric inhibitor.

Membrane-bound pyrophosphatases are homodimeric integral membrane proteins that hydrolyze pyrophosphate into orthophosphates, coupled to the active transport of protons or sodium ions across membranes. They are important in the life cycle of bacteria, archaea, plants, and parasitic protists, but no homologous proteins exist in vertebrates, making them a promising drug target. Here, we report the first nonphosphorus allosteric inhibitor of the thermophilic bacterium Thermotoga maritima membrane-bound pyrophosphatase and its bound structure together with the substrate analog imidodiphosphate. The unit cell contains two protein homodimers, each binding a single inhibitor dimer near the exit channel, creating a hydrophobic clamp that inhibits the movement of ß-strand 1-2 during pumping, and thus prevents the hydrophobic gate from opening. This asymmetry of inhibitor binding with respect to each homodimer provides the first clear structural demonstration of asymmetry in the catalytic cycle of membrane-bound pyrophosphatases.

Cartography of rhodopsin-like G protein-coupled receptors across vertebrate genomes.

We conduct a cartography of rhodopsin-like non-olfactory G protein-coupled receptors in the Ensembl database. The most recent genomic data (releases 90-92, 90 vertebrate genomes) are analyzed through the online interface and receptors mapped on phylogenetic guide trees that were constructed based on a set of ~14.000 amino acid sequences. This snapshot of genomic data suggest vertebrate genomes to harbour 142 clades of GPCRs without human orthologues. Among those, 69 have not to our knowledge been mentioned or studied previously in the literature, of which 28 are distant from existing receptors and likely new orphans. These newly identified receptors are candidates for more focused evolutionary studies such as chromosomal mapping as well for in-depth pharmacological characterization. Interestingly, we also show that 37 of the 72 human orphan (or recently deorphanized) receptors included in this study cluster into nineteen closely related groups, which implies that there are less ligands to be identified than previously anticipated. Altogether, this work has significant implications when discussing nomenclature issues for GPCRs.

Structure-Activity Relationships of 1-Benzoylazulenes at the OX1 and OX2 Orexin Receptors.

We previously demonstrated the potential of di- or trisubstituted azulenes as ligands (potentiators, weak agonists, and antagonists) of the orexin receptors. In this study we investigated 27 1-benzoylazulene derivatives, uncovering seven potentiators of the orexin response on OX1 and two weak dual orexin receptor agonists. For potentiators, replacement of the azulene scaffold by indole retained the activity of four out of six compounds. The structure-activity relationships for agonism and potentiation can be summarized into a bicyclic aromatic ring system substituted with two hydrogen-bond acceptors (1-position, benzoyl; 6-position, carboxyl/ester) within 7-8 Šof each other; a third acceptor at the 3-position is also well tolerated. The same pharmacophoric signature is found in the preferred conformations of the orexin receptor agonist Nag26 from molecular dynamics simulations. Subtle changes switch the activity between weak agonism and potentiation, suggesting overlapping binding sites.

Determinants of Orexin Receptor Binding and Activation-A Molecular Dynamics Study.

We assess the stability of two previously suggested binding modes for the neuropeptide orexin-A in the OX2 receptor through extensive molecular dynamics simulations. As the activation determinants of the receptor remain unknown, we simulated an unliganded receptor and two small-molecular ligands, the antagonist suvorexant and the agonist Nag26 for comparison. Each system was simulated in pure POPC membrane as well as in the 25% cholesterol-POPC membrane. In total, we carried out 36 µs of simulations. Through this set of simulations, we report a stable binding mode for the C-terminus of orexin-A. In addition, we suggest interactions that would promote orexin receptor activation, as well as others that would stabilize the inactive state.