GPCR structure, function, drug discovery and crystallography: report from Academia-Industry International Conference (UK Royal Society) Chicheley Hall, 1-2 September 2014

G-protein coupled receptors (GPCRs) are the targets of over half of all prescribed drugs today. The UniProt database has records for about 800 proteins classified as GPCRs, but drugs have only been developed against 50 of these. Thus, there is huge potential in terms of the number of targets for new therapies to be designed. Several breakthroughs in GPCRs biased pharmacology, structural biology, modelling and scoring have resulted in a resurgence of interest in GPCRs as drug targets. Therefore, an international conference, sponsored by the Royal Society, with world-renowned researchers from industry and academia was recently held to discuss recent progress and highlight key areas of future research needed to accelerate GPCR drug discovery. Several key points emerged. Firstly, structures for all three major classes of GPCRs have now been solved and there is increasing coverage across the GPCR phylogenetic tree. This is likely to be substantially enhanced with data from x-ray free electron sources as they move beyond proof of concept. Secondly, the concept of biased signalling or functional selectivity is likely to be prevalent in many GPCRs, and this presents exciting new opportunities for selectivity and the control of side effects, especially when combined with increasing data regarding allosteric modulation. Thirdly, there will almost certainly be some GPCRs that will remain difficult targets because they exhibit complex ligand dependencies and have many metastable states rendering them difficult to resolve by crystallographic methods. Subtle effects within the packing of the transmembrane helices are likely to mask and contribute to this aspect, which may play a role in species dependent behaviour. This is particularly important because it has ramifications for how we interpret pre-clinical data. In summary, collaborative efforts between industry and academia have delivered significant progress in terms of structure and understanding of GPCRs and will be essential for resolving problems associated with the more difficult targets in the future.

Stoichiometric Molecularly Imprinted Polymers for the Recognition of Anti-Cancer Pro-drug Tegafur

Molecularly Imprinted Polymers (MIPs) targeting tegafur, an anti-cancer 5-fluorouracil pro-drug, have been prepared by stoichiometric imprinting using 2,6-bis(acrylamido)pyridine (BAAPy) as the functional monomer. Solution association between tegafur and BAAPy was studied by 1H NMR titration, which confirmed the formation of 1:1 complexes with an affinity constant of 574±15 M-1 ¬in CDCl3. Evaluation of the synthesised materials by HPLC and equilibrium rebinding experiments revealed high selectivity of the imprinted polymer for the pro-drug versus 5-fluorouracil and other competing analytes, with maximum imprinting factors of 25.3 and a binding capacity of 45.1 μmol g-1. The synthesised imprinted polymer was employed in solid-phase extraction of the pro-drug using an optimised protocol that included a simple wash with the porogen used in the preparation of the material. Tegafur recoveries of up to 96% were achieved from aqueous samples and 92% from urine samples spiked with the template and three competing analytes. The results demonstrate the potential of the prepared polymers in the pre-concentration of tegafur from biological samples, which could be an invaluable tool in the monitoring of patient compliance and drug uptake and excretion.

Medicinal chemistry approaches to malaria drug discovery

Tese de doutoramento, Farmácia (Química Farmacêutica e Terapêutica), Universidade de Lisboa, Faculdade de Farmácia, 2016

Use of cellular impedance to characterize ligand functional selectivity at G protein-coupled receptors

Les récepteurs couplés aux protéines G (RCPGs) représentent la plus grande famille de cibles thérapeutiques pour le traitement d’une panoplie de pathologies humaines. Bien que plusieurs décennies de recherche aient permis de façonner nos connaissances sur ces protéines membranaires, notre compréhension des déterminants moléculaires de leur activité signalétique reste encore limitée. De ces domaines de recherche, une avancée récente a mis à jour un nouveau phénomène, appelé sélectivité fonctionnelle des ligands, qui a bouleversé les paradigmes décrivant leu fonctionnement de ces récepteurs. Ce concept émane d’observations montrant que l’activité pharmacologique de certains ligands n’est pas nécessairement conservée sur tout le répertoire signalétiques connu du récepteur et peu se restreindre à l'activation sélective d’un sous-groupe de voies de signalisation.Ce nouveau modèle pharmacologique de l'activation des RCPG ouvre de nouvelles possibilités pour la découverte de médicaments plus efficace et sûr, ciblant les RCPGs. En effet, il permet la conception de molécules modulant spécifiquement les voies signalétiques d’intérêt thérapeutique, sans engager les autres voies qui pourraient mener à des effets secondaires indésirables ou de la tolérance. Cette thèse décrit l'utilisation d'une nouvelle approche sans marquage, basée sur la mesure du changement l'impédance cellulaire. Par la mesure des changements cellulaires, comme la morphologie, l’adhésion et/ou la redistribution des macromolécules, cette approche permet de mesurer de façon simultanée l'activité de plusieurs voies de signalisation impliqués dans ces réponses. Utilisant le récepteur β2-adrénergique (β2AR) comme modèle, nous avons démontré que les variations dans l’impédance cellulaire étaient directement liées à l’activation de multiples voies de signalisation suite à la stimulation du récepteur par son ligand. L’agoniste type du β2AR, l’isoprotérénol, s’est avéré induire une réponse d’impédance dose-dépendante constituée, dans le temps, de plusieurs caractéristiques distinctes pouvant être bloquées de façon compétitive par l’antagoniste ICI118,551 Par l’utilisation d’inhibiteurs sélectifs, nous avons été en mesure de déterminer la contribution de plusieurs voies signalétiques canoniques, comme les voies dépendantes de Gs et Gi, la production d’AMPc et l’activation de ERK1/2, sur ces changements. De plus, la dissection de la réponse d’impédance a permis d’identifier une nouvelle voie de mobilisation du Ca2+ contribuant à la réponse globale des changements initiés par la stimulation du β2AR. Dans une autre étude, nous avons rapporté que la réponse calcique induite par le β2AR serait attribuable à une transactivation Gs-dépendant du récepteur purinergique P2Y11, lui-même couplé à la protéine Gq. La mesure d’impédance permettant de distinguer et de décrire une pléiade d’activités signalétiques, nous avons émis l’hypothèse que des ligands arborant des profils signalétiques différents généreraient des réponses d’impédance distinctes. Le criblage d’une librairie de ligands spécifiques au β2AR a révélé une grande variété de signatures d’impédance. Grâce au développement d’une approche computationnelle innovatrice, nous avons été en mesure de regrouper ces signatures en cinq classes de composés, un regroupement qui s’est avéré hautement corrélé avec le profil signalétique des différents ligands. Nous avons ensuite combiné le criblage de composés par impédance avec l’utilisation d’inhibiteurs sélectifs de voies signalétiques afin d’augmenter la résolution du regroupement. En évaluant l’impact d’une voie signalétique donnée sur la signature d’impédance, nous avons été en mesure de révéler une plus grande variété de textures parmi les ligands. De plus, cette méthode s’est avérée efficace pour prédire le profil signalétique d’une librairie de composés non caractérisés, ciblant le β2AR. Ces travaux ont mené à l’élaboration d’une méthode permettant d’exprimer visuellement la sélectivité fonctionnelle de ligands et ont révélé de nouvelles classes de composés pour ce récepteur. Ces nouvelles classes de composés ont ensuite été testées sur des cardiomyocytes humains, confirmant que les composés regroupés dans différentes classes produisent des effets distincts sur la contractilité de ces cellules. Globalement, ces travaux démontrent la pertinence de l’utilisation de l’impédance cellulaire pour une évaluation précise des différences fonctionnelles parmi les composés ciblant les RCPGs. En fournissant une représentation pluridimensionnelle de la signalisation émanant des RCPGs à l’aide d’un seul essai ne requérant pas de marquage, les signatures d’impédance représentent une stratégie simple et innovante pour l’évaluation de la fonctionnalité sélective des ligands. Cette méthode pourrait être d’une grande utilité dans le processus de découverte de nouveaux médicaments.

Beta-blocker selectivity and airways obstruction

In this Drug Points article the author describes a case of unrecognised airways disease where prescribing timolol resulted in shortness of breath and comments on the issues it raises.

Structure-Based Approach for the Study of Estrogen Receptor Binding Affinity and Subtype Selectivity

Estrogens exert important physiological effects through the modulation of two human estrogen receptor (hER) subtypes, alpa (hER alpha) and beta (hER beta). Because the levels and relative proportion of hER alpha and hER beta differ significantly in different target cells, selective hER ligands could target specific tissues or pathways regulated by one receptor subtype without affecting the other. To understand the structural and chemical basis by which small molecule modulators are able to discriminate between the two subtypes, we have applied three-dimensional target-based approaches employing a series of potent hER-ligands. Comparative molecular field analysis (CoMFA) studies were applied to a data set of 81 hER modulators, for which binding affinity values were collected for both hER alpha and hER beta. Significant statistical coefficients were obtained (hER alpha, q(2) = 0.76; hER beta, q(2) = 0.70), indicating the internal consistency of the models. The generated models were validated using external test sets, and the predicted values were in good agreement with the experimental results. Five hER crystal structures were used in GRID/PCA investigations to generate molecular interaction fields (MIF) maps. hER alpha and hER beta were separated using one factor. The resulting 3D information was integrated with the aim of revealing the most relevant structural features involved in hER subtype selectivity. The final QSAR and GRID/PCA models and the information gathered from 3D contour maps should be useful for the design or novel hER modulators with improved selectivity.

Gaining ligand selectivity in thyroid hormone receptors via entropy

Nuclear receptors are important targets for pharmaceuticals, but similarities between family members cause difficulties in obtaining highly selective compounds. Synthetic ligands that are selective for thyroid hormone (TH) receptor beta (TR beta) vs. TR alpha reduce cholesterol and fat without effects on heart rate; thus, it is important to understand TR beta-selective binding. Binding of 3 selective ligands (GC-1, KB141, and GC-24) is characterized at the atomic level; preferential binding depends on a nonconserved residue (Asn-331 beta) in the TR beta ligand-binding cavity (LBC), and GC-24 gains extra selectivity from insertion of a bulky side group into an extension of the LBC that only opens up with this ligand. Here we report that the natural TH 3,5,3`-triodothyroacetic acid (Triac) exhibits a previously unrecognized mechanism of TR beta selectivity. TR x-ray structures reveal better fit of ligand with the TR alpha LBC. The TR beta LBC, however, expands relative to TR alpha in the presence of Triac (549 angstrom(3) vs. 461 angstrom(3)), and molecular dynamics simulations reveal that water occupies the extra space. Increased solvation compensates for weaker interactions of ligand with TR beta and permits greater flexibility of the Triac carboxylate group in TR beta than in TR alpha. We propose that this effect results in lower entropic restraint and decreases free energy of interactions between Triac and TR beta, explaining subtype-selective binding. Similar effects could potentially be exploited in nuclear receptor drug design.

Estrogen receptor: Structural differences and potential implications on selectivity examined by the GRID/CPCA approach

Selective Estrogen Receptor Modulators ( SERMs) have been developed, but the selectivity towards the subtypes ( ER or ER is not well understood. Based on three-dimensional structural properties of ligand binding domains, a model that takes into account this aspect was developed via molecular interaction fields and consensus principal component analysis (GRID/CPCA).

2D QSAR and similarity studies on cruzain inhibitors aimed at improving selectivity over cathepsin L

Hologram quantitative structure-activity relationships (HQSAR) were applied to a data set of 41 cruzain inhibitors. The best HQSAR model (Q(2) = 0.77; R-2 = 0.90) employing Surflex-Sim, as training and test sets generator, was obtained using atoms, bonds, and connections as fragment distinctions and 4-7 as fragment size. This model was then used to predict the potencies of 12 test set compounds, giving satisfactory predictive R-2 value of 0,88. The contribution maps obtained from the best HQSAR model are in agreement with the biological activities of the study compounds. The Trypanosoma cruzi cruzain shares high similarity with the mammalian homolog cathepsin L. The selectivity toward cruzam was checked by a database of 123 compounds, which corresponds to the 41 cruzain inhibitors used in the HQSAR model development plus 82 cathepsin L inhibitors. We screened these compounds by ROCS (Rapid Overlay of Chemical Structures), a Gaussian-shape volume overlap filter that can rapidly identify shapes that match the query molecule. Remarkably, ROCS was able to rank the first 37 hits as being only cruzain inhibitors. In addition, the area under the curve (AUC) obtained with ROCS was 0.96, indicating that the method was very efficient to distinguishing between cruzain and cathepsin L inhibitors. (c) 2007 Elsevier Ltd. All rights reserved.

Thermodynamics of lipophilic drug binding to intestinal fatty acid binding protein and permeation across membranes

Intestinal fatty acid binding protein (I-FABP) is present at high levels in the absorptive cells of the intestine (enterocytes), where it plays a role in the intracellular solubilization of fatty acids (FA). However, I-FABP has also been shown to bind to a range of non-FA ligands, including some lipophilic drug molecules. Thus, in addition to its central role in FA trafficking, I-FABP potentially serves as an important intracellular carrier of lipophilic drugs. In this study we provide a detailed thermodynamic analysis of the binding and stability properties of I-FABP in complex with a series of fibrate and fenamate drugs to provide an insight into the forces driving drug binding to I-FABP. Drug binding and selectivity for I-FABP are driven by the interplay of protein−ligand interactions and solvent processes. The Gibbs free energies (ΔG°) determined from dissociation constants at 25 °C ranged from −6.2 to −10 kcal/mol. The reaction energetics indicate that drug binding to I-FABP is an enthalpy−entropy driven process. The relationship between I-FABP stability and drug binding affinity was examined by pulse proteolysis. There is a strong coupling between drug binding and I-FABP stability. The effect of an I-FABP protein sink on the kinetics and thermodynamics of tolfenamic acid permeation across an artificial phospholipid membrane were investigated. I-FABP significantly decreased the energy barrier for desorption of tolfenamic acid from the membrane into the acceptor compartment. Taken together, these data suggest that the formation of stable drug−I-FABP complexes is thermodynamically viable under conditions simulating the reactant concentrations likely observed in vivo and maybe a significant biochemical process that serves as a driving force for passive intestinal absorption of lipophilic drugs.

Chemiluminescence detection of MDMA in street drug samples using tris (2,2'-bipyridine) ruthenium (III)

Tris(2,2'-bipyridine)ruthenium(II) chemiluminescence was investigated for the detection of 3,4-methylenedioxymethamphetamine (MDMA) and several related compounds in street drug samples. Optimization using flow injection analysis showed that the selectivity of the reagent can be targeted towards the detection of secondary amines by altering the pH of the reaction environment. The greater selectivity of this mode of detection, compared to UV-absorbance, reduces the probability of false positive results from interfering compounds. The detection limit for MDMA under these conditions was 0.48 μM. A HPLC method incorporating post-column tris(2,2'-bipyridine)ruthenium(II) chemiluminescence detection was applied to the determination of MDMA in five street drug samples. The results obtained were in good agreement with quantification performed using traditional UV-absorbance detection, which demonstrates the viability of this method for confirmatory analysis of drug samples. This is the first report of tris(2,2'-bipyridine)ruthenium(II) chemiluminescence for the detection of MDMA and related amphetamine derivatives.

AN INVESTIGATION OF THE ELECTRONIC-STRUCTURE OF THE ANTITUMOR DRUG ELLIPTICINE AND ITS DERIVATIVES .1. GEOMETRICAL AM1 STUDY

Ellipticine and its derivatives are a class of molecules that show antitumor and cytotoxic activity with a multimodal mechanism of action. In this paper we report a preliminary Austin Method One (AM1) study of ellipticine and some molecules derived from it. We have observed a relationship between charge density distribution and biological selectivity. A mechanism that could improve cytotoxic activity is proposed.

Electrochemical determination of antimalarial drug amodiaquine in maternal milk using a hemin-based electrode

Voltammetric analysis of amodiaquine using a hemin biosensor revealed a well-defined peak at 0.14 V (vs. Ag/AgCl), corresponding to the oxidation of amodiaquine at pH 7.0. The electrodic behavior indicated that the oxidation process was irreversible, and that it was controlled by diffusion. In addition to advantages such as high selectivity and sensitivity, the method developed could be used for the analysis of breast milk containing amodiaquine without any need for prior sample treatment, an important consideration in routine analysis laboratories. Measurements of the drug contained in breast milk were used to validate the technique. The detection limit for standard solutions was 3.30 mg L-1, and the quantification limit was 11.0 mg L-1. ©The Electrochemical Society.

Crystal structure of Jararacussin-I: The highly negatively charged catalytic interface contributes to macromolecular selectivity in snake venom thrombin-like enzymes

Snake venom serine proteinases (SVSPs) are hemostatically active toxins that perturb the maintenance and regulation of both the blood coagulation cascade and fibrinolytic feedback system at specific points, and hence, are widely used as tools in pharmacological and clinical diagnosis. The crystal structure of a thrombin-like enzyme (TLE) from Bothrops jararacussu venom (Jararacussin-I) was determined at 2.48 Å resolution. This is the first crystal structure of a TLE and allows structural comparisons with both the Agkistrodon contortrix contortrix Protein C Activator and the Trimeresurus stejnegeri plasminogen activator. Despite the highly conserved overall fold, significant differences in the amino acid compositions and three-dimensional conformations of the loops surrounding the active site significantly alter the molecular topography and charge distribution profile of the catalytic interface. In contrast to other SVSPs, the catalytic interface of Jararacussin-I is highly negatively charged, which contributes to its unique macromolecular selectivity. © 2012 The Protein Society.

Glycosomal Targets for Anti-Trypanosomatid Drug Discovery

Glycosomes are peroxisome-related organelles found in all kinetoplastid protists, including the human pathogenic species of the family Trypanosomatidae: Trypanosoma brucei, Trypanosoma cruzi and Leishmania spp. Glycosomes are unique in containing the majority of the glycolytic/gluconeogenic enzymes, but they also possess enzymes of several other important catabolic and anabolic pathways. The different metabolic processes are connected by shared co-factors and some metabolic intermediates, and their relative importance differs between the parasites or their distinct life-cycle stages, dependent on the environmental conditions encountered. By genetic or chemical means, a variety of glycosomal enzymes participating in different processes have been validated as drug targets. For several of these enzymes, as well as others that are likely crucial for proliferation, viability or virulence of the parasites, inhibitors have been obtained by different approaches such as compound libraries screening or design and synthesis. The efficacy and selectivity of some initially obtained inhibitors of parasite enzymes were further optimized by structure-activity relationship analysis, using available protein crystal structures. Several of the inhibitors cause growth inhibition of the clinically relevant stages of one or more parasitic trypanosomatid species and in some cases exert therapeutic effects in infected animals. The integrity of glycosomes and proper compartmentalization of at least several matrix enzymes is also crucial for the viability of the parasites. Therefore, proteins involved in the assembly of the organelles and transmembrane passage of substrates and products of glycosomal metabolism offer also promise as drug targets. Natural products with trypanocidal activity by affecting glycosomal integrity have been reported.