Two new atom centered fragment descriptors and scoring function enhance classification of antibacterial activity

Classification of pharmacologic activity of a chemical compound is an essential step in any drug discovery process. We develop two new atom-centered fragment descriptors (vertex indices) - one based solely on topological considerations without discriminating atomor bond types, and another based on topological and electronic features. We also assess their usefulness by devising a method to rank and classify molecules with regard to their antibacterial activity. Classification performances of our method are found to be superior compared to two previous studies on large heterogeneous data sets for hit finding and hit-to-lead studies even though we use much fewer parameters. It is found that for hit finding studies topological features (simple graph) alone provide significant discriminating power, and for hit-to-lead process small but consistent improvement can be made by additionally including electronic features (colored graph). Our approach is simple, interpretable, and suitable for design of molecules as we do not use any physicochemical properties. The singular use of vertex index as descriptor, novel range based feature extraction, and rigorous statistical validation are the key elements of this study.

Quantitative relationships between molecular structures, environmental temperatures and octanol-air partition coefficients of PCDD/Fs

A concise quantitative model that incorporates information on both environmental temperature M and molecular structures, for logarithm of octanol-air partition coefficient (K-OA) to base 10 (logK(OA)) of PCDDs, was developed. Partial least squares (PLS) analysis together with 14 quantum chemical descriptors were used to develop the quantitative relationships between structures, environmental temperatures and properties (QRSETP) model. It has been validated that the obtained QRSETP model can be used to predict logK(OA) of other PCDDs. Molecular size, environmental temperature (T), q(+) (the most positive net atomic charge on hydrogen or chlorine atoms in PCDD molecules) and E-LUMO (the energy of the lowest unoccupied molecular orbital) are main factors governing logK(OA) of PCDD/Fs under study. The intermolecular dispersive interactions and thus the size of the molecules play a leading role in governing logK(OA). The more chlorines in PCDD molecules, the greater the logK(OA) values. Increasing E-LUMO values of the molecules leads to decreasing logK(OA) values, implying possible intermolecular interactions between the molecules under study and octanol molecules. Greater q(+) values results in greater intermolecular electrostatic repulsive interactions between PCDD and octanol molecules and smaller logK(OA) values. (C) 2002 Elsevier Science B.V. All rights reserved.

Peripheries of molecular projection in three-dimensional space and studies of quantitative structure-activity relationships

A new index, i.e., the periphery representation of the projection of a molecule from 3D space to a 2D plane is described. The results, correlation with toxicity of substituted nitrobenzenes, obtained by using periphery descriptors are much better than that obtained by using the areas (i.e., shadows) of projections of the compounds. Even better results were achieved by using the combination of periphery descriptors and the projections areas as well as the indicated variable K reflecting the action of group NO position on the benzene ring.

Quantitative structure-hydrophobicity relationships of molecular fragments and beyond

Quantitative structure-property relationship (QSPR) models were firstly established for the hydrophobic substituent constant (πX) using the theoretical descriptors derived solely from electrostatic potentials (EPSs) at the substituent atoms. The descriptors introduced are found to be related to hydrogen-bond basicity, hydrogen-bond acidity, cavity, or dipolarity/polarizability terms in linear solvation energy relationship, which endows the models good interpretability. The predictive capabilities of the models constructed were also verified by rigorous Monte Carlo cross-validation. Then, eight groups of meta- or para- disubstituted benzenes and one group of substituted pyridines were investigated. QSPR models for individual systems were achieved with the ESP-derived descriptors. Additionally, two QSPR models were also established for Rekker's fragment constants (foct), which is a secondary-treatment quantity and reflects average contribution of the fragment to logP. It has been demonstrated that the descriptors derived from ESPs at the fragments, can be well used to quantitatively express the relationship between fragment structures and their hydrophobic properties, regardless of the attached parent structure or the valence state. Finally, the relations of Hammett σ constant and ESP quantities were explored. It implies that σ and π, which are essential in classic QSAR and represent different type of contributions to biological activities, are also complementary in interaction site.

Development, implementation and application of electronic structural descriptors to the analysis of the chemical bonding, aromaticity and chemical reactivity

En la literatura sobre mecànica quàntica és freqüent trobar descriptors basats en la densitat de parells o la densitat electrònica, amb un èxit divers segons les aplicacions que atenyin. Per tal de que tingui sentit químic un descriptor ha de donar la definició d'un àtom en una molècula, o ésser capaç d'identificar regions de l'espai molecular associades amb algun concepte químic (com pot ser un parell solitari o zona d'enllaç, entre d'altres). En aquesta línia, s'han proposat diversos esquemes de partició: la teoria d'àtoms en molècules (AIM), la funció de localització electrònica (ELF), les cel·les de Voroni, els àtoms de Hirshfeld, els àtoms difusos, etc. L'objectiu d'aquesta tesi és explorar descriptors de la densitat basats en particions de l'espai molecular del tipus AIM, ELF o àtoms difusos, analitzar els descriptors existents amb diferents nivells de teoria, proposar nous descriptors d'aromaticitat, així com estudiar l'habilitat de totes aquestes eines per discernir entre diferents mecanismes de reacció.

Caracterização morfológica, citogenética e molecular de uma população de tangerineiras híbridas de 'Clementina fina' (Citrus clementina Hort. ex Tan.) e 'Montenegrina' (Citrus deliciosa Ten.)

A produção citrícola se encontra dispersa por todos os continentes e no Brasil, os citros são a produção frutícola de maior volume de produção. A produção de citros de mesa, como as tangerinas, possibilita ao produtor obter maior valor pelo seu produto. O mercado consumidor é ávido por novas variedades e para tanto, um programa de melhoramento deve estar sempre em busca de genótipos que atendam ao mercado consumidor, bem como a cadeia produtiva. Na Estação Experimental Agronômica da Universidade Federal do Rio Grande do Sul, está localizada uma população de tangerineiras híbridas oriundas do cruzamento da tangerineira ‘Clementina Fina’ (Citrus clementina Hort. ex Tan.) e ‘Montenegrina’ (Citrus deliciosa Ten.) a qual foi caracterizada neste estudo, avaliando-se características morfológicas de acordo com os descritores propostos pelo International Board for Plant Genetic Resources, além da identificação da época de maturação, viabilidade de pólen, número cromossômico e caracterização molecular, utilizando marcadores do tipo microssatélites. Através da análise morfológica foi possível distinguir todas as 96 plantas avaliadas, porém não foi possível agrupar a F1 em grupos distintos de cada um dos genitores. A época de maturação de frutos das plantas se concentra entre a primeira quinzena de abril até a primeira quinzena de agosto. Todas as plantas analisadas apresentaram um alto grau de viabilidade de pólen, variando entre 79,04 e 98,08 %. Todas as plantas avaliadas são diplóides com um número cromossômico de 2n=18. Utilizando 12 pares de primers de microssatélites foi possível diferenciar 90 acessos do estudo, e agrupar a F1 em indivíduos mais próximos do genitor feminino e do genitor masculino. O PIC (Conteúdo de Informação de Polimorfismo) dos primers variou de 0,27 a 0,65. Não foi possível estabelecer uma relação entre a caracterização utilizando marcadores morfológicos e a caracterização utilizando marcadores moleculares.

Molecular physicochemical parameters predicting antioxidant activity of Brazilian natural products

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

Developing descriptors to predict mechanical properties of nanotubes

Descriptors and quantitative structure property relationships (QSPR) were investigated for mechanical property prediction of carbon nanotubes (CNTs). 78 molecular dynamics (MD) simulations were carried out, and 20 descriptors were calculated to build quantitative structure property relationships (QSPRs) for Young's modulus and Poisson's ratio in two separate analyses: vacancy only and vacancy plus methyl functionalization. In the first analysis, C N2/CT (number of non-sp2 hybridized carbons per the total carbons) and chiral angle were identified as critical descriptors for both Young's modulus and Poisson's ratio. Further analysis and literature findings indicate the effect of chiral angle is negligible at larger CNT radii for both properties. Raman spectroscopy can be used to measure CN2/C T, providing a direct link between experimental and computational results. Poisson's ratio approaches two different limiting values as CNT radii increases: 0.23-0.25 for chiral and armchair CNTs and 0.10 for zigzag CNTs (surface defects <3%). In the second analysis, the critical descriptors were CN2/CT, chiral angle, and MN/CT (number of methyl groups per total carbons). These results imply new types of defects can be represented as a new descriptor in QSPR models. Finally, results are qualified and quantified against experimental data. © 2013 American Chemical Society.

Caracterizações morfológica e molecular de acessos de pimenta (Capsicum chinense Jaqc.)

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Protease-inhibiting, molecular modeling and antimicrobial activities of extracts and constituents from Helichrysum foetidum and Helichrysum mechowianum (compositae)

Background Helichrysum species are used extensively for stress-related ailments and as dressings for wounds normally encountered in circumcision rites, bruises, cuts and sores. It has been reported that Helichysum species are used to relief abdominal pain, heart burn, cough, cold, wounds, female sterility, menstrual pain. Results From the extracts of Helichrysum foetidum (L.) Moench, six known compounds were isolated and identified. They were 7, 4′-dihydroxy-5-methoxy-flavanone (1), 6′-methoxy-2′,4, 4′-trihydroxychalcone (2), 6′-methoxy-2′,4-dihydroxychalcone -4′-O-β-D-glucoside (3), apigenin (4), apigenin-7-O-β-D-glucoside (5), kaur-16-en-18-oic acid (6) while two known compounds 3,5,7-trihydroxy-8-methoxyflavone (12), 4,5-dicaffeoyl quinic acid (13) together with a mixture of phytosterol were isolated from the methanol extract of Helichrysum mechowianum Klatt. All the compounds were characterized by spectroscopic and mass spectrometric methods, and by comparison with literature data. Both extracts and all the isolates were screened for the protease inhibition, antibacterial and antifungal activities. In addition, the phytochemical profiles of both species were investigated by ESI-MS experiments. Conclusions These results showed that the protease inhibition assay of H. foetidum could be mainly attributed to the constituents of flavonoids glycosides (3, 5) while the compound (13) from H. mechowianum contributes to the stomach protecting effects. In addition, among the antibacterial and antifungal activities of all the isolates, compound (6) was found to possess a potent inhibitor effect against the tested microorganisms. The heterogeneity of the genus is also reflected in its phytochemical diversity. The differential bioactivities and determined constituents support the traditional use of the species. Molecular modelling was carried out by computing selected descriptors related to drug absorption, distribution, metabolism, excretion and toxicity (ADMET).

A novel set of ss-N-biaryl ether sulfonamide hydroxamates as potential MMPs inhibitors: Molecular dynamics simulations and molecular properties evaluation

Matrix metalloproteinases (MMPs) constitute a family of zinc-dependent proteases involved in the extracellular matrix degradation. MMP-2 and MMP9 are overexpressed in several human cancer types, including melanoma, thus the development of new compounds to inhibit MMPs' activity is desirable. Molecular dynamic simulation and molecular properties calculations were performed on a set of novel beta-N-biaryl ether sulfonamide-based hydroxamates, reported as MMP-2 and MMP-9 inhibitors, for providing data to develop an exploratory analysis. Thermodynamic, electronic, and steric descriptors have significantly discriminated highly active from moderately and less active inhibitors of MMP-2 whereas apparent partition coefficient at pH 1.5 was also significant for the MMP-9 data set. Compound 47 was considered an outlier in all analysis, indicating the presence of a bulky substituent group in R3 is crucial to this set of inhibitors for the establishment of molecular interactions with the S1 subsite of both enzymes, but there is a limit. (C) 2012 Wiley Periodicals, Inc.

Molecular Surface Features in Modeling the HIV-1 RT Inhibitory Activity of 2-(2,6-Disubstituted phenyl)-3-(substituted pyrimidin-2-yl)-thiazolidin-4-ones

The human immunodeficiency virus-1 reverse transcriptase inhibitory activity of 2-(2,6-disubstituted phenyl)-3-(substituted pyrimidin-2-yl)-thiazolidin-4-ones have been analyzed using combinatorial protocol in multiple linear regression (CP-MLR) with several electronic and molecular surface area features of the compounds obtained from Molecular Operating Environment (MOE) software. The study has indicated the role of different charged molecular surface areas in modeling the inhibitory activity of the compounds. The derived models collectively suggested that the compounds should be compact without bulky substitutions on its peripheries for better HIV-1 RT inhibitory activity. It also emphasized the necessity of hydrophobicity and compact structural features for their activity. The scope of the descriptors identified for these analogues have been verified by extending the dataset with different 2-(disubstituted phenyl)-3-(substituted pyridin-2-yl)-thiazolidin-4-ones. The joint analysis of extended dataset highlighted the information content of identified descriptors in modeling the HIV-1 RT inhibitory activity of the compounds.

CP-MLR Directed QSAR Studies on the Antimycobacterial Activity of Functionalised Alkenols - Topological Descriptors in Modeling the Activity

The antimycobacterial activity of nitro/ acetamido alkenol derivatives and chloro/ amino alkenol derivatives has been analyzed through combinatorial protocol in multiple linear regression (CP-MLR) using different topological descriptors obtained from Dragon software. Among the topological descriptor classes considered in the study, the activity is correlated with simple topological descriptors (TOPO) and more complex 2D autocorrelation descriptors (2DAUTO). In model building the descriptors from other classes, that is, empirical, constitutional, molecular walk counts, modified Burden eigenvalues and Galvez topological charge indices have made secondary contribution in association with TOPO and / or 2DAUTO classes. The structure-activity correlations obtained with the TOPO descriptors suggest that less branched and saturated structural templates would be better for the activity. For both the series of compounds, in 2DAUTO the activity has been correlated to the descriptors having mass, volume and/ or polarizability as weighting component. In these two series of compounds, however, the regression coefficients of the descriptors have opposite arithmetic signs with respect to one another. Outwardly these two series of compounds appear very similar. But in terms of activity they belong to different segments of descriptor-activity profiles. This difference in the activity of these two series of compounds may be mainly due to the spacing difference between the C1 (also C6) substituents and rest of the functional groups in them.

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.

Flexibility of short-strand RNA in aqueous solution as revealed by molecular dynamics simulation:are A′-RNA and A-RNA distinct conformational structures?

We use molecular dynamics simulations to compare the conformational structure and dynamics of a 21-base pair RNA sequence initially constructed according to the canonical A-RNA and A'-RNA forms in the presence of counterions and explicit water. Our study aims to add a dynamical perspective to the solid-state structural information that has been derived from X-ray data for these two characteristic forms of RNA. Analysis of the three main structural descriptors commonly used to differentiate between the two forms of RNA namely major groove width, inclination and the number of base pairs in a helical twist over a 30 ns simulation period reveals a flexible structure in aqueous solution with fluctuations in the values of these structural parameters encompassing the range between the two crystal forms and more. This provides evidence to suggest that the identification of distinct A-RNA and A'-RNA structures, while relevant in the crystalline form, may not be generally relevant in the context of RNA in the aqueous phase. The apparent structural flexibility observed in our simulations is likely to bear ramifications for the interactions of RNA with biological molecules (e.g. proteins) and non-biological molecules (e.g. non-viral gene delivery vectors). © CSIRO 2009.