Funcions densitat i semblança molecular quàntica: nous desenvolupaments i aplicacions

La present tesi, tot i que emmarcada dins de la teoria de les Mesures Semblança Molecular Quántica (MQSM), es deriva en tres àmbits clarament definits: - La creació de Contorns Moleculars de IsoDensitat Electrònica (MIDCOs, de l'anglès Molecular IsoDensity COntours) a partir de densitats electròniques ajustades. - El desenvolupament d'un mètode de sobreposició molecular, alternatiu a la regla de la màxima semblança. - Relacions Quantitatives Estructura-Activitat (QSAR, de l'anglès Quantitative Structure-Activity Relationships). L'objectiu en el camp dels MIDCOs és l'aplicació de funcions densitat ajustades, ideades inicialment per a abaratir els càlculs de MQSM, per a l'obtenció de MIDCOs. Així, es realitza un estudi gràfic comparatiu entre diferents funcions densitat ajustades a diferents bases amb densitats obtingudes de càlculs duts a terme a nivells ab initio. D'aquesta manera, l'analogia visual entre les funcions ajustades i les ab initio obtinguda en el ventall de representacions de densitat obtingudes, i juntament amb els valors de les mesures de semblança obtinguts prèviament, totalment comparables, fonamenta l'ús d'aquestes funcions ajustades. Més enllà del propòsit inicial, es van realitzar dos estudis complementaris a la simple representació de densitats, i són l'anàlisi de curvatura i l'extensió a macromolècules. La primera observació correspon a comprovar no només la semblança dels MIDCOs, sinó la coherència del seu comportament a nivell de curvatura, podent-se així observar punts d'inflexió en la representació de densitats i veure gràficament aquelles zones on la densitat és còncava o convexa. Aquest primer estudi revela que tant les densitats ajustades com les calculades a nivell ab initio es comporten de manera totalment anàloga. En la segona part d'aquest treball es va poder estendre el mètode a molècules més grans, de fins uns 2500 àtoms. Finalment, s'aplica part de la filosofia del MEDLA. Sabent que la densitat electrònica decau ràpidament al allunyar-se dels nuclis, el càlcul d'aquesta pot ser obviat a distàncies grans d'aquests. D'aquesta manera es va proposar particionar l'espai, i calcular tan sols les funcions ajustades de cada àtom tan sols en una regió petita, envoltant l'àtom en qüestió. Duent a terme aquest procés, es disminueix el temps de càlcul i el procés esdevé lineal amb nombre d'àtoms presents en la molècula tractada. En el tema dedicat a la sobreposició molecular es tracta la creació d'un algorisme, així com la seva implementació en forma de programa, batejat Topo-Geometrical Superposition Algorithm (TGSA), d'un mètode que proporcionés aquells alineaments que coincideixen amb la intuïció química. El resultat és un programa informàtic, codificat en Fortran 90, el qual alinea les molècules per parelles considerant tan sols nombres i distàncies atòmiques. La total absència de paràmetres teòrics permet desenvolupar un mètode de sobreposició molecular general, que proporcioni una sobreposició intuïtiva, i també de forma rellevant, de manera ràpida i amb poca intervenció de l'usuari. L'ús màxim del TGSA s'ha dedicat a calcular semblances per al seu ús posterior en QSAR, les quals majoritàriament no corresponen al valor que s'obtindria d'emprar la regla de la màxima semblança, sobretot si hi ha àtoms pesats en joc. Finalment, en l'últim tema, dedicat a la Semblança Quàntica en el marc del QSAR, es tracten tres aspectes diferents: - Ús de matrius de semblança. Aquí intervé l'anomenada matriu de semblança, calculada a partir de les semblances per parelles d'entre un conjunt de molècules. Aquesta matriu és emprada posteriorment, degudament tractada, com a font de descriptors moleculars per a estudis QSAR. Dins d'aquest àmbit s'han fet diversos estudis de correlació d'interès farmacològic, toxicològic, així com de diverses propietats físiques. - Aplicació de l'energia d'interacció electró-electró, assimilat com a una forma d'autosemblança. Aquesta modesta contribució consisteix breument en prendre el valor d'aquesta magnitud, i per analogia amb la notació de l'autosemblança molecular quàntica, assimilar-la com a cas particular de d'aquesta mesura. Aquesta energia d'interacció s'obté fàcilment a partir de programari mecanoquàntic, i esdevé ideal per a fer un primer estudi preliminar de correlació, on s'utilitza aquesta magnitud com a únic descriptor. - Càlcul d'autosemblances, on la densitat ha estat modificada per a augmentar el paper d'un substituent. Treballs previs amb densitats de fragments, tot i donar molt bons resultats, manquen de cert rigor conceptual en aïllar un fragment, suposadament responsable de l'activitat molecular, de la totalitat de l'estructura molecular, tot i que les densitats associades a aquest fragment ja difereixen degut a pertànyer a esquelets amb diferents substitucions. Un procediment per a omplir aquest buit que deixa la simple separació del fragment, considerant així la totalitat de la molècula (calcular-ne l'autosemblança), però evitant al mateix temps valors d'autosemblança no desitjats provocats per àtoms pesats, és l'ús de densitats de Forats de fermi, els quals es troben definits al voltant del fragment d'interès. Aquest procediment modifica la densitat de manera que es troba majoritàriament concentrada a la regió d'interès, però alhora permet obtenir una funció densitat, la qual es comporta matemàticament igual que la densitat electrònica regular, podent-se així incorporar dins del marc de la semblança molecular. Les autosemblances calculades amb aquesta metodologia han portat a bones correlacions amb àcids aromàtics substituïts, podent així donar una explicació al seu comportament. Des d'un altre punt de vista, també s'han fet contribucions conceptuals. S'ha implementat una nova mesura de semblança, la d'energia cinètica, la qual consisteix en prendre la recentment desenvolupada funció densitat d'energia cinètica, la qual al comportar-se matemàticament igual a les densitats electròniques regulars, s'ha incorporat en el marc de la semblança. A partir d'aquesta mesura s'han obtingut models QSAR satisfactoris per diferents conjunts moleculars. Dins de l'aspecte del tractament de les matrius de semblança s'ha implementat l'anomenada transformació estocàstica com a alternativa a l'ús de l'índex Carbó. Aquesta transformació de la matriu de semblança permet obtenir una nova matriu no simètrica, la qual pot ser posteriorment tractada per a construir models QSAR.

Synteza, badania strukturalne i spektroskopowe oraz aktywność przeciwnowotworowa i przeciwdrobnoustrojowa nowych pochodnych salinomycyny

Wydział Chemii: Pracownia Chemii Bioorganicznej

The extrapolation problem and how population modeling can help

We argue that population modeling can add value to ecological risk assessment by reducing uncertainty when extrapolating from ecotoxicological observations to relevant ecological effects. We review other methods of extrapolation, ranging from application factors to species sensitivity distributions to suborganismal (biomarker and "-omics'') responses to quantitative structure activity relationships and model ecosystems, drawing attention to the limitations of each. We suggest a simple classification of population models and critically examine each model in an extrapolation context. We conclude that population models have the potential for adding value to ecological risk assessment by incorporating better understanding of the links between individual responses and population size and structure and by incorporating greater levels of ecological complexity. A number of issues, however, need to be addressed before such models are likely to become more widely used. In a science context, these involve challenges in parameterization, questions about appropriate levels of complexity, issues concerning how specific or general the models need to be, and the extent to which interactions through competition and trophic relationships can be easily incorporated.

QSAR studies of multidentate nitrogen ligands used in lanthanide and actinide extraction processes

Quantitative structure activity relationships (QSARs) have been developed to optimise the choice of nitrogen heterocyclic molecules that can be used to separate the minor actinides such as americium(III) from europium(III) in the aqueous PUREX raffinate of nuclear waste. Experimental data on distribution coefficients and separation factors (SFs) for 47 such ligands have been obtained and show SF values ranging from 0.61 to 100. The ligands were divided into a training set of 36 molecules to develop the QSAR and a test set of 11 molecules to validate the QSAR. Over 1500 molecular descriptors were calculated for each heterocycle and the Genetic Algorithm was used to select the most appropriate for use in multiple regression equations. Equations were developed fitting the separation factors to 6-8 molecular descriptors which gave r(2) values of >0.8 for the training set and values of >0.7 for the test set, thus showing good predictive quality. The descriptors used in the equations were primarily electronic and steric. These equations can be used to predict the separation factors of nitrogen heterocycles not yet synthesised and/or tested and hence obtain the most efficient ligands for lanthanide and actinide separation. (C) 2003 Elsevier B.V. All rights reserved.

Classical and hologram QSAR studies on a series of inhibitors of trypanosomatid Glyceraldehyde-3-Phosphate Dehydrogenase

Leishmaniasis and trypanosomiasis are major causes of morbidity and mortality in both tropical and subtropical regions of the world. The current available drugs are limited, ineffective, and require long treatment regimens. Due to the high dependence of trypanosomatids on glycolysis as a source of energy, some glycolytic enzymes have been identified as attractive targets for drug design. In the present work, classical Two-Dimensional Quantitative Structure -Activity Relationships (2D QSAR) and Hologram QSAR (HQSAR) studies were performed on a series of adenosine derivatives as inhibitors of Leishmania mexicana Glyceraldehyde-3-Phosphate Dehydrogenase (LmGAPDH). Significant correlation coefficients (classical QSAR, r(2)=0.83 and q(2) =0.81; HQSAR, r(2)=0.91 and q(2) =0.86) were obtained for the 56 training set compounds, indicating the potential of the models for untested compounds. The models were then externally validated using a test set of 14 structurally related compounds and the predicted values were in good agreement with the experimental results (classical QSAR, r(pred)(2) = 0.94; HQSAR, r(pred)(2) = 0.92).

Discovery of New Inhibitors of Schistosoma mansoni PNP by Pharmacophore-Based Virtual Screening

Schistosomiasis is considered the second most important tropical parasitic disease, with severe socioeconomic consequences for millions of people worldwide. Schistosoma monsoni, one of the causative agents of human schistosomiasis, is unable to synthesize purine nucleotides de novo, which makes the enzymes of the purine salvage pathway important targets for antischistosomal drug development. In the present work, we describe the development of a pharmacophore model for ligands of S. mansoni purine nucleoside phosphorylase (SmPNP) as well as a pharmacophore-based virtual screening approach, which resulted in the identification of three thioxothiazolidinones (1-3) with substantial in vitro inhibitory activity against SmPNP. Synthesis, biochemical evaluation, and structure activity relationship investigations led to the successful development of a small set of thioxothiazolidinone derivatives harboring a novel chemical scaffold as new competitive inhibitors of SmPNP at the low-micromolar range. Seven compounds were identified with IC(50) values below 100 mu M. The most potent inhibitors 7, 10, and 17 with 1050 of 2, 18, and 38 mu M, respectively, could represent new potential lead compounds for further development of the therapy of schistosomiasis.

Conformational Properties of Angiotensin II and Its Active and Inactive TOAC-Labeled Analogs in the Presence of Micelles. Electron Paramagnetic Resonance, Fluorescence, and Circular Dichroism Studies

The interaction between angiotensin II (AII, DRVYIHPF) and its analogs carrying 2,2,6,6-tetramethylpiperidine-1-oxyl-4-amino-4-carboxylic acid (TOAC) and detergents-negatively charged sodium dodecyl sulfate (SDS) and zwitterionic N-hexadecyl-N,N-dimethyl-3-ammonio-1-propanesulfonate (HPS)-was examined by means of EPR, CD, and fluorescence. EPR spectra of partially active TOAC(1)-AII and inactive TOAC(3)-AII in aqueous solution indicated fast tumbling, the freedom of motion being greater at the N-terminus. Line broadening occurred upon interaction with micelles. Below SDS critical micelle concentration, broader lines indicated complex formation with tighter molecular packing than in micelles. Small changes in hyperfine splittings evinced TOAC location at the micelle-water interface. The interaction with anionic micelles was more effective than with zwitterionic micelles. Peptide-micelle interaction caused fluorescence increase. The TOAC-promoted intramolecular fluorescence quenching was more, pronounced for TOAC(3)-AII because of the proximity between the nitroxide and Tyr(4). CD spectra showed that although both AII and TOAC(1)-AII presented flexible conformations in water, TOAC(3)-AII displayed conformational restriction because of the TOAC-imposed bend (Schreier et al., Biopolymers 2004, 74, 389). In HPS, conformational changes were observed for the labeled peptides at neutral and basic pH. In SDS, all peptides underwent pH-dependent conformational changes. Although the spectra suggested similar folds for All and TOAC(1)-AII, different conformations were acquired by TOAC(3)-AII. The membrane environment has been hypothesized to shift conformational equilibria so as to stabilize the receptor-bound conformation of ligands. The fact that TOAC(3)-AII is unable to acquire conformations similar to those of native AII and partially active TOAC(1)-AII is probably the explanation for its lack of biological activity. (C) 2009 Wiley Periodicals, Inc. Biopolymers (Pept Sci) 92: 525-537, 2009.

Integration of Ligand- and Target-Based Virtual Screening for the Discovery of Cruzain Inhibitors

A myriad of methods are available for virtual screening of small organic compound databases. In this study we have successfully applied a quantitative model of consensus measurements, using a combination of 3D similarity searches (ROCS and EON), Hologram Quantitative Structure Activity Relationships (HQSAR) and docking (FRED, FlexX, Glide and AutoDock Vina), to retrieve cruzain inhibitors from collected databases. All methods were assessed individually and then combined in a Ligand-Based Virtual Screening (LBVS) and Target-Based Virtual Screening (TBVS) consensus scoring, using Receiving Operating Characteristic (ROC) curves to evaluate their performance. Three consensus strategies were used: scaled-rank-by-number, rank-by-rank and rank-by-vote, with the most thriving the scaled-rank-by-number strategy, considering that the stiff ROC curve appeared to be satisfactory in every way to indicate a higher enrichment power at early retrieval of active compounds from the database. The ligand-based method provided access to a robust and predictive HQSAR model that was developed to show superior discrimination between active and inactive compounds, which was also better than ROCS and EON procedures. Overall, the integration of fast computational techniques based on ligand and target structures resulted in a more efficient retrieval of cruzain inhibitors with desired pharmacological profiles that may be useful to advance the discovery of new trypanocidal agents.

Interaction of giant extracellular Glossoscolex paulistus hemoglobin (HbGp) with zwitterionic surfactant N-hexadecyl-N,N-dimethyl-3-ammonio-1-propanesulfonate (HPS): Effects of oligomeric dissociation

The present work focuses on the interaction between the zwitterionic surfactant N-hexadecyl-N,N-dimethyl-3-ammonio-1-propanesulfonate (HPS) and the giant extracellular hemoglobin of Glossoscolex paulistus (HbGp). Electronic optical absorption, fluorescence emission and circular dichroism spectroscopy techniques, together with Gel-filtration chromatography, were used in order to evaluate the oligomeric dissociation as well as the autoxidation of HbGp as a function of the interaction with HPS. A peculiar behavior was observed for the HPS-HbGp interaction: a complex ferric species formation equilibrium was promoted, as a consequence of the autoxidation and oligomeric dissociation processes. At pH 7.0, HPS is more effective up to 1 mM while at pH 9.0 the surfactant effect is more intense above 1 mM. Furthermore, the interaction of HPS with HbGp was clearly less intense than the interaction of this hemoglobin with cationic (CTAC) and anionic (SDS) surfactants. Probably, this lower interaction with HPS is due to two factors: (i) the lower electrostatic attraction between the HPS surfactant and the protein surface ionic sites when compared to the electrostatic interaction between HbGp and cationic and anionic surfactants, and (ii) the low cmc of HPS, which probably reduces the interaction of the surfactant in the monomeric form with the protein. The present work emphasizes the importance of the electrostatic contribution in the interaction between ionic surfactants and HbGp. Furthermore, in the whole HPS concentration range used in this study, no folding and autoxidation decrease induced by this surfactant were observed. This is quite different from the literature data on the interaction between surfactants and tetrameric hemoglobins, that supports the occurrence of this behavior for the intracellular hemoglobins at low surfactant concentration range. Spectroscopic data are discussed and compared with the literature in order to improve the understanding of hemoglobin-surfactant interaction as well as the acid isoelectric point (pI) influence of the giant extracellular hemoglobins on their structure-activity relationship. (c) 2007 Elsevier B.V. All rights reserved.

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.

Antitumoral, mutagenic and (anti)estrogenic activities of tingenone and pristimerin

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

The reactivity of oytho-methoxy-substituted catechol radicals with sulfhydryl groups: Contribution for the comprehension of the mechanism of inhibition of NADPH oxidase by apocynin

Redox processes are involved in the mechanism of action of NADPH oxidase inhibitors such as diphenyleneiodonium and apocynin. Here, we studied the structure-activity relationship for apocynin and analogous ortho-methoxy-substituted catechols as inhibitors of the NADPH oxidase in neutrophils and their reactivity with peroxidase. Aiming to alter the reduction potential, the ortho-methoxy-catechol moiety was kept constant and the substituents at para position related to the hydroxyl group were varied. Two series of compounds were employed: methoxy-catechols bearing electron-withdrawing groups (MC-W) such as apocynin, vanillin, 4-nitroguaiacol, 4-cyanoguaiacol, and methoxy-catechol bearing electron-donating groups (MC-D) such as 4-methylguaiacol and 4-ethylguaiacol. We found that MC-D were weaker inhibitors compared to MD-W. Furthermore, the radicals generated by oxidation of MC-W via MPO/H(2)O(2), but not for MC-D, were able to oxidize glutathione (GSH) as verified by the formation of thiyl radicals, depletion of GSH, and recycling of the ortho-methoxy-catechols during their oxidations. The capacity of oxidizing sulfhydryl (SH) groups was also verified when ovalbumin was incubated with MC-W, but not for MC-D. Since the effect of apocynin has been correlated with inactivation of the cytosolic fractions of the NADPH oxidase complex and its oxidation during the inhibitory process develops a special role in this process, we suggest that the close relationship between the reactivity of the radicals of MC-W compounds with thiol groups and their efficacy as NADPH oxidase inhibitor could be the chemical pathway behind the mechanism of action of apocynin and should be taken into account in the design of new and specific NADPH oxidase inhibitors. (c) 2007 Elsevier B.V. All rights reserved.

Investigating the effect of different positioning of lysine residues along the peptide chain of mastoparans for their secondary structures and biological activities

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

Antinociceptive profile of 2,3,6-trisubstituted piperidine alkaloids: 3-O-acetyl-spectaline and semi-synthetic derivatives of (-)-spectaline

In early studies, we have reported the antinociceptive profile of (-)-spectaline, a piperidine alkaloid from Cassia spectabilis. The present study describes the synthesis, the antinociceptive and anti-inflammatory activities of a series of 2,3,6-trialkyl-piperidine alkaloids: the natural (-)-3-O-acetyl-spectaline (LASSBio-755) and ten semi-synthetic spectaline derivatives. Structure-activity relationship (SARs) studies were performed. The structures of all synthesized derivatives were confirmed by means of nuclear magnetic resonance. Compounds were evaluated for their analgesic (acetic acid-induced mouse abdominal constrictions, hot-plate test, formalin-induced pain test) and some of them for the anti-inflammatory activities (carrageenan-induced rat paw edema test). The pharmacological results showed that several of the new compounds given orally at a dose of 100 mu mol/kg significantly inhibited the acetic acid-induced abdominal constrictions, but they were less active than (-)-spectaline. LASSBio-755 and LASSBio-776 were the most actives with 37% and 31.7% of inhibition. In the formalin-induced pain only LASSBio-776 was able to inhibit by 34.4% the paw licking response of the inflammatory phase, (-)-spectaline and LASSBio-755 did show any activity. In the carrageenan-induced rat paw edema, only (-)-spectaline exhibited an anti-inflammatory profile, showing an ED(50) value of 56.6 mu mol/kg. Our results suggest different mechanisms of action for the analgesic activity observed for LASSBio-776 (3-O-Bocspectaline), LASSBio-755 (3-O-acetyl-spectaline) and (-)-spectaline (LASSBio-754). The antinociceptive profile of some of the semi-synthetic spectaline derivatives extends our research concerning the chemical and pharmacological optimization of isolated natural products in the search of new drug candidates from brazilian biodiversity.

Molecular orbital calculation for the model compounds of kainoid amino acids, agonists of excitatory amino acid receptors. Does the kainoid C4-substituent directly interact with the receptors?

Kainoid amino acids are agonists of the AMPA/kainate receptors and exhibit highly potent neuroexcitatory activity. From the results of extensive structure-activity relationship studies, we previously postulated that the C4-substituent of the kainoid amino acids interacts with an allosteric site of the glutamate receptor with electron-donating character. In order to investigate the mode of action in more detail, molecular orbital calculation for model compounds of the kainoid were performed. The results indicated that the HOMO energy level of the C4-substituent is involved in the potent neuroexcitatory activity, thus supporting our hypothesis. (C) 2002 Elsevier B.V. Ltd. All rights reserved.