Combining Free Access Methodologies to Build 2D. QSAR Models: A Case Study for Piperazinenitroimidazole

Transcription

1 Research Article imedpub Journals Journal of Informatics and Data Mining Abstract Combining ree Access Methodologies to Build D QSAR Models: A Case Study for Piperazine-nitroimidazole Analogues against S. aureus Quantitative Structure-Activity Relationships (QSAR) approaches involve the search for mathematic models that correlates a set of biological properties (or activities) from a series of analogues with their molecular features, also called molecular descriptors. These studies seek to understand how molecular features, previously identified, can influence the biological responses, allowing the design and future synthesis of potentially most active compounds. In this work we present the use of free and open source softwares to data processing in a D QSAR perspective, using a set of metronidazole-piperazine analogs with activity against S. aureus, a relatively common species of bacteria that has represented a serious public health problem, due their mechanism of resistance to drugs. The combination of different methodologies conduced to a robust and predictive QSAR model. These strategies can be used in programs of medicinal organic synthesis and drug design. Keywords: QSAR, ree softwares, Piperazine-nitroimidazole analogues, S. aureus Edilson B. Alencar ilho and Aline A. Santos Collegiate of Pharmaceutical Sciences, ederal University of San rancisco Valley, Petrolina, Pernambuco, Brazil, Corresponding author: Edilson B Alencar ilho edilson.beserra@univasf.edu.br Collegiate of Pharmaceutical Sciences, ederal University of San rancisco Valley, Petrolina, Pernambuco, Brazil, Tel: Citation: ilho EBA, Santos AA. Combining ree Access Methodologies to Build D QSAR Models: A Case Study for Piperazinenitroimidazole Analogues against S. aureus. J Inform Data Min., 1:1. Received: September 11, 015; Accepted: October 4, 015; Published: October 9, 015 Introduction In the context of bioactive compounds, Quantitative Structure- Activity Relationships (QSAR) approaches involves the search for mathematic models that correlates a set of biological properties (or activities) from a series of analogues with their molecular features, also called molecular descriptors. These studies seek to understand how molecular features previously identified can influence the biological responses, allowing the design and future synthesis of potentially most active compounds. QSAR s has in general the following steps: Selection of analogues compounds with particular activity (important to a better comprehensive structure-activity relationships); Determination of molecular descriptors by experimental or computational molecular modelling methods; Application of a variable selection methodology to descriptors; Construction and validation of mathematic models by statistical/chemometric tools; Interpretation of final model. Multivariate regression is the core of QSAR analysis. It is based on obtaining a model Y=bX+e, where Y is a column matrix with values of molecular property or activity (dependent variable), X is an I x J matrix with values of J molecular descriptors of I compounds. [1] The J x 1 vector b contains the calibration regression coefficients and e represents the I x 1 error vector. [1] The equation can be obtained by Multiple Linear Regression (MLR), [] but when there is a great number of variables it is not adequate. Thus, techniques as Principal Component Regression (PCR) and Partial Least Squares (PLS), based on Principal Component Analysis (PCA), can be used [3-5]. The kind of molecular descriptors used in QSAR modelling provides the different nomenclatures as D, 3D, 4D QSAR among others. The most cited are 3D and 4D methodologies and they are based on alignment of a set of compounds in a common spatial arrangement, defined by a common portion of the molecules, which are involved in a three-dimensional grid to generation of descriptors via interaction potentials with selected probes (specific functional groups models, charged or not). 4D differentiates by the uses of conformational sampling profile of each compound [6]. On the other hand, D methodologies use general molecular or atomic properties from quantum chemical or molecular mechanic calculations, experimental parameters (as LogP, pka), fragmental index and other numerical values, as long as they do 1

2 not represent interaction potentials with probes. In this context, E-DRAGON online plataform represents a free remote version of DRAGON software, which is dedicated to calculation of molecular descriptors of various kinds. [7-9] These descriptors can be used to evaluate molecular structure-activity relationships, similarity analysis and screening of molecule databases. [7,8,10] Initially, molecular geometries have to be converted to appropriate format, usually sdf, smiles or mol formats. Then the file can be loaded on the page, providing more than 1,660 descriptors, divided into 0 logical blocks. [9,10] Several topological and geometrical descriptors can be calculated, besides atom type, functional group and fragment counts. [9,10] Variable selection is other important step in QSAR analysis. Ordered predictors Selection (OPS) is an algorithm developed and published in 009 by Teófilo and coworkers, and can be downloaded by previous registration in the software page [1,11]. This methodology is based on sort the variables from an initial informative vector, followed by a systematic investigation of Partial Least Squares (PLS) regression models with the aim of finding the most relevant set of variables by comparing the crossvalidation parameters of the models systematically obtained [1]. Considering our continuous work in Cheminformatics area, this paper presents a D QSAR study of a set of compounds with bactericidal activity, using free and open source softwares. The adopted methodologies can be used in other datasets to optimize the synthesis of bioactive compounds, in drug design programs. Methodology Because of our interest on antimicrobial and antiparasitic compounds, specially those containing nitroimidazole moiety, the set of compounds was selected from a paper of Wang and coworkers [1] (igure 1) and corresponds to a series of metronidazole derivatives containing piperazine skeleton with antibacterial activity against Staphylococcus aureus. This is a relatively common bacteria which can cause various types of infections, from a simple infection (pimples, boils) to more serious health problems (meningitis, pneumonia, endocarditis). A serious aspect of this microorganism is the acquired resistance of some strains to certain antibiotics such as vancomycin, methicillin, among others [1]. Thus, the search for new agents for the treatment of S. aureus infections is an attractive and important field of research. The original activity was considered as a Minimum Inhibitory Concentration (MIC) in µg/ml, by MTT method. [1] or QSAR purpose, the values were converted to molar units and expressed as LogMIC (pmic) (igure 1). or calculation of descriptors, molecules were submitted to a geometry optimization using PM3 semiempirical method in Arguslab 4.0 free software [13]. The optimized geometries were converted to sdf format using free OpenBabel, [14] an appropriate extension to load the structures on the E-Dragon platform. This platform is easy to use and quickly generates a txt extension with a matrix of 1,666 D descriptors, as described on introduction. These descriptors were organized in an Excel worksheet J x I and then pasted in a text file. We opted to not choose training and test sets [15]. All the 15 molecules were used in QSAR analysis considering the small number of compounds with activity available in the literature, so as not to lose important information in the modeling [15]. In this way, this study represents an initial direction for optimization of metronidazole-piperazine analogues. Another txt file contain a column with pmic values were saved. The two text files were opened in QSAR Modelling free software [16] to run OPS algorithm. In this software there is an option to autoscaling variables, important when we work with different units of measurement [4]. The parameters for the OPS procedure were: Window of 10 initial variables, among the most correlated; increment of 5 variables until a maximum value of 30%. SPRESS was chosen as the sort parameter, the model with the lowest value corresponding to the selected variables. SPREES represents the standard deviation of PRESS, which is inserted in the following root [1]: SPRESS = pred exp ( yi ) (1) I NVL 1 whereicorresponds to the number of samples used and NVL pred represents the number of latent variables for PLS model; y i is exp the predicted value and y i the experimental value of biological activity. It is important to consider that the variable selection process with QSAR modeling can be done with more than one type of validation parameter, beyond SPRESS, as the coefficient of determination of leave-one-out cross validation,q [17]. During this process, depending on the number of initially selected descriptors, the generated matrix may be subjected to new selections, up to a minimal set of descriptors that meet the model validation parameters and have the physical sense in the structure-activity relationships. Thus, the development of a QSAR model also depends of a little patience and experience of modeler to identify descriptors that meet the criteria of robustness and model predictability, as well enable consistent interpretation of it. If the ratio between the numbers of samples/descriptors is <5, a multivariate regression model as PCR or PLS has to be used to generate the coefficients of final regression model. They are implemented on QSAR modeling software. If >5, MLR technique can be used, with the advantage that preserves the original variables in final model, more easily to interpret. The two principal quality parameters are the coefficient of determination of calibration (R ) and coefficient of determination of leave-one-out cross validation (Q ). In R evaluation, all the samples are considered to construct the model before prediction (Equation ). In Leave-one-out, one sample is taken at a time and a regression model is constructed, which is used to predict the activity of the sample taken. In these determinations all compounds are used. The values predicted for all samples are used in calculation of quality parameter Q (Equation 3). Literature reports that acceptable values of these parameters are: Q > 0.5; R > 0.6. [17,0] exp cal i i exp mean i R = 1 ( y y ) ( y ) ()

3 Q 1n n 3n 4n igure 1 Table 1 Biological activity against S. aureus in pmic and E-Dragon molecular descriptors selected after successive runs of OPS algorithm. Composto pmic G(N N) GATS8e RD100e H0e 1n n n n n n n n n n n n n n n = 1 H 3CO N 5n 6n 7n OCH 3 8n NO Compounds obtained from paper of Wang and coworkers (014). ( y y ) exp pred i i exp mean i ( y ) N 3C H 3 C R = H 3CO N N R *cal represents the predicted values using model of calibration 9n CH 3 10n 11n 1n 13n 14n 15n (3) constructed with all samples; mean represents the mean of experimental values of dependent variable. Lack of correlation between independent variables in a QSAR is also important. In general is usually correlations >0.6 are avoided [18]. The free software Build QSAR [19,0] is an way to calculate correlation coefficients between descriptors, as well generate MLR or PCR models. Considering the number of final selected variables, in this work we opted for MLR modeling using this program. Results and Discussion: ully optimized geometries of compounds are presented in igure. Molecular descriptors generated by E-Dragon plataform and selected using OPS algorithm can be visualized in Table 1 and igure. After variable selection, MLR model was constructed using Build QSAR. The regression vector obtained in this step is: pmic= 0.749(G(N...N)) (GATS8e) (RD100e) 4.360(H0e) (4) The coefficients of determination were R = 0.87 and Q = 0.77, indicating that a robust and predictive model was obtained. The standard error of prevision was s = 0.5. Value of SPRESS = 0.68 and the correlation matrix can be visualized in (Table ). Together, this results show a good QSAR model for this class of new anti S. aureus agents. In addition to allow the prediction of compounds activity, in a QSAR model is important to know the molecular descriptors used, analyzing the particular characteristics, seeking to rationalize the physical sense for the design of more active molecules. However, it is important to note that the overall effect on a molecule is given by the sum of the contributions of each selected descriptor. Therefore, we can analyze can analyze the influence of the descriptor in the activity individually but not neglect the concomitant influence of other descriptors. G(N...N): This descriptor is related as the sum of the geometric distance between nitrogen atoms [1]. The greater its numerical value, the lower the activity of molecules. Observing the structures, we conclude that molecules with axial nitroimidazole, which leads to smaller distances N...N, are more active. In other words, most compact conformational geometries showed to be important to biological activity; GATS8e: Related to the "autocorrelation Geary" at a distance of 8 bonds, weighted by electronegativity [1]. Thus, greater electronegativity difference between two atoms at a distance of 8 bonds, leads to greater the value of the descriptor. The higher values of this descriptor, the lower activity. Electronegative atoms in para (4) position of phenyl ring leads to a greater value of the descriptor, decreasing the activity (compounds n, 3n and 15n are the less active, with presence of halogens on this position);. RD100e: Related to the Radial Distribution unction weighted by electronegativity and radius 10.0 A [1]. Probability of finding an atom in a radius of spherical volume 10. Higher number of atoms in a molecular volume increases its value. Considering regression vector, higher numerical value of this descriptor, increases the 3

4 igure 1n 4n n 5n 7n 8n 9n 10n 11n 1n 13n 14n 15n 3D optimized geometries of piperazine-metronidazole derivatives. Table Correlation matrix between selected descriptors. Descriptor G(N...N) GATS8e RD100e H0e G(N...N) GATS8e RD100e H0e activity. Molecules containing phenyl groups and presence of more than one phenyl ring have a higher value of this descriptor. This may be related to π π stacking interactions, which have fundamental importance in drug-target pharmacodynamics. This 3n 6n can be visualized on compound 13n, the most active of series, which has two phenyl groups without substituents. On the other hand, compound 15n, even having two phenyl groups has one of the lowest activities. Within the discussions about GATS8e, this molecule has two halogens in para position, what may explain the decreased effect on activity. H0e: H autocorrelation without distance, weighted by electronegativity. Denotes the influence of electronegative atoms in molecular form [1]. The more electronegative atoms in the molecules, the higher the value of this descriptor. Regression vector shows that the higher values of this descriptor, the lower activity. Molecules containing more electronegative atoms tend to be less active. In addition to information of GATS8e descriptor, these electronegative atoms can be attributed to halogens (n, 3n, 7n, 15n). Compounds 8n and 13n are the most actives and have no chlorine or fluorine electronegative atoms. Conclusions In accordance with obtained results, combination of free software s presented here conduced to a robust and predictive QSAR model. These data processing strategies can be used for drug design and lead compounds optimization, especially in academical level, assisting the efforts of the chemical synthesis and pharmacology. In relation to the explored data set, four descriptors were very important to explain the biological activity of the compounds and the obtained equation can be used to make predictions of other piperazine-metronidazole analogues ideally designed based on these characteristics. Molecules with a more "compact" conformational behavior, with phenyl and non-electronegative substituents (halogens) especially that not occupy the para positions, should lead to more active analogues against S. aureus. In addition, the generated equation can be used to predict the activity of new analogues compounds, saving synthetic efforts. Acknowledgments The authors are grateful to Research Support oundation of the State of Pernambuco (ACEPE), National Counsel of Technological and Scientific Development (CNPq) and ederal University of San rancisco Valley (UNIVAS), for financial support. 4

5 References 1 Teófilo R, Martins J P, erreira M M C (009) Sorting variables by using informative vectors as a strategy for feature selection in multivariate regression, Journal of Chemometrics, 3: Uyanic G K, Güler N (013) A study on multiple linear regression analysis.4 th International Conference on New Horizons in Education. Procedia - Social and Behavioral Sciences,106: Wold S, Esbensen K, Geladi P (1987) Principal Component Analysis. Chemometrics and Intelligent Laboratory Systems, : Beebe K R, Pell R J, Seasholtz M B (1998)Chemometrics: A Pratical Guide. Wiley, New York. 5 Wold S, Sjöström M, Eriksson L, (001) PLS-regression: a basic tool of chemometrics. Chemometrics and Intelligent Laboratory Systems, 58: Martins J P, Barbosa E G, Pasqualoto K M, erreira M M C (009) LQTA-QSAR: A New 4D-QSAR Methodology, Journal Chemical Information and Modeling, 49: Todeschini R, Consonni V (009) Molecular Descriptors for Chemoinformatic.WILEY-VCH,Weinheim, Germany. 8 Tetko I V (005) Computing chemistry on the web, Drug Discovery Today, 10: Virtual Computational Chemistry Laboratory. 10 Tetko I V, Gasteiger J, Todeschini R, Mauri A, Livingstone D, et al. (005) Virtual computational chemistry laboratory- design and description, Journal of Computational Aided Molecular Design, 19: Laboratório de Quimiometria Teórica e Aplicada. 1 Wang S,Yin Y, Qiao, Wu X, Sha S, et al. (014) Synthesis, molecular docking and biological evaluation ofmetronidazole derivatives containing piperazine skeleton aspotential antibacterial agents. Bioorganic & Medicinal Chemistry, : Thompson M A (014) ArgusLab 4.0; Planaria Software LLC: Seattle, WA. 14 O'Boyle N M, Banck M, James C A, Morley C, Vandermeersch T et al. (011) Open Babel: An open chemical toolbox. Journal of Cheminformatics, 3: Silva J B P, Ramos M N, Neto BB, Melo S, alcão E P S (008) Catanho, M.T.J.A, Quantitative Structure-Activity Relationships (QSAR) of 4-Amino-,6-diarylpyrimidine-5-carbonitriles with Antiinflammatory Activity, Journal of the Brazilian Chemical Society, 19: Martins J P A, erreira M M C (013) QSAR modeling: um novo pacote computacional open source para gerar e validar modelos QSAR. Química Nova, 36: Kiralj R, erreira M M C (009) Basic validation procedures for regression models in QSAR and QSPR studies: theory and application. Journal of the Brazilian Chemical Society 0: erreira M M C, Montanari C A, Gaudio A C (00) Seleção de variáveis em QSAR. Química Nova,5: De Oliveira D B, Gaudio A C (001) Build QSAR: A New Computer Program for QSAR Analysis. Quantitative Structure-Activity Relationships, 19: Alencar ilho E B, Weber K C, Vasconcellos M L A A (014) Selection of D/3D molecular descriptors and QSAR modeling of aromatic Morita Baylis Hillman adducts with leishmanicidal activities, Medicinal Chemistry Research,3: Strand Life Science P. Ltda, Theory Descriptors (007) 5