key: cord-0841615-ohxojj4s
authors: Ahmed, Atteeque; Saeed, Aamer; Ejaz, Syeda Abida; Aziz, Mubashir; Hashmi, Muhammad Zaffar; Channar, Pervaiz Ali; Abbas, Qamar; Raza, Hussain; Shafiq, Zahid; El-Seedi, Hesham R.
title: Novel adamantyl clubbed iminothiazolidinones as promising elastase inhibitors: design, synthesis, molecular docking, ADMET and DFT studies
date: 2022-04-19
journal: RSC advances
DOI: 10.1039/d1ra09318e
sha: 0703d76f1b0433912007c79a88a0ac20834fc1ae
doc_id: 841615
cord_uid: ohxojj4s

Porcine Pancreatic Elastase (PPE) is a serine protease that is homologous to trypsin and chymotrypsin that are involved in various pathologies like inflammatory disease, Chronic Obstructive Pulmonary Disease (COPD), acute respiratory distress syndrome, cystic fibrosis, and atherosclerosis. PPE if remained uninhibited would lead to digestion of important connective tissue. We developed new structurally diverse series of adamantyl-iminothiazolidinone hybrids to divulge elastase inhibition assay. To identify potent derivatives, in silico screening was conducted and in vitro studies disclosed that the compounds 5a, 5f, 5g, and 5h showed excellent binding energies and low IC(50) values. In silico studies including molecular docking, DFT studies (using the B3LYP/SVP basis set in the gas phase) drug likeness scores and molecular dynamic simulation studies were conducted to evaluate protein–ligand interactions and to determine the stability of top ranked conformation. In silico studies further supported the results of in vitro experiments and suggest these derivatives as novel inhibitors of elastase enzyme.

Elastase is a pancreatic enzyme that is resistant to degradation in the gastrointestinal tract. It is made by pancreatic acinar cells. Its measurement in faeces has been widely used to determine pancreatic exocrine adequacy because to its simplicity; nonetheless, it is insensitive to mild to moderate pancreatic insufficiency. 1,2 It is initially synthesized as an inactive zymogen before being activated by trypsin in the duodenum. Reduced faecal elastase concentrations have been seen in various gastrointestinal illnesses (e.g., duodenal enteropathy), in circumstances where elastase is destroyed by bacteria (e.g., inammatory bowel disease), and as a consequence of diarrhea's dilution effect. 3 Because the test looks for human elastase, it does not react with the pancreatic enzyme replacement drug derived from swine that is used to treat cystic brosis. Elastin is also present in the skin's bottom layer, alongside collagen, as well as nerves, blood arteries, lymph nodes, and muscles. 4 Collagen is a fundamental component of the extracellular matrix (ECM) and has a direct impact on the skin. The preservation of collagen structure is linked to the intrinsic and photoaging processes of the skin. 5, 6 Elastase are serine proteases that hydrolyze the amides and esters present in many proteins, including elastin. Elastic bers are composed of elastin, a protein found in the majority of connective tissues and tendons. Because elastase destroys elastic bers and leads to skin inammation and other skin-related issues, decreasing elastase activity may be an effective method of avoiding skin ageing. 7 Adamantyl moiety is a non-polar cage structure that provides lipophilic character to the molecule and improves the barrier crossing ability of a molecule to the biological system. In most cases, the adamantyl group tunes the hydrophobicity of the molecule by enhancing the penetrating ability and therapeutic activity of several trial compounds through various mechanisms. [8] [9] [10] The steric factor of the adamantyl group promotes the stability and pharmacologic half-life of compounds. From the last few decades, adamantyl appending compounds as potential pharmacophore, attracted the attention of several synthetic chemists and pharmacists. 11 As a result, adamantyl analogs passed from the intensive biological investigation. 12 To date, several marketed drugs incorporating with adamantyl scaffold are now available. Adamantyl-containing compounds exhibit a broad spectrum of biological activities such as antiviral, anti-HIV, antifungal, anti-diabetic, antibacterial, 13 anti-inammatory, 14 antimicrobial, hypoglycemic activity, treatment of COVID-19, 15, 16 and treatment for neurological disease. 17, 18 Acyl thioureas signify a class of pharmaceutically important compounds to develop the architect of a potential drug. These molecules have been explored for many pharmaceutical activities. [19] [20] [21] and in agriculture as an insecticide and weed killer. In addition, they are the precursor for the synthesis of a number of heterocyclic compounds and metal complexes 22 Acyl thiourea have also been applied as ion sensors and for the estimation of metal ions. [23] [24] [25] These electron-rich molecules act as ligand and organocatalysts in various asymmetric synthesis [26] [27] [28] Heterocyclic skeletons containing ve-member rings with two heteroatoms are promising core structures, commonly encountered in organic and medicinal chemistry. Thiazolidinone motifs are biologically signicant ve-membered heterocycle that is very common substructures in various pharmacological active molecules and known to exhibit a wide range of biological activities such as anticancer, antifungal, anti-inammatory, antimicrobial, antidiabetic, inhibiting neuraminidase of inuenza virus, and anti-HIV, 29 antischistosomal activity. Among these are inhibitors for necroptosis, selective GSK-3b inhibitor, anticancer activity, 30 aldose reductase inhibitor, 31 potent antiproliferative agent and inhibitor for non-membrane protein tyrosine phosphate. The aforesaid facts inspired us to synthesize novel ethyl 3-adamantanyl-2-((2-methyl benzoyl)imino)-4-oxothiazolidin-5-ylidene acetates (5a-j) for the appraisal of elastase inhibition assay. The structures of some pharmacologically active compounds and designed molecules are shown in Fig. 1 . 

The target compounds were synthesized according to synthetic route described in Scheme 1. Suitably substituted benzoic acids were converted to corresponding acid chlorides (1a-j) using dichloromethane as solvent. The acid chloride was converted to isothiocyanate by adding potassium thiocyanate in acetone. Reaction of the latter in situ with an equimolar amount of amantadine (2) afforded the acyl thioureas (3a-j). Ethyl 4ethoxypent-4-en-2-ynoate with (3a-j) at room temperature in dry ethanol furnished the desired compounds (5a-j) in excellent yields and high purity. The complete structures of synthesized derivatives are given in (Fig. 2) .

The important SAR features can be deduced from this study. The inhibitory potential of the synthesized amantadine derivatives was investigated using elastase from porcine pancreas. The relationships between different substituted and unsubstituted patterns were analyzed on their ability to inhibit elastase. Such type of analysis assists to infer the SAR among the synthesized screened derivatives. Oleanolic acid is a pentacyclic triterpenoid and has been used as stander. In the current study, amantadine thiazolidinone derivatives were synthesized to nd out more effective and selective antagonists of elastase. Among the synthesized analogues, compounds 5a, 5f, 5g, and 5h were showed high potency against elastase. The presence of unsubstituted phenyl group in 5g, provide effectiveness to develop high potency against elastase with IC 50 of 0.124 AE 0.022 mM (Table 1) which is 50 times more potent as compared to oleanolic acid value of 5.996 AE 0.882 mM. With respect to structural comparison compounds with a methyl 5a, ouro 5f, chloro 5h, and bromo 5d substituted phenyl group exhibit excellent inhibitory activity. This is due to positive mesomeric effect of these substituents. Interestingly the compound 5g doesn't have any substitution on benzoylimino group of a compound, but it showed excellent IC 50 value and good binding energy. It can be perceived that absence of mesomeric effect might increase activity against PPE enzyme.

An inhibition kinetic study was carried out in order to better understand the mode of action of the synthetic compound against elastase. Based on our IC 50 results, we choose our most potent compound 5g for further investigation to determine the type of inhibition and the constant of inhibition. 1/V versus substrate N-succinyl-Ala-Ala-p-nitroanilide 1/[S] in the presence of different inhibitor concentrations was plotted using a Lineweaver-Burk plot in which the results were a series of straight lines; however, the result of a Lineweaver-Burk plot of compound 5g revealed that V max remained constant without signicantly affecting the slopes. While K m increases with increasing concentration, V max remains constant with only a statistically insignicant difference. A competitive inhibition model was used to explain this behavior (Fig. 3A) , and a second plot (Fig. 3B ) of slope against concentration of 5g revealed the E I dissociation constant. The inhibitor concentration of 5g was plotted against the slope of the graph, and the calculated K i value was found to be 0.6 mM.

rstly, using the B3LYP/SVP basis set in the gas phase, and their values are shown in Table 2 .

The geometries of selected compounds were optimized to lowest energy gradient and no imaginary frequencies were observed which demonstrate that all structures were true local minima. The optimized structures of selected compounds are shown in Fig. 4 .

The molecular orbital (MO) analysis is of prime importance in quantum chemistry. It has been broadly dening chemical behavior. The highest occupied molecular orbital (HOMO) and lower unoccupied molecular orbital (LUMO) are the chief molecular orbitals in a compound. They are used to explain chemical properties such as reactivity, stability, and kinetics. The FMO orbitals of synthesized compounds are shown in Fig. 5 .

The hardness (h) dened the extent to which a molecule is hard or so. The soer the molecule, the better will be the reactivity.

Electronegativity (X) characterizes the power to attract electrons. All these properties were estimated at. The compound 5c showing smallest HOMO-LUMO energy gap value (0.064 eV) predicts a good chemical reactivity. The compound 5j is showing a small hardness value (0.071), high polarizability thus leading to be the so molecule among all compounds. The higher the electronegativity value in the case of compound 5d showing (0.173) value exposes that 5d has good electron attracting power and acts as a better electrophile among all compounds. While compounds 5f and 5a also showed good reactivity aer 5c, with the energy gap values 0.067 eV and 0.069 eV, respectively. In the same manner, just aer 5j the compounds 5b, 5g and 5d showing high polarizability with the values of 0.072, 0.076 and 0.078, respectively. Energetic parameters of compounds (5a-j) are showing in Table 3 below. Molecular docking studies. Porcine Pancreatic Elastase (PPE) is a serine protease that is homologous to trypsin and chymotrypsin. It is composed of single polypeptide chain of 240 amino acids that are linked with four disulphide linkages. Moreover, it is complexed with Ca 2+ atom which is necessary for its stability. 32 PPE consist of four subsites, among these most of inhibitors bind to S1 and S2 active sites. S1 site is consider as a primary specicity active site (Fig. 6 ). PPE plays pivotal role in onset and progression of inammatory disease like cancer and COPD. 33 PPE if remained uncontrolled would lead to digestion of important connective tissue. Keeping in view the importance of PPE, in silico studies were conducted to support the result of in vitro assays. Crystallographic structure of PPE was downloaded from Protein data bank. (PDB ID: 1BMA; resolution 1.99). The important amino acids residues of catalytic site were selected and docked with potent thiazolidinone derivatives. The important amino acid residues of PPE protein are as follows; THR221, SER203, HIS60, SER222, PHE223, VAL224, SER225, TRP179, CYS199, GLY201, GLY198, THR152, LEU227, CYS229, ALA104 THR152, GLY228, ARG226, VAL103, SER203, HIS60, and GLN200. Synthesized derivatives were selected on the basis of IC 50 values obtained during in vitro assays. Based on IC 50 value, 5a, 5f, 5g and 5h was selected for in silico studies as these derivatives showed excellent inhibitory potential during in vitro assays.

Synthesized derivatives were evaluated for binding affinities and binding scores with amino acid residues of active pocket. Binding affinities were determined on the basis of predicted inhibitory constant value (k i ). Most of the compounds exhibited potent binding scores and demonstrated good binding affinities. The best conformation of protein-ligand was selected and analyzed further for bonding and non-bonding interactions. Docking scores of potent compounds with their IC 50 values are given in Table 4 .

The detailed 3D and 2D binding interactions of derivative 5a within active pocket of PPE enzyme is shown in Fig. 7 . The amino acid residues involved in bonding and non-bonding interactions with 5a were THR152, LEU227, CYS229, GLY228, SER225, ARG226, TRP179, VAL103, PHE223, VAL224, SER203, HIS60, SER222 and GLN200. Briey substituted benzoyl imino ring of parent compound was making strong pi-sulfur linkage with C229. Compound 5a showed good binding energy and IC 50 value which might be due to substitution of methyl group on benzoylimino ring. It was observed that compound 5a was producing hydrogen bonding of moderate strength with GLN200 with bond length of 3.1Å. Presence of methyl group induce positive mesomeric effect by donating electrons. Moreover, acetate group of molecule was making alkyl and pi-alkyl interactions with PHE223 and VAL103. Amantadine ring was making van der Waals interactions with VAL222 and SER225. The detailed 3D and 2D binding interactions of derivative 5f within active pocket of PPE enzyme is shown in Fig. 8 . The amino acid residues involved in bonding and non-bonding interactions with 5f TRP179, SER225, VAL224, PHE223, GLN200, CYS199, SER203, SER222 and ARG226. Briey ouro substituted benzoylimino ring of parent compound was making strong p-p stacked interaction with PHE223. Compound 5f showed good binding energy and IC 50 value which might be due to substitution of uorine on benzoyl imino ring. Fluorine induce positive mesomeric effect by donating electrons. Moreover, amantadine ring was making alkyl interaction with VAL103. In addition, single hydrogen-acceptor bonding was observed with ARG226 with bond length of 2.9Å. The van der Waals interactions were exhibited by oxothiazolidinylidene ring with TRP179 and CYS199.

The detailed 3D and 2D binding interactions of derivative 5g within active pocket of PPE enzyme is shown in Fig. 9 . The amino acid residues involved in bonding and non-bonding interactions with 5g were PHE223, ARG226, VAL103, ASP63, THR100, HIS60, ASP202, SER222, SER203, THR221, GLN200, CYS199 and VAL224. Briey, compound 5g was making two hydrogen bonds with VAL224 and GLN200 with bond length of 2.95Å and 3.21Å respectively. Interestingly, compound 5g doesn't bear any substitution on benzoylimino group, but it showed excellent in vitro IC 50 value and exhibited good binding energy during molecular docking studies. It can be predicted that presence of stabilizing hydrogen bonds with good cutoff value might have increased the activity against PPE enzyme. Moreover, central amantadine ring was making hydrophobic interaction with SER222. Other amino acids involved in van der Waals interactions were SER203 and SER103.

The detailed 3D and 2D binding interactions of derivative 5h within active pocket of PPE enzyme is shown in Fig. 10 . The amino acid residues involved in bonding and non-bonding interactions with 5h were SER225, GLN200, SER203, SER222, HIS60, VAL224, VAL103, PHE223, TRP179 and ARG226. Briey, 5h compound is appeared to be most potent derivative which showed highest binding energy of À28.82 kJ mol À1 and excellent in vitro activity (IC 50 value). It might be due to substitution of chlorine on benzoyl imino ring of compound. Substituted chlorine was making strong p-alkyl interaction with TRP179. Also, benzoyl imino ring was making p-p stacked interaction with PHE223. Moreover, it was observed that amantadine ring was playing pivotal role in making hydrophobic interactions with amino acid residues f active site. Amantadine ring of 5h compound was producing p and p-alkyl interaction with VAL103. Whereas acetate group of 5h compound was involved in p-alkyl interaction with VAL224. Furthermore, compound 5h was producing hydrogen bond with SER225 having bond length of 2.3Å. The amino acid residues involved in van der Waals interaction with 5h compound were GLN200, SER203, SER222 and HIS60.

The co-crystal ligand was re-docked into active pocket of elastase enzyme which was aimed to validate the docking protocol and evaluation of binding score and binding affinity. It was observed that co-crystal ligand was producing important bonding and nonbonding interactions with targeted protein.

The amino acid residues which were involved in important stabilizing interactions were THR221, SER203, HIS60, SER222, PHE223, VAL224, SER225, TRP179, GLN200, CYS199, GLY201, GLY198, THR152, LEU227, CYS229, ALA104 and VAL103. Cocrystal ligand was involved in important stabilizing hydrogen bonding with GLN200, VAL224, PHE223 and HIS60. The docking score was found to be À28.03 kJ mol À1 . In comparison to synthesized derivatives, it was observed that compound 5h was producing better docking scores than co-crystal ligand. It can be predicted that synthesized derivatives are better inhibitors of PPE enzyme than co-crystal ligand. The 3D interactions of cocrystal ligand with targeted protein is shown in Fig. 11 .

Molecular dynamic simulations: analysis of RMSD and RMSF. The root mean square deviation (RMSD) is a metric that indicate the average change in displacement of atoms with respect to frame. The value of RMSD can be obtained by analysis of MD trajectory with respect to time. Y-Axis of a graph depicts the evolution of protein RMSD and x-axis indicate the course of time in nanosecond. The RMSD of ligand, targeted protein and protein ligand complex is shown in Fig. 12 .

The RMSD of ligand is plotted by aligning protein-ligand on the reference protein backbone and RMSD of ligand heavy atoms were calculated. Fig. 12 depicts the evolution of RMSD of protein, ligand and protein-ligand complex. The RMSD value of protein ligand shows that complex attained stability initially up to 30 ns and aer that it showed uctuations in RMSD up to 3 angstrom. Aer 60 ns of simulation time, complex again stability and remained around 2 angstrom which is perfectly acceptable. Aer 80 ns, complex RMSD dropped to 1.2 angstrom and remained stable around 1.5 angstrom. In case of targeted protein, RMSD value was quite stable and it attained stability aer 10 ns and remained stable around 1 angstrom which is perfect. Green line is indicating evolution of ligand RMSD which is showing uctuations within 2 angstrom and aer 80 ns of simulation time, it got stable at 2 angstrom.

The local changes around protein chain are characterize by root mean square uctuation (RMSF) values. Fig. 13 depicting the RMSF value of residues index of targeted protein.

Peaks sections are showing those residues which uctuate more during MD simulation. Usually a residue belonging to N and C terminal tends to uctuate more whereas more compact regions like secondary structures part are more rigids and show less uctuation. The stability of ligand complexed with protein are indicated by low RMSF values of amino acid residues of binding site. MD simulations demonstrated total 13 hydrogen bonds were produced between ligand and protein. Hydrogen bonds formed during simulations are tabulated in Table 5 .

ADMET properties of compounds (5a-j). Prediction of ADMET properties is prime step in the process of drug discovery and drug development. In silico ADMET evaluation models were developed to assist the scientist in drug design and lead optimization. Physicochemical properties, absorption, distribution, metabolism, toxicity, excretion, and medicinal properties of adamantyl hybrid iminothiazolidinones derivatives are calculated through the online web server ADMET lab 2.0. It is an integrated online platform for accurate and comprehensive prediction of ADMET properties. ADMET factors that were considered are physicochemical properties, blood-brain barrier (BBB), Caco-2 permeability, volume of distribution (VD), Pglycoprotein (PGP) substrate, P-glycoprotein (PGP) substrate, plasma protein binding, human intestinal absorption (HIA), MDCK permeability, clearance (CL), half-life (T 1/2 ), eye corrosion, eye irritation, respiratory toxicity, AMES toxicity, carcinogenicity, and synthetic accessibility score. Total ten adamantyl hybrid iminothiazolidinones derivative were subjected ADMET study using online web server ADMETlab 2.0. It was observed that physicochemical properties of all compounds were meeting the criteria of drug like rule i.e., Lipinski rule of ve (Table 6) .

Moreover, it was observed that all compounds showed a positive value of Caco-2 permeability. Caco-2 permeability is human colon epithelial cell line which corresponds to better human intestinal absorption. A drug having a Caco-2+ value corresponds to high intestinal absorption. The calculated value for human intestinal absorption (HIA) showed that all compounds have probability of good absorption from intestinal membrane. Furthermore, the calculated value for human intestinal absorption was better than standard oleanolic acid for compound 5h and 5i. A compound having positive value of blood-brain barrier have better lipophilicity prole and it can absorb readily from plasma membranes. Calculated values for blood-brain barrier (BBB) and blood placental barriers (BPB) for all compounds were found to be better than standard. It was noticed that the compounds 5a, 5d, 5h and 5i showed better lipophilicity prole and values was observed as BBB++. In terms of Plasma glycoprotein (PGP) substrate and PGP-inhibitor, it was observed that the output value of all compounds had the probability of being PGP-substrate, and PGP-inhibitor. Plasma protein binding (PPB) is an important factor to determine safety of drugs, as drug having high value of PPB (>90%) have narrow therapeutic index whereas drugs having low value of PPB are much safer. In the present study, compounds 5b, 5c, 5f, 5g, and 5j showed low PPB values, which means that these compounds have a broad therapeutic index. In terms of carcinogenicity, compounds 5f, 5h, and 5i were found to be non-carcinogenic whereas all other compounds have the probability of being carcinogenic. AMES toxicity prole depicted that compounds have probability of being toxic. Synthetic accessibility score (SA score) is designed to estimate ease of synthesis of drug-like molecules, and it was observed that all compounds have positive value of SA score. Overall, all compounds showed better ADMET prole, all values are tabulated in Table 7 given below.

Solvent and reagents were purchased from Sigma-Aldrich and dried before use. The melting point was nding a digital Gallenkamp (SANYO) model MPD BM 3 apparatus and was uncorrected. 1 H NMR and 13 C NMR spectra were determined in CDCl 3 or acetone-d 6 solutions at 300 MHz and 75.4 MHz, respectively using a Bruker AM-300 spectrophotometer. FTIR spectra were recorded using an FTS 3000 MX spectrophotometer. The chemical shis were measured in ppm units and coupling constants (J) in hertz.

Substituted aryl carboxylic acid (3.97 mmol, 1.2 mol eq.) was reacted with thionyl chloride (4.30 mmol, 1.3 mol eq.) under inert conditions at reux in dichloromethane. When the reaction was completed, the dichloromethane was removed under reduced pressure. The resulting acid chloride was converted to isothiocyanate by adding potassium thiocyanate (4.96 mmol, 1.5 mol eq.) in acetone. Amantadine (3.31 mmol, 1 mol eq.) was added to the reaction mixture and heated to 40 C for 1.5 hours. The reaction mixture was pure in crushed ice the crude product precipitated out and recrystallized in ethanol. The obtained acyl thioureas (1.27 mmol, 1 mol eq.) was further reacted with ethyl 4-ethoxypent-4-en-2-ynoate (2.54 mmol, 2 mol eq.) at room temperature. The nal product precipitated out during reaction. The product was ltered, dried, recrystallized in ethanol. Other compounds were prepared using the same method. SER39-side 5h-side 10 GLY38-main 5h-side 11

5h-side ASN153-side 12

5h-side SER225-main 13 THR44-side 5h-side In vitro methodology Elastase inhibition assay. The inhibition activity of elastase (elastase from porcine pancreas) was determined by following the methods of ref. 34 and 35 with minor modications as discussed in our previously reported article. 36 By measuring the absorbance of the solution at 410 nm, the amount of released pnitroaniline, which had been hydrolyzed by the substrate (Nsuccinyl-Ala-Ala-Ala-p-nitroanilide), could be used to determine the elastase inhibition activity. It was necessary to prepare a 0.8 mM solution of the N-succinyl-Ala-Ala-Ala-p-nitroanilide in a 0.2 M Tris-HCl buffer in order to conduct the experiment (pH 8.0). The buffer solution (130 mL) was added to the test sample (10 mL) in a 96-well microplate, and the mixture was mixed thoroughly. The microplate was pre-incubated for 10 minutes at 25 C before the addition of an commercially available elastase stock solution (10 mL, 0.0375 unit per mL) was performed. Following the addition of the enzyme, the microplate was maintained at 25 C for 30 minutes, and the absorbance at 410 nm was measured using a microplate reader (SpectraMax ABS, USA). The IC 50 values were determined through the series of triplicate experiments with eight different concentrations as a serial dilution using default protocol in the SoMax pro soware. Using the following formula, it was possible to calculate the elastase inhibition activities:

Elastase inhibition activity (%) ¼ (OD control À OD sample  100)/ OD control where OD control and OD sample represents the optical densities in the absence and presence of sample, respectively. Oleanolic acid was used as the standard inhibitor for elastase. The IC 50 of the derivatives was calculated by assaying the eight different concentrations of the tested derivatives (serial dilutions). The observed values were analyzed through non-linear calculations using graph pad prism soware 5 as discussed in previously reported article. 35 Protocol for kinetics

To determine the mode of inhibition, a kinetic analysis was performed using our previously published method 34 in the following manner: the compound 5g was chosen based on its high IC 50 value and high potency. It was possible to study the kinetics of compound 5g by varying the concentration of Nsuccinyl-Ala-Ala-Ala-p-nitroanilide in the presence of various concentrations of compound 5g (0.00, 0.062, 0.124, and 0.248 mM). Briey, the N-succinyl-Ala-Ala-Ala-p-nitroanilide concentration was changed from 2, 1, 0.5, 0.25, 0.125, and 0.0625 mM for kinetics studies, and the rest of the procedure was the same as described in the elastase inhibition assay protocol. The maximum initial velocities were determined from the initial linear portion of absorbance up to 10 minutes aer the addition of the enzyme at a rate of one minute intervals until the enzyme was removed. The type of inhibition on the enzyme was determined by plotting the inverse of velocities (1/V) against the inverse of substrate concentration 1/[S] mM À1 on a Lineweaver-Burk plot. SoMaxPro was used to process the data aer it was collected. An analysis of a secondary plot of 1/V against inhibitor concentration revealed that the E I dissociation constant K i was observed.

Quantum chemistry through density functional theory calculation. The density functional theory (DFT) is one of widely used method to determine the electron density and energetic parameter of the compounds. It is used to determine the structure of atoms, molecules, crystals, surfaces, and their interactions. These calculations are done by using the Gaussian 09W program on the Linux based workstation. 37 B3LYP method was used to determine vibrational wavenumbers with an SVP basis set. It is known that the B3LYP functional assists a good description of harmonic vibrational wavenumbers for small and medium-sized compounds. The optimized structural parameters are further use for calculations. The output check les were analyzed by GuassView 6.0. 38 Molecular docking studies. Molecular docking studies were performed using Molecular Operating Environment (MOE) 2015.10 soware package. 39 The docking studies were conducted to evaluate and investigate binding energies and binding affinities of synthesized derivatives within activation loop of targeted protein. The crystallographic structure of targeted protein porcine pancreatic elastase (PDB ID ¼ 1 BMA; resolution 1.99 angstrom) was obtained from protein data bank (https://www.rcsb.org/). 40 All the steps for protein preparation that are pre requisite for docking procedure were performed. Initially co-crystal ligand and het atoms were removed moreover, water molecules were also removed as they interfere with ligand molecules during docking calculations. The addition of polar hydrogen atoms with their standard 3D geometry, addition of partial Gasteiger charges using MMFF94x as a force eld was undertaken. Furthermore, xation of errors in connections, repairing missing residues, 3D protonation, and xation of potential and energy minimization was also performed using protein preparation wizard of MOE soware. The selection of same active site amino acid residues of co crystal inhibitor was carried out using site nder utility of MOE 2015.10. Dummies atoms were created at the site of selected amino acid residues. The structures synthesized thiazolidinone derivatives (5a-j) were generated using Chemdraw ultra 12.0. The generated structures were optimized to steepest energy gradient prior to docking with targeted protein. The ligand database comprising of synthesized thiazolidinone derivatives were docked into selected amino acid residues of active site. Docking protocol was as follows, placement method was set to triangular matcher, scoring methodology was set to London dG, renement methodology was set to rigid receptor and GBVI/WSA dG was the scoring function for selection of best pose from 100 poses for each compound. Aer completion of docking protocol, best poses for each compound were evaluated for their ttings into active pocket and best protein-ligand conformations were selected on the basis of docking score for further analysis. The analysis and visualization of best docked pose was carried out using PyMol. 41 The putative 2D ligand interactions were created using LigPlot + 2.2. 42 It is run from intuitive java interface which provide detailed 2D interactions of proteinligand complex. Molecular docking protocol was validated by redocking co-crystal ligand benzyl methyl aminimide inhibitor with amino acid residues of active site.

ADMET properties. Determination of biochemical processes from drug administration to elimination plays a vital role in lead optimization. An ideal drug candidate should be administered and absorbed into the systemic circulation and must be non-toxic and eliminate without affecting the biological activity. These processes are seeming to be distinct, but they are closely inter-related, so the determination of ADMET properties have prime importance in the drug discovery process. Since traditional methods were time-consuming and much research till 1990 went into vain due to the appearance of undesirable effects in the middle of the drug discovery process. So, the determination of ADMET properties is of prime importance to rolling out undesired effects of a drug candidate at the initial stage of the drug discovery process. Efficient and reliable online prediction models developed for in silico determination of ADMET properties. ADMET properties like absorption, distribution, metabolism, excretion, toxicity and physicochemical properties of eight compounds were done using online in silico prediction model ADMET lab 2.0. 43 All compounds were converted to SMI format using ChemDraw Ultra and these SMI structures were loaded into online server ADMET lab 2.0. It also facilitates the manually generating the desired structures using JMSE editor.

Molecular dynamic simulations. The molecular dynamic study of top ranked conformation was performed using Nano scale molecular dynamic (NAMD) soware 44 on CUDAaccelerated GPU system having 16 core processor and 64 GB RAM memory. Visualization of MD simulations was carried out using VMD soware. 45 MD simulations were carried out in order to determine binding interactions and stability of proteinligand complex under accelerated conditions. Top ranked protein-ligand complex was selected and topology les were generated for protein and ligand using CHARMM36 forceeld. 46 The system was neutralized by addition of NaCl charges. The system was minimized to steepest-energy gradient to remove any close contacts between atoms. The system was equilibrated in NVT ensemble for 500 000 steps followed by equilibration in NPT ensemble for additional 500 000 steps. Aerward, simulation was performed for 100 ns under periodic boundaries condition. 47 The calculation of binding energy, van der Waals and electrostatic interactions were carried out using PME method. 48 

A structurally diverse series of adamantyl-iminothiazolidinone conjugates was synthesized, characterized, and tested for in vitro elastase inhibition assay and subjected to in silico ADMET prediction. The inhibition studies revealed compounds 5a, 5f, 5g, and 5h to show signicant activity. Compound 5g with no substitution on benzoyl ring having IC 50 value (IC 50 0.124 mM) elastase was 50 times more potent as compared to oleanolic acid having IC 50 5.996 mM. Kinetic studies showed competitive mode of inhibition and excellent binding energies. Overall, all compounds showed better ADMET prole, indicating their potential to serve as lead compounds in drug discovery. The comprehensive in silico studies including molecular docking studies, drug likeness properties and stability of top ranked protein-ligand complex was conducted to evaluate and support in vitro studies. All these studies suggested potent inhibitory potential of synthesized derivatives.

The authors declare no conict of interests.

82 (s, 1H, C]CH), 4.31 (q, 2H

3s,5s,7s)-adamantan-1-yl)-2-((3-methylbenzoyl) imino)-4-oxothiazolidin-5-ylidene)acetate (5b

300 MHz) d ppm, 8.21 (s, 1H, Ar-H), 8.10 (d, 1H, J ¼ 7.8 Ar-H), 7.58 (d, 1H, J ¼ 7.8 Ar-H), 7.45 (t, 1H, J ¼ 7.8 Ar-H)

White solid: yield: 82%; R f : 0.58

41 (t, 1H, J ¼ 8.1 Hz, Ar-H), 7.16 (dd, 1H, J ¼ 7.8, 2.1 Hz, Ar-H), 6.85 (s, 1H

bromobenzoyl) imino)-4-oxothiazolidin-5-ylidene) acetate (5d)

CDCl 3 , 300 MHz) d ppm, 8.07 (d, 2H

3s,5s,7s)-adamantan-1-yl)-2-((3-chlorobenzoyl) imino)-4-oxothiazolidin-5-ylidene) acetate (5e

CDCl 3 , 300 MHz) d ppm, 8.21 (s, 1H, Ar-H) 8.08 (d, 1H

87 (s, 1H, C]CH), 4.32 (q, 2H, J ¼ 7.2 Hz, CH 2 ), 2.77 (s, 6H, CH 2 , adamantyl), 2.28 (s, 3H, CH adamantyl), 1.82-1.80 (q, 6H

Ar-C)

White crystalline solid: yield: 87%; R f : 0.48; mp 231-233 C; IR: 2980, 2907, 2850 (aliphatic C-H stretching

C-S) cm À1 . 1 H NMR (300 MHz, CDCl 3 , 300 MHz) d ppm

26 (s, 3H, CH adamantyl), 1.86-1.74 (q, 6H, J ¼ 12.3 Hz, CH 2 , adamantyl), 1.35 (t, 3H, J ¼ 7.2 Hz, CH 3 ), 13 C NMR (75 MHz, CDCl 3 , 300 MHz) d ppm, 175.2 (C]O of benzoyl)

77 (d, 6H, J ¼ 2.4 Hz, CH 2 , adamantyl), 2.26 (s, 3H, CH adamantyl), 1.80 (q, 6H, J ¼ 12 Hz, CH 2 adamantyl), 1.35 (t, 3H, J ¼ 7.2 Hz, CH 3 ). 13 C NMR (75 MHz, CDCl 3 , 300 MHz) d ppm

CDCl 3 , 300 MHz) d ppm, 8.15 (dd, 2H, J ¼ 6.9 Hz, J ¼ 1.8 Hz Ar-H), 7.49 (dd, 2H, J ¼ 6.6 Hz, J ¼ 1.8 Hz Ar-H)

CH 2 adamantyl), 1.35 (t, 3H, J ¼ 7.2 Hz, CH 3 ). 13 C NMR (75 MHz, CDCl 3 , 300 MHz) d ppm, 175.5 (C]O of benzoyl)

87 (s, 1H, C]CH), 4.31 (q, 2H, J ¼ 9 Hz, CH 2 ), 2.71 (d, 6H, J ¼ 2.4 Hz, CH 2 , adamantyl), 2.22 (s, 3H, CH adamantyl), 1.75 (q, 6H, J ¼ 12 Hz, CH 2 adamantyl), 1.35 (t, 3H, J ¼ 7.2 Hz, CH 3 ). 13 C NMR (75 MHz, CDCl 3 , 300 MHz) d ppm

Ar-H), 6.83 (s, 1H, C]CH), 4.30 (q, 2H, J ¼ 9 Hz, CH 2 ), 3.90 (s, 3H, O-CH 3 ), 2.76 (d, 6H, J ¼ 2.4 Hz, CH 2 , adamantyl), 2.26 (s, 3H, CH adamantyl), 1.80 (q, 6H, J ¼ 12 Hz, CH 2 adamantyl), 1.34 (t, 3H, J ¼ 7.2 Hz, CH 3 ). 13 C NMR (75 MHz, CDCl 3 , 300 MHz) d ppm

Ar-C)

Large scale simulations of particulate systems using the PME method

The authors gratefully acknowledge the Quaid-i-Azam University for environment conducive for research.