Molecular quantum similarity in QSAR: applications in computer-aided molecular design

La present tesi està centrada en l'ús de la Teoria de Semblança Quàntica per a calcular descriptors moleculars. Aquests descriptors s'utilitzen com a paràmetres estructurals per a derivar correlacions entre l'estructura i la funció o activitat experimental per a un conjunt de compostos. Els estudis de Relacions Quantitatives Estructura-Activitat són d'especial interès per al disseny racional de molècules assistit per ordinador i, en particular, per al disseny de fàrmacs. Aquesta memòria consta de quatre parts diferenciades. En els dos primers blocs es revisen els fonaments de la teoria de semblança quàntica, així com l'aproximació topològica basada en la teoria de grafs. Ambdues teories es fan servir per a calcular els descriptors moleculars. En el segon bloc, s'ha de remarcar la programació i implementació de programari per a calcular els anomenats índexs topològics de semblança quàntica. La tercera secció detalla les bases de les Relacions Quantitatives Estructura-Activitat i, finalment, el darrer apartat recull els resultats d'aplicació obtinguts per a diferents sistemes biològics.

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.

Role of Halogen Bonds in Thyroid Hormone Receptor Selectivity: Pharmacophore-Based 3D-QSSR Studies

Most physiological effects of thyroid hormones are mediated by the two thyroid hormone receptor subtypes, TR alpha and TR beta. Several pharmacological effects mediated by TR beta might be beneficial in important medical conditions such as obesity, hypercholesterolemia and diabetes, and selective TR beta activation may elicit these effects while maintaining an acceptable safety profile, To understand the molecular determinants of affinity and subtype selectivity of TR ligands, we have successfully employed a ligand- and structure-guided pharmacophore-based approach to obtain the molecular alignment of a large series of thyromimetics. Statistically reliable three-dimensional quantitative structure-activity relationship (3D-QSAR) and three-dimensional quantitative structure-selectivity relationship (3D-QSSR) models were obtained using the comparative molecular field analysis (CoMFA) method, and the visual analyses of the contour maps drew attention to a number of possible opportunities for the development of analogs with improved affinity and selectivity. Furthermore, the 3D-QSSR analysis allowed the identification of a novel and previously unmentioned halogen bond, bringing new insights to the mechanism of activity and selectivity of thyromimetics.

Fragment-Based QSAR and Molecular Modeling Studies on a Series of Discodermolide Analogs as Microtubule-Stabilizing Anticancer Agents

Inhibition of microtubule function is an attractive rational approach to anticancer therapy. Although taxanes are the most prominent among the microtubule-stabilizers, their clinical toxicity, poor pharmacokinetic properties, and resistance have stimulated the search for new antitumor agents having the same mechanism of action. Discodermolide is an example of nontaxane natural product that has the same mechanism of action, demonstrating superior antitumor efficacy and therapeutic index. The extraordinary chemical and biological properties have qualified discodermolide as a lead structure for the design of novel anticancer agents with optimized therapeutic properties. In the present work, we have employed a specialized fragment-based method to develop robust quantitative structure - activity relationship models for a series of synthetic discodermolide analogs. The generated molecular recognition patterns were combined with three-dimensional molecular modeling studies as a fundamental step on the path to understanding the molecular basis of drug-receptor interactions within this important series of potent antitumoral agents.

Two-dimensional QSAR studies on arylpiperazines as high-affinity 5-HT(1A) receptor ligands

5-HT(1A) receptor plays an important role in the delayed onset of antidepressant action of a class of selective serotonin reuptake inhibitors. Moreover, 5-HT(1A) receptor levels have been shown to be altered in patients suffering from major depression. In this work, hologram quantitative structure-activity relationship (HQSAR) studies were performed on a series of arylpiperazine compounds presenting affinity to the 5-HT(1A) receptor. The models were constructed with a training set of 70 compounds. The most significant HQSAR model (q(2) = 0.81, r(2) = 0.96) was generated using atoms, bonds, connections, chirality, and donor and acceptor as fragment distinction, with fragment size of 6-9. Predictions for an external test set containing 20 compounds are in good agreement with experimental results showing the robustness of the model. Additionally, useful information can be obtained from the 2D contribution maps.

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).

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.

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.

A quantum chemical study on a set of non-imidazole H(3) antihistamine molecules

Molecular orbital calculations were carried out on a set of 28 non-imidazole H(3) antihistamine compounds using the Hartree-Fock method in order to investigate the possible relationships between electronic structural properties and binding affinity for H3 receptors (pK(i)). It was observed that the frontier effective-for-reaction molecular orbital (FERMO) energies were better correlated with pK(i) values than highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) energy values. Exploratory data analysis through hierarchical cluster (HCA) and principal component analysis (PCA) showed a separation of the compounds in two sets, one grouping the molecules with high pK(i) values, the other gathering low pK(i) value compounds. This separation was obtained with the use of the following descriptors: FERMO energies (epsilon(FERMO)), charges derived from the electrostatic potential on the nitrogen atom (N(1)), electronic density indexes for FERMO on the N(1) atom (Sigma((FERMO))c(i)(2)). and electrophilicity (omega`). These electronic descriptors were used to construct a quantitative structure-activity relationship (QSAR) model through the partial least-squares (PLS) method with three principal components. This model generated Q(2) = 0.88 and R(2) = 0.927 values obtained from a training set and external validation of 23 and 5 molecules, respectively. After the analysis of the PLS regression equation and the values for the selected electronic descriptors, it is suggested that high values of FERMO energies and of Sigma((FERMO))c(i)(2), together with low values of electrophilicity and pronounced negative charges on N(1) appear as desirable properties for the conception of new molecules which might have high binding affinity. 2010 Elsevier Inc. All rights reserved.

A partial least squares regression study with antioxidant flavonoid compounds

The quantitative structure-activity relationship of a set of 19 flavonoid compounds presenting antioxidant activity was studied by means of PLS (Partial Least Squares) regression. The optimization of the structures and calculation of electronic properties were done by using the semiempirical method AMI. A reliable model (r(2) = 0.806 and q(2) = 0.730) was obtained and from this model it was possible to consider some aspects of the structure of the flavonoid compounds studied that are related with their free radical scavenging ability. The quality of the PLS model obtained in this work indicates that it can be used in order to design new flavonoid compounds that present ability to scavenge free radicals.

The Effects of N Helix Deletion and Mutant F29W on the Ca2+ Binding and Functional Properties of Chicken Skeletal Muscle Troponin C

To assess the structural and functional significance of the N helix (residues 3-13) of avian recombinant troponin C (rTnC), we have constructed NHdel, in which residues 1-11 have been deleted, both in rTnC and in the spectral probe mutant F29W (Pearlstone, J. R., Borgford, T., Chandra, M., Oikawa, K., Kay, C. M., Herzberg, O., Moult, J., Herklotz, A., Reinach, F. C., and Smillie, L.B. (1992) Biochemistry 31, 6545-6553). Comparison of the far- and near-UV CD spectra (±Ca2+) of F29W and F29W/ NHdel and titration of the Ca2+-induced ellipticity and fluorescence changes indicates that the deletion has little effect on the global fold of the molecule but reduces the Ca2+ affinity of the N domain, but not the C domain, by 1.6-1.8-fold. Comparisons of the mutants NHdel, F29W, and F29W/NHdel with rTnC have been made using several functional assays. In reconstituted troponin-tropomyosin actomyosin subfragment 1 and myofibrillar ATPase systems, both F29W and NHdel have significantly reduced Ca2+-activated enzymic activities. These effects are cumulative in the double mutant F29W/ NHdel. On the other hand, maximal isometric tension development in Ca2+-activated reconstituted skinned fibers is not affected with F29W and NHdel, although the Ca2+ sensitivity of NHdel in this system is markedly reduced. We conclude that both mutations, NHdel and F29W, are functionally deleterious, possibly affecting interactions of the N domain with troponin I and/or T.

Esparfloxacino e as fluorquinolonas: Uma revisao de sua classificacao, atividade antibacteriana, mecanismo de acao, caracteristicas farmacocineticas, tolerabilidade e perspectivas farmaceuticas

A revision was accomplished in the literature on the quinolones, antibacterial class that presents wide action spectrum, focusing, mainly, the sparfloxacin, third generation fluorquinolone which has potente activity against Gram positive organisms, including Streptococcus pneumoniae and methicillin-resistant Staphylococcus aureus (MRSA), Gram negative organisms, Legionelia spp, Mycoplasma spp, Chlamydia spp and Mycobacterium spp, including multidrug resistant organisms.

Conserved sequences in the β subunit of archaeal and eukaryal translation initiation factor 2 (elF2), absent from elF5, mediate interaction with elF2γ

The eukaryotic translation initiation factor 2 (eIF2) binds the methionyl-initiator tRNA in a GTP-dependent mode. This complex associates with the 40 S ribosomal particle, which then, with the aid of other factors, binds to the 5' end of the mRNA and migrates to the first AUG codon, where eIF5 promotes GTP hydrolysis, followed by the formation of the 80 S ribosome. Here we provide a comparative sequence analysis of the β subunit of eIF2 and its archaeal counterpart (aIF2β). aIF2β differs from eIF2β in not possessing an N-terminal extension implicated in binding RNA, eIF5 and eIF2B. The remaining sequences are highly conserved, and are shared with eIF5. Previously isolated mutations in the yeast eIF2β, which allow initiation of translation at UUG codons due to the uncovering of an intrinsic GTPase activity in eIF2, involve residues that are conserved in aIF2β, but not in eIF5. We show that the sequence of eIF2B homologous to aIF2β is sufficient for binding eIF2γ, the only subunit with which it interacts, and comprises, at the most, 78 residues, eIF5 does not interact with eIF2γ, despite its similarity with eIF2β, probably because of a gap in homology in this region. These observations have implications for the evolution of the mechanism of translation initiation.

In vitro cytotoxicity of some natural and semi-synthetic isocoumarins from Paepalanthus bromelioides

Numerous natural compounds have a potential for therapeutic applications, but may have to be chemically modified to alter toxic side effects. We investigated structural parameters that could affect the cytotoxicity of isocoumarins similar to 9,10-dihydroxy-5,7-dimethoxy-1H-naphtho(2,3c)pyran-1-one (paepalantine 1). Paepalantine 1 has antimicrobial activity, as well as significant in vitro cytotoxic effects in the McCoy cell line. Two other natural and two semi-synthetic isocoumarins with similar structures obtained from the capitula of Paepalanthus bromelioides were tested on the same cell line by the neutral red assay. Substitution of the 9 and/or 10-OH group made these compounds less cytotoxic.