Exploring 2-Tetradecanoylimino-3-aryl-4-methyl-1,3-thiazolines Derivatives as Alkaline Phosphatase Inhibitors: Biochemical Evaluation and Computational Analysis

The current study focused on the laboratory approach in conjunction with computational methods for the synthesis and bioactivity assessment of unique 2-tetradecanoylimino-3-aryl-4-methyl-1,3-thiazolines (2a–2k). Processes included cyclizing 1-aroyl-3-arylthioureas with propan-2-one in the presence of trimethylamine and bromine. By using spectroscopic techniques and elemental analyses, structures were elucidated. To assess the electronic properties, density functional theory (DFT) calculations were made, while binding interactions of synthesized derivatives were studied by the molecular docking tool. Promising results were found during the evaluation of bioactivity of synthesized compounds against alkaline phosphatase. The drug likeliness score, an indicator used for any chemical entity posing as a drug, was within acceptable limits. The data suggested that most of the derivatives were potent inhibitors of alkaline phosphatase, which in turn may act as lead molecules to synthesize derivatives having desired pharmacological profiles for the treatment of specific diseases associated with abnormal levels of ALPs.


Introduction
Alkaline phosphatases (ALPs) are the isozymes concentrated in the outer layer of the cell membrane. The clinical relevance of ALPs in various diseases, especially hepatic disorders, makes them an important area of research among medicinal chemists. Elevated ALP levels are considered in the diagnostic procedure of a wide variety of bile duct and liver-related morbidities, while abnormally low levels are presented in Wilson disease [1]. Alkaline phosphatase, a subclass of the superfamily of hydrolase enzymes, acts as a catalyst in the hydrolysis of various phosphates, pyrophosphate, sulphate, and sulphonate esters. Four alkaline phosphatase isozymes are expressed in the human body, including three tissue-specific ALPs and one non-tissue-specific ALP (TNAP). Tissue-specific ALPs are localized in placenta, germ cells, and intestine while tissue-nonspecific ALPs are concentrated in bone, kidney, liver, and the central nervous system [2][3][4][5][6][7]. Human serum also contains ALPs and has been used in the diagnosis of liver disorders, bone diseases, and multiple carcinomas [7][8][9]. Human intestinal ALPs (h-IAPs) are located at the intestinal brush border. Important roles of h-IAPs include -HCO 3 secretion to maintain pH, intestinal lipopolysaccharides detoxification, and the regulation of lipid absorption [10,11]. Decreased circulating level of ALPs results in inflammation due to weakened effects of endotoxin detoxification [12], while upregulation occurs in inflammatory bowel disease. Limited numbers of h-IAPs are known, owing to the fact that the structural similarity among TNAP and IAP is higher [13]. Due to wider implications of alkaline phosphatase, the search for a selective and potent inhibitor of ALPs is of greatest importance. A literature review indicated the efficient potential of thiazole derivatives as alkaline phosphatase inhibitors [14,15].
In the search for potent, selective, and drug-like molecules of intestinal alkaline phosphatase inhibitor intended to be used in human populations, the synthesis of 2tetradecanoylimino-3-aryl-4-methyl-1,3-thiazolines derivatives was carried out.
In the present study, new 2-tetradecanoylimino-3-aryl-4-methyl-1,3-thiazolines were synthesized by considering the advantage of avoidance of unwanted side products and overall high yields. With the acetone cyclization of corresponding 1-tetradecanoyl-3arylthiourea taking place, the reaction takes place in the presence of bromine and triethyl amine.
To gain insights into geometric and electronic behaviors, density functional theory (DFT) studies were performed. In vitro free radical scavenging activity and alkaline phosphatase inhibition activity were determined. To study the binding interactions of synthesized derivatives with alkaline phosphatase, molecular docking was performed. Drug likeliness, an important characteristic of any chemical entity to be purposed as a medicinal agent, was calculated using an online tool. In the search for potent, selective, and drug-like molecules of intestinal alkaline phosphatase inhibitor intended to be used in human populations, the synthesis of 2-tetradecanoylimino-3-aryl-4-methyl-1,3-thiazolines derivatives was carried out.
In the present study, new 2-tetradecanoylimino-3-aryl-4-methyl-1,3-thiazolines were synthesized by considering the advantage of avoidance of unwanted side products and overall high yields. With the acetone cyclization of corresponding 1-tetradecanoyl-3-arylthiourea taking place, the reaction takes place in the presence of bromine and triethyl amine.
To gain insights into geometric and electronic behaviors, density functional theory (DFT) studies were performed. In vitro free radical scavenging activity and alkaline phosphatase inhibition activity were determined. To study the binding interactions of synthesized derivatives with alkaline phosphatase, molecular docking was performed. Drug likeliness, an important characteristic of any chemical entity to be purposed as a medicinal agent, was calculated using an online tool.  Under a nitrogen gas atmosphere, bromine solution (0.1 mL) in dry acetone (10 mL) was added dropwise to a 1-tetradecanoyl-3-arylthiourea (0.5 g) solution under stirring in 20 mL dry acetone that contained 0.3 mL triethylamine. The next step was the overnight stirring of the solution at room temperature after the final reactant was added. To track the reaction's progress, the thin layer chromatographic technique was employed. To remove crude solids, the reaction mixture was concentrated and filtered. By recrystallizing with ethanol, derivatives of 2-tetradecanoylimino-3-aryl-4-methyl-1,3-thiazoline were collected. The scheme of synthesis in given in Figure 2.

Alkaline Phosphatase Assay
As already described by Jamshed et al. in 2011, the spectrophotometric technique was used for the assessment of calf intestinal alkaline phosphatase (CIALP) activity [41]. The reaction mixture included 50 mM Tris-HCl buffer with 0.1 mM ZnCl 2 and 5 mM MgCl 2 (pH 9.5) and test substance (0.1 mM with final DMSO 1% (v/v). Then 0.025 U/mL CIALP was added to the pre-incubated mixture. Mixing was conducted for 10 min. The reaction mixture was then added with 10 µL of the substrate, having 0.5 mM p-NPP (para nitrophenylphosphate disodium salt). It was once more incubated for 30 min at 37 • C. A change in absorbance of the released p-nitrophenolate was observed at 405 nm with a 96-well microplate reader (OPTI MAX, Milton Freewater, OR, USA). Triplicate independent repetition of experiments was done. As a reference inhibitor of calf ALP, KH 2 PO 4 was used.

Free Radical Scavenging Assay
By modifying the previously described method [42] the 2,2-diphenyl-1 picrylhydrazyl (DPPH) assay was used to assess the radical scavenging activity. The reaction mixture contained 20 µL of test compounds at progressively higher concentrations, 100 µL of DPPH (150 µM), and 200 µL of methanol to adjust the volume in each well. At room temperature, the reaction mixture was let to sit for another 30 min. The reference inhibitor utilized was ascorbic acid (also known as vitamin C). Measurements were taken using an OptiMax tunable microplate reader for the wavelength of 517 nm comparing reaction rates and calculating the % inhibition brought on by the presence of the tested inhibitors. The independent repetition of experiments was done in a triplicate manner.

Density Functional Theory (DFT) Calculations
Density functional theory studies are a computational quantum mechanical modelling technique used for the prediction of the electronic properties of chemical compounds. For quantum mechanical calculations, frontier molecular orbitals (HOMO, LUMO), the molecular electrostatic map, and electron density were measured. DFT studies were conducted using Gaussian 16 on a workstation having an AMD Ryzen 9 5950x @ 16 core processor with 64 GB of installed RAM memory (Santa Clara, CA, USA). For representing the quantum mechanics region, the density functional theory (DFT) method is critical; all calculations were done using the B3LYP/6-311G basis set. Visualization of calculations were performed using Gauss view 6 [43,44].

Molecular Docking Studies
The Protein Data Bank (PDB) https://www.rcsb.org/ was used to retrieve the threedimensional (3D) structure of human placental alkaline phosphatase, having the PDB ID: 1ZED, accessed on 18 September 2022. Utilizing the Discovery Studio Visualizer 2021 software, the chosen protein structure was minimized.
The structures of all synthesized derivatives and the reference compound p-nitrophenylphosphonate were drawn using ChemDraw 12.0 followed by energy minimization by using Chem3D Pro 12.0. In addition to an active pocket, alkaline phosphatases have peripheral allosteric binding sites, which play roles in enzyme regulation. Keeping in mind the enzyme structural studies, dimensions of p-nitrophenyl-phosphonate (PNP) were used for molecular docking studies. The reason for the selection of p-nitrophenyl-phosphonate as the reference ligand is its binding capability with placental alkaline phosphatase [45].
The dimensions of PNP 902, presented as a ligand of interest, were taken, i.e., x = +57.430500, y = 28.435571, and z = −2.404500 with size values of X = 40, Y = 40, and Z = 40. The value of 8 was taken as the default exhaustiveness. Synthesized ligands (2a-2k) were docked using AutoDock Vina [46] against alkaline phosphatase. The predicted docked complexes were analyzed using the values of binding energy (Kcal/mol). Ligand-protein interactions were visualized in three dimensions (3D) and two dimensions (2D) using Discovery Studio 21 [47].

Chemistry
In accordance with reported method [48], a series of 2-tetradecanoylimino-3-aryl-4methyl-1,3-thiazolines derivatives was synthesized by reacting 1-tetradecanoyl-3-arylthioureas (1a-1k) with bromine in dry acetone in the presence of triethyl amine to afford the cyclized 1,3-thiazolines derivatives (2a-2k), as delineated in Figure 2. The synthesized compounds were characterized by FT-IR and NMR spectroscopy. FT-IR spectra of synthesized compounds exhibited absorptions for the C-H aromatic at 3020 cm −1 , C-H thiazoline at 2917 cm −1 , and C=O at 1642 cm −1 . The 1H NMR spectra of compound (2a-2k) contained the characteristic one proton singlet of the thiazoline ring at δ 6.40 ppm. The alkyl protons gave rise to signals at δ 2.50-0.86, whereas the aromatic protons appeared in the range δ 7.53-7.21 ppm. In 13C NMR spectra, the signal for carbonyl carbon appeared at δ 182.1 ppm, and for imine carbon of the thiazoline ring, it was observed at δ 169.

Alkaline Phosphatase Assay
All the synthesized compounds (2a-2k) were evaluated for their inhibition potential to calf intestinal alkaline phosphatase (CIAP). The standard used was KH 2 PO 4 , and the results of evaluation as half maximal inhibitory concentration (IC 50 ) are summarized in Table 1.

Structure-Activity Relationship
Remarkably, most of the synthesized compounds exhibited many-fold inhibition potential towards alkaline phosphatase compared to the standard. Compounds 2a, 2k, and 2c showed better and more efficient inhibition potential, as represented by their respective IC 50 values. Compound 2h, with two nitro atoms at the 2-and 4-positions of the phenyl ring, was found to be the least potent due to their electron withdrawing nature. Compound 2k, containing a 4-methyl group, was found to be moderately potent. Compounds with cyano and carboxylic and sulfonic acid substitutions at the 4-position showed various degrees of inhibition. The derivative 2f, having a chlorine atom at the 3-position, was found to be the most potent among the series, which may be attributed to the electron-withdrawing inductive effect of the chlorine atom, which perturbs the electron density and possibly results in the enhancement of alkaline phosphatase activity. In general, compounds with an electron donating nature showed better inhibition of CIAP as compared to those with an electron withdrawing nature.

Free Radical Scavenging Activity
Radical scavenging activities of compounds (2a-2k) were compared, with vitamin C used as a reference. Compound 2a exhibited maximum potency, comparable to the standard ascorbic acid. Compounds 2a, 2b, 2d, and 2i also showed significant radical scavenging potency. The remaining compounds were quite inactive. % age free radical scavenging activity is graphically represented in Figure 3.

Structure-Activity Relationship
Remarkably, most of the synthesized compounds exhibited many-fold inhibition potential towards alkaline phosphatase compared to the standard. Compounds 2a, 2k, and 2c showed better and more efficient inhibition potential, as represented by their respective IC50 values. Compound 2h, with two nitro atoms at the 2-and 4-positions of the phenyl ring, was found to be the least potent due to their electron withdrawing nature. Compound 2k, containing a 4-methyl group, was found to be moderately potent. Compounds with cyano and carboxylic and sulfonic acid substitutions at the 4-position showed various degrees of inhibition. The derivative 2f, having a chlorine atom at the 3-position, was found to be the most potent among the series, which may be attributed to the electronwithdrawing inductive effect of the chlorine atom, which perturbs the electron density and possibly results in the enhancement of alkaline phosphatase activity. In general, compounds with an electron donating nature showed better inhibition of CIAP as compared to those with an electron withdrawing nature.

Free Radical Scavenging Activity
Radical scavenging activities of compounds (2a-2k) were compared, with vitamin C used as a reference. Compound 2a exhibited maximum potency, comparable to the standard ascorbic acid. Compounds 2a, 2b, 2d, and 2i also showed significant radical scavenging potency. The remaining compounds were quite inactive. % age free radical scavenging activity is graphically represented in Figure 3.

Density Functional Theory (DFTs) Calculations
DFT calculations of 2-tetradecanoylimino-3-aryl-4-methyl-1,3-thiazolines derivatives (2a-2k) including optimization, global reactivity descriptors, and FMOs analysis were made using the B3LYP/6-311G basis set. All derivatives were optimized to the steepest energy gradient. Optimized structures of the most potent compounds are given in Figure 4, while other structures are given in the supplementary information ( Figure S9 from Supplementary Materials).  Optimization energy, polarizability, and dipole moment play significant roles in determining the chemical reactivity of a compound. The analysis of the HOMO and LUMO energies facilitated the understanding of the various molecular properties as well as electrical and optical properties of the molecule. Moreover, the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) energy gap also depict the chemical reactivity of compounds. Compounds having small HOMO and LUMO energy gaps are referred to as reactive molecules, whereas large energy gaps correspond to stability and less reactivity. Table 2 shows the optimization energies, dipole moments, polarizability, HOMO and LUMO energies, and their energy gaps. Compound 2d, which has the highest energy gap (∆Egap = 0.182) tended to be a stable molecule. Similarly, compound 2h had the lowest energy gap (∆Egap = 0.092), justifying its highest reactivity. The molecule with the highest HOMO energy, 2a (EHOMO = −0.226), turned out to be the best electron donor, while compound 2k, which had the lowest ELUMO of −0.024 was shown to be the best electron acceptor.  Optimization energy, polarizability, and dipole moment play significant roles in determining the chemical reactivity of a compound. The analysis of the HOMO and LUMO energies facilitated the understanding of the various molecular properties as well as electrical and optical properties of the molecule. Moreover, the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) energy gap also depict the chemical reactivity of compounds. Compounds having small HOMO and LUMO energy gaps are referred to as reactive molecules, whereas large energy gaps correspond to stability and less reactivity. Table 2 shows the optimization energies, dipole moments, polarizability, HOMO and LUMO energies, and their energy gaps. Compound 2d, which has the highest energy gap (∆Egap = 0.182) tended to be a stable molecule. Similarly, compound 2h had the lowest energy gap (∆Egap = 0.092), justifying its highest reactivity. The molecule with the highest HOMO energy, 2a (EHOMO = −0.226), turned out to be the best electron donor, while compound 2k, which had the lowest ELUMO of −0.024 was shown to be the best electron acceptor.
HOMO and LUMO orbitals of the most potent compounds are visualized in Figure 5, while HOMO and LUMO orbitals of the remaining compounds are shown in Figure S10 in the Supplementary Materials.
Global chemical reactivity descriptors By using the highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) energy values, we evaluated the following parameters by using their respective formulas [49]. Global reactivity descriptor's values are given in Table 3.
Compound 2h had the highest chemical potential and electronegativity. Similarly, compound 2d was the harder compound, whereas compound 2h was the softer compound. Compound 2a had the highest electrophilicity index.

Molecular Docking Studies Molecular Docking Analysis
Synthesized compounds (2a-2k) were docked against alkaline phosphatase for the prediction of the best conformational position in the active pocket of the target protein. On the basis of energy values (Kcal/mol) and hydrogen/hydrophobic interactions patterns, docked complexes were analyzed. Binding energy values of synthesized compounds and reference ligand are given in Table 4. Due to the fact that all synthesized compounds had the same basic chemical nucleus, significant energy difference are not shown. HOMO and LUMO orbitals of the most potent compounds are visualized in Figure 5, while HOMO and LUMO orbitals of the remaining compounds are shown in Figure S10 in the Supplementary Materials. Global chemical reactivity descriptors By using the highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) energy values, we evaluated the following parameters by using their respective formulas [49]. Global reactivity descriptor's values are given in Table 3.
Compound 2h had the highest chemical potential and electronegativity. Similarly, compound 2d was the harder compound, whereas compound 2h was the softer compound. Compound 2a had the highest electrophilicity index.  2a, 2c, 2f, 2k). Red, yellow, and blue atoms indicate O, S, and N atoms.

Binding Interactions
The 3D and 2D interactions of the most potent ligand, 2f, are shown in Figure 6.    Amino acid residues present in the active pocket that are responsible for binding interactions are: PHE208, THR212, PRO213, ASP214, PRO215, ASN232, VAL234, GLN235, LEU238, TYR246, VAL247, and TRP248. Figure 6 shows that all ligands were contained in the active pocket of alkaline phosphatase. Amino acid residue TRP248 was involved in three types of interactions, i.e., conventional hydrogen bond with carbonyl oxygen, Pi-Pi T-shaped interaction with a thiazole ring, and alkyl interaction with an aliphatic side chain. TYR246 contributed in the Amino acid residues present in the active pocket that are responsible for binding interactions are: PHE208, THR212, PRO213, ASP214, PRO215, ASN232, VAL234, GLN235, LEU238, TYR246, VAL247, and TRP248. Figure 6 shows that all ligands were contained in the active pocket of alkaline phosphatase. Amino acid residue TRP248 was involved in three types of interactions, i.e., conventional hydrogen bond with carbonyl oxygen, Pi-Pi T-shaped interaction with a thiazole ring, and alkyl interaction with an aliphatic side chain. TYR246 contributed in the formation of the Pi sulphur bond with a sulphur atom of the thiazole ring and two alkyl interactions with the aliphatic side chain. Amino acid PHE208 formed Pi Sigma interactions with the methyl group attached with thiazole ring.VAL234 and LEU238 were involved in alkyl interactions with the side chain, while PRO213 and PRO215 took part in Pi alkyl interactions with phenyl ring and thiazole rings, respectively. In vitro studies and molecular docking studies suggested the potency of derivative 2f as an alkaline phosphatase inhibitor.
The 2D and 3D interactions of human placental ALP with reference ligand para nitrophenyl phosphonate are given in Supplementary File Figure S11 from Supplementary Materials.

Chemo-Informatics Analysis
The expected attributes including polar surface area (PSA), density, molar volume, logP, molecular weight, and RO5 were predicted to justify their drug-like behavior [50]. According to research results, for the assessment of drug absorption in the drug development process, calculation of PSA is crucial [51]. Lipophilicity and molar refractivity play important roles in binding to receptors, bioavailability, and cellular absorption. According to earlier research, the usual ranges of molecular weight of 160 to 480 g/mol, molar refractivity of 40 to 130 cm 3 and polar surface area (PSA) of 89 are considered optimal [52]. For the total number of atoms in the drug-like molecule, the permitted range is between 20 and 70 [53]. The results shown in Table 5 indicate that all the predicted chemo-informatic properties of synthesized compounds lied within the standard ranges except logP. Furthermore, compounds 2a-2k were validated by Lipinski's rule. Our results justified that, except for logP values, all other parameters were followed, so the Lipinski rule was fulfilled. Lipinski violation should be employed for the compounds with poor absorption. However, large number of drug molecules are present that violate RO5 [54].

Conclusions
In the present study, 11 synthetic derivatives of 2-tetradecanoylimino-3-aryl-4-methyl-1,3-thiazolines were synthesized and characterized. Lab data of the alkaline phosphatase inhibition assay revealed the higher potency of 2f, 2a, 2e and 2k. IC 50 values were correlated with molecular docking studies, which verified the lab results. Important information about electronic properties and global reactivity descriptors were obtained through DFT calculations. Finally, chemo-informatic analyses were made to justify their behavior as druglike molecules. Results were found within the acceptable limits. Thus, these compounds are an important addition among ALPIs. However further studies are required for their use as medicinal agents for the treatment of diseases linked with alkaline phosphatase abnormalities.