QSAR Studies on Thiazolidines: A Biologically Privileged Scaffold

A large number of drugs and biologically relevant molecules contain heterocyclic systems. Often the presence of hetero atoms or groupings imparts preferential specificities in their biological responses. Amongst the heterocyclic systems, thiazolidine is a biologically important scaffold known to be associated with several biological activities. Some of the prominent biological responses attributed to this skeleton are antiviral, antibacterial, antifungal, antihistaminic, hypoglycemic, anti-inflammatory activities. This diversity in the biological response profiles of thiazolidine has attracted the attention of many researchers to explore this skeleton to its multiple potential against several activities. Many of these synthetic and biological explorations have been subsequently analyzed in detailed quantitative structure-activity relationship (QSAR) studies to correlate the respective structural features and physicochemical properties with the activities to identify the important structural components in deciding their activity behavior. In this, drugs or any biologically active molecules may be viewed as structural frames consisting of strategically positioned functional groups that will interact effectively with the complementary groups/sites of the receptor. With this in focus, the present article reviews the QSAR studies of diverse biological activities of the thiazolidines published during the past decade.

MECHANISMS OF SPECIES DIFFERENCES IN BRAIN ACETYLCHOLINESTERASE SENSITIVITY TO ORGANOPHOSPHATE INHIBITION

In vitro incubation of acetylcholinesterase from brain tissue of several species with organophosphate compounds indicated that the concentrations required to inhibit 50% of acetylcholinesterase activity (IC(,50)) differed from species to species for the same compound (Murphy, et al., 1968; Andersen, et al., 1972, 1977 and 1978).^ The hypothesis that non-specific binding proteins (Lauwerys and Murphy, 1969a,b) exerts a protective effect on acetylcholinesterase, and thus cause the differences observed in IC(,50) studies was tested by a ('3)H-DFP binding experiment. It was found that differences in the amount of non-specific binding protein cannot explain the observed differences observed in IC(,50) studies.^ An alternative hypothesis, that acetylcholinesterase from different species have different affinities for binding and/or different rates of phosphorylation by organophosphate insecticides was tested by determining the apparent affinity constant (k(,a)) and apparent rate of phosphorylation (k(,p)). Kinetic studies indicated that acetylcholinesterases from different species have different sensitivities to inhibition by organophosphate insecticides, and the differences are due to different affinities for binding and/or different rates of phosphorylation by the same organophosphate compound.^ Studies of the spontaneous reactivation of acetylcholinesterase after inhibition by organophosphate insecticides also indicated that acetylcholinesterases from different species have different rates and extents of spontaneous reactivation. This further substantiates the hypothesis that acetylcholinesterases from different species have different kinetic characteristics with respect to organophosphate insecticides inhibition.^ Eleven paraoxon analogs were synthesized for a quantitative structure-activity relationship study. It was found that the electron-withdrawing power ((sigma)) and hydrophobicity ((PARAGR)) of the substituent are important in determining the anti-cholinesterase activity of paraoxon analogs. Thus, predictions of species differences in acetylcholinesterase sensitivities to paraoxon analogs can be made if the physicochemical parameters ((sigma) and (PARAGR)) of the substituents are known.^ In another approach, i.e. enzyme modeling, the sensitivity of rat brain acetylcholinesterase to organophosphate insecticides was used as the independent variable to predict the sensitivities of acetylcholinesterases from other species to the same compound. Regression equations were derived for each species based on nineteen organophosphate insecticides studied. It was found, that in addition to paraoxon analogs, this method is also applicable to other organophosphate compounds with wide variations in structure. Thus, the sensitivities of acetylcholinesterases from other species can also be predicted from the sensitivity of rat brain acetylcholinesterase. ^

QSAR e dinâmica molecular no estudo de sistemas biomoleculares: predição da atividade biológica de antagonistas do receptor sigma-1 e simulações de bicamadas lipídicas

Diferentes abordagens teóricas têm sido utilizadas em estudos de sistemas biomoleculares com o objetivo de contribuir com o tratamento de diversas doenças. Para a dor neuropática, por exemplo, o estudo de compostos que interagem com o receptor sigma-1 (Sig-1R) pode elucidar os principais fatores associados à atividade biológica dos mesmos. Nesse propósito, estudos de Relações Quantitativas Estrutura-Atividade (QSAR) utilizando os métodos de regressão por Mínimos Quadrados Parciais (PLS) e Rede Neural Artificial (ANN) foram aplicados a 64 antagonistas do Sig-1R pertencentes à classe de 1-arilpirazóis. Modelos PLS e ANN foram utilizados com o objetivo de descrever comportamentos lineares e não lineares, respectivamente, entre um conjunto de descritores e a atividade biológica dos compostos selecionados. O modelo PLS foi obtido com 51 compostos no conjunto treinamento e 13 compostos no conjunto teste (r² = 0,768, q² = 0,684 e r²teste = 0,785). Testes de leave-N-out, randomização da atividade biológica e detecção de outliers confirmaram a robustez e estabilidade dos modelos e mostraram que os mesmos não foram obtidos por correlações ao acaso. Modelos também foram gerados a partir da Rede Neural Artificial Perceptron de Multicamadas (MLP-ANN), sendo que a arquitetura 6-12-1, treinada com as funções de transferência tansig-tansig, apresentou a melhor resposta para a predição da atividade biológica dos compostos (r²treinamento = 0,891, r²validação = 0,852 e r²teste = 0,793). Outra abordagem foi utilizada para simular o ambiente de membranas sinápticas utilizando bicamadas lipídicas compostas por POPC, DOPE, POPS e colesterol. Os estudos de dinâmica molecular desenvolvidos mostraram que altas concentrações de colesterol induzem redução da área por lipídeo e difusão lateral e aumento na espessura da membrana e nos valores de parâmetro de ordem causados pelo ordenamento das cadeias acil dos fosfolipídeos. As bicamadas lipídicas obtidas podem ser usadas para simular interações entre lipídeos e pequenas moléculas ou proteínas contribuindo para as pesquisas associadas a doenças como Alzheimer e Parkinson. As abordagens usadas nessa tese são essenciais para o desenvolvimento de novas pesquisas em Química Medicinal Computacional.

Molecular size as the main determinant of solute maximum flux across the skin

One of the most important determinants of dermatological and systemic penetration after topical application is the delivery or flux of solutes into or through the skin. The maximum dose of solute able to be delivered over a given period of time and area of application is defined by its maximum flux (J(max), mol per cm(2) per h) from a given vehicle. In this work, J(max) values from aqueous solution across human skin were acquired or estimated from experimental data and correlated with solute physicochemical properties. Whereas epidermal permeability coefficients (k(p)) are optimally correlated to solute octanol-water partition coefficient (K-ow) and molecular weight (MW) was found to be the dominant determinant of J(max) for this literature data set: log J(max)=-3.90-0.0190MW (n=87, r(2)=0.847, p

Empirical prediction of peptide octanol-water partition coefficients

Peptides are of great therapeutic potential as vaccines and drugs. Knowledge of physicochemical descriptors, including the partition coefficient P (commonly expressed in logarithm form: logP), is useful for screening out unsuitable molecules and also for the development of predictive Quantitative Structure-Activity Relationships (QSARs). In this paper we develop a new approach to the prediction of LogP values for peptides based on an empirical relationship between global molecular properties and measured physical properties. Our method was successful in terms of peptide prediction (total r2 = 0.641). The final model consisted of 5 physicochemical descriptors (molecular weight, number of single bonds, 2D-VDW volume, 2D-VSA hydrophobic and 2D-VSA polar). The approach is peptide specific and its predictive accuracy was high. Overall, 67% of the peptides were able to be predicted within +/-0.5 log units from the experimental values. Our method thus represents a novel prediction method with proven predictive ability.

On the hydrophobicity of peptides:comparing empirical predictions of peptide log P values

Peptides are of great therapeutic potential as vaccines and drugs. Knowledge of physicochemical descriptors, including the partition coefficient logP, is useful for the development of predictive Quantitative Structure-Activity Relationships (QSARs). We have investigated the accuracy of available programs for the prediction of logP values for peptides with known experimental values obtained from the literature. Eight prediction programs were tested, of which seven programs were fragment-based methods: XLogP, LogKow, PLogP, ACDLogP, AlogP, Interactive Analysis's LogP and MlogP; and one program used a whole molecule approach: QikProp. The predictive accuracy of the programs was assessed using r(2) values, with ALogP being the most effective (r( 2) = 0.822) and MLogP the least (r(2) = 0.090). We also examined three distinct types of peptide structure: blocked, unblocked, and cyclic. For each study (all peptides, blocked, unblocked and cyclic peptides) the performance of programs rated from best to worse is as follows: all peptides - ALogP, QikProp, PLogP, XLogP, IALogP, LogKow, ACDLogP, and MlogP; blocked peptides - PLogP, XLogP, ACDLogP, IALogP, LogKow, QikProp, ALogP, and MLogP; unblocked peptides - QikProp, IALogP, ALogP, ACDLogP, MLogP, XLogP, LogKow and PLogP; cyclic peptides - LogKow, ALogP, XLogP, MLogP, QikProp, ACDLogP, IALogP. In summary, all programs gave better predictions for blocked peptides, while, in general, logP values for cyclic peptides were under-predicted and those of unblocked peptides were over-predicted.

The use of computational QSAR analysis in the toxicological evaluation of a series of 2-pyridylcarboxamidrazone candidate anti-tuberculosis compounds

A series of N1-benzylidene pyridine-2-carboxamidrazone anti-tuberculosis compounds has been evaluated for their cytotoxicity using human mononuclear leucocytes (MNL) as target cells. All eight compounds were significantly more toxic than dimethyl sulphoxide control and isoniazid (INH) with the exception of a 4-methoxy-3-(2-phenylethyloxy) derivative, which was not significantly different in toxicity compared with INH. The most toxic agent was an ethoxy derivative, followed by 3-nitro, 4-methoxy, dimethylpropyl, 4-methylbenzyloxy, 3-methoxy-4-(-2-phenylethyloxy) and 4-benzyloxy in rank order. In comparison with the effect of selected carboxamidrazone agents on cells alone, the presence of either N-acetyl cysteine (NAC) or glutathione caused a significant reduction in the toxicity of INH, as well as on the 4-benzyloxy derivative, although both increased the toxicity of a 4-N,N-dimethylamino-1-naphthylidene and a 2-t-butylthio derivative. The derivatives from this and three previous studies were subjected to computational analysis in order to derive equations designed to establish quantitative structure activity relationships for these agents. Twenty-five compounds were thus resolved into two groups (1 and 2), which on analysis yielded equations with r2 values in the range 0.65-0.92. Group 1 shares a common mode of toxicity related to hydrophobicity, where cytotoxicity peaked at logP of 3.2, while Group 2 toxicity was strongly related to ionisation potential. The presence of thiols such as NAC and GSH both promoted and attenuated toxicity in selected compounds from Group 1, suggesting that secondary mechanisms of toxicity were operating. These studies will facilitate the design of future low toxicity high activity anti-tubercular carboxamidrazone agents. © 2003 Elsevier Science B.V. All rights reserved.

Étude de la relation entre les conformations et la signalisation des 7TMRs

L’interaction d’un ligand avec un récepteur à sept domaines transmembranaires (7TMR) couplé aux protéines G, mène à l’adoption de différentes conformations par le récepteur. Ces diverses conformations pourraient expliquer l’activation différentielle des voies de signalisation. Or, le lien entre la conformation et l’activité du récepteur n’est pas tout à fait claire. Selon les modèles classiques pharmacologiques, comme le modèle du complexe ternaire, il n’existe qu’un nombre limité de conformations qu’un récepteur peut adopter. Afin d’établir un lien entre la structure et la fonction des récepteurs, nous avons choisi dans un premier temps, le récepteur de chimiokine CXCR4 comme récepteur modèle. Ce dernier, est une cible thérapeutique prometteuse, impliqué dans l’entrée du VIH-1 dans les cellules cibles et dans la dissémination de métastases cancéreuses. Grâce au transfert d’énergie par résonance de bioluminescence (BRET) nous pouvons détecter les changements conformationnels des homodimères constitutifs de CXCR4 dans les cellules vivantes. En conséquence, nous avons mesuré les conformations de mutants de CXCR4 dont les mutations affecteraient sa fonction. Nous montrons que la capacité des mutants à activer la protéine Galphai est altérée suite au traitement avec l’agoniste SDF-1. Notamment, ces mutations altèrent la conformation du récepteur à l’état basal ainsi que la réponse conformationnelle induite suite au traitement avec l’agoniste SDF-1, l’agoniste partiel AMD3100 ou l’agoniste inverse TC14012. Ainsi, différentes conformations de CXCR4 peuvent donner lieu à une activation similaire de la protéine G, ce qui implique une flexibilité des récepteurs actifs qui ne peut pas être expliquée par le modèle du complexe ternaire (Berchiche et al. 2007). Également, nous nous sommes intéressés au récepteur de chimiokine CCR2, exprimé à la surface des cellules immunitaires. Il joue un rôle important dans l’inflammation et dans des pathologies inflammatoires telles que l’asthme. CCR2 forme des homodimères constitutifs et possède différents ligands naturels dont la redondance fonctionnelle a été suggérée. Nous avons étudié le lien entre les conformations et les activations d’effecteurs (fonctions) de CCR2. Notre hypothèse est que les différents ligands naturels induisent différentes conformations du récepteur menant à différentes fonctions. Nous montrons que les réponses de CCR2 aux différents ligands ne sont pas redondantes au niveau pharmacologique et que les chimiokines CCL8, CCL7 et CCL13 (MCP-2 à MCP-4) sont des agonistes partiels de CCR2, du moins dans les systèmes que nous avons étudiés. Ainsi, l’absence de redondance fonctionnelle parmi les chimiokines liant le même récepteur, ne résulterait pas de mécanismes complexes de régulation in vivo, mais ferait partie de leurs propriétés pharmacologiques intrinsèques (Berchiche et al. 2011). Enfin, nous nous sommes intéressés au récepteur de chimiokine CXCR7. Récemment identifié, CXCR7 est le deuxième récepteur cible de la chimiokine SDF-1. Cette chimiokine a été considérée comme étant capable d’interagir uniquement avec le récepteur CXCR4. Notamment, CXCR4 et CXCR7 possèdent un patron d’expression semblable dans les tissus. Nous avons évalué l’effet de l’AMD3100, ligand synthétique de CXCR4, sur la conformation et la signalisation de CXCR7. Nos résultats montrent qu’AMD3100, tout comme SDF-1, lie CXCR7 et augmente la liaison de SDF-1 à CXCR7. Grâce au BRET, nous montrons aussi qu’AMD3100 seul est un agoniste de CXCR7 et qu’il est un modulateur allostérique positif de la liaison de SDF-1 à CXCR7. Aussi, nous montrons pour la première fois le recrutement de la beta-arrestine 2 à CXCR7 en réponse à un agoniste. L’AMD3100 est un ligand de CXCR4 et de CXCR7 avec des effets opposés, ce qui appelle à la prudence lors de l’utilisation de cette molécule pour l’étude des voies de signalisation impliquant SDF-1 (Kalatskaya et al. 2009). En conclusion, nos travaux amènent des évidences qu’il existe plusieurs conformations actives des récepteurs et appuient les modèles de structure-activité des récepteurs qui prennent en considération leur flexibilité conformationnelle.

Conception, synthèse et diversification de l'azaGly-Pro : un mime de tour beta, outil d'étude structure-activité : application au développement d'un nouvel agent tocolytique

Les accouchements prématurés constituent un problème médical majeur en constante augmentation et ce, malgré tous les efforts mis en œuvre afin de contrer le déclenchement des contractions avant terme. Cette thèse relate du ''design'' rationnel d'un nouvel agent thérapeutique (i.e., tocolytique) qui serait capable de 1) arrêter les contractions, et 2) prolonger la gestation. Pour ce faire, une nouvelle cible, la prostaglandine F2α et son récepteur ont été sélectionnés et le peptidomimétisme a été choisi afin de résoudre cette problématique. L'introduction contient un historique rapide de la conception à la synthèse (''drug design'') du peptide parent, le PDC113, premier peptide a avoir démontré des aptitudes tocolytiques suffisantes pour faire du peptidomimétisme. La deuxième partie de l'introduction présente les concepts du peptidomimétisme appliqués au PDC113 qui ont permis d'accéder au PDC113.824, inhibiteur allostérique du récepteur de la prostaglandine F2α, et explique comment ce mime nous a permis d'élucider les mécanismes de signalisation intracellulaire impliqués dans la contraction musculaire lisse. Cette thèse présente la conception, la synthèse et l'étude structure-activité de mimes de repliement de tour β au sein du mime peptidique original (PDC113.824) dans lequel nous avons remplacé l'azabicycloalkane central (l'indolizidin-2-one) par une série d'autres azabicycloalcanes connus et des acides aza-aminés dont nous avons élaboré la synthèse. Dans un premier temps, une nouvelle stratégie de synthèse en solution de l'aza-glycyl-proline à partir de la diphényle hydrazone et du chloroformate de p-nitrophényle a été réalisée. Cette stratégie a permis d'éliminer les réactions secondaires de cyclisation intramoléculaires communément obtenues lors de l'introduction d'acides aza-aminés avec les protections traditionnelles de type carbamate en présence de phosgène, mais aussi de faciliter l'accès en une étape à des dérivés peptidiques du type aza-glycyle. L'élongation de l'aza-glycyl-proline en solution nous a permis d'accéder à un nouveau mime tetrapeptidique du Smac, un activateur potentiel de l'apoptose au sein de cellules cancéreuses. Par la suite, nous avons développé une stratégie de diversification sélective de l'azote α du résidu azaglycine en utilisant différents types d'halogénures d'alkyle en présence de tert-butoxyde de potassium. Afin de valider le protocole d'alkylation de l'aza-dipeptide, différents halogénures d'alkyle ont été testés. Nous avons également démontré l'utilité des aza-dipeptides résultants en tant que ''building block'' afin d'accéder à une variété d'azapeptides. En effet, l'aza-dipeptide a été déprotégée sélectivement soit en N-terminal soit en C-terminal, respectivement. D'autre part, la libération de l'amine de l'ester méthylique de l'aza-alkylglycyl-proline a conduit à une catégorie de composés à potentiel thérapeutique, les azadicétopipérazines (aza-DKP) par cyclisation intramoléculaire. Enfin, notre intérêt quant au développement d'un nouvel agent tocolytique nous a amené à développer une nouvelle voie de synthèse en solution du PDC113.824 permettant ainsi d'élucider les voies de signalisation intracellulaires du récepteur de la prostaglandine F2α. Afin de valider l'importance de la stéréochimie et d'étudier la relation structure/ activité du mime, nous avons remplacé l'indolizidin-2-one (I2aa) centrale du PDC113.824 par une série d'autres azabicycloalcanes et azadipeptides. Les azabicycloalcanes D-I2aa, quinolizidinone, et indolizidin-9-one ont été synthétisés et incorporés au sein du dit peptide ne donnant aucune activité ni in vitro ni ex vivo, validant ainsi l'importance du tour β de type II' pour le maintien de l'activité biologique du PDC113.824. Finalement, l'insertion d'une série de dérivés aza(alkyl)glycyl-prolyles a mené à de nouveaux inhibiteurs allostériques du récepteur de la PGF2α, l'un contenant l'azaglycine et l'autre, l'azaphénylalanine. Cette thèse a ainsi contribué, grâce à la conception et l'application de nouvelles méthodes de synthèse d'aza-peptides, au développement de nouveaux composés à potentiel thérapeutique afin d'inhiber le travail prématuré.

Possible pre-dissociation of diorganotin dihalide complexes: relationship between antitumour activity and structure

An examination of crystallographic data has indicated that the structure/activity relationship for diorganotin dihalide complexes is different from that of other metal dihalides, in that the SnN bond lengths appear to determine the antitumour activity.

Use of self-organizing maps and molecular descriptors to predict the cytotoxic activity of sesquiterpene lactones

Some sesquiterpene lactones (SLs) are the active compounds of a great number of traditionally medicinal plants from the Asteraceae family and possess considerable cytotoxic activity. Several studies in vitro have shown the inhibitory activity against cells derived from human carcinoma of the nasopharynx (KB). Chemical studies showed that the cytotoxic activity is due to the reaction of alpha,beta-unsaturated carbonyl structures of the SLs with thiols, such as cysteine. These studies support the view that SLs inhibit tumour growth by selective alkylation of growth-regulatory biological macromolecules, such as key enzymes, which control cell division, thereby inhibiting a variety of cellular functions, which directs the cells into apoptosis. In this study we investigated a set of 55 different sesquiterpene lactones, represented by 5 skeletons (22 germacranolides, 6 elemanolides, 2 eudesmanolides, 16 guaianolides and nor-derivatives and 9 pseudoguaianolides), in respect to their cytotoxic properties. The experimental results and 3D molecular descriptors were submitted to Kohonen self-organizing map (SOM) to classify (training set) and predict (test set) the cytotoxic activity. From the obtained results, it was concluded that only the geometrical descriptors showed satisfactory values. The Kohonen map obtained after training set using 25 geometrical descriptors shows a very significant match, mainly among the inactive compounds (similar to 84%). Analyzing both groups, the percentage seen is high (83%). The test set shows the highest match, where 89% of the substances had their cytotoxic activity correctly predicted. From these results, important properties for the inhibition potency are discussed for the whole dataset and for subsets of the different structural skeletons. (C) 2008 Elsevier Masson SAS. All rights reserved.

In Vitro Phytotoxicity and Antioxidant Activity of Selected Flavonoids

The knowledge of flavonoids involved in plant-plant interactions and their mechanisms of action are poor and, moreover, the structural characteristics required for these biological activities are scarcely known. The objective of this work was to study the possible in vitro phytotoxic effects of 27 flavonoids on the germination and early radical growth of Raphanus sativus L. and Lepidium sativum L., with the aim to evaluate the possible structure/activity relationship. Moreover, the antioxidant activity of the same compounds was also evaluated. Generally, in response to various tested flavonoids, germination was only slightly affected, whereas significant differences were observed in the activity of the various tested flavonoids against radical elongation. DPPH test confirms the antioxidant activity of luteolin, quercetin, catechol, morin, and catechin. The biological activity recorded is discussed in relation to the structure of compounds and their capability to interact with cell structures and physiology. No correlation was found between phytotoxic and antioxidant activities.

Mechanism of Action and Relationship Between Structure and Biological Activity of Ctx-Ha: A New Ceratotoxin-like Peptide from Hypsiboas albopunctatus

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

New synthetic bile acid analogue agonists of FXR and TGR5 receptors: Analytical methodologies for the study of their physico-chemical properties, pharmacokinetic activity and metabolism.

This thesis reports an integrated analytical approach for the study of physicochemical and biological properties of new synthetic bile acid (BA) analogues agonists of FXR and TGR5 receptors. Structure-activity data were compared with those previous obtained using the same experimental protocols on synthetic and natural occurring BA. The new synthetic BA analogues are classified in different groups according also to their potency as a FXR and TGR5 agonists: unconjugated and steroid modified BA and side chain modified BA including taurine or glycine conjugates and pseudo-conjugates (sulphonate and sulphate analogues). In order to investigate the relationship between structure and activity the synthetic analogues where admitted to a physicochemical characterization and to a preliminary screening for their pharmacokinetic and metabolism using a bile fistula rat model. Sensitive and accurate analytical methods have been developed for the quali-quantitative analysis of BA in biological fluids and sample used for physicochemical studies. Combined High Performance Liquid Chromatography Electrospray tandem mass spectrometry with efficient chromatographic separation of all studied BA and their metabolites have been optimized and validated. Analytical strategies for the identification of the BA and their minor metabolites have been developed. Taurine and glycine conjugates were identified in MS/MS by monitoring the specific ion transitions in multiple reaction monitoring (MRM) mode while all other metabolites (sulphate, glucuronic acid, dehydroxylated, decarboxylated or oxo) were monitored in a selected-ion reaction (SIR) mode with a negative ESI interface by the following ions. Accurate and precise data where achieved regarding the main physicochemical properties including solubility, detergency, lipophilicity and albumin binding . These studies have shown that minor structural modification greatly affect the pharmacokinetics and metabolism of the new analogues in respect to the natural BA and on turn their site of action, particularly where their receptor are located in the enterohepatic circulation.