QUANTITATIVE STRUCTURE ACTIVITY RELATIONSHIP STUDIES OF SULFADIAZINE- SALICYLALDEHYDE SCHIFF BASE AS AN ANTI-MYCOBACTERIAL AGENTS
HTML Full TextQUANTITATIVE STRUCTURE ACTIVITY RELATIONSHIP STUDIES OF SULFADIAZINE- SALICYLALDEHYDE SCHIFF BASE AS AN ANTI-MYCOBACTERIAL AGENTS
Prabha Thangavelu *, Madhavan Ramakrishnan, Kowsalya Mani, Kirubanithi Somasundram, Dinesh Srinivas and Sivakumar Thangavel
Department of Pharmaceutical Chemistry, Nandha College of Pharmacy, Koorapalayam Pirivu, Pitchandam Palayam Post, Erode - 638052, Tamilnadu, India.
ABSTRACT: Quantitative structure activity relationship is among the most widely used computational technology for analogue-based drug design. A molecular modelling approach using sulfadiazine- salicylaldehyde Schiff base analogue as anti-mycobacterial activity from recently reported literature were taken and was designed using MOE 2009.10. Sulfanilamide sulfa drugs have been shown to inhibit dihydropteroate synthetase. In order to develop a pharmacophoric model for this inhibition, quantitative structure-activity relationship (QSAR) for sulfa drugs active against DHPS have been studied. Multiple linear regression analysis was performed to derive quantitative structure activity relationship models which were further evaluated internally as well as externally for the prediction of activity. Accurate minimum inhibitory concentrations (µM) were collected for 19 analogs, and other parameters, such as partition coefficients, molar refractivity, dipole moment, electrostatic interactions, lowest unoccupied molecular orbital (LUMO), highest unoccupied molecular orbital (HOMO), topological polar surface area were calculated. A data set of 19 analogs, all having a common sulphadiazine moiety, provided a cross-validated correlation coefficient (r2) value of 0.5998 and root mean square error value of 0.3829. The resulting QSAR pharmacophore model generated from the present study should be useful in the design of a similar group of more potent substituted compounds of Schiff base of sulpha drugs analogues as anti-mycobacterial agents.
Keywords: |
QSAR, Sulphadiazine- salicylaldehyde Schiff base, Molecular descriptors and Anti-mycobacterial activity
INTRODUCTION: Tuberculosis (TB), an airborne infectious disease caused predominantly by Mycobacterium tuberculosis (MTB), is a global health problem and a leading cause of death among adults in the developing world 1. TB is a global emergency and one of the most important killing infectious diseases worldwide today.
The recent statistics reveal that one third of the world’s population is currently infected, and more than 5,000 people die from TB every day.
As per the WHO data for Southeast Asia and Africa, there is an alarming rise in new TB cases. Indeed, with the enormous global burden of TB and the alarming rise in the number of clinical isolates displaying drug resistance. The increasing global TB burden its health and socio-economic impacts, the presence of drug resistance, etc. leads to the search for newer therapeutic agents with divergent and unique structure and with a mechanism of action possibly different from that of existing drugs are urgently required to tackle the menace 2, 3, 4, 5.
An MTB strain is limiting the choice of anti-tubercular antibiotics to toxic, injectable second-line antibiotics. Therefore, there is an interest in search of alternative antibiotics used as anti tubercular agents. Among such antibiotics, sulfonamides are structural analogues of para-aminobenzoic acid (PABA), inhibiting the dihydro-pteroate synthase (DHPS) and blocking the folic acid synthesis pathway of MTB 6. Therefore, the aim of our paper publication is related to build a Quantitative structure - activity relationship (QSAR) models using multiple regression methods for synthesized sulfadiazine- salicylaldehyde Schiff base to explore the substitutional requirement which is essential for the improved anti-mycobacterial activity.
QSAR research has been considered a major tool in drug discovery to explore the ligand- receptor/ enzyme interactions, especially when the structural details of the target are not known. QSAR can aid in identifying functional groups with their various structural parameters used to increase the bioactivity leads designing new structures with the variation in enhanced bioactivity depends on changes in chemical structure 7. QSAR is an effective way for optimizing or correlating the biological activity within a congeneric sequence with certain structural parameters or with functional determinants, such as lipophilicity, polarizability, or electronic and steric properties 2, 8. Once a correlation between structure and activity/ property is found, any number of compounds including those not synthesized, yet can readily be screened for the selection of structures with desired properties 9.
Thus, in this article, we have selected a series of Schiff bases derived from the sulfa drug sulfadiazine with various salicylaldehydes 10. The author of this journal paper done work on Schiff base derived sulfa drug and they were characterized and evaluated against gram-positive and gram-negative bacteria, yeasts, moulds, Mycobacterium tuberculosis, nontuberculous mycobacteria and their cytotoxicity. Some derivatives showed more promising possession for these activities than original compound such as sulfadiazine; hence provide the hopeful hits for further antimicrobial drug development 10.
So, in this work we attempt to establish a quantitative structure-activity relationship between antimycobacterial activities of a series of 19 bioactive molecules derived from sulfadiazine- salicylaldehyde Schiff base was performed by using multiple linear regression analysis. There is no QSAR studies have been conducted out on these analogues so far. Moreover, the principal components analysis and contingency analysis has served to classify the compounds according to their activities and the variability of the descriptors. It allows reducing the number of variables not significant and making the information less redundant. The multiple linear regressions have served to select the descriptors used as the input parameters for the multiples non-linear regression and cross-validation to validate models used with the process leave-one-out (LOO) 11. This study has urged new efforts to find new anti-tubercular drug candidates with novel modes of action, develop pipelines for drug discovery and development and, in particular, try to find new regimens that can considerably shorten the duration of effective therapy which would improve patient compliance and slow down the emergence of drug resistant strains.
MATERIALS AND METHODS: MIC is defined as the lowest concentration of compound that inhibited visible growth. MIC values were manually converted into pMIC (predicted MIC; pMIC = -log MIC). For the development of QSAR models have been considered by using software MOE 2009.10. Biological activity was selected as the dependent variable and remaining molecular descriptors variables were selected as independent variables. The dataset molecules are validated by both internal and external validation procedures.
Dataset: The chemical structures of data set of 19 sulfadiazine- salicylaldehyde Schiff base analogue that already been synthesized and evaluated for their mycobacterium inhibitory activity in-vitro are represented in Fig. 1 10. The micromolar concentration of MIC (µM) was converted to logarithmic values for performing the QSAR study. The Biological data used in this study were a minimum inhibitory concentration of sulfadiazine- salicylaldehyde Schiff base has an activity against MTB Table 1.
FIG. 1: STRUCTURE OF SULPHADIAZINE SALICYLALDEHYDE SCHIFF BASE
TABLE 1: ANTI-MYCOBACTERIAL ACTIVITY OF SULPHADIAZINE-SALICYLALDEHYDE SCHIFF BASE DERIVATIVES
S. no. | Code | R | Experimental
MIC (µM) |
pMIC
(-log MIC) |
1 | 1a | H | 16 | 4.7959 |
2 | 1b | 5-F | 125 | 3.9031 |
3 | 1c | 5-Cl | 125 | 3.9031 |
4 | 1d | 5-Br | 16 | 4.7959 |
5 | 1e | 5-I | 62.5 | 4.2041 |
6 | 1f | 5-NO2 | 250 | 3.6021 |
7 | 1g | 5-CH3 | 125 | 3.9031 |
8 | 1h | 5-CH3O | 32 | 4.4949 |
9 | 1i | 5-OH | 125 | 3.9031 |
10 | 1j | 5-tert-Bu | 32 | 4.4949 |
11 | 1k | 6-Cl | 125 | 3.9031 |
12 | 1l | 3-Cl | 32 | 4.4949 |
13 | 1m | 3,5-Cl2 | 125 | 3.9031 |
14 | 1n | 3-Br-5-Cl | 125 | 3.9031 |
15 | 1o | 3-I-5-Cl | 125 | 3.9031 |
16 | 1p | 3,5-I2 | 250 | 3.6021 |
17 | 1q | - | 125 | 3.9031 |
18 | SDZ | - | 32 | 4.4949 |
19 | INH | - | 0.5 | 6.3010 |
SDZ: Sulphadiazine; INH: Isoniazide
QSAR Study: QSAR is an application of combinatorial chemistry to analyze experimental data and build numerical models of the data for prediction and interpretation. QSAR analysis was performed for 19 sulfadiazine- salicylaldehyde Schiff base compounds and their biological activity corresponding to MTB was collected from the reported literature 10 and analysis was done by using the QSAR module by using MOE 2009.10 software 12. The QSAR model was constructed for the above compounds with selected molecular descriptors viz. partition coefficients and molar refractivity, dipole moment, electrostatic interactions, lowest unoccupied molecular orbital (LUMO), highest unoccupied molecular orbital (HOMO) and the topological polar surface area from the QSAR descriptor panel.
The experimentally measured pMIC (µM) for 19 compounds of sulfadiazine- salicylaldehyde Schiff base against MTB were entered manually as activity fields. The regression analysis was performed; the RMSE and R2 values were derived from the QSAR fit. $PRED descriptor is the dependent variable for the activity field pMIC generated by QSAR validate. A QSAR correlation plot was generated by plotting the values of pMIC on the x-axis and the predicted values ($PRED) on the y-axis for all the above compounds. Furthermore, in order to identify active compounds in the data set “QSAR Contingency” application was implemented and followed by principal component analysis (PCA) by plotting a 3D graphical scatter plot with three vector values 13.
RESULTS AND DISCUSSION: In this study, a data set of 19 compounds was subjected to linear free energy regression analysis for model generation. Preliminary analysis was carried out in terms of correlation analysis. The correlations coefficients of different molecular descriptors with anti-mycobacterial activity are presented in Table 2. Some statistically significant QSAR models were chosen for discussion. The derived models in QSAR from multiple linear regression with forwarding stepwise shows a good correlation between biological activity and parameters.
The descriptors show the positive correlation among the parameters selected for the derived QSAR model. The positive coefficients suggest that inclusion of such carbon atoms in the molecules lead to increased biological activity. The computed molecular descriptors reported and predicted biological activities (pMIC) of the ligand dataset were shown in Table 2.
In order to validate the anti-mycobacterial activity, the QSAR model was constructed for these 19 compounds with experimentally derived pMIC (µM) values, system predicted pMIC, $PRED and all other standard molecular descriptors and reported in Table 1 and 2. Furthermore, the regression graph was plotted for pMIC vs. $PRED. Fig. 2. A cross-validated correlation coefficient (r2) value is 0.5998 and root mean square error value is 0.3829. The Z-Score was predicted for all the compounds, wherein a Z-score value of 2.5 and above can be considered to indicate the molecules which are outliers to the fit showed in Table 2. All the compounds showed significant Z-scores except compound 10 had a moderate value of 2.0571 Table 2. The Principal component analysis was carried out using three eigenvalues, viz. PCA1, PCA2 and PCA3, and a 3D graphical plot was created, which included 98% of the variance and presented in Fig. 3. In the plot, all the values were found to lie in the range of 3 to +3 and are shown with different colored spots which correspond to the pMIC of compounds Fig. 3. The result of contingency analysis showed the selected data set for QSAR study was good and all are within the standard limit as suggested. The useful value of C (contingency coefficient) is above 0. 6, V (cramer's V) is above 0.2, U (entropic uncertainty) is above 0.2 and R (linear Correlation) is above 0.2 Table 3.
TABLE 2: MOLECULAR DESCRIPTORS, EXPERIMENTAL AND PREDICTED pMIC (µM) VALUES OF 19 SULPHADIAZINE- SALICYLALDEHYDE SCHIFF BASE ON QSAR STUDY
S. no. | log P (o/w) | DM | EELE | HOMO | LUMO | MR | TPSA | pMIC | $PRED | $RES | $Z-SCORE |
1 | 1.9680 | 5.5421 | -673190.18 | -9.3618 | -1.1313 | 9.5620 | 104.5400 | 4.7959 | 4.3288 | 0.4671 | 1.2199 |
2 | 2.1580 | 5.8215 | -729813.56 | -9.2605 | -1.3275 | 9.6268 | 104.5400 | 3.9031 | 4.1519 | -0.2488 | 0.6496 |
3 | 2.5970 | 4.3427 | -720246.06 | -9.3652 | -1.1689 | 10.0508 | 104.5400 | 3.9031 | 4.1817 | -0.2786 | 0.7277 |
4 | 2.8030 | 3.9223 | -717033.93 | -9.4064 | -1.2076 | 10.3077 | 104.5400 | 4.7959 | 4.1918 | 0.6041 | 1.5777 |
5 | 3.1950 | 4.2123 | -717859.43 | -9.4686 | -1.1325 | 10.8967 | 104.5400 | 4.2041 | 4.1892 | 0.0149 | 0.0389 |
6 | 1.9400 | 2.1130 | -816182.75 | -9.9899 | -1.3927 | 10.0555 | 150.3600 | 3.6021 | 3.8820 | -0.2799 | 0.7310 |
7 | 2.3030 | 8.2412 | -725464.18 | -9.1281 | -0.9054 | 10.0179 | 104.5400 | 3.9031 | 4.1654 | -0.2623 | 0.6851 |
8 | 1.9610 | 5.0467 | -774665.87 | -8.9640 | -1.1874 | 10.2083 | 113.7700 | 4.4949 | 4.0117 | 0.4832 | 1.2618 |
9 | 1.6970 | 6.7929 | -729117.18 | -8.9668 | -1.3101 | 9.6923 | 124.7700 | 3.9031 | 4.1540 | -0.2509 | 0.6553 |
10 | 3.5040 | 8.1242 | -872140.50 | -9.1249 | -1.1238 | 11.3742 | 104.5400 | 4.4949 | 3.7072 | 0.7877 | 2.0571 |
11 | 2.5580 | 9.5335 | -731670.25 | -9.4407 | -1.0136 | 10.0617 | 104.5400 | 3.9031 | 4.1461 | -0.2430 | 0.6345 |
12 | 2.5580 | 6.5893 | -720365.75 | -9.4431 | -1.0608 | 10.0617 | 104.5400 | 4.4949 | 4.1814 | 0.3135 | 0.8188 |
13 | 3.2240 | 6.5733 | -773697.62 | -9.3811 | -1.4083 | 10.5247 | 104.5400 | 3.9031 | 4.0147 | -0.1116 | 0.2916 |
14 | 3.4300 | 5.3825 | -769072.87 | -9.4816 | -1.1832 | 10.8056 | 104.5400 | 3.9031 | 4.0292 | -0.1261 | 0.3293 |
15 | 3.8220 | 6.2276 | -770192.50 | -9.4114 | -1.4355 | 11.3947 | 104.5400 | 3.9031 | 4.0257 | -0.1226 | 0.3202 |
16 | 4.4200 | 6.2836 | -769047.93 | -9.5111 | -1.3875 | 12.2023 | 104.5400 | 3.6021 | 4.0293 | -0.4272 | 1.1156 |
17 | 1.1560 | 6.6635 | -771228.06 | -10.0097 | -1.2203 | 9.2580 | 139.3600 | 3.9031 | 4.0225 | -0.1194 | 0.3117 |
18 | -0.0260 | 7.1634 | -415397.40 | -9.1805 | -0.5092 | 6.4432 | 97.9700 | 4.4949 | 5.1342 | -0.6393 | 1.6695 |
19 | -0.7970 | 1.5982 | -182473.57 | -10.3111 | -0.5516 | 3.7316 | 68.0100 | 6.3010 | 5.8619 | 0.4391 | 1.1467 |
logP (o/w): Log octanal/water partition coefficient, DM: Dipole moment, EELE: electrostatic interactions, LUMO: lowest unoccupied molecular orbital, HOMO: highest unoccupied molecular orbital, MR: Molecular refractivity, TPSA: Topological polar surface area, pMIC: experimentally calculated log MIC, $PRED, $RES and $Z-SCORE: predicted pMIC , residual and Z-score values.
FIG. 2: FITNESS GRAPH BETWEEN OBSERVED AND PREDICTED BINDING AFFINITY FOR DATA SET MOLECULES. Linear correlation graph of experimentally measured pMIC values of data set of 19 compounds against MTB. The linearity of QSAR model was shown with error value (RMSE) and correlation actor (R2)
FIG. 3: PRINCIPAL COMPONENT ANALYSIS. Principal Component Analysis (PCA) plot of all these compounds, 3D graphical plot was constructed with three eigenvectors PCA1, PCA2 and PCA3 in the range of 3 to +3 and coordinate value of each compound is represented with a colored spot.
TABLE 3: QSAR CONTINGENCY ANALYSIS (DATABASE FIELDS AND SUGGESTED DESCRIPTORS FOR QSAR)
S. no. | C | V | U | R | Field |
1 | 0.83271 | 0.43412 | 0.48112 | 0.17498 | AM1_dipole |
2 | 0.79487 | 0.37816 | 0.42087 | 0.59978 | AM1_Eele |
3 | 0.83786 | 0.44308 | 0.46851 | 0.10267 | AM1_HOMO |
4 | 0.83530 | 0.43859 | 0.48789 | 0.44269 | AM1_LUMO |
5 | 0.84000 | 0.44691 | 0.41928 | 0.36180 | logP(o/w) |
6 | 0.82910 | 0.42808 | 0.43736 | 0.54755 | mr |
7 | 0.81562 | 0.40694 | 0.43214 | 0.43110 | TPSA |
8 | 0.79487 | 0.37816 | 0.42087 | 0.59978 | $PRED |
9 | 0.87678 | 0.52631 | 0.68071 | 0.40022 | $RES |
10 | 0.86207 | 0.49105 | 0.60277 | 0.21975 | $Z-SCORE |
11 | 0.84454 | 0.45527 | 0.49842 | 0.25462 | $XPRED |
12 | 0.86046 | 0.48750 | 0.63876 | 0.65453 | $XRES |
13 | 0.86046 | 0.48750 | 0.63876 | 0.68296 | $XZ-SCORE |
14 | 0.87678 | 0.52631 | 0.68071 | 0.40022 | Residues |
CONCLUSION: A cross-validated correlation coefficient (r2) value is 0.5998 and root mean square error value is 0.3829. Moreover, all the compounds showed significant fitness score ($Z-scores) and lie within the limit. The result of principal component analysis and contingency analysis showed the selected data set for QSAR study was good and all are within the standard limit as suggested in the literature.
The QSAR studies on these compounds have shown that the predicted pMIC values of the compounds have an acceptable correlation with the experimental values from the generated regression, principal component analysis (PCA) plots and the contingency analysis for QSAR.
Quantitative structure activity relationship studies revealed that the anti-mycobacterial activities of these derivatives against the test microorganisms are mainly governed by log P value, dipole moment, HOMO, LUMO, and the molar refractivity, a steric parameter. Thus a proper substitution of the group on aromatic ring probably improves the effectiveness of these derivatives as anti-mycobacterial agents.
In conclusion, QSAR pharmacophore model generated from the present study should be useful in the designing a similar group of more potent substituted compounds of Schiff base of sulpha drugs analogues as anti-mycobacterial agents.
ACKNOWLEDGEMENT: The authors are highly thankful to the Principal and Management of Nandha College of Pharmacy, Erode, Tamil Nadu for providing the necessary facilities to carry out this research.
CONFLICT OF INTEREST: The authors declare no conflicts of interest.
REFERENCES:
- Piyush K, Pinki R and Pooja S: Medicinal Chemistry and Development of Antitubercular Drugs. Int J Pharm Sci and Res. 2013; 4: 2080-2093.
- Vikramjeet J, Balasubramanian N, Munish A, Dharmarajan S, Perumal Y, Erik DC, Christophe P and Jan B: Isonicotinic acid hydrazide derivatives: synthesis, antimicrobial activity, and QSAR studies. Med Chem Res. 2012; 21: 1451-1470.
- Elham P, Hadi K and Nahid SR: A study on antitubercular and antimicrobial activity of isoniazid derivative. Zahedan J Res Med Sci. 2015; 17: 7-10.
- Jerad SA, Geetha LS, Shiva KV, Madhuraj M, Durga V and Surya PR: Design, synthesis, characterization, and biological evaluation of some novel antitubercular agents targeting thymidylate synthetase x, enoyl-ACP reductase, and l, d-transpeptidase. J Pharm Res. 2018; 12: 3.
- Krátký M, Vinšová J, Novotná E, Mandíková J, Wsól V, Trejtnar F, Ulmann V, Stolaříková J, Fernandes S, Bhat S and Liu JO: Salicylanilide derivatives block Mycobacterium tuberculosis through inhibition of isocitrate lyase and methionine aminopeptidase. Tuberculosis (Edinb) 2012; 92: 434-439.
- Sirwan MA and Michel D: In-vitro susceptibility of Mycobacterium tuberculosis to Trimethoprim and Sulfonamides in France. Antimicrob Agents Chemother 2013; 57: 6370-6371.
- Corwin H and Toshio F: A Method for the Correlation of Biological Activity and Chemical Structure. J Am Chem Soc. 1964; 86: 1616–1626.
- Sivaprakasam P, Xiea A and Doerksen RJ: Probing the physicochemical and structural requirements for glycogen synthase kinase-3a inhibition: 2D-QSAR for 3-anilino-4-phenylmaleimides. Bio Med Chem. 2006; 14: 8210-8218.
- Sharma MC, Kohli DV, Chaturvedi SC and Smita S: Molecular modelling studies of some substitued 2-butylbenzimidazoles angiotensin II receptor antagonists as antihypertensive agents. Digest J Nanomat Biostruc. 2009; 4: 843-856.
- Martin K, Magdaléna D, Jiˇrí J, Klára K, František T, Jiˇrina S and Jarmila V: Sulfadiazine Salicylaldehyde-Based Schiff Bases: Synthesis, Antimicrobial Activity and Cytotoxicity. Molecules 2017; 22: 1573.
- Boukarai Y, Khalil F and Bouachrine M: QSAR study of flavonoid derivatives as in-vitro inhibitors agents of Aldose Reductase (ALR2) enzyme for diabetic complications. J Materials Envir Sci. 2017; 8: 1533-1545.
- Molecular Operating Environment (MOE), 2009.10; Chemical Computing Group ULC, 1010 Sherbooke St. West, Suite #910, Montreal, QC, Canada, H3A 2R7, 2018.
- Tirumalasetty MB, Sivarathri SR, Baki VB, Savita D, Aluru R, Thirunavakkarasu S and Wudayagiri R: Molecular docking, molecular dynamics simulation, biological evaluation and 2D QSAR analysis of flavonoids from Syzygium alternifolium as potent anti-Helicobacter pylori agents. The Royal Soc Chem. 2017; 7: 18277-18292.
How to cite this article:
Thangavelu P, Ramakrishnan M, Mani K, Somasundram K, Srinivas D and Thangavel S: Quantitative structure activity relationship studies of sulfadiazine- salicylaldehyde schiff base as an anti-mycobacterial agents. Int J Pharm Sci & Res 2018; 9(12): 5204-09. doi: 10.13040/IJPSR.0975-8232.9(12).5204-09.
All © 2013 are reserved by International Journal of Pharmaceutical Sciences and Research. This Journal licensed under a Creative Commons Attribution-NonCommercial-ShareAlike 3.0 Unported License.
Article Information
19
5204-5209
478
928
English
IJPSR
P. Thangavelu *, M. Ramakrishnan, K. Mani, K. Somasundram, D. Srinivas and S. Thangavel
Nandha College of Pharmacy, Koorapalayam Pirivu, Pitchandam Palayam Post, Erode, Tamil Nadu, India.
drtpappa@yahoo.com
13 April, 2018
13 June, 2018
26 June, 2018
10.13040/IJPSR.0975-8232.9(12).5204-09
01 December, 2018