Structure-activity relationships from in vitro efficacies of the thiazolide series against the intracellular apicomplexan protozoan Neospora caninum

Nitazoxanide (2-acetolyloxy-N-(5-nitro 2-thiazolyl) benzamide; NTZ) represents the parent compound of a novel class of broad-spectrum anti-parasitic compounds named thiazolides. NTZ is active against a wide variety of intestinal and tissue-dwelling helminths, protozoa, enteric bacteria and a number of viruses infecting animals and humans. While potent, this poses a problem in practice, since this obvious non-selectivity can lead to undesired side effects in both humans and animals. In this study, we used real time PCR to determine the in vitro activities of 29 different thiazolides (NTZ-derivatives), which carry distinct modifications on both the thiazole- and the benzene moieties, against the tachyzoite stage of the intracellular protozoan Neospora caninum. The goal was to identify a highly active compound lacking the undesirable nitro group, which would have a more specific applicability, such as in food animals. By applying self-organizing molecular field analysis (SOMFA), these data were used to develop a predictive model for future drug design. SOMFA performs self-alignment of the molecules, and takes into account the steric and electrostatic properties, in order to determine 3D-quantitative structure activity relationship models. The best model was obtained by overlay of the thiazole moieties. Plotting of predicted versus experimentally determined activity produced an r2 value of 0.8052 and cross-validation using the "leave one out" methodology resulted in a q2 value of 0.7987. A master grid map showed that large steric groups at the R2 position, the nitrogen of the amide bond and position Y could greatly reduce activity, and the presence of large steric groups placed at positions X, R4 and surrounding the oxygen atom of the amide bond, may increase the activity of thiazolides against Neospora caninum tachyzoites. The model obtained here will be an important predictive tool for future development of this important class of drugs.

Structure-activity relationships for a new family of sulfonylurea herbicides

Acetohydroxyacid synthase (AHAS; EC 2.2.1.6) catalyzes the first common step in branched-chain amino acid biosynthesis. The enzyme is inhibited by several chemical classes of compounds and this inhibition is the basis of action of the sulfonylurea and imidazolinone herbicides. The commercial sulfonylureas contain a pyrimidine or a triazine ring that is substituted at both meta positions, thus obeying the initial rules proposed by Levitt. Here we assess the activity of 69 monosubstituted sulfonylurea analogs and related compounds as inhibitors of pure recombinant Arabidopsis thaliana AHAS and show that disubstitution is not absolutely essential as exemplified by our novel herbicide, monosulfuron (2-nitro-N-(4'-methyl-pyrimidin-2'-yl) phenyl-sulfonylurea), which has a pyrimidine ring with a single meta substituent. A subset of these compounds was tested for herbicidal activity and it was shown that their effect in vivo correlates well with their potency in vitro as AHAS inhibitors. Three-dimensional quantitative structure-activity relationships were developed using comparative molecular field analysis and comparative molecular similarity indices analysis. For the latter, the best result was obtained when steric, electrostatic, hydrophobic and H-bond acceptor factors were taken into consideration. The resulting fields were mapped on to the published crystal structure of the yeast enzyme and it was shown that the steric and hydrophobic fields are in good agreement with sulfonylurea-AHAS interaction geometry.

Toxicidade de fármacos nitrofurânicos

During the structural designing of new drugs, it is possible predict the influence of specific chemical groups on pharmacological activity. Among these, the nitro group has potential antiparasitic activity, being present in many antimicrobial drugs, such as metronidazole, nitrofurazone, furazolidone, oxamniquine and chloramphenicol. Also, the introduction of the nitro group into a molecule can modify the physicochemical and electronic properties of the substance. Besides antimicrobial drugs, this group is also found in other drug classes, such as antiulcer, anti-inflamatory and anxiolytic. However, the use of the nitro group in drug design has encountered restrictions, due to the associated toxicity. This article is a review of the toxicity of nitrofuran compounds, as well the possible mechanisms involved and the strategy of latentiation by molecular modification to decrease their toxicity.

Relation between antioxidant activity and electronic structure of Phenols

We describe a new physicochemical descriptor of the antioxidant activity of phenols, the energy difference between the two highest occupied molecular orbitals, which we believe will improve quantitative structure-activity relationship studies about these compounds. (C) 2003 Wiley Periodicals, Inc.

Development, validation and uncertainty analysis of quantitative structure and activity relationship models for Log P of disinfection by -products

Hydrophobicity as measured by Log P is an important molecular property related to toxicity and carcinogenicity. With increasing public health concerns for the effects of Disinfection By-Products (DBPs), there are considerable benefits in developing Quantitative Structure and Activity Relationship (QSAR) models capable of accurately predicting Log P. In this research, Log P values of 173 DBP compounds in 6 functional classes were used to develop QSAR models, by applying 3 molecular descriptors, namely, Energy of the Lowest Unoccupied Molecular Orbital (ELUMO), Number of Chlorine (NCl) and Number of Carbon (NC) by Multiple Linear Regression (MLR) analysis. The QSAR models developed were validated based on the Organization for Economic Co-operation and Development (OECD) principles. The model Applicability Domain (AD) and mechanistic interpretation were explored. Considering the very complex nature of DBPs, the established QSAR models performed very well with respect to goodness-of-fit, robustness and predictability. The predicted values of Log P of DBPs by the QSAR models were found to be significant with a correlation coefficient R2 from 81% to 98%. The Leverage Approach by Williams Plot was applied to detect and remove outliers, consequently increasing R 2 by approximately 2% to 13% for different DBP classes. The developed QSAR models were statistically validated for their predictive power by the Leave-One-Out (LOO) and Leave-Many-Out (LMO) cross validation methods. Finally, Monte Carlo simulation was used to assess the variations and inherent uncertainties in the QSAR models of Log P and determine the most influential parameters in connection with Log P prediction. The developed QSAR models in this dissertation will have a broad applicability domain because the research data set covered six out of eight common DBP classes, including halogenated alkane, halogenated alkene, halogenated aromatic, halogenated aldehyde, halogenated ketone, and halogenated carboxylic acid, which have been brought to the attention of regulatory agencies in recent years. Furthermore, the QSAR models are suitable to be used for prediction of similar DBP compounds within the same applicability domain. The selection and integration of various methodologies developed in this research may also benefit future research in similar fields.

Photic driving in the electroencephalogram of children and adolescents: harmonic structure and relation to the resting state

In order to identify latent bioelectrical oscillators, 15 normal subjects (aged 9-17 years, 8 males, 7 females) were subjected to intermittent photic stimulation. The EEG amplitude spectra corresponding to the 11 fixed frequencies of stimulation presented (3-24 Hz) were combined to form "profiles" of the driving reaction in the right occipital area. The driving response varied with frequency, and was demonstrable in 70-100% of cases (using as criterion peak amplitudes 20% larger than those of the neighbors). The strongest responses were observed at the frequency closest to the alpha peak of the resting EEG. A secondary profile maximum was in the theta band. In 10 subjects, this maximum exceeded half the alpha peak (with an average of 72.4% of the alpha peak), while in the resting spectra, theta amplitudes were much lower than the alpha maxima. This responsiveness in theta activity seems to be characteristic of prepubertal and pubertal subjects. The profiles and resting EEG spectra showed a highly significant Pearson's correlation, with the peak in the theta band of the profiles being the main difference observed between them. The correlation coefficient was significantly correlated with the ratio of the maxima in the theta and alpha bands (R = -0.77, P<0.001). The correlation coefficient between profile and resting spectrum may be a useful indicator in screening methods used to reveal latent cerebral oscillators. Profiles for the second and third harmonics were correlated with those of the first harmonic (fundamental frequency), when considering the corresponding EEG frequencies. Peak frequencies in all three profiles were close to those of the individual's background alpha rhythm, and peak amplitudes in higher harmonics were not much lower than those of the fundamental frequency (mean values of 84 and 63%, for second and third harmonics, respectively).

Understanding the bioactivity of plant tannins: developments in analysis methods and structure-activity studies

Tannins, typically segregated into two major groups, the hydrolyzable tannins (HTs) and the proanthocyanidins (PAs), are plant polyphenolic secondary metabolites found throughout the plant kingdom. On one hand, tannins may cause harmful nutritional effects on herbivores, for example insects, and hence they work as plants’ defense against plant-eating animals. On the other hand, they may affect positively some herbivores, such as mammals, for example by their antioxidant, antimicrobial, anti-inflammatory or anticarcinogenic activities. This thesis focuses on understanding the bioactivity of plant tannins, their anthelmintic properties and the tools used for the qualitative and quantitative analysis of this endless source of structural diversity. The first part of the experimental work focused on the development of ultra-high performance liquid chromatography−tandem mass spectrometry (UHPLC-MS/MS) based methods for the rapid fingerprint analysis of bioactive polyphenols, especially tannins. In the second part of the experimental work the in vitro activity of isolated and purified HTs and their hydrolysis product, gallic acid, was tested against egg hatching and larval motility of two larval developmental stages, L1 and L2, of a common ruminant gastrointestinal parasite, Haemonchus contortus. The results indicated clear relationships between the HT structure and the anthelmintic activity. The activity of the studied compounds depended on many structural features, including size, functional groups present in the structure, and the structural rigidness. To further understand tannin bioactivity on a molecular level, the interaction between bovine serum albumin (BSA), and seven HTs and epigallocatechin gallate was examined. The objective was to define the effect of pH on the formation on tannin–protein complexes and to evaluate the stability of the formed complexes by gel electrophoresis and MALDI-TOF-MS. The results indicated that more basic pH values had a stabilizing effect on the tannin–protein complexes and that the tannin oxidative activity was directly linked with their tendency to form covalently stabilized complexes with BSA at increased pH.

The correlation between electronic structure and antimalarial activity of tetrahydropyridines

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Structure-activity relationship of mastoparan analogs: effects of the number and positioning of Lys residues on secondary structure, interaction with membrane-mimetic systems and biological activity

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Computational analyses on the structure-function relation in ion channels

Ion channels are protein molecules, embedded in the lipid bilayer of the cell membranes. They act as powerful sensing elements switching chemicalphysical stimuli into ion-fluxes. At a glance, ion channels are water-filled pores, which can open and close in response to different stimuli (gating), and one once open select the permeating ion species (selectivity). They play a crucial role in several physiological functions, like nerve transmission, muscular contraction, and secretion. Besides, ion channels can be used in technological applications for different purpose (sensing of organic molecules, DNA sequencing). As a result, there is remarkable interest in understanding the molecular determinants of the channel functioning. Nowadays, both the functional and the structural characteristics of ion channels can be experimentally solved. The purpose of this thesis was to investigate the structure-function relation in ion channels, by computational techniques. Most of the analyses focused on the mechanisms of ion conduction, and the numerical methodologies to compute the channel conductance. The standard techniques for atomistic simulation of complex molecular systems (Molecular Dynamics) cannot be routinely used to calculate ion fluxes in membrane channels, because of the high computational resources needed. The main step forward of the PhD research activity was the development of a computational algorithm for the calculation of ion fluxes in protein channels. The algorithm - based on the electrodiffusion theory - is computational inexpensive, and was used for an extensive analysis on the molecular determinants of the channel conductance. The first record of ion-fluxes through a single protein channel dates back to 1976, and since then measuring the single channel conductance has become a standard experimental procedure. Chapter 1 introduces ion channels, and the experimental techniques used to measure the channel currents. The abundance of functional data (channel currents) does not match with an equal abundance of structural data. The bacterial potassium channel KcsA was the first selective ion channels to be experimentally solved (1998), and after KcsA the structures of four different potassium channels were revealed. These experimental data inspired a new era in ion channel modeling. Once the atomic structures of channels are known, it is possible to define mathematical models based on physical descriptions of the molecular systems. These physically based models can provide an atomic description of ion channel functioning, and predict the effect of structural changes. Chapter 2 introduces the computation methods used throughout the thesis to model ion channels functioning at the atomic level. In Chapter 3 and Chapter 4 the ion conduction through potassium channels is analyzed, by an approach based on the Poisson-Nernst-Planck electrodiffusion theory. In the electrodiffusion theory ion conduction is modeled by the drift-diffusion equations, thus describing the ion distributions by continuum functions. The numerical solver of the Poisson- Nernst-Planck equations was tested in the KcsA potassium channel (Chapter 3), and then used to analyze how the atomic structure of the intracellular vestibule of potassium channels affects the conductance (Chapter 4). As a major result, a correlation between the channel conductance and the potassium concentration in the intracellular vestibule emerged. The atomic structure of the channel modulates the potassium concentration in the vestibule, thus its conductance. This mechanism explains the phenotype of the BK potassium channels, a sub-family of potassium channels with high single channel conductance. The functional role of the intracellular vestibule is also the subject of Chapter 5, where the affinity of the potassium channels hEag1 (involved in tumour-cell proliferation) and hErg (important in the cardiac cycle) for several pharmaceutical drugs was compared. Both experimental measurements and molecular modeling were used in order to identify differences in the blocking mechanism of the two channels, which could be exploited in the synthesis of selective blockers. The experimental data pointed out the different role of residue mutations in the blockage of hEag1 and hErg, and the molecular modeling provided a possible explanation based on different binding sites in the intracellular vestibule. Modeling ion channels at the molecular levels relates the functioning of a channel to its atomic structure (Chapters 3-5), and can also be useful to predict the structure of ion channels (Chapter 6-7). In Chapter 6 the structure of the KcsA potassium channel depleted from potassium ions is analyzed by molecular dynamics simulations. Recently, a surprisingly high osmotic permeability of the KcsA channel was experimentally measured. All the available crystallographic structure of KcsA refers to a channel occupied by potassium ions. To conduct water molecules potassium ions must be expelled from KcsA. The structure of the potassium-depleted KcsA channel and the mechanism of water permeation are still unknown, and have been investigated by numerical simulations. Molecular dynamics of KcsA identified a possible atomic structure of the potassium-depleted KcsA channel, and a mechanism for water permeation. The depletion from potassium ions is an extreme situation for potassium channels, unlikely in physiological conditions. However, the simulation of such an extreme condition could help to identify the structural conformations, so the functional states, accessible to potassium ion channels. The last chapter of the thesis deals with the atomic structure of the !- Hemolysin channel. !-Hemolysin is the major determinant of the Staphylococcus Aureus toxicity, and is also the prototype channel for a possible usage in technological applications. The atomic structure of !- Hemolysin was revealed by X-Ray crystallography, but several experimental evidences suggest the presence of an alternative atomic structure. This alternative structure was predicted, combining experimental measurements of single channel currents and numerical simulations. This thesis is organized in two parts, in the first part an overview on ion channels and on the numerical methods adopted throughout the thesis is provided, while the second part describes the research projects tackled in the course of the PhD programme. The aim of the research activity was to relate the functional characteristics of ion channels to their atomic structure. In presenting the different research projects, the role of numerical simulations to analyze the structure-function relation in ion channels is highlighted.

Structure-function relation in retinoid regulated gene expression

Retinoids are known to inhibit proliferation of and induce terminal differentiation of many normal and transformed cells. It has been postulated that retinoids exert their effect by altering gene expression. HL-60 cells and macrophages both respond to retinoic acid action by the rapid induction of the enzyme tissue transglutaminase. The induction has been shown to be due to increased transcription of the transglutaminase gene. The first part of the dissertation studied the structure-function relationship of retinoid-regulated transglutaminase induction, differentiation and proliferation in HL-60 cells using retinoid analogs. The results indicated strict structural constraints and a strong structure-function correlation between transglutaminase induction and differentiation; those retinoids that induced transglutaminase also induced differentiation, those analogs that did not induce transglutaminase could not induce differentiation. The ability of the retinoids to induce transglutaminase in HL-60 cells was paralleled in macrophages. However, the antiproliferative effect of the retinoids displayed less stringent structural constraints than their differentiation- and transglutaminase-inducing properties. Specifically all the retinoids were able to inhibit proliferation to varying extents. It is concluded that the induction of transglutaminase and of differentiation by retinoids is mediated by receptors. While receptor mediation cannot be entirely ruled out, with the current data no definitive statement can be made about the antiproliferative activity of retinoids. Also, the concordance in the ability of the retinoids to induce transglutaminase and the ability to induce differentiation of HL-60 cells suggests that the former is an early response of the cells to retinoids and differentiation a later consequence on the same pathway. Using the induction of transglutaminase as an index of the direct, or primary, effect of retinoids on gene expression, the second part of the dissertation investigates, by 2D gel electrophoresis, the alteration in the rates of synthesis of other proteins in macrophages and HL-60 cells in response to short incubations with retinoic acid. Any changes in parallel with transglutaminase were taken to indicate proteins directly under the control of retinoic acid. It is concluded that retinoic acid regulates the expression of a circumscribed set of genes in a cell-specific manner. The results support the hypothesis that retinoids exert their multiple effects on myeloid cells, in part, by receptor-mediated alternations in gene expression. ^

Trypanocidal structure-activity relationship for cis- and trans-methylpluviatolide

The trypanocidal activity of racemic mixtures of cis- and trans-methylpluviatolides was evaluated in vitro against trypomastigote forms of two strains of Trypanosoma cruzi, and in the enzymatic assay of T. cruzi gGAPDH. The cytotoxicity of the compounds was assessed by the MTT method using LLC-MK2 cells. The effect of the compounds on peroxide and NO production were also investigated. The mixture of the trans stereoisomers displayed trypanocidal activity (IC(50) similar to 89.3 mu M). Therefore, it was separated by chiral HPLC, furnishing the and (+) (-)-enantiomers. Only the (-)-enantiomer was active against the parasite (IC(50) similar to 18.7 mu M). Despite being inactive, the (+)-enantiomer acted as an antagonistic competitor. Trans-methylpluviatolide displayed low toxicity for LLC-MK(2) cells, with an IC(50) of 6.53 mM. Furthermore, methylpluviatolide neither inhibited gGAPDH activity nor hindered peroxide and NO production at the evaluated concentrations. (C) 2008 Elsevier Ltd. All rights reserved.

Schistosomicidal and trypanocidal structure-activity relationships for (+/-)-licarin A and its (-)- and (+)-enantiomers

(+/-)-Licarin A (1) was obtained by oxidative coupling, and its enantiomers, (-)-licarin A (2) and (+)-licarin A (3), were resolved by chiral HPLC. Schistosomicidal and trypanocidal activities of these compounds were evaluated in vitro against Schistosoma mansoni adult worms and trypomastigote forms of Trypanosoma cruzi. The racemic mixture (1) displayed significant schistosomicidal activity with an LC(50) value of 53.57 mu M and moderate trypanocidal activity with an IC(50) value of 127.17 mu M. On the other hand, the (-)-enantiomer (2), displaying a LC(50) value of 91.71 mu M, was more active against S. mansoni than the (+)-enantiomer (3), which did not show activity. For the trypanocidal assay, enantiomer 2 showed more significant activity (IC(50) of 23.46 mu M) than enantiomer 3, which showed an IC(50) value of 87.73 mu M. Therefore, these results suggest that (+/-)-licarin A (1) and (-)-licarin A (2) are promising compounds that could be used for the development of schistosomicidal and trypanocidal agents. (C) 2011 Elsevier Ltd. All rights reserved.

Structure-activity studies of conantokins as human N-methyl-D-aspartate receptor modulators

The activities of conantokin-G (con-G), conantokin-T (con-T), and several novel analogues have been studied using polyamine enhancement of [H-3]MK-801 binding to human glutamate-N-methyl-D-aspartate (NMDA) receptors, and their structures have been examined using CD and H-1 NMR spectroscopy. The potencies of con-G[A7], con-G, and con-T as noncompetitive inhibitors of spermine-enhanced [H-3]MK-801 binding to NMDA receptor obtained from human brain tissue are similar to those obtained using rat brain tissue. The secondary structure and activity of con-G are found to be highly sensitive to amino acid substitution and modification. NMR chemical shift data indicate that con-G, con-G[D8,D17], and con-G[A7] have similar conformations in the presence of Ca2+. This consists of a helix for residues 2-16, which is kinked in the vicinity of Gla10. This is confirmed by 3D structure calculations on con-G[A7]. Restraining this helix in a linear form (i.e., con-G[A7,E10-K13]) results in a minor reduction in potency. Incorporation of a 7-10 salt-bridge replacement (con-G[K7-E10]) prevents helix formation in aqueous solution and produces a peptide with low potency. Peptides with the Leu5-Tyr5 substitution also have low potencies (con-G[Y5,A7] and con-G[Y5,K7]) indicating that Leu5 in con-G is important for full antagonist behavior. We have also shown that the Gla-Ala7 substitution increases potency, whereas the Gla-Lys7 substitution has no effect. Con-G and con-G[K7] both exhibit selectivity between NMDA subtypes from mid-frontal and superior temporal gyri, but not between sensorimotor and mid-frontal gyri. Asn8 and/or Asn17 appear to be important for the ability of con-G to function as an inhibitor of polyamine-stimulated [3H]MK-801 binding, but not in maintaining secondary structure. The presence of Ca2+ does not increase the potencies of con-G and con-T for NMDA receptors but does stabilize the helical structures of con-G, con-G[D8,D17], and, to a lesser extent, con-G[A7]. The NMR data support the existence of at least two independent Ca2+-chelating sites in con-G, one involving Gla7 and possibly Gla3 and the other likely to involve Gla10 and/or Gla14.