[1, 8]-Naphthyridine derivatives as dual inhibitor of alkaline phosphatase and carbonic anhydrase

[1,8]-Naphthyridine derivatives have been reported to possess important biological activities and may serve as attractive pharmacophores in the drug discovery process. [1,8]-Naphthyridine derivatives (1a-1l) were evaluated for inhibitory potential for isozymes of carbonic anhydrase (CA) and alkaline phosphatase (ALP). CAs have been reported to carry out reversible hydration of CO2 into HCO3−, secretion of electrolytes, acid–base regulation, bone resorption, calcification, and biosynthetic reactions. Whereas ALPs hydrolyze monophosphate esters with the release of inorganic phosphate and play an important role in bone mineralization. Both enzymes have been found to be over-expressed and raised functional activities in patients suffering from rheumatoid arthritis. The discovery of dual inhibitors of these enzymes may provide a synergistic effect to cure bone disorders such as rheumatoid arthritis and ankylosing spondylitis. Among the test compounds, the most potent inhibitors for CA-II, CA-IX, and CA-XII were 1e, 1g, and 1a with IC50 values of 0.44 ± 0.19, 0.11 ± 0.03 and 0.32 ± 0.07 µM, respectively. [1,8]-Naphthyridine derivatives (1a–1l) were approximately 4 folds more potent than standard CA inhibitor acetazolamide. While in the case of ALPs, the most potent compounds for b-TNAP and c-IAP were 1b and 1e with IC50 values of 0.122 ± 0.06 and 0.107 ± 0.02 µM, respectively. Thus, synthesized derivatives proved to be 100 to 800 times more potent as compared to standard inhibitors of b-TNAP and c-IAP (Levamisole and L-phenyl alanine, respectively). In addition, selectivity and dual inhibition of [1,8]-Naphthyridine derivatives confer precedence over known inhibitors. Molecular docking and molecular simulation studies were also conducted in the present studies to define the type of interactions between potential inhibitors and enzyme active sites. Supplementary Information The online version contains supplementary material available at 10.1186/s13065-023-01052-8.


Introduction
Carbonic anhydrases (CAs; EC: 4.2.1.1)are a family of pervasive Zn 2+ metalloenzymes found in almost all living organisms, from prokaryotic cells to more complex eukaryotes.CAs catalyze the reaction of reversible hydration of CO 2 via two-step reactions to yield HCO 3 − and a proton [1].To date, 16 distinct mammalian CA isozymes having specific activity, physiological roles, kinetic properties, tissue specificity, sub-cellular localization, and susceptibility to various inhibitors have been identified.Carbonic anhydrase I, II, III, VII, and XIII are cytosolic isozymes; CA-IV, IX, XII, and XIV are trans-membrane; CA-VA, and VB are present in mitochondrial.CA-VI is a secretory protein commonly found in saliva and breast milk [2].Cell surface CA-IX and XII isoforms have high expression in cancerous cells and are responsible for cell adhesion, proliferation, and mutagenesis.CA-IX has minimal expression in normal tissues, whereas it was found to be over-expressed on the cell surface of solid tumors.High expression of CA-IX and CA-XII in hypoxic tumors, coupled with the poor prognosis and the aggressive phenotype, make them candidate targets for cancer therapy [3].
Alkaline phosphatases (APs, EC: 3.1.3.1),belong to the family of ecto-nucleotidases, are metalloenzymes having two Zn 2+ and one Mg 2+ ions present in active sites for optimal activities of enzymes [4].Alkaline phosphatases are classified into two categories; tissue-specific APs and tissue non-specific alkaline phosphatase (TNAP).Tissue-specific APs are further subdivided into placental (PLAP), intestinal (IAP), and germ cell (GCAP) [5].APs catalyze nucleotides and release inorganic phosphate (Pi).APs can hydrolyze a range of substrates, i.e., glucosephosphates, phosphatides, and inorganic phosphates [6].APs are assumed to be imperative in differentiating osteoblasts adipocytes and produce adenosine by sequential de-phosphorylation [6][7][8].TNAP is predominantly found in bone and other mineralizing tissues as it regulates mineralization and bone formation [9].Decreased levels of TNAP lead to skeletal hypo-mineralization [10].Overexpression of TNAP causes mineral deposition in soft tissues and also causes osteoarthritis [11].IAP regulates the accumulation of lipids and is involved in adipogenesis.Moreover, it protects against bacterial lipopolysaccharides (LPS) [12].In the duodenum, the gastric proton is neutralized with secreted bicarbonate to generate carbon dioxide and water, whereas IAP and cytosolic CAs are responsible for maintaining homeostasis [13].Selective inhibition of IAP inhibits bacterial growth [14].
Alkaline phosphatase and carbonic anhydrases are involved in mineral deposition in bone which may lead to rheumatoid arthritis [15].It has been reported that CA-IX and CA-XII over-express in the inflamed synovium of patients having juvenile idiopathic arthritis [16,17].The raised activity of serum alkaline phosphatase has also been observed in patients affected with rheumatoid arthritis [18,19].Inhibition of both enzymes may have a synergistic effect in treating this pathological condition.Designing and synthesis of compounds that inhibit both enzymes can be very fruitful for treating arthritis [20].Figure 1 shows that CA and ALP are involved in bone mineralization.
Naphthyridine derivatives are reported to exhibit antimicrobial, antiviral, anti-inflammatory, immune-modulatory, and anti-cancerous activities [22,23].Nalidixic acid was reported as the first antimicrobial agent in 1962.According to the literature [1,8]-Naphthyridine derivatives are potential antihypertensive molecules.Some heterocyclic condensed [1,8]-Naphthyridine derivatives were also reported as antimicrobial agents [24,25].In short, the spectrum of activities of these derivatives includes antiviral, anticancer, and anti-inflammatory activities [25,26].These derivatives are also reported to show immune-modulatory and anti-cancerous activity [27].In the current studies, the role of [1,8]-Naphthyridine derivatives was investigated for inhibitory potential against CA-II, IX, XII, and also for b-TNAP and c-IAP.The experimental results were subjected to molecular docking and MD simulations to study the possible interactions among active site amino acids and inhibitors.
Data obtained from the inhibition of CA-II suggests that one of the most potent compound was 1e with IC 50 (µM) ± SEM value of 0.44 ± 0.19.The potency of this series was due to the substitution of R group.The most potent inhibitor of CA-II among the screened compounds, was found to be 1g that possess a methoxy substitution at meta position of benzene ring, Ig has IC 50 value of 0.10 ± 0.04 µM.In compound 1e butane was attached at pyrrolo [1,8]-Naphthyridine ring resulting in maximum inhibition.While in the case of compounds 1d and 1i the toluene was present which significantly decreased its inhibitory action with loss of activity, the substitution of anisole in 1f was also causing a decrease in inhibition as in 1b.Inhibition of fluorobenzene-substituted compounds was dependent on fluorine.In compound 1c fluorine at para-position was exhibiting more inhibition with IC 50 ± SEM value of 0.90 ± 0.5 µM, than ortho-position in compound 1j IC 50 ± SEM value of 2.14 ± 0.3 µM.There was no significant change with the substitution of benzene or tert-butylbenzene.Only three compounds 1d, 1f, and 1i have shown less than 50% inhibition due to substitution of toluene and methoxy phenol.While in the case of CA IX, it was found that the compound 1g was the most potent compound having an IC 50 (µM) ± SEM value of 0.11 ± 0.03 µM.Compound 1g exhibited its activity due to the presence of OMe group at meta position of the ring.Further substitution of this R effected the inhibition.The compound 1c has fluorobenzene with fluorine at position para position.If this fluorine is transferred from para to ortho position as in case of compound 1j, the IC 50 is elevated to 1.65 ± 0.98 µM.Substitution of anisole with fluorobenzene had not much effect on its inhibition, while substitution of thiophene or tert-butylbenzene at the same position decreased inhibition drastically.Interestingly substituting pyrrolo [1,8]-Naphthyridine ring with butane in 1e also showed effective inhibition with an IC 50 (µM) ± SEM value of 1.61 ± 0.26 µM.Only two compounds 1b and 1k have shown inhibition of less than 50%.
All compounds were potent inhibitors of CA-XII.The most potent was 1a as in the case of CA-XII with IC 50 (µM) ± SEM value of 0.32 ± 0.07 µM.Compound 1i is inactive against CA-II, while it shows significant activity against the other two isozymes.Similarly, compound 1d is inactive against CAII while showing marked activity against CA-IX and CA-XII (Table 2).
Compound 1b was most potent in the case of b-TNAP with IC 50 ± SEM of 0.12 ± 0.06 µM, where Bu is substituted at para position of the ring.1b was found to be 160 times more potent as compared to the standard inhibitor of TNAP, Levamisole (IC 50 = 19.21± 0.001 µM).Compound 1f exhibited a slight increase in IC 50 because of substitution of methoxy phenyl at same position.IC 50 was increased to 0.14 ± 0.12 µM.Compounds 1l, 1e and 1k showed relatively close IC 50 values of 0.24 ± 0.06, 0.32 ± 0.11, and 0.36 ± 0.07 µM, respectively.IC 50 value of compound 1k is due to the substitution of the furan ring.In compound 1e free pentane is substituted at pyrrolo [1,8]-Naphthyridine ring which is responsible for its activity.In compound 1l activity is attributed to cyclohexane substituted at pyrrolo [1,8]-Naphthyridine ring.Substitution with fluorine as in the case of 1j, toluene in 1i and anisole in 1h resulted in loss of activity.
Compound 1e expressed the highest potency against c-IAP as butane is fused with a naphthyridine ring with an IC 50 xvalue of 0.10 ± 0.02 µM and found to be 800 times more potent as compared to L-phenyl alanine, a standard inhibitor for IAP.Any other substitution resulted in the elevation of IC 50 value as shown in compound 1c substitution of fluorine at benzene ring increased IC 50 value to 5.80 ± 1.90 µM.Substitution of anisole resulted in a loss of activity in 1h.In compound 1h anisole is substituted which resulted in a complete loss of activity.Activity in compound 1d is due to the substitution of toluene at pyrrolo [1,8]-Naphthyridine ring, having an IC 50 of 0.30 ± 0.20 µM.Compound 1g exhibited an IC 50 value of 0.35 ± 0.21 µM due to substitution of methoxy phenyl at benzene ring attached to pyrrolo [1,8]-Naphthyridine nucleus.

Density functional theory (DFT) calculations
For the last 30 years, among most popular techniques for doing ab initio calculations on the structure of atoms, molecules, crystals, surfaces, and their interactions is density functional theory (DFT).These calculations are used to understand the electronic behavior of a compound which helps in predicting the chemical reactivity.All the DFT calculations of this study [1,8]-Naphthyridine derivatives (1a-1l) were performed using B3LYP functional theory and 6-31G basis set.The 6-31G is a medium basis set most commonly used for compounds with atoms up to Argon [41].Using Gaussian software, a number of parameters, including the dipole moment, molecular polarizability, and optimization energy, were determined for all the derivatives [42].The optimization energy of a compound is the amount of energy required to get the most stable configuration, i.e. lower the optimization energy of a compound higher its stability, which may be represented by its optimized structures.The optimized structures of all the derivatives are given in Fig. 2.

Table 2
In vitro inhibition assay for compounds 1a-l against b-TNAP and c-IAP.Levamisole and L-phenyl alanine were used as standard

HOMO and LUMO orbitals
The energy of any compound for its highest occupied and lowest unoccupied molecular orbitals, i.e., HOMO and LUMO, can be used to determine the range of chemical properties related to its electronic distribution [43,44].Any molecule with a high HOMO value appears to be an excellent electron donor; electron acceptors have higher LUMO values.Determining local reactivity at various sites in a compound was made more accessible with the help of the overall FMO calculation.A compound's HOMO/LUMO energy difference significantly indicates its reactivity profile.Any compound with a more significant energy gap tends to be less reactive, as indicated by the high hardness value of the compound.According to the DFT results, all 12 compounds were highly reactive, evident from their lowest average energy gap, i.e., ΔEgap = − 0.136, which is almost common for all derivatives.The ΔEgap results for all the other compounds are positive and are given in the Table 3, along with other parameters.The dipole moment of any substance is related to the charge segregation in a molecule and helps to find its active sites.The higher value of dipole moment speaks for its more polar character with maximum reactive points.According to the DFT results, the derivative 1g has the highest dipole moment value, i.e., 3.07D, which is consistent with both in vitro and in silico studies.The lowest dipole moment value of 1j is also relevant to other studies.In the same high polarizability value of compound 1g augments the DFT findings.The structure of the HOMO and LUMO orbitals for all the compounds are shown in Fig. 3 below, while the polarizability, dipole moment, energy, and relative energy gap are calculated and presented in the Table 3.

Chemical descriptors
Various parameters used to assess a compound's reactivity include chemical hardness and softness, electrophilicity index, electronegativity, and electronic chemical potential [45].The chemical formula (ELUMO-EHOMO)/2) can be used to determine the chemical hardness of a chemical system based on its stability and reactivity.The high value of hardness for a compound is related to its least reactivity and high stability.The ability of an atom in a molecule to draw protons towards itself is known as electronegativity, which is quantified by the formula X = -(EHOMO + ELUMO)/2).The electrophilicity index measures a molecule's capacity for accepting electrons by utilizing its electronic chemical potential and chemical hardness.Table 4 contains the values of several chemical descriptors for all substances that were determined to be similar.

Molecular docking studies against carbonic anhydrase II (3K34)
The in vitro analysis revealed that 1g was found to be the most potent inhibitor of CA II.The nitrogen heteroatom of pyrrolo-naphthyridine ring of 1g formed a metallic linkage with Zn 2+ metal, and the same pyrrolo-naphthyridine ring was involved in hydrogen bonding with His94, His119, and Thr199, while forming π-π linkages with Val143 and Leu198.The compound 1g also exhibited van der Waals interactions with Gln92, His94, Val121, Phe131, Val143, Leu198, Thr200, and Trp209 residues (Fig. 5).

Molecular docking studies against carbonic anhydrase-XII (5MSA)
The enzyme inhibition assay revealed 1a and 1i to be the most potent inhibitors of the carbonic anhydrase XII enzyme.The nitrogen heteroatom of pyrrolonaphthyridine ring of 1a formed a metallic linkage with Zn 2+ metal and forming π-π linkages with His91, while π-alkyl bonding with Val141 and Leu197 residues.Pyrrolo-naphthyridine ring established a hydrogen bonding with Thr198 residue.The compound 1a also exhibited van der Waals interactions with His91, His93, Leu197, Thr199, Pro201 and Trp208 (Fig. 7).
The nitrogen heteroatom of pyrrolo-naphthyridine ring of 1j formed a metallic linkage with Zn 2+ metal and forming π-π linkages with His91, while π-alkyl bonding with Leu197 residue.Pyrrolo-naphthyridine ring was also involved in the formation of hydrogen bonding with Thr198 residue.Fluorophenyl moiety of the compound 1j was involved in the formation of hydrogen bonding with Gln89.The compound 1j also expressed van der Waals interactions with His91, Val141, Leu197, and Thr199 residues (Fig. 9).

Molecular docking studies against human intestinal alkaline phosphatase (h-IAP)
Initially, the docking methodology was validated by docking the L-phenylalanine as standard inhibitor of IAP into the active pocket of h-IAP.
It was observed during the molecular docking that the test compounds imitated the interactions of L-phenylalanine.In vitro analysis found 1a, 1d, 1e, 1g, 1j, and 1k as effective inhibitors of the h-IAP protein.The nitrogen heteroatom of pyrrolo-naphthyridine ring of 1a formed a  conventional hydrogen bonding with His317 residue, and the pyrrolo-naphthyridine ring itself was involved in the formation π-cationic bonding with Arg166 residue of the protein.The compound 1a also exhibited van der Waals interactions with Asp42, Ser92, Arg166, Asp316, His317, His437, and the two Zn 2+ metals found in the active pocket of the enzyme (Fig. 10).
he nitrogen heteroatom of pyrrolo-naphthyridine ring of 1d formed hydrogen bonding with His317 and His153 residues, and the pyrrolo-naphthyridine ring itself was involved in the formation π-cationic bonding with Arg166, Asp316, and Zn 2+ metal of the active pocket.His432 formed π-π interaction with the phenyl ring of compound 1d.The compound 1d also exhibited van der Waals interactions with Asp42, Ser92, Phe107, Arg166, His317, and His432 protein residues (Fig. 11).

Molecular docking studies of test compounds against human tissue non-specific alkaline phosphatase (h-TNAP)
The binding residues for test compounds in h-TNAP were validated with its known ligand levamisole.It was observed during the molecular docking that the test In vitro analysis found 1a, 1b, 1e, 1f, 1k, and 1l as compelling inhibitors of the h-TNAP protein.
The pyrrolo-naphthyridine ring of 1f was involved in forming π-π linkages with His321 and His434 residues while π-cationic linkages with Asp277 and Arg318 residues.The methoxy moiety of methoxyphenyl ring formed a metallic linkage with Zn 2+ metal, and the ring itself was involved in forming π-π linkage with His324 and π-cationic linkage with Asp320 residues.The compound 1f also exhibited van der Waals interactions with His154, His321, His324, and Glu325 residues (Fig. 13).

HYDE assessment of selective and potent test compounds against CA-II, CA-IX, CA-XII, h-IAP and h-TNAP
The HYDE affinity assessment tool of LeadIT was put into effect for the top 300 ranking docked conformations within the active sites of crystal structures of CA II, CA IX, CA XII, h-IAP and h-TANP.This affinityserved in the selection of appropriate binding mode of the most potent compound.The FlexX score of the selective compounds and their binding free energy ΔG are given in Table 5.

Molecular dynamic simulations
After carrying out the docking studies, the binding results of the most potent derivatives were validated by MD simulation to augment the research findings further.The protein-ligand complexes of both proteins i.e.Carbonic Anhydrase (II and IX) and Alkaline Phosphatase with the most potent [1,8]-Naphthyridine derivative for both groups of proteins i.e.1g was simulated for 100 ns in aqueous conditions by NAMD software.The Root Mean Square Deviation (RMSD) value and Root Mean Square Fluctuation (RMSF) value were used to analyze the conformational stability of proteins and their complexes.The RMSD graph details any structural variation caused by the ligandprotein interaction under a simulated environment.The results of MD simulation were assessed from the RMSD plots of both proteins as well as their derivatives.These RMSD plots were compared to evaluate the stability of the protein complexes.Further analysis of the MD simulation was done by plotting RMSF graphs, which is used to detect any changes in the target protein's C and N terminal lobe amino acid residues during the trajectory (see Additional file 1).
All the results of MD simulations were evaluated using RMSD graphs.The RMSD value less than 2 is considered good and indicates the stability of the ligand-protein complex.

Carbonic anhydrase II with 1g
According to simulation results, the docking complex of carbonic anhydrase II with 1g was found stable throughout the simulation time.The average RMSD value of protein is less than 2, which indicates the higher stability of the protein backbone, while the RMSD value of the protein-ligand complex was also smooth and within the acceptable range indicating its stability with an average RMSD value of 2 which is slightly higher than protein RMSD.In the same way, RMSF of the protein chain was also calculated, showing the slight fluctuation of Amino Acid Residues returning to its initial position with an average RMSF value of 1.5 Å. Results are elaborated in the following Fig. 14.

Carbonic anhydrase IX with 1g
In the same way docking complex of Carbonic anhydrase IX was also found stable throughout the simulation trajectory.The average RMSD value for the protein-ligand complex is slightly higher than 2 but is comparable to the protein RMSD graph indicating the stronger and stable association of the protein-ligand complex.The RMSF graph was also found stable with sight fluctuations which can be attributed to the flexibility of protein.The RMSD and RMSF graphs for carbonic anhydrase IX with 1g are shown in the figure below.Results of molecular simulation are presented in Fig. 15.

Alkaline phosphatase with 1g
The MD simulation results of protein alkaline phosphatase with 1g also indicate stable interaction with an average RMSD value of 2.2 and are comparable to the protein RMSD value.The average RMSF value of 1.5 indicates the stability of protein residues during the trajectory (Fig. 16).
When cells gained about 80-90% confluency, they were washed twice with PBS and transfected with 0.2 µg plasmid solution.Cells were incubated for 6 h in 5% CO 2 incubator, followed by the addition of fresh growth media.Cells were again incubated for 24 h.Hygromycin B (0.2 mg/mL) was used as a selection antibiotic for transfected cells.Cells were harvested after five days and lysed.Then desired proteins were purified via Ni-NTA resin using imidazole as affinity.SDS PAGE was carried out to confirm the molecular weight of purified protein.

Inhibition assay for carbonic anhydrase
The CA isozymes inhibition assay was performed according to the reported method after some modifications [30].
Assay buffer was prepared that contained 50 mM tris-sulfate, 0.1 mM ZnCl 2, and pH was adjusted at 7.6.Initially, 60 µL assay buffer was added in each well of 96 wells plate, 10 µL test compound (final concentration 100 µM), and followed by the addition of 10 µL working enzyme solution.
The reaction mixture was incubated for 10 min at 37 °C.Absorbance was measured as pre-read at 348 nm using a microplate reader (FLUOstar Omega, Germany).Then 20 µL of the substrate (6 mM) was added to the reaction mixture to initiate the reaction.The reaction mixture was incubated at 37 °C for 30 min, and absorbance was measured again.The activity was compared with a negative control, a reaction mixture without any inhibitor.Results were reported as the mean of triplicate experiments (± SEM).The % inhibition was calculated by following the formula.

Alkaline phosphatase inhibition assay
[1, 8]-Naphthyridine derivatives were screened for inhibitor potential towards b-TNAP and c-IAP.The assay was based on previously reported protocol [31].The assay volume was 100 µL, and the final concentration of each tested compound was 100 µM per well.

Density functional theory (DFTs)
For all the DFT calculations, Guassian09 programme [32] was used with the B3LYP functional scheme and 3-21G basis set to carry out the density functional theory (DFTs) calculations [33].The Computational calculation of electronic structure for atoms and molecules makes use of this efficient theory.The chemical structures of all the chosen derivatives were drawn in Chemdraw 12.0 to facilitate DFT experiments.ChemDraw 3D Pro was used in a similar manner for superficial energy minimization, and files were saved in the common Sybyl mol2 format.
Frontier molecular orbital (FMO), global, and local reactivity descriptors were obtained, together with the optimized geometric parameters.Fchk files were analyzed with Guass View 6 [34].

Selection of the protein structures and preparation of ligands
Crystallographic structures of CA-II, CA-IX, and CA-XII with PDB IDs 3K34, 6G9U, and 5MSA, respectively, were procured from the Protein Data Bank.The x-ray crystal structures of h-TNAP and h-IAP are currently not present in the Protein Data Bank; homology-modeled already reported structures from our group have resorted to docking.The crystal structures were prepared utilizing the Protein Preparation Wizard, which was implemented within the MOE software, employing the default settings.Hydrogen atoms were incorporated into the molecular operating environment (MOE) utilizing the protonate 3D protocol.Subsequently, the protein structures underwent minimization employing the MMFF94x force field until reaching an RMSD gradient of 0.1 kcal•mol −1 Å −1 [35].Structures of selected compounds and their stereochemistry was patched by energy minimization.h-TNAP and h-IAP models did not have any incorporated ligands, so the validation of these enzymes models was done with the positive standards used in the biological assay.Inbuilt MOE site finder was applied for the selection of possible binding sites of the target protein.Zn 2+ was kept in the center of enzyme-selected active sites.After initial validation, molecular docking of the targeted compounds was performed.

Docking of ligands in targeted proteins
The docking studies were executed by docking software LeadIT [36].Software default parameters were used to perform the molecular docking studies of the reference standards that were used in in vitro assay as well as for selected compounds.For each ligand, the most promising docked pose was selected and was further assessed through the HYDE assessment tool.The Discovery Studio Visualizer (v19.1.0.18287) was used to represent 3D interaction of selected docked poses [37].

Molecular dynamics simulation studies
The molecular dynamic simulation studies were performed using Nanoscale molecular dynamic (NAMD) software on a CUDA-accelerated GPU machine with a 16-core CPU and 64 GB RAM memory with top-ranked conformation.The visualization of results was done using VMD software [38].The MD simulation is used to ascertain the binding interactions with targeted proteins as well as stability of the protein-ligand complex under accelerated conditions.The top-ranked protein-ligand complex was selected, and topology files for the protein and the ligand were created using the CHARMM36 forcefield [39].The required NaCl charges were added to neutralize the system.To eliminate any close atomic interactions, the system was reduced to the sharpest energy gradient.The system was equilibrated in an NVT ensemble for 500,000 steps, followed by another 500,000 steps in an NPT ensemble.After that, a simulation with specific periodic boundaries was run for 100 ns [40].

Conclusion
In the summary, [1,8]-Naphthyridine derivatives (1a-1l) were found to be potent inhibitors of carbonic anhydrase as well as alkaline phosphatase isozymes.Among the test compounds, the most potent inhibitors for CA-II, CA-IX, and CA-XII were 1e, 1g, and 1a with IC 50 values of 0.44 ± 0.19, 0.11 ± 0.03 and 0.32 ± 0.07 µM, respectively.While in the case of ALPs, the most potent compounds for b-TNAP and c-IAP were 1b and 1e with IC 50 values of 0.122 ± 0.06 and 0.107 ± 0.02 µM, respectively.In silico studies, including Density Function Theory, Molecular docking, and MD simulations, confirmed strong binding interactions exist between active sites of enzymes and potential inhibitors.These dual inhibitors of both the enzymes may prove to be very efficient in the treatment of various bone disorders, especially rheumatoid arthritis.Based on this study, more potent molecules can be designed as dual inhibitors of both enzymes for further investigation in the future.

Fig. 1
Fig.1Mineral deposition by two enzymes.a Carbonic anhydrase (CA) and a alkaline phosphatase (ALP), has been reported to be done in three steps.Step 1: CA catalyzes the reaction that leads to the formation of calcium carbonate bio-seed.Step 2: Non-enzymatic reaction, that transforms carbonate-phosphate by using the orthophosphate generated as a result of ALP-mediated hydrolysis of polyphosphate.Step 3: Formation of amorphous calcium carbonate as well as hydroxyapatite crystals[15,21]

a n = 3
results were obtained from three individual experiments b Percentage of inhibition measured at 100 µM final concentration of selected compounds Compounds codes b-TNAP c-IAP IC 50 (µM) ± SEM a / % Inhibition b

Fig. 14 AFig. 15 AFig. 16 A
Fig. 14 A RMSD plots of Protein and Protein-ligand complex of Carbonic Anhydrase II.B RMSF plots of protein-ligand complex

Table 3
Optimization energies, HOMO and LUMO energies, and energy gap for compounds (

Table 5
HYDE assessment of the top-ranking poses