Designing of Novel 6(H)-1,3,4-Thiadiazine Derivatives as MMP12 Inhibitors: A MLR and Docking Approac...

Ajeet  , Laxmi Tripathi, Praveen Kumar

  Open Access OPEN ACCESS  Peer Reviewed PEER-REVIEWED

Designing of Novel 6(H)-1,3,4-Thiadiazine Derivatives as MMP12 Inhibitors: A MLR and Docking Approach

Ajeet  1, Laxmi Tripathi1, Praveen Kumar1,

1Department of Pharmaceutical Chemistry, S. D. College of Pharmacy and Vocational Studies, Muzaffarnagar, India


Here 6(H)-1,3,4-thiadiazine analogues have been used to correlate the inhibiting constant with the eccentric connectivity index (ECI), fragment complexity (FC), McGowan volume (MG) and topological polar surface area (TPSA) for studying the quantitative structure activity relationship (QSAR). Correlation may be an adequate predictive model which can help to provide guidance in designing and subsequently yielding greatly specific compounds that may have reduced side effects and improved pharmacological activities. The literature survey revealed that there are so many models available for MMP-1, MMP-2, MMP-3, MMP-7, MMP-8, MMP-9 and MMP-13, but QSAR model in respect of MMP-12 inhibition was necessity. We have used multiple linear regression (MLR) for developing QSAR model. For the validation of the developed QSAR model, statistical analysis such as cross validation test, standard deviation, quality factor, fischers test, root mean square deviation (RMSD), variance; and internal validation such as Y-randomization test have been performed and all the tests validated this QSAR model with fraction of variance r2 = 0.9364 and LOO-CV variance q2= 0.9146. Thirteen novel 6(H)-1,3,4-thiadiazine analogues have been designed and their inhibiting constant have been calculated with the developed QSAR model. It was found that the calculated inhibiting constant of these analogues were within the same range as of the training set. Further, these 6(H)-1,3,4-thiadiazine analogues have been docked with the catalytic domain of human matrix metalloproteinase (MMP12) which shows better docking score as compared to the (N-hydroxy-2-(N-hydroxyethyl)biphenyl-4-ylsulfonamido) acetamide, a MMP12 inhibitor. The results suggested that the designed novel 6(H)-1,3,4-thiadiazine analogues could be developed as a good MMP12 inhibitors.

At a glance: Figures

Cite this article:

  •  , Ajeet, Laxmi Tripathi, and Praveen Kumar. "Designing of Novel 6(H)-1,3,4-Thiadiazine Derivatives as MMP12 Inhibitors: A MLR and Docking Approach." American Journal of Pharmacological Sciences 1.2 (2013): 29-34.
  •  , A. , Tripathi, L. , & Kumar, P. (2013). Designing of Novel 6(H)-1,3,4-Thiadiazine Derivatives as MMP12 Inhibitors: A MLR and Docking Approach. American Journal of Pharmacological Sciences, 1(2), 29-34.
  •  , Ajeet, Laxmi Tripathi, and Praveen Kumar. "Designing of Novel 6(H)-1,3,4-Thiadiazine Derivatives as MMP12 Inhibitors: A MLR and Docking Approach." American Journal of Pharmacological Sciences 1, no. 2 (2013): 29-34.

Import into BibTeX Import into EndNote Import into RefMan Import into RefWorks

1. Introduction

Matrix Metalloproteinase (MMPs) is a family of calcium dependent zinc containing endopeptidase, which is responsible for the remodeling of tissue and extracellular matrix degradation, including gelatin, matrix glycoprotein, proteoglycan etc. These are excreted by a variety of connective tissue and some pro-inflammatory cells like lymphocytes, osteoblasts, macrophages, endothelial cells, neutrophils, and fibroblasts. MMPs have been considered as a target for cancer therapy and a number of natural and synthetic MMPs inhibitors have been identified as cytostatic and anti-angiogenic agents. MMPs consists of four domains, these are hinge region, N-terminal pro-domain, catalytic domain, and C-terminal hemopexin-like domain. It may be responsible for interaction with tissue inhibitors of MMPs as well as for the macromolecular substrate recognition [1]. Matrix metalloproteinase-12 (MMP-12) also known as macrophage metalloelastase (MME) or macrophage elastase (ME) is an enzyme that in humans is encoded by the MMP12 gene. The literature review revealed that inhibitory activity of various compounds with the MMP-1, MMP-2, MMP-3, MMP-7, MMP-8, MMP-9 and MMP-13 has been correlated with hydrophobicity, Kier’s first order valence molecular connectivity index, electrotopological state indices of the atom and polarizability. Hydrophobicity is the most important descriptor and have been found to describe the inhibitory activity of 6(H)-1,3,4-thiadiazine derivatives [2]. Further, it was concluded that development of a QSAR model dealing with MMP-12 inhibition was necessary.

The QSAR studies are perfect tool for understanding the drug design process in terms of their chemical-pharmacological activity interaction. QSAR studies can focus on mechanism of action of ligands with human, bacteria, virus, membranes, enzymes etc. The QSAR methodology comprises of computationally derived descriptors to correlate with pharmacological activities. These descriptors are principally of four types such as electronic, stearic, hydrophobic and topological indices [3]. The descriptors used for developing the QSAR model are eccentric connectivity index (ECI) [4], fragment complexity (FC) [5], McGowan volume (MG) [6] and topological polar surface area (TPSA) [7]. Rational drug design helps to facilitate and fasten the drug designing process, which involves various methods to identify novel compound, out of them one method is the docking of molecules with the receptor [8].

In the present study, we developed a QSAR model on a series of 6(H)-1,3,4-thiadiazine analogues with respect to their MMP12 inhibition. Further, these 6(H)-1,3,4-thiadiazine analogues have been docked with the catalytic domain of MMP12.

Figure 1.Structure of 6(H)-1,3,4-thiadiazine used for model development

2. Materials and Method

All the bioactivity values and information about 2D structure of 6(H)-1,3,4-thiadiazine derivatives (Figure 1) were taken from literature [1]. IC50 is referred as the molar concentration of a compound that inhibits 50% growth of bacteria [3]; log1/Ki (Ki is inhibition constant) is subsequent variable that comprises the bioactivity parameter for the QSAR model. In order to calculate the 2D molecular descriptors, PaDEL descriptor software [9] which incorporate CDK library for descriptor calculation have been used after optimitizing the 6(H)-1,3,4-thiadiazine derivatives. For the development of QSAR model, MLR [3] has been employed and all were validated through statistics. The novel 6(H)-1,3,4-thiadiazine molecules have been designed and optimized through ChemDraw Ultra 7.0.and their bioactivity values have been calculated from developed QSAR model by putting the descriptor values in the QSAR equation. The docking studies [8] of novel 6(H)-1,3,4-thiadiazine have been performed with the Molegro Virtual Docker (Free trial version).

2.1. Statistical Parameters
2.1.1. Fraction of Variance (R2)

The value of fraction of variance may vary between 0 (means model without explanatory power) and 1 (means perfect model). QSAR model having r2 > 0.6 will only be considered for validation [3].

2.1.2. Cross-validation Test (q2)

A QSAR model must have q2> 0.5 for the predictive ability [3].

2.1.3. Standard Deviation (s)

The smaller s value is always required for the predictive QSAR model.

2.1.4. r2-q2 < 0.3

The difference between r2 and q2 should never be exceeding by 0.3. A large difference suggests the following: presence of outliers, over-fitted model, and presence of irrelevant variables in data [3].

2.1.5. Quality Factor (Q)

Overfitting and chance correlation, due to excess number of descriptors, can be detected by Q value. Positive value for this QSAR model suggests its high predictive power and lack of overfitting.

2.1.6. Fischer Statistics (F)

The F value of QSAR model was compared with their literature value at 95% level.

2.2. Model Validation

The QSAR model validation was carried with statistical analysis and with internal validation.

2.3. Docking Studies

The designed novel 6(H)-1,3,4-thiadiazine derivatives were docked into catalytic domain of humanMMP12. For this study, X-ray crystal structure of human MMP12 was taken from protein data bank with PDB id 3N2V.

2.3.1. Determination of Binding Site

The catalytic sites of human MMP12 have been found by DETECT CAVITIES under the PREPARATION option of Molegro Virtual Docker (Free trial version) which detects 5 cavities by default with a 10 Å resolution.

3. Results and Discussion

From the data (Table 1), QSAR equation 1 have been developed and 95% confidence intervals are given in parantheses as follows

3.1. Validation of QSAR Model

A quantitative assessment of model robustness has been performed through model validation. All the statistical results of model validation have been given (Table 2).

Table 1. Descriptors used to derive QSAR equation along with bioactivities for the inhibition by 6(H)-1,3,4-thiadiazine analogues with substituents

3.2. Internal Validation

Y-Randomization Test:, This technique is used to establish the QSAR model robustness. For this test, the dependent variable vector is randomly shuffled, and a new QSAR model is developed using the unchanged independent variable. This process was repeated for five times. The statistical data of r2 for five runs are given (Table 3). The values r2<0.6in Y-randomization test confirm the robustness of this QSAR model [3].

Table 3 Results of internal validation: Y-randomization test (5 runs)

According to the developed QSAR model, the 6(H)-1,3,4-thiadiazine derivatives must have less ECI for enhanced MMP12 inhibitory activity at X, X1 andY substituents. A negative coefficient of FC containing X, X1 and Ysubstituents favors the MMP12 inhibitory activity of 6(H)-1,3,4-thiadiazine derivatives. As far as effects of the TPSA on the MMP12 inhibitory activity of 6(H)-1,3,4-thiadiazinederivatives concerned, the developed QSAR model suggested that an decrement in TPSA at substituents X, X1 and Y will definitely be favorable to the activity; along with positive McGowan volume, as discussed by Verma and Hansch, 2010 [3]. A comparison (MLR plots) of observed values and predicted values of log1/Kifor 6(H)-1,3,4-thiadiazine derivatives used for development of QSAR equation is shown in Figure 2 and multiple linear graph is shown in Figure 3.

Figure 2. Multiple linear regression plot for QSAR study
Figure 3. Multiple linear graph between Number of data points and bioactivities
3.3. Designing and Optimization of Novel 6(H)-1,3,4-thiadiazine Derivatives

The novel 6(H)-1,3,4-thiadiazine derivatives were designed and their energy minimization for highest stability was performed. Details of substituents, calculated descriptor values and predicted bioactivity through derived QSAR model have been given (Table 4).

3.4. Docking Results
3.4.1. Binding Site Analysis

The experimental analysis of binding site shows that Leu 214, Tyr 240, Glu 219, Leu 181, Ala 182, Ile 180, His 228, Pro 238, Phe 237, Val 235, His 218 could be the catalytic site residue present in the structure of human MMP12.

Table 4. Results of QSAR studies on novel 6(H)-1,3,4-thiadiazine derivatives

3.4.2. Docking Studies of Novel 6(H)-1,3,4-thiadiazine Derivatives with Human MMP12

Docking studies showed that derivatives THZ2-14 were docked in the same binding site as docked by (N-hydroxy-2-(N-hydroxyethyl)biphenyl-4-ylsulfonamido) acetamide; a pre-existing human MMP12 inhibitor. The protein-ligand interaction energy of novel inhibitors were given by Pose organizer option in Molegro virtual docker (Free trial version) for 5 different poses of novel inhibitors. The best pose interaction energy was found to be as -179.676, -180.116, -180.162, -189.003, -182.567, -186.026, -175.835, -175.277, -174.593, -176.54, -117.702, -186.236 and -190.092. Here, negative values for interaction energy would reflect the positive docking approach. Number of hydrogen bonds and other binding details with cavities (Table 5) and docking images (Table 6) are given.

Table 5. Docking results of novel 6(H)-1,3,4-thiadiazine derivatives

Table 6. Docked photographs of novel 6(H)-1,3,4-thiadiazine derivatives

3.5. Comparison of Docking Results with Pre-existing Human MMP-12 Inhibitor (N-Hydroxy-2-(N-Hydroxyethyl)Biphenyl-4-Ylsulfonamido)Acetamide

Docking analysis of (N-hydroxy-2-(N-hydroxyethyl) biphenyl-4-ylsulfonamido)acetamide with the MMP-12 revealed that interacting residues (amino acids) were His 222, His 228, Pro 238, Phe 237, Val 235, His 218, Thr 239, His 183, Glu 219, Ala 182, Thr 215, Tyr 240 and Leu 214 (Figure 4). These results were similar to the docking analysis of novel designed 6(H)-1,3,4-thiadiazine derivatives.

Figure 4. Docking image of (N-hydroxy-2-(N-hydroxyethyl)biphenyl-4-ylsulfonamido) acetamide with the MMP-12 showing interacting amino acids

On docking analysis, the docking images of all novel 6(H)-1,3,4-thiadiazine derivatives were found to be superimposed with the (N-hydroxy-2-(N-hydroxyethyl)biphenyl-4-ylsulfonamido) acetamide, a pre-existing MMP-12 inhibitor (Figure 5).

Figure 5. Superimposed docking images of novel 6(H)-1,3,4-thiadiazine derivatives (showing with wireframe model) with the pre-existing ligand (N-hydroxy-2-(N-hydroxyethyl)biphenyl-4-ylsulfonamido) acetamide with the MMP-12 showing interacting amino acids (showing with ball and stick model)

4. Conclusion

Computational study of novel 6(H)-1,3,4-thiadiazine derivatives proved them potential human MMP12 inhibitors. Although a systemic biochemical study is necessary to confirm the findings. On comparing the chemical structure of 6(H)-1,3,4-thiadiazine derivatives with (N-hydroxy-2-(N-hydroxyethyl) biphenyl-4-ylsulfonamido) acetamide, a pre-existing human MMP12 inhibitor; it is concluded that a bicyclic system, a sulfonamido linkage and a acetamide functionality are essential pharmacophoric requirements in designing of human MMP12 inhibitors.


[1]  Verma, R.P.,Hansch, C., “Matrix metalloproteinases (MMPs): Chemical–biological functions and (Q)SARs,” Bioorg. Med. Chem, 15, 2223-2268, 2007.
In article      CrossRefPubMed
[2]  Belaaouaj A, Shipley J.M., Kobayashi D.K., Zimonjic D.B., Popescu N., Silverman G.A., Shapiro S.D., “Human macrophage metalloelastase. Genomic organization, chromosomal location, gene linkage, and tissue-specific expression,” J. Biol. Chem, 270 (24), 14568–14575, 1995.
In article      PubMed
[3]  Verma, R.P.,Hansch, C., “QSAR modeling of taxane analogues against colon cancer,” Eur. J. Med. Chem. 45, 1470-1477, 2010.
In article      CrossRefPubMed
[4]  Sardana, S.,Madan, A.K., “Topological models for prediction of antihypertensive activity of substituted benzylimidazoles,” J. Mol.Struct, 638, 41-49, 2003.
In article      CrossRef
[5]  Gregg, S.,Eiso, A.B., Jan, S., “Integration of fragment screening and library design,” Drug Discovery Today, 12(23/24), 1032-1039, 2007.
In article      
[6]  Abraham, H. M., Ibrahim, A.,Zissimos, M. A., ”Determination of sets of solute descriptors from chromatographic measurements.,” Jour. Chromat, 1037, 29-47, 2004.
In article      CrossRef
[7]  Hou, T.J.,Xu X.J., “ADME Evaluation in Drug Discovery. 3. Modeling Blood-Brain Barrier Partitioning Using Simple Molecular Descriptors,” J. Chem. Inf. Comput. Sci, 43, 2137-2152, 2003.
In article      CrossRefPubMed
[8]  Sharma, J.,Ramanathan, K.,Rao, S., “Identification of Potential Inhibitors against Acetylcholinesterase Associated With Alzheimer's Diseases: A Molecular Docking Approach,” J. Comput. Method Mol. Design, 1(1), 44-51, 2011.
In article      
[9]  Yap, C.W., 2011. “PaDEL-descriptor: An open source software to calculate molecular descriptors and fingerprints,” J. Comput. Chem. 32(7), 1466-1474, 2011.
In article      
  • CiteULikeCiteULike
  • Digg ThisDigg
  • MendeleyMendeley
  • RedditReddit
  • Google+Google+
  • StumbleUponStumbleUpon
  • Add to DeliciousDelicious
  • FacebookFacebook
  • TwitterTwitter
  • LinkedInLinkedIn