Gangireddy Sujeevan Reddyab,
Sharda Shuklaab,
Harshavardhan Bhuktarab,
Kazi Amirul Hossaina,
Rebecca Kristina Edwina,
Varadaraj Bhat Giliyarub,
Parimal Misraa and
Manojit Pal*a
aDr. Reddy's Institute of Life Sciences, University of Hyderabad Campus, Gachibowli, Hyderabad 500 046, India. E-mail: manojitpal@rediffmail.com
bManipal College of Pharmaceutical Sciences, Manipal Academy of Higher Education, Madhav Nagar, Manipal 576 104, Karnataka, India
First published on 21st September 2022
In view of the reported chorismate mutase (CM or MtbCM) inhibitory activities of 3-indolylmethyl substituted (pyrazolo/benzo)triazinone derivatives the structurally similar 3-(benzofuran-2-ylmethyl) substituted (pyrazolo/benzo)triazinones were designed and evaluated in silico against CM. The docking of target molecules was performed at the interface site of MtbCM (PDB: 2FP2). All the best ranked molecules participated in a strong H-bonding with the ILE67 of the B chain at the backbone position in addition to several hydrophobic/van der Waals interactions with the hydrophobic residues. Based on encouraging docking results, the one-pot synthesis of newly designed benzofuran derivatives was carried out using tandem Pd/Cu-catalyzed Sonogashira cross-coupling followed by intramolecular cyclization of 2-iodophenols with appropriate terminal alkynes. A range of novel 3-(benzofuran-2-ylmethyl) substituted (pyrazolo/benzo)triazinone derivatives were prepared in high (>80%) yields. Three molecules i.e. 3h, 3i and 3m that participated in good interaction with CM in silico showed encouraging (64–65%) inhibition at 30 μM in vitro. An SAR within this class of molecules suggested that the benzotriazinone series in general was better than the pyrazolotriazinone series. Based on molecular docking in silico, CM inhibition in vitro and computational ADME prediction the benzofuran derivatives 3i and 3m seemed to be of further medicinal interest in the context of discovery and development of new anti-tubercular agents.
Fig. 1 Reported inhibitors of CM based on (pyrazolo/benzo)triazinone derivatives A and the proposed template B. |
Moreover, CM inhibitory potential of benzofuran derivatives have been documented recently.10 Herein, we report the in silico studies, chemical synthesis and in vitro CM inhibitory potential of a series compounds based on B (Fig. 1).
Since the binding energy of our studied ligands at interface site was lower than the other possible sites of MtbCM, hence we thought that our molecules were preferably stabilized at the interface site of MtbCM (which is a homodimer). Nevertheless, the docking of our molecules at the interface site of MtbCM (PDB: 2FP2) was performed using the open-source tool AutoDock Vina.13 Additionally, the MarvinSketch14 and AutoDock tool15 were also used for conducting the docking studies and the docking score of the best pose of each ligand along with that of the reported compound8 A-1 is presented in Table 1. To check the reproducibility of docking results, the docking of each individual ligand was performed at least for five times when the maximum difference in the scores obtained was found to be within the range of ±0.2. The validation of docking protocol was carried out by re-docking the co-crystal ligand (TSA) and calculating the RMSD difference between co-crystal one and docked one in PyMOL.16 The maximum RMSD being 1.713 (<2) clearly confirmed the correctness of the docking protocol employed.
a Figure within the bracket represents AutoDock Vina score (kcal mol−1). |
---|
Based on docking scores (Table 1) the compound 3h, 3i and 3m (score ∼ −10 kcal mol−1) appeared as potential hits in the in silico studies and was comparable to the known compound A-1. Indeed, the compound 3h (Fig. 2) formed a strong H-bond through the ketone carbonyl group with the ILE67 of B chain at the backbone position. Additionally, the molecule also participated in the hydrophobic/van der Waals interactions with the hydrophobic residues (LEU65, PRO66, TYR110, PHE113, LEU65 of both A and B chain) and hydrophobic regions (polar and charged residues e.g. GLU68 of both A and B chain and GLN64, LYS117 of A chain). The compound 3i (Fig. 3) showed H-bonding interaction through its nitro group with the backbone atom of ILE67 of B chain in addition to a π–π interaction with TYR110 of B chain (which could be the reason of its slightly better docking score than compound 3h). Similar hydrophobic contacts mainly with LEU56, PRO66, PHE113, GLU68, etc. of both the chain were also found in this case. Like compound 3i, the involvement of nitro group in the H-bond interaction was noted for compound 3m (Fig. 4) in addition to its participation in the usual hydrophobic interactions.
Fig. 2 2D and 3D interaction (H-bond is shown in yellow colour) of compound 3h with the interface residues of MtbCM (PDB: 2FP2) [prepared using Maestro visualizer (Schrödinger, LLC)]. |
Fig. 3 2D and 3D interaction (H-bond and π–π interaction shown in yellow and cyan colour, respectively) of compound 3i with the interface residues of MtbCM (PDB: 2FP2). |
Fig. 4 2D and 3D interaction (H-bond and π–π interaction shown in yellow and cyan colour, respectively) of compound 3m with the interface residues of MtbCM (PDB: 2FP2). |
Nevertheless, in the light of promising overall outcome of the in silico study it was become necessary to correlate these observations with the data generated through an in vitro assay. Hence, the chemical synthesis as well as some preliminary in vitro evaluation of compounds presented in Table 1 was undertaken and the results are presented in the following sections. Notably, the in silico studies, synthesis and in vitro evaluation of 3-(benzofuran-2-ylmethyl) substituted (pyrazolo/benzo)triazinone derivatives as potential inhibitors of chorismate mutase (CM) was not known previously.
Among the several possible strategies considered for the synthesis of targeted compounds based on B (Fig. 1), we envisioned that the construction of the benzofuran ring could be a convenient as well as effective approach for serving the purpose. While numerous approaches are known for the construction of 2-substituted benzofuran ring,17 the one-pot method18–35 involving the tandem transition-metal-catalyzed Sonogashira cross-coupling36 followed by intramolecular cyclization of 2-halogenated phenols with terminal alkynes gained considerable interest due to its convenience, functional group tolerance and environmental advantages. Indeed, one of us has been involved in exploring this strategy since 1992 (ref. 37) and its subsequent applications in the preparation of potential bioactive compounds.38–41 In further continuation of this efforts we decided to adopt a similar one-pot synthetic strategy in the current study for accessing the targeted 3-(benzofuran-2-ylmethyl) substituted (pyrazolo/benzo)triazinone derivatives. Thus, the Pd/Cu-catalyzed coupling of the required terminal alkyne (1) [e.g. 5-methyl-3-(prop-2-ynyl)-7-propyl-3H-pyrazolo[4,3-d][1,2,3]triazin-4(5H)-one (1a) or 3-(prop-2-yn-1-yl)benzo[d][1,2,3]triazin-4(3H)-one (1b) or 6-bromo-3-(prop-2-yn-1-yl)benzo[d][1,2,3]triazin-4(3H)-one (1c)] obtained via the reported method8 with appropriate 2-iodophenols (2) afforded the expected products (3) in good to excellent yields (Scheme 1).
Scheme 1 Pd/Cu-catalyzed synthesis of 3-(benzofuran-2-ylmethyl) substituted (pyrazolo/benzo)triazinone derivatives. |
We used a similar reaction conditions reported by us previously20,37a with some suitable modifications as the use of terminal alkyne 1 for the coupling with 2-iodophenols was not explored earlier. For example, a longer reaction time i.e. 16–24 h was required in the previously reported method when the reaction was performed using 2–3.5 mol% of (Ph3P)2PdCl2, 3–5 mol% of CuI and 2 equiv. of Et3N in DMF. In the current study, we have observed that the reaction time could be decreased significantly with the marginal increase in the catalyst loading such as the use of 5 mol% of (Ph3P)2PdCl2 and CuI. Nevertheless, a range of reaction conditions by changing the catalyst, base and solvent was examined for the coupling of alkyne 1a with 2-iodophenol 2a (Table 1). Among the combination of various Pd-catalyst, ligand PPh3 and co-catalyst CuI tested (entries 1–7, Table 2), the combination of (Ph3P)2PdCl2 and CuI (entry 7, Table 2) was found to be most effective for affording the desired product 3a in high yield. Further, the use of a solvent such as DMF (entries 1–4, Table 2) or EtOH (entries 5 and 6, Table 2) was counterproductive and was not necessary in the present case because the excess of Et3N could play the role of base as well as solvent (entry 7, Table 2). Indeed, the reaction was completed within 2 h when performed in Et3N alone whereas product yield was not increased when longer reaction time was used (entry 7 vs. 8, Table 2). However, in view of the known toxicity of Et3N an attempt was made to reduce its volume by conducting the reaction using Et3N as a base, and dry DMF as a solvent in the presence of Pd(PPh3)2Cl2 and CuI as catalysts for 4 h. The desired product 3a was obtained in 55% yield in this case that was less than that of entry 7 of Table 2. The reaction did not proceed in the absence of co-catalyst CuI (entry 9, Table 2) or the catalyst Pd(PPh3)2Cl2 (entry 10, Table 2). Thus, the condition of entry 7 of Table 2 was used for the preparation of various other analogues of 3a.
Entry | Catalyst | Base | Solvent | Time (h) | % Yieldb |
---|---|---|---|---|---|
a Reactions were carried out using alkyne 1a (1 mmol), 2a (1.2 mmol), catalyst and base (2.0 mmol) in a solvent (5 mL) under nitrogen at 60 °C.b Isolated yield.c 5 mol% of each catalyst used.d 5 mol% of each catalyst and 10 mol% of PPh3 used. | |||||
1 | 10% Pd/C/CuI | K2CO3 | DMF | 6 | 0c |
2 | 10% Pd/C/PPh3/CuI | K2CO3 | DMF | 6 | 15d |
3 | Pd(OAc)2/PPh3/CuI | K2CO3 | DMF | 6 | 20d |
4 | Pd(OAc)2/PPh3/CuI | Et3N | DMF | 6 | 30d |
5 | Pd(OAc)2/PPh3/CuI | Et3N | EtOH | 6 | 20d |
6 | Pd(PPh3)2Cl2/CuI | Et3N | EtOH | 6 | 40c |
7 | Pd(PPh3)2Cl2/CuI | Et3N | Et3N | 2 | 85c |
8 | Pd(PPh3)2Cl2/CuI | Et3N | Et3N | 6 | 83c |
9 | Pd(PPh3)2Cl2 | Et3N | Et3N | 6 | 0c |
10 | CuI | Et3N | Et3N | 6 | 0c |
Three different terminal alkynes (1a–c) were employed for the reaction with a variety of 2-iodophenols (2a–f) in the presence of (Ph3P)2PdCl2 and CuI in Et3N alone at 60 °C (Table 1). The 2-iodophenols may contain the electron withdrawing group like CH3CO or NO2, halogen such as Cl or Br and electron donating group like CH3. The reactions proceeded smoothly in all these cases and reached to the completion within 2–3 hours affording a range of desired products (3) in high (>80%) yields. All the 3-(benzofuran-2-ylmethyl) substituted (pyrazolo/benzo)triazinone derivatives (3a-n) obtained were characterized with the help of 1H as well as 13C NMR and MS spectral data. Briefly, the presence of benzofuran ring in these compounds was indicated by the appearance of C-3 proton in 1H NMR spectra as a singlet in the range δ 6.9–7.2 and C-3 carbon in 13C NMR at 106 ppm. The linker “CH2” moiety appeared near δ 4.5 and 45 ppm in the 1H and 13C NMR spectra, respectively. Further the signal near 155 ppm in 13C NMR spectra was due to the “CO” carbon. These are shown in the selected 1H and 13C NMR data of two representative compounds 3c and 3i (Fig. 5). The compound 3c was further characterized by the presence of NCH3 group at δ 4.17 and the n-propyl protons at δ 2.93 (t, CH2), 1.86–1.73 (m, CH2) and 0.94 (t, CH3) in 1H NMR spectra.
From the view point of reaction mechanism,20,22–24,29,32 the involvement of two key steps e.g. Sonogashira coupling followed by intramolecular cyclization in the current approach are depicted in Scheme 2. Thus, the active Pd(0) species generated from Pd(II) participated in oxidative addition to 2 to afford the intermediate E-1 that on trans-organometallation with copper acetylide (generated in situ from compound 1 and CuI) followed by reductive elimination of Pd(0) furnished the internal alkyne E-3 via E-2 thereby completing the first catalytic cycle. In the second catalytic cycle the interaction of Cu(I) species with the alkyne moiety of E-3 (seemed to be aided by the proximate carbonyl group) formed E-4 that on intramolecular ring closure afforded E-5 and then E-6 and finally the product 3 with the regeneration of CuI to complete the second catalytic cycle.
Scheme 2 The plausible reaction mechanism for coupling-cyclization of fused 1,2,3-triazin-4-one based alkyne (1) with o-iodophenols (2) under Pd/Cu-catalysis. |
Next, the CM inhibitory potential of the synthesized compounds 3 was assessed in vitro42,43 with the help of an assay that measures the catalytic activity of CM during the conversion of substrate chorismate to the product prephenate. The known compound A-1 (Table 1) and a reported inhibitor44 C [i.e. 4-(3,4-dimethoxyphenethylamino)-3-nitro-5-sulfamoylbenzoic acid] was used as the reference compound in this assay. At the concentration of 30 μM the compound 3h, 3i and 3m showed 63.75 ± 1.09, 65.24 ± 0.98 and 64.90 ± 1.05% inhibition respectively that was comparable to 69.67 ± 0.57% inhibition of the reference compound A-1 (53.23 ± 0.13% inhibition at 10 μM was observed for the known inhibitor C). The other compounds that showed significant inhibition of CM included 3k (63.66 ± 2.07), 3n (63.75 ± 1.12) and 3e (61.28 ± 2.03), 3g (60.27 ± 3.22) and 3j (60.09 ± 2.35). Among rest of the compounds 3a and 3d appeared to be the least active (>45% inhibition) that was reflected in the results of docking studies. Indeed, the typical H-bonding interaction (as observed in case of active molecules) was missing in case of compound 3a (Fig. 6) due to the lack of any group at C-5 position of its benzofuran ring. Further, though the molecule participated in the hydrophobic interactions but the number of such contacts was lower than that of active molecules. On the other hand, the compound 3d (Fig. 7) though participated in a H-bond interaction through one of the trazinone nitrogen with the TYR110 of B chain however fewer hydrophobic contacts (with PRO66, ILE67, GLU68 of both A and B chain, and LYS60, TRP61, GLN64, PHE113 of A chain) were noted in this case.
Fig. 6 2D and 3D interaction of compound 3a with the interface residues of MtbCM (PDB: 2FP2). |
Fig. 7 2D and 3D interaction of compound 3d with the interface residues of MtbCM (PDB: 2FP2). |
Nevertheless, in general both in silico and in vitro evaluation suggested that the compounds belonging to the benzotriazinone series emerged as better active molecules than that from pyrazolotriazinone series. Indeed, all the best active molecules i.e. 3h, 3i and 3m emerged from the benzotriazinone series. Thus, an overview of SAR (Structure–Activity-Relationship) noted for the current series of benzofuran derivatives is presented in Fig. 8. Briefly, the “NO2” or “CH3CO” group at C-5 of the benzofuran ring was favorable for CM inhibition while other substituents like Cl, Br or Me or no substituent at this position reduced the activity. A “H” atom at C-7 position was generally better than the “Cl” or “Br” group at the same position in terms of activity. The inhibitory activity was influenced by the aryl/heteroaryl ring fused with the triazinone moiety and as mentioned above, the benzene ring rather than the pyrazole moiety was found to be beneficial for the CM inhibition of the current series of compounds.
In view of usefulness of the early detection of ADME (absorption, distribution, metabolism, and excretion) or pharmacokinetic properties of NCEs the computational ADME prediction of compounds 3h, 3i and 3m along with the reference compound A-1 was carried out using SwissADME web-tool.45 It is evident from Table 3 that desirable ADME was predicted for compound 3i and 3m with high GI absorption, no BBB (Blood Brain Barrier) penetration and no P-gp substrate potential. Indeed, the superior bioavailability score over the reference molecule A-1 was predicted for these compounds with no violation of Lipinski's criteria in contrary to A-1. However, BBB penetration was predicted for compound 3h. Thus, the benzofuran derivatives 3i and 3m seemed to be of further medicinal interest for the discovery and development of new and promising inhibitors of CM.
Properties | Molecules | ||
---|---|---|---|
A-1 | 3i | 3m | |
a logP: lipophilicity.b logS (ESOL): water solubility, calculated by ESOL method which Quantitative Structure–Property Relationship (QSPR) based model.c GI: gastrointestinal.d BBB: blood brain barrier.e P-gp: permeability glycoprotein.f MW: molecular weight. | |||
(i) Physicochemical | |||
Molecular weight (g mol−1) | 511.00 | 322.28 | 401.17 |
Consensus logPa | 4.41 | 2.34 | 2.88 |
logS (ESOL)b | −6.41 (poorly soluble) | −4.32 (moderately soluble) | −5.22 (moderately soluble) |
(ii) Pharmacokinetics | |||
GIc absorption | High | High | High |
BBBd permeation | No | No | No |
P-gpe substrate | No | No | No |
(iii) Druglikenss | |||
Lipinski rule | 1 violation: MWf > 500 | No violation | No violation |
Veber rule | No violation | No violation | No violation |
Bioavailability score | 0.17 | 0.55 | 0.55 |
Footnote |
† Electronic supplementary information (ESI) available: Experimental procedures and spectral data (file type: word). Copies of spectra of all target compounds synthesized. See https://doi.org/10.1039/d2ra05255e |
This journal is © The Royal Society of Chemistry 2022 |