Exploration of newly synthesized azo-thiohydantoins as the potential alkaline phosphatase inhibitors via advanced biochemical characterization and molecular modeling approaches

In the current study, Azo-Thiohydantoins derivatives were synthesized and characterized by using various spectroscopic techniques including FTIR, 1H-NMR, 13C-NMR, elemental and HRMS analysis. The compounds were evaluated for alkaline phosphatase activity and it was observed that among all the synthesized compounds, derivative 7e exhibited substantial inhibitory activity (IC50 = 0.308 ± 0.065 µM), surpassing the standard inhibitor (L–Phenyl alanine, IC50 = 80.2 ± 1.1 µM). Along with this, these derivatives were comprehensively examined regarding the electronic properties and reactivity of the synthesized compounds using Density Functional Theory (DFT) calculations, where the results were found very promising and the synthesized compound were found stable. After that, SwissADME evaluations highlighted compounds for their favorable physicochemical properties, including solubility and drug-likeness. Molecular docking exhibited the strong binding affinities of 7f and 7e derivatives with intestinal alkaline phosphatase (IAP), further supported by Molecular Dynamics (MD) simulations. This comprehensive integration of experimental and computational approaches sheds the light on the potential therapeutic applications of the synthesized compounds. By providing a detailed investigation of these aspects, this research opens the avenues for the development of novel pharmacologically active compounds with diverse applications. Supplementary Information The online version contains supplementary material available at 10.1186/s13065-024-01149-8.

Exploration of newly synthesized azothiohydantoins as the potential alkaline phosphatase inhibitors via advanced biochemical characterization and molecular modeling approaches Introduction Alkaline phosphatase (ALP) is a group of isoenzymes that catalyze the hydrolysis of phosphate monoesters at basic pH [1].It is an ubiquitous membrane-bound glycoprotein located on the outer layer of the cell membrane [2].ALP is found in various tissues throughout the human body, including the liver, bone, placenta, kidneys, and intestines.The enzyme is involved in several biological processes, including bone mineralization, cellular signal transmission, and phospholipid metabolism.The reference intervals for ALP levels in the blood of adults and females at 37 °C can be found in clinical chemistry reference ranges [3].ALP tests are commonly used to diagnose liver damage or bone disorders.
Among the types of ALP, the tissue-nonspecific alkaline phosphatase (TNAP) is an enzyme that is found in various tissues throughout the human body [4], including the liver, bone, placenta, kidneys, and intestines.It plays a significant role in regulating mineralization of bone tissue, as well as phospholipid metabolism and cellular signal transmission.Any mutation in the TNAP gene can lead to a rare disorder called hypophosphatasia [5], which is characterized by poor mineralization of bones and teeth, as well as other health complications.Additionally, TNAP is a marker for liver damage, and ALP tests are commonly used to assess liver function.TNAP has also been implicated in regulating purinergic signaling, which plays a role in immune response and inflammation.
Another type i.e., intestinal alkaline phosphatase (IAP) is a brush-border enzyme that is expressed throughout the gastrointestinal tract and is involved in several biological processes [6], including gut mucosal defense, regulation of gut homeostasis, and modulation of inflammatory mediators.IAP is secreted by enterocytes and is present on the apical surface [7].It plays an important role in detoxifying bacterial lipopolysaccharides (LPS), which are the major components of the bacterial cell wall and can cause inflammation and tissue damage.IAP has been shown to be decreased during conditions that commonly affect the gut [8], i.e. inflammatory bowel disease (IBD) and type-2 diabetes.It appears that IAP supplementation may have therapeutic potential in the treatment of these conditions by reducing inflammation and protecting the gut barrier [9].In addition, evidence to reveal that IAP may play a role in the regulation of mineral metabolism, although more research is needed in this area.From above it is evident that the optimum level of IAP is associated with the normal bone functioning and other cellular pathways [1].The role of IAP in disease conditions and the significance of their inhibitors have led to the discovery of numerous compounds with nitrogen-containing fused heterocyclic skeletons.These compounds have shown promising inhibitory values and many of them have also demonstrated anticancer activity [10].However, the analysis of the biological data from reported derivatives indicates that there is a lack of selectivity in the mode of inhibition within the same group [11].
The results of previous studies have generated interest in aromatic heterocycles with two nitrogen atoms in their rings.These compounds have been found to exhibit a diverse range of biological activities, including anticancer, analgesic, anti-inflammatory, antimicrobial, antiviral, anticonvulsant, antihistaminic, and anti-HIV properties [12][13][14].The significance of hydantoins inhibitor in medicinal chemistry has been emphasized by its recent success [15].Thiohydantoins are chemical compounds that are similar to hydantoins, but with sulfur atoms replacing one or both of the oxygen atoms in the carbonyl groups, as shown in Fig. 1 [16].Thiohydantoins and their derivatives have been the subject of extensive research for over 145 years, focusing on their chemistry and properties.The presence of the hydantoin moiety is significant as it serves as a crucial pharmacophore found in numerous biologically active compounds [17].
The current study is focused on synthesizing new derivatives of azo-thiohydantoins, taking into consideration the need for IAP inhibitors and the significance of azo and thiohydantoins moieties.The synthesized derivatives were investigated to determine their potential as inhibitors of the targeted enzyme, specifically IAP.The study utilized multiple computational methods, such as Molecular docking, Density functional theory (DFT) calculations, and molecular dynamic simulations, to support the obtained results.These approaches were employed to enhance the understanding of how these derivatives interact with the IAP active pocket.The ADMET profile of all the compounds was calculated in order to evaluate their safety level.

Material and method
Solvents were pre-treated and distilled as per conventional methods before usage.Reagents were procured from Sigma Aldrich and were used as received, without additional purification.The Stuart SMP3 apparatus was employed to determine the melting point.NMR spectra were taken on a Bruker 300, with 1 H-NMR at 300 MHz and 13 C-NMR at 75.5 MHz.Chemical shifts are presented in ppm, using tetramethylsilane or the solvent's residual resonance as a reference standard.Reactions' progress was tracked using thin layer chromatography (TLC).For TLC, we utilized aluminum sheets covered with silica gel F254 from Merck, and compounds were visualized under UV light at 254 and 360 nm wavelengths.

Synthesis protocol
To a stirred solution of azo-thiosemicarbazone (6) (1 mmol) in dry ethanol (20 mL), were added successively the Ethyl chloroacetate (1 mmol) and fused sodium acetate (3 mmol).The reaction mixture was refluxed for 6 h and the progress of the reaction was monitored by TLC.On completion the precipitates formed were cooled and poured into water, filtered and dried.Recrystallized from acetone to afford the purified product which were Azo-Thiohydantoin derivatives (7) as given in Scheme 1.

In-vitro studies Alkaline phosphatase inhibition assay
The measurement of intestinal alkaline phosphatase (IAP) activity was conducted using a spectrophotometric method, as described earlier [31,32].The assay utilized a reaction mixture consisting of 50 mM Tris-HCl buffer, which included 5 mM MgCl 2 , 0.1 mM ZnCl 2 , and had a pH of 9.5.Additionally, 10 µL of the compound was added to the mixture.The mixture was subjected to preincubation with 5 µL of IAP (0.025 U/mL) for a duration of 10 min.The reaction was started by introducing 10 µL of the substrate, which was 0.5 mM p-NPP (para nitrophenyl phosphate disodium salt).Following the addition, the assay underwent an incubation period of 30 min at a temperature of 37 °C.The change in absorbance resulting from the release of p-nitrophenolate was monitored at a wavelength of 405 nm using a Thermo Scientific™ Multiskan GO™ microplate reader.The experiment was repeated three times, with each repetition being done in triplicate.The reference inhibitor for CIAP was L-Phenyl alanine.The calculation of inhibition activity and IC 50 values were carried out as discussed earlier [33].

Computational studies Density functional theory
In this study, Density Functional Theory (DFT) was employed using the B3LYP method to accurately determine the electronic density of the molecules.The Gaussian 9w software was utilized to perform comprehensive geometry optimization of the compounds [34].Additionally, frequency calculations including Frontier Molecular Orbital (FMO) analysis, and global reactivity descriptors were also using the 6-31G basis set [34,35].This basis set consists of primitive and Gaussian type orbitals, which are advantageous for accurate FMO analysis.For visualization of the generated files, GaussView 6 software was employed [36].

Physicochemical and pharmacokinetic characteristics
The physicochemical and pharmacokinetic attributes of the synthesized derivatives were ascertained through SwissADME (http://www.swissadme.ch/).This computational tool facilitated a comprehensive assessment of several pharmacokinetic parameters, notably absorption, distribution, metabolism, and excretion (ADME) [37].SwissADME utilize an array of sophisticated algorithms and predictive models to prognosticate both drug-likeness and potential toxicity [38].The algorithms undergirding the bioavailability radar chart are underpinned by advanced machine learning and statistical Fig. 1 Synthesis of Azo-Thiohydantoin derivatives methodologies, calibrated against vast molecular datasets with delineated properties.Within the confines of the BOILED-Egg model, salient ADME properties, such as blood-brain barrier (BBB) permeation, passive human gastrointestinal absorption (HIA), and designation as either substrate or non-substrate for permeability glycoprotein, were distinctly identified [39].

Molecular Docking Analysis
A molecular docking analysis was conducted to examine the binding interaction between the synthesized derivatives and the targeted protein.The research was conducted by using the homology modeling in which the tissue non-specific ALP was used as the template and model was prepared using MODELLAR as described in our earlier publication.The structure of the model was validated via Ramachandran plot [40].Prior to conducting docking analysis, the protein structure of IAP was prepared using MGL tools This involved the removal of heteroatoms and water molecules, followed by the addition of polar hydrogen atoms and Kollman charges.Following the rendering of the protein structure, corrections were made to account for any missing residues [41].The 3D structure of the synthesized derivatives was determined using ChemDraw 3D [42].Energy minimization was performed to obtain the most stable arrangement of atoms.The dimensions of grid box for active pocket were fixed as follows: (x; 28.093476, y; 19.735517; z; 24.598425) for IAP.Analysis using docking was conducted on the target protein utilizing AutoDock 4.2, with the default genetic algorithm serving as the scoring function.The RMSD value along with re-docking of cocrystal ligand into active pocket was used to validate the docking process.
The newly synthesized derivatives were then docked within the active pocket of target protein, and total 100 different configurations were observed.The most stable configuration with least binding energy was retrieved and analyzed in 2D and 3D forms to elaborate the interactions formed between synthesized compounds and targeted protein.The outcomes of this study are to pave off the way for the designing of novel compounds possessing better binding affinities to the targeted protein.

Molecular dynamic simulations
Molecular Dynamic (MD) simulations provide a computational approach for studying the dynamic behavior of molecular interactions within docked complexes [43].The Nanoscale Molecular Dynamics (NAMD) software version 1.9.3 was utilized in this investigation to make predictions about the dynamic interactions occurring within a protein-ligand complex.The complex with the highest binding energy, as determined through molecular docking results, was selected for the MD simulations.
The topology files necessary for the analysis were generated using the CHARM-GUI software [44].The study employed the CHARMM36 force field parameters to accurately represent the ligand and protein molecules in a system with periodic boundary conditions.In order to replicate the conditions found in the body, the TIP3P water system was chosen as the solvent.Additionally, counter ions (NaCl, 0.15 M) were added to balance the charges in the simulation box.The system underwent equilibration in both the NVT and NPT ensembles using the Berendson thermostat and barostat, respectively.This process occurred over a timestep of 2 fs and at a temperature of 300 K before the main simulation was started.MD simulations were conducted for a duration of 100 nanoseconds.To assess the stability and flexibility of the complex during the simulation, we utilized the Visual Molecular Dynamics (VMD) software version 1.9.3.The analysis involved computing the Root Mean Square Deviation (RMSD) and the Root Mean Square Fluctuation (RMSF) using the Ewald summation method [45].,782 (C-Cl str.arom.), 1 H-NMR (300 MHZ, DMSO-d6); δ (ppm) 12.36(s, 1H, OH), 11.631(s, 1H, NH), 8.72 (1H, s), 8.21 (1H, dd, J = 1.9, 0.5 Hz), 7.94 (1H, dd, J = 8.3, 1.9 Hz), 7.41 (1H, dd, J = 7.9, 1.2 Hz), 7.16 (1H, dd, J = 8.1, 7.9 Hz), 6.85 (dd, J = 8.3, 0.5 Hz), 6.83-6.89(

Alkaline phosphatase inhibition assay
The derivative 7e, and 7f showed comparable enzyme inhibition activity i.e. 0.308 ± 0.065, and 0.836 ± 0.097 µM, respectively.The structure activity relationship of these derivatives displayed the combine effect of both mesomeric and inductive effect of various substituted groups i.e., carbonyl and Nitro group on parent ring structure.Carbonyl and nitro group functional groups are strong electron withdrawing groups due to greater electronegativity and impart specific physicochemical properties to synthesized derivatives leading to enhanced anti-enzymatic potential.
The 7e derivatives exhibited greater inhibitory activity as compared to the standard inhibitor i.e., L-Phenyl alanine IC 50 = 80.2 ± 1.1 µM.The substitutional effect on parent molecule was also confirmed in case of 7a, 7b, 7c, and 7d derivative with further reduction of enzyme inhibitory activity.In case of 7a compound the phenyl ring without any other substitution showed weak mesomeric effect that resulted in least inhibitory effect i.e. 4.847 ± 0.985.IC 50 values of synthesized derivatives 7a-7f are given in Table 1.

Density functional theory calculations
The study focused on conducting a thorough analysis of the energies of the Highest Occupied Molecular Orbital (HOMO) and Lowest Unoccupied Molecular Orbital (LUMO) using Density Functional Theory (DFT) calculations.The DFT calculations utilized the B3LYP method along with widely recognized 631-G basis set, renowned for its ability to accurately depict molecular properties.Table 2 provides the HOMO-LUMO energies, which are crucial for comprehending the electronic structure and reactivity of the compounds being studied.Table 2 presents a concise overview of the significant findings derived from our study, which specifically examined the electronic characteristics of the compounds being analyzed.
The Table 2 contains the optimization energies, as well as the energies of the Highest Occupied Molecular Orbital (HOMO) and Lowest Unoccupied Molecular Orbital (LUMO).Moreover, it provides the computed difference between the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) energies, which is a vital factor in comprehending the electronic characteristics and reactivity of the molecules, especially in the gaseous state.The optimized structures and highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) orbitals are depicted in Fig. 2.
A molecule with a larger HOMO-LUMO energy gap is considered to be "hard." This means that it possesses higher kinetic stability and is less prone to undergo chemical reactions.The larger gap indicates that there is a greater energy barrier between the occupied and unoccupied orbitals, which hinders electron transfer and reactivity.Additionally, a "hard" molecule tends to have lower polarization, as the electron density is more evenly distributed.On the other hand, a molecule with a smaller HOMO-LUMO energy gap is referred to as "soft." Such molecules exhibit lower stability and are more reactive.The smaller gap indicates a lower energy barrier, allowing for easier electron transfer and facilitating chemical reactions.These molecules tend to be more polarized, with uneven distribution of electron density.
By analyzing the HOMO-LUMO energy band gap, we can gain insights into the stability, reactivity, and polarization tendencies of a molecule, which are important factors in understanding its chemical behavior and interactions with other substances.

Global chemical reactivity descriptors
By using HOMO LUMO energy values, we evaluated the following parameters by using their respective formulas: Hardness: η = 1/2(ELUMO -EHOMO); Hardness essentially measures the resistance of a chemical system to change.In simpler terms, a molecule's 'hardness' indicates its reactivity and stability.A larger gap between HOMO and LUMO energies means greater stability and lower reactivity.Softness: S = 1/2η; Softness indicates how easily a molecule can donate or accept electrons, essentially a measure of its flexibility or adaptability in reactions.Electronegativity: χ = -1/2(ELUMO + EHOMO); Electronegativity represents the tendency of a molecule to attract electrons.Chemical potential: µ = -χ; This represents the potential energy change of a system when an additional electron is added.Electrophilicity index: ω = µ/2η gauges a molecule's likelihood to accept electrons.Higher values suggest the molecule is a good electron acceptor, making it more electrophilic.Table 3 provides an overview of the Global Reactivity Descriptors for compounds 7a-7f.These descriptors offer valuable information about the reactivity and stability of each compound.In particular, compound 7f and 7e stands out with an Electronegativity value of -0.170 eV and − 0.161 eV respectively.This value signifies its strong electronattracting capability, indicating that it has a high power to draw electrons towards itself.Consequently, compound 7f is classified as the "softest" compound among the studied compounds.This means that it exhibits lower stability, higher reactivity, and a higher tendency to undergo chemical reactions.Furthermore, compound 7f also demonstrates a high electrophilicity index, which measures the molecule's ability to accept electrons through its chemical potential.This index indicates that compound 7f has a pronounced capability to accept electrons, making it an attractive candidate for electron-transfer reactions.
In summary, the analysis of Global Reactivity Descriptors in Table 3 reveals that compound 7f exhibits the characteristics of a "soft" molecule, with high electrophilicity and a strong ability to attract electrons.These findings provide valuable insights into the reactivity and chemical behavior of compound 7f compared to the other compounds studied.

Physicochemical and pharmacokinetic characteristics
The synthesized derivatives (7a-7f) underwent evaluation for their physicochemical properties using the Swis-sADME tool.Among these, the derivatives 7a, 7c, and 7d demonstrated high oral bioavailability.The comprehensive physicochemical and pharmacokinetic profiles of these derivatives are described in Fig. 3 and supplementary file (Table S1).
In terms of solubility, all derivatives displayed effective dissolution.However, 7a, 7b, 7e, and 7f emerged as the most soluble.The compounds 7a,7c and 7d displayed a more pronounced drug likeness characteristic in comparison to all synthesized compounds.Derivatives 7a, 7c, and 7d exhibited high passive human gastrointestinal absorption, while 7b, 7e, and 7f showed low levels.
The graphical representation, illustrated by the boiled egg diagram as shown in Fig. 4, aids in the visualization of drug likeness properties, especially in relation to the blood-brain barrier (BBB) permeation.Within this diagram the yolk designated as compounds with a potential to permeate the BBB.Notably, none of the derivatives permeated the blood-brain barrier (BBB).While the white region denoted as compounds unlikely to cross the Fig. 3 Structures of synthesized derivatives and radar charts describing their physicochemical and pharmacokinetic properties BBB and are not P-glycoprotein (P-gp) substrates.The analysis showed that the derivatives 7b, 7e, and 7f were identified as P-gp substrates.The external region of the egg identified as the compounds anticipated to be P-gp substrates, implying they might exhibit reduced absorption and distribution due to potential active efflux from cells.
In the boiled egg diagram, the x-axis conventionally illustrates the topological polar surface area (TPSA), while the y-axis depicts the Wildman-Crippen LogP (WLOGP), serving as a metric for lipophilicity based on molecular structure.While the boiled egg diagram offers a streamlined visualization of the intricate factors influencing ADME, it remains invaluable for preliminary evaluations in drug discovery endeavors.Additionally, SwissTarget Prediction tool was utilized to predict the potential protein targets ranked by the probability score as shown in Fig. 5.The higher the score, the more likely the molecule binds to that protein target.However, the results were further validated with experimental data for a conclusive understanding.

Molecular docking analysis
Molecular docking is an important computational tool for determination of interactions between ligand and protein at molecular level.In current study we were aimed to predict the binding scores of synthesized derivatives with homology model of targeted protein, namely intestinal alkaline phosphatase (IAP).Homology models was retrieved from our previously reported study [1].
In recent docking experiments, the interaction of compounds 7f and 7e with the active site of IAP demonstrated notable binding affinities, with binding energies of -7.8 kcal/mol and − 7.9 kcal/mol respectively as given in Table 4. Upon in-silico evaluation, it was noted that compound 7f possessed the most favorable binding energy, whereas the 7e derivative showcased the second highest.Despite the subtle structural variations between the derivatives 7e and 7f, there is a pronounced disparity in their binding energies, which is further elaborated by in-vitro analyses.
The two-dimensional and three-dimensional binding mode of 7f derivative with IAP is shown in Fig. 6.Analysis revealed that the amino acid residues SER92, HIS432, and HIS317 engage in conventional hydrogen bonding with the 7f derivative.Concurrently, the THR431 residue establishes a carbon-hydrogen bond with the 7f.It's imperative to highlight that the metal ions, specifically Mg532 and Zn530, play a pivotal role in this binding interaction.While Mg532 functions as a metal acceptor, Zn530 contributes through a pi-cation interaction.This intricate binding interactions attributed the potential significance and specificity of these interactions in the context of IAP's active site.
In Fig. 7, a detailed representation of both the twodimensional and three-dimensional binding modes of the 7e derivative with IAP is presented.Upon rigorous   analysis, it becomes evident that certain amino acid residues, namely ARG314, GLU321, TYR276, and GLN108, engage in conventional hydrogen bonding interactions with compound 7e.This underlines the pivotal roles these specific residues might play in facilitating the binding and stabilization of the 7e derivative within IAP's active site.Moreover, other residues such as HIS432, HIS317, HIS320, and ASP316, along with the metal ion Zn530, are observed to participate in pi-anion interactions.These interactions are characterized by the overlap between the pi-system of an aromatic residue and an anionic moiety, which results in a stabilizing force.Such interactions emphasize the intricate nature of the binding, where multiple residues and forces cooperate to accommodate the ligand, in this case, the 7e derivative, within the protein's binding pocket.2D and 3D binding modes of 7a-7d are given in supplementary file (Figure S4-S7).These results were further validated through invitro analysis as described next.
Figure 8 provides a comprehensive depiction of the binding modes between L-Phenyl alanine and IAP, showcasing both two-dimensional and three-dimensional perspectives.Through thorough examination, it is apparent that specific amino acid residues-namely GLY434, ASP436, THR82 and TYR83-establish the conventional hydrogen bond with the L-phenylalnine.Additionally, TYR83 participate in pi-pi stacked interaction with the standard inhibitor compound, L-Phenyl alanine.

Molecular dynamics simulations
Molecular Dynamics (MD) simulations were executed to investigate the conformational dynamics of both the apo Fig. 7 The predicted 2D and 3D binding pose of compound 7e within the active pocket of IAP Fig. 6 The predicted 2D and 3D binding pose of compound 7f within the active pocket of IAP protein and its complex with the ligand.The stability of these structures was quantitatively examined using Root Mean Square Deviation (RMSD) and Root Mean Square Fluctuation (RMSF) over the simulation duration.
From the RMSD plot (Fig. 9 (a)), it is discernible that the protein-ligand complex offers superior stability in comparison to the isolated apo protein and ligand.Specifically, Fig. 9(a) elucidates that the complex retains consistent conformational steadiness throughout the simulation, as evinced by an RMSD of 1.88Å.Intriguingly, the apo protein and ligand (7e) independently showcased robust stability, evidenced by RMSD values of 1.70Å and 1.60Å, respectively.
RMSF calculations, depicted in Fig. 9(b), reveal an average fluctuation of 12.14Å, emphasizing the dynamic flexibility of the amino acid residues, a feature pivotal for potent protein-ligand interactions.Notably, the ligand (7e) remained anchored within the protein's active site during the simulation.The pronounced peaks in Fig. 9(b) signify residues undergoing significant vibrational oscillations, while the troughs pinpoint residues fostering robust hydrophobic and hydrophilic interactions with the ligand.
The Solvent Accessible Surface Area (SASA) plot (Fig. 9(c)) underscores segments of the protein exhibiting heightened solvent exposure, suggesting these regions possess inherent flexibility and are less compact.Such domains potentially partake in the dynamic functional activity and molecular interplay of the protein.Conversely, the troughs in the SASA plot allude to internally situated residues, integral to the molecule's structural stability and folding.
Radius of Gyration (Rg) analysis, represented in Fig. 9(d), provides a lens into the compactness of the protein's conformation.Periods with lower Rg values correlate with instances where the protein assumes a more condensed, globular structure, indicative of its stabilized states.By juxtaposing insights from both SASA and Rg plots, one can derive a comprehensive understanding of the evolving conformation and resilience of the apo protein across various conditions.Together, these analytical tools offer invaluable perspectives into the nuanced structure-function dynamics of biomolecules, facilitating advanced studies in protein stability and mechanics.
Electrostatic interactions between ligand 7e and the IAP protein are intricately shaped by the amino acid residues responsible for hydrogen bond formation.An in-depth analysis of Table 5; Fig. 10 highlights the pivotal roles of specific IAP residues in fostering a strong hydrogen bond network with the ligand.
Notably, residues SER429, GLU321, HSD432, and SER428 exhibit a pronounced propensity to engage in hydrogen bonding with the ligand's side atoms.This examination into hydrogen bond dynamics reveals not only the system's stability nuances but also provides insights into the complex dance between the ligand and protein, deepening our grasp on the underlying forces dictating their synergy.

Conclusion
In conclusion, the study comprehensively examined the electronic properties and reactivity of the synthesized compounds through Density Functional Theory (DFT) calculations.The HOMO-LUMO energy band gap analysis demonstrated its pivotal role in assessing molecule stability and reactivity.SwissADME evaluations highlighted compounds 7a, 7c, and 7d for their favorable physicochemical properties, including solubility and This comprehensive integration of computational and experimental approaches sheds light on the potential therapeutic applications of the synthesized compounds.The study underscores the importance of understanding electronic properties, binding mechanisms, and structure-activity relationships in drug discovery and design.By providing a detailed investigation of these aspects, this research opens avenues for the development of novel pharmacologically active compounds with diverse applications.

Fig. 5
Fig. 5 Predicted classes of proteins using Swiss ADME's target prediction tool

Fig. 8
Fig.8The predicted 2D and 3D binding pose of compound L-Phenyl alanine within the active pocket of IAP

Table 1
Intestinal Alkaline phosphatase inhibitory activity of 7a-7f derivativesValues are presented as Mean ± SEM (Standard error of the mean)

Table 2
Optimization energies, HOMO and LUMO energies, and their gap estimated in gas

Table 3
An overview of the Global Reactivity Descriptors for synthesized compounds (7a-7f)

Table 4
Molecular docking scores of compounds with Intestinal Alkaline Phosphatase (IAP)

Table 5
Hydrogen Bond Donor and Acceptor Groups of 7e and IAP protein