Computational Analysis and Biological Activities of Oxyresveratrol Analogues, the Putative Cyclooxygenase-2 Inhibitors

Polyphenols are a large family of naturally occurring phytochemicals. Herein, oxyresveratrol was isolated from ethanolic crude extracts of Artocarpus lacucha Buch.-Ham., and chemically modified to derive its lipophilic analogues. Biological screening assays showed their inhibitory potency against cyclooxygenase-2 (COX-2) with very low cytotoxicity to the MRC-5 normal cell lines. At the catalytic site of COX-2, docking protocols with ChemPLP, GoldScore and AutoDock scoring functions were carried out to reveal hydrogen bonding interactions with key polar contacts and hydrophobic pi-interactions. For more accurate binding energetics, COX-2/ligand complexes at the binding region were computed in vacuo and implicit aqueous solvation using M06-2X density functional with 6-31G+(d,p) basis set. Our computational results confirmed that dihydrooxyresveratrol (4) is the putative inhibitor of human COX-2 with the highest inhibitory activity (IC50 of 11.50 ± 1.54 µM) among studied non-fluorinated analogues for further lead optimization. Selective substitution of fluorine provides a stronger binding affinity; however, lowering the cytotoxicity of a fluorinated analogue to a normal cell is challenging. The consensus among biological activities, ChemPLP docking score and the binding energies computed at the quantum mechanical level is obviously helpful for identification of oxyresveratrol analogues as a putative anti-inflammatory agent.


Introduction
Inflammation is the intrinsic immune response triggered by injury and infection. At the sites of tissue damage, the biosynthesis of prostaglandins in such inflamed tissue is significantly increased [1]. Prostaglandins-the lipid prostanoids derived from arachidonic acid-play a key role in innate immunity and inflammatory response. One of the major pathways of prostaglandin regulation in vivo involves the catalytic process of cyclooxygenase (COX) enzymes, of which three isoforms have been identified and are referred to as COX-1, COX-2 and the COX-1 splice variant, COX-3. COX isozymes catalyse the oxygenation of polyunsaturated fatty acids at the cyclooxygenase active site, and reduce the produced hydroperoxide to an alcohol at the peroxidase active site [2]. The basis of their difference is hydroperoxide-dependent activation; COX-2 approximately requires lower hydroperoxide level by a 10-fold amount [3]. COX-1 and COX-2 isoforms share approximately 60% sequence homology of amino acids with similar shape of the catalytic region [4,5]. However, the substitution of isoleucine numbered 523 in COX-1 with valine in COX-2 enlarges the volume of the COX-2 active site by approximately 25% [4,5]. It is worth noting that this valine residue at the pocket of COX-2 in this work is numbered 509 according to the originally retrieved protein structure of 3LN1 as a PDB entry code. Other a housekeeping enzyme to maintain homeostasis such as gastric epithelial cytopro kidney function and platelet aggregation; on the other hand, COX-2 is invo inflammation [7]. Therefore, COX-2 has become a target of anti-inflam therapeutics. Deleterious side effects can occur as a consequence of drug admin used to treat inflammatory conditions if the administered drugs synerg suppresssed the functions of COX-1.
The design of novel inhibitors to selectively bind to a specific COX isoform w harmful adverse effects has been a challenge. Most of the FDA-approved nons anti-inflammatory drugs (NSAIDs) such as diclofenac, ibuprofen, indom naproxen, oxaprozin and piroxicam are nonselective, inhibiting both COX-1 and isoforms, and are associated with well-known side effects. Recent progress in de coxibs as a COX-2 selective NSAID has been addressed in six broad categories [8, celecoxib remains legally available for use after rofecoxib and valdecoxib were wit due to their cardiotoxicity [10]. Gastrointestinal side effects and cardiovascular thr risk of celecoxib and other COX-2-selective drugs have also been reported [11,1 the use of celecoxib in patients with cardiovascular problems should be avoided an alternative to coxibs, medicinal plants exhibit potent anti-inflammatory proper due to a variety of their phytochemicals including phenolics [14][15][16][17], terpeniods alkaloids [20] and flavonoids [17,21,22].
Phenolic compounds are a major group of secondary plant metabolites and strong antioxidative, anticarcinogenic and anticyclooxigenase properties [23]. S based phenolic compounds are widely found in nature with a diverse array of ph logical and biochemical properties [24][25][26]. Trans-stilbenoids-the polypheno bene-are more stable and biologically relevant mediators with acidic and amp properties [8]. Oxyresveratrol (trans-2,3′,4,5′-tetrahydroxystilbene (1), Figure 1) is ral stilbenoid found in the root, stem and heartwood of plants including Artocarp ophyllus L. (Moraceae) [27], Scirpus maritimus (Cyperaceae) [28] and Morus alba L ceae) [29]. It has been brought to our attention due to a broad spectrum of biolo tivities such as antioxidant [30], anti-tyrosinase [31], antiviral [32], neuroprotec and anti-inflammatory activities [33,34]. Interestingly, the anti-inflammatory ac oxyresveratrol has been mentioned via various pathways such as the COX enzym oxide synthase (iNOS) and NF-B activation [33]. It has been reported that oxyres reduced not only paw edema induced by carrageenan in rats [35], but also expression levels of COX-2 isoform [36]. It can inhibit a 5-lipoxygenase (5-LO) e which induced inflammatory mediators named leukotrienes (LTs) for asthma with the half-maximal inhibitory concentration (IC50) in µM [34]. Despite the biological activity of COX inhibitors reported in the literature [37], insufficient metabolic stability and a high degree of non-specific binding have been the issues for consideration in lead optimization. Incorporation of fluorine into stilbene-like compounds is one effective way to increase their lipophilicity, resulting in an increase in metabolic stability and membrane permeation. Electron withdrawal by fluorine introduces a strong polarization of the C-F bond, which pronounces the electronic effects for drug-target interactions as in the case of celecoxib-the fluorine-containing clinical NSAID and a selective COX-2 inhibitor [38]. Thus, we introduced fluorine atom(s) to the double bond of an oxyresveratrol scaffold, yielding the fluorinated oxyresveratrol analogues.
To explore ligand-protein interactions, molecular docking is one of the most used approaches in structure-based drug design, which requires the three-dimensional structure of a macromolecular target [39,40]. The docking process involves: (i) sampling or generating a set of all possible ligand conformations and orientations; and (ii) evaluating the binding affinity of each ligand pose placed at the active site of its target. Therefore, top-ranked poses of docked ligand(s) depend on the quality of the search algorithm and scoring function. Different sampling algorithms such as shape matching, systematic, stochastic and simulation methods can be applied for a conformational search [39,41,42]. Then, the binding affinity of each generated pose is estimated and ranked by the scoring functions, which are traditionally classified into the force field-based, knowledge-based, and empirical scoring functions [39,41,42]. Effective scoring functions have been a challenge in molecular docking and often limit their reliability [43,44]. Therefore, quantum mechanical (QM) methods have been increasingly used in drug design due to a demand for the accurate estimation of binding affinities [45][46][47]. Some QM computed parameters, such as a dipole moment and atomic charges, have been found to be useful for improving the accuracy of a docking protocol [43,48].
In this work, we exploit the structure-based computational approaches to investigate the binding interactions of oxyresveratrol and its analogues at the catalytic site of COX-2. Molecular docking procedures with different scoring functions are used to predict the binding modes and affinities of each COX-2/ligand complex. We perform energy calculations on the optimized COX-2 complexes with the best-docked pose using M06-2X density functional and 6-31G(d) basis set in both vacuo and implicit aqueous solvation. The computed binding energetics of complexes are discussed and compared with the cytotoxicity and inhibitory potency screened by using the human COX-2 assay kit. This computational analysis is aimed at providing a more accurate description of the binding affinity of studied oxyresveratrol analogues.

Binding Affinity of COX-2/Celecoxib Complex
To understand the nature of the catalytic pocket of COX-2, we performed a short molecular dynamics simulation of the complex of COX-2/co-crystallized celecoxib (PDB entry: 3LN1, Figure S1, in the Supplementary Materials) to relax the complex geometry. A long and hydrophobic channel of the cyclooxygenase active site is applicable to NSAIDs for binding [2,49]. Figure 2 illustrated the molecular interactions of a relaxed pose of co-crystallized celecoxib, the known selective COX-2 inhibitor. Its trifluoromethyl group was surrounded by the methyl groups of Val102, Leu345 and Leu517 at the binding entrance formed by Arg106 and Tyr341. These two amino acids appeared to reveal the hydrogen bonds and/or charged interactions with the carboxyl groups of non-selective NSAIDs [11,50,51]. Unlike traditional NSAIDs, celecoxib placed its sulfamoylphenyl moiety in a side pocket located at the polar His75, Glu178 and Arg499 residues. This pocket, containing the hydrogen bond acceptor, is crucially responsible for the COX-2 selectivity [5,52,53]. A similar binding mode could also be observed for other coxibs and non-selective NSAIDs [50].
Hydrogen bonding between the sulfonyl group and amino groups Phe504 and Arg499 could, respectively, be observed in 2.46 Å and 2.19 Å as being as strong as that found between the sulfonamide group of celecoxib and the carbonyl side chain of Gln178 in 2.06 Å and Leu338 in 1.97 Å. This agreed with the previous study [51]. We note that CH-pi interactions of this sulfamoylphenyl ring with Val509 (in 3.87 Å) and Ser339 (in 2.26 Å) were also helpful for conformational stabilization (Figure 2b). Other types of noncovalent interactions of celebcoxib occurred at the binding site and are provided in Table S1. The entrance of the COX-2 channel was where the hydrogen bonding occurred between the fluoromethyl group of celecoxib and Arg106 in 2.33 Å and 2.98 Å, as well as between the pyrazole nitrogen of celecoxib and Tyr341 in 3.24 Å. Regarding the per- residue decomposition of binding energy reported in the previous study, it was reported that Arg106, Leu338, Ser339, Arg499, Phe504, and Val509 are the crucial amino acids and commonly contribute to the binding energy of COX-2/celecoxib complex [54,55]. Hydrogen bonding between the sulfonyl group and amino groups Phe504 and Arg499 could, respectively, be observed in 2.46 Å and 2.19 Å as being as strong as that found between the sulfonamide group of celecoxib and the carbonyl side chain of Gln178 in 2.06 Å and Leu338 in 1.97 Å. This agreed with the previous study [51]. We note that CH-pi interactions of this sulfamoylphenyl ring with Val509 (in 3.87 Å) and Ser339 (in 2.26 Å) were also helpful for conformational stabilization (Figure 2b). Other types of noncovalent interactions of celebcoxib occurred at the binding site and are provided in Table S1. The entrance of the COX-2 channel was where the hydrogen bonding occurred between the fluoromethyl group of celecoxib and Arg106 in 2.33 Å and 2.98 Å, as well as between the pyrazole nitrogen of celecoxib and Tyr341 in 3.24 Å. Regarding the perresidue decomposition of binding energy reported in the previous study, it was reported that Arg106, Leu338, Ser339, Arg499, Phe504, and Val509 are the crucial amino acids and commonly contribute to the binding energy of COX-2/celecoxib complex [54,55].

Validation of Docking Protocols
To validate the reliability of three scoring functions, we docked the relaxed pose of celecoxib at the binding site of COX-2 and investigated the top-scored pose ranked by AutoDock, ChemPLP and GoldScore functions. Scoring function in AutoDock4 is a semiempirical free energy force field [56], whereas GoldScore [57] is the empirical fitness function, originally provided with GOLD Suite [58]. The ChemPLP [59] empirical function uses multiple linear potentials and derives hydrogen bonding terms from the total free energy change [60]. Pose prediction by these scoring functions was examined by the fact that the higher the score of GoldScore and ChemPLP functions (or the more negative value in the estimated binding affinity from AutoDock), the better pose is suggested.
Docking results revealed that the pyrazole and sulfomoylphenyl rings of the best-

Validation of Docking Protocols
To validate the reliability of three scoring functions, we docked the relaxed pose of celecoxib at the binding site of COX-2 and investigated the top-scored pose ranked by AutoDock, ChemPLP and GoldScore functions. Scoring function in AutoDock4 is a semi-empirical free energy force field [56], whereas GoldScore [57] is the empirical fitness function, originally provided with GOLD Suite [58]. The ChemPLP [59] empirical function uses multiple linear potentials and derives hydrogen bonding terms from the total free energy change [60]. Pose prediction by these scoring functions was examined by the fact that the higher the score of GoldScore and ChemPLP functions (or the more negative value in the estimated binding affinity from AutoDock), the better pose is suggested.
Docking results revealed that the pyrazole and sulfomoylphenyl rings of the bestdocked poses were nicely overlaid on those rings of the relaxed pose of co-crystallized celecoxib ( Figure 2); therefore, hydrogen bonding interactions with the key amino acids were still preserved. The binding affinity of the best-docked pose of celecoxib estimated by AutoDock, was −11.24 kcal/mol (Table 1), agreed well with the previously reported values ranging from −10.01 kcal/mol [61] to −11.00 kcal/mol [62,63]. The fitness score given by ChemPLP and GoldScore for the best-docked poses of celecoxib was, respectively, 95.56 and 79.07. AutoDock, ChemPLP and GoldScore scoring functions can reproduce the native pose of celecoxib with the RMSD value of 0.68, 0.55 and 0.56 Å, respectively (Figures 2c and S2).

Molecular Docking of Oxyresveratrol Analogues
First, we explore the binding modes of our stilbenoid compounds 1-3 and dihydrostilbenoid compounds 4-6 by using molecular docking techniques. The best-docked pose of all six compounds was bound in the COX-2 pocket at the same position of the native pose of co-crystallized celecoxib regardless of the applied scoring functions. Hydrogen bonding between Arg499 and the methoxy or acetoxy group of C3 on ring B of oxyresveratrol analogues provides the binding stability, in a similar fashion to that which occurred with sulfonamide of celecoxib ( Figure 3). This key interaction was also observed in other coxibs such as etoricoxib and valdecoxib and chromone derivatives bound to COX-2 with high affinity [63].
AutoDock ranked the docked poses by the estimated binding energy. In our study, acetoxylated analogues (3 and 6) of stilbenoid and dihydrooxyresveratrol were suggested to be the strongest-bound inhibitors for COX-2 with the binding energy of −9.21 and −9.26 kcal/mol compared with −11.24 kcal/mol of celecoxib (Table 1). Figure 3a,b show the binding modes of the best-docked pose of compounds 3 and 6 overlaid with that of celecoxib. Ring B of both compounds was aligned in the disposition of sulfamoylphenyl of celecoxib, while their ring A slipped nearby the position of the pyrazole ring of celecoxib.
With a similar ranking trend to AutoDock, compound 3 and compound 6 were also the top-ranked ligands suggested by GoldScore fitness function with the docked score of 80.20 and 78.62 regarding that of 79.07 for celecoxib (Table 1). Hydrogen-bonding interactions were observed between Arg499 with acetoxy groups on ring B of the best-docked pose of compounds 3 ( Figure 3c) and 6 ( Figure 3d). The acetoxy groups on ring B of compound 3 were able to form hydrogen bonds with Arg106 at the binding entrance, however with Tyr341 at the binding entrance in case of compound 6.
ChemPLP fitness function-the effective docking score of the GOLD docking Suite-predicted the highest score of 83.45 for acetoxylated analogue (6) of dihydrooxyresveratrol, against that of 95.56 for celecoxib. Both carboxylate groups on ring A of compound 6 can form hydrogen bonds with Arg499 ( Figure 3f) in a similar fashion to sulfonamide of celecoxib. This was no longer observed for dihydrooxyresveratrol (4), the second-best ligand suggested by ChemPLP fitness function, with a score of 74.09 (Table 1). Instead, there was a hydrogen bond between the hydroxyl group at C1 of ring A of compound 4 and Ser339 (Figure 3e). ChemPLP also suggested that oxyresveratrol (1), was the third-best ligand with the score of 68.34 (Table 1). However, its binding mode was different from the top two ligands, i.e., ring B was interchanged in the side pocket of COX-2 so that one of its hydroxyl groups on this ring formed hydrogen bonds with Gln178 and Leu338. Molecules 2022, 27, x 9 of 30  (4) and (f) compound 6 ranked by ChemPLP-were superimposed onto the relaxed pose of co-crystallized celecoxib (green). Some residues are omitted for clarity.  (4) and (f) compound 6 ranked by ChemPLP-were superimposed onto the relaxed pose of co-crystallized celecoxib (green). Some residues are omitted for clarity.

Biological Activity of Oxyresveratrol Analogues
Oxyresveratrol was successfully extracted and isolated from Artocarpus lacucha Buch.-Ham. It was chemically modified to obtain the methoxylated and acetoxylated analogues including those in a series of hydrogenated form with excellent yields.
The in vitro COX-2 inhibitor screening assay and cytotoxicity (CC 50 ) in normal MRC-5 cell lines strongly indicates that dihydrooxyresveratrol (4) and its parent oxyresveratrol (1) are the potent lead among other compounds studied with the IC 50 values of 11.50 ± 1.54 and 14.50 ± 2.04 (Table 1). Compared with the IC 50 value of 0.09 ± 0.01 µM for celecoxib, oxyresveratrol (1) and its stilbenoid derivatives with methoxy (compound 2) and acetoxy (compound 3) substituents showed their IC 50 values of 14.50 ± 2.04, 25.00 ± 2.34, and 23.00 ± 2.02 µM, respectively. Upon hydrogenation, the IC 50 values for dihydrostilbenoids 4, 5, and 6 were, respectively, 11.50 ± 1.54, 18.10± 2.07 and 19.00 ± 2.01 µM, lower than that of their corresponding parent stilbenoids. Therefore, we note the high potency of dihydrostilbenoids for COX-2 inhibition. The most potent compound in the series is dihydrooxyresveratrol (4), a reductive form of oxyresveratrol (1), due to its character of being a good hydrogen donor with some degree of hydrogenation.
On the MRC-5 cell lines which are normal human fibroblasts, the cytotoxic effect of all analogues was evaluated and compared with the CC 50 value of >100 µM for doxorubicin as a positive control. The naturally occurring oxyresveratrol (1) and its hydrogenated derivative (compound 4) demonstrated the significant differences in cytotoxicity with CC 50 values of 57.48 ± 2.34 and 106.02 ± 3.86 µM. In comparison with these polyhydroxy parent compounds 1 and 4, the methoxy-containing stilbenoid (2) with CC 50 value of 118.21 ± 2.69 µM, the methoxy-containing dihydrostilbenoid (5) with CC 50 value of 130.32 ± 3.04 µM, and acetoxy-containing dihydrostilbenoid (6) with CC 50 value of 120.51 ± 2.37 µM are high in cytotoxicity, although they are highly lipophilic and suitable for further lead optimization. Considering the cytotoxicity of acetoxy-containing stilbenoid (3), its CC 50 value was 77.66 ± 3.07 µM, which was lower than that of doxorubicin. This is probably due to a slow hydrolysis occurring in the cells to release hydroxylated metabolites. The acetoxylated derivatives do not only present some extent of biological activities, but are also a worthy choice for undergoing the structural modification.

Rescoring Oxyresveratrol Analogues by Quantum Mechanical (QM) Binding Energy
Regarding our biological and docking studies, it is difficult to come to a consensus on the structural-activity relationship of oxyresveratrol analogues. However, the binding energy (BE) computed using M06-2X density functional [64] and 6-31G+(d,p) basis set is able to re-score the oxyresveratrol analogues bound to COX-2 considering their biological potency. This computed BE represents the binding affinity of COX-2/ligand complexes, derived from the difference of energy between the bound state of a complex optimized at a M02-2X/6-31G(d) level of theory and the unbound states of both the receptor and ligand according to Equation (2) in the section of computational methods. Among the six nonfluorinated compounds previously discussed, dihydrooxyresveratrol (4) was suggested to be the strongest bound ligand in a pocket of COX-2 with the BE of −85.29 kcal/mol in vacuo and of −47.16 kcal/mol in an implicit aqueous solvation (Table 1 and Figure 4). This agreed well with its lowest IC 50 value of 11.50 µM (Table 1). Nine kcal/mol was compensated by counterpoise (CP) correction for the basis set superposition error on the computed BE of the noncovalently bound complex; therefore, the corrected IEs for the complex of COX-2/compound 4 in vacuo and implicit aqueous solvation were −75.98 and −37.86 kcal/mol, respectively.
Interactions of oxyresveratrol analogues with solvated water molecules and key residues of the amino acid at a COX-2 pocket play a role in their binding affinity (Figures 4, 5 and S3). The hydroxyl group at C5 on ring B of compound 4 interacted within 1.42 Å with Ser339, whereas that at C3 of the same ring interacted with Glu510 in 1.49 Å (Figure 5c) and resulted in a slight shift of ring B of compound 4 from the sulfamoylphenyl motif of celecoxib (Figure 4a). It is worth noting that hydrogen bonding interactions within 2 Å between the hydroxyl group positioned at C3 of ring B and Arg499, as well as that at C1 of ring A and Tyr371, occurred not only in the case of compound 4, but also in oxyresveratrol (1) and acetoxylated compounds 3 and 6 ( Figure 5). In addition to key interactions with Arg499 and Tyr371, T-shaped pi-pi interactions of aryl ring B of these compounds with His75 and Tyr34 could be observed within 5 Å. For oxyresveratrol (1), aromatic pi interactions between its ring A with an alkyl chain of Leu338, Ala513 and Val335 were also observed as well as hydrogen bonding between its hydroxyl groups at C5 of ring B with Leu338 in 1.50 Å (Figure 5a). For acetoxylated analogue (3) of oxyresveratrol, the hydrogen bond of its acetoxy group at C3 of ring A with Ser516 could be observed in 2.80 Å (Figure 5b), whereas this occurred with Ser516 in 2.87 Å in the case of the acetoxylated analogue (6) of dihydrooxyresveratrol (Figure 5d). Interactions of oxyresveratrol analogues with solvated water molecules and key residues of the amino acid at a COX-2 pocket play a role in their binding affinity (Figures 4b  and 5). The hydroxyl group at C5 on ring B of compound 4 interacted within 1.42 Å with Ser339, whereas that at C3 of the same ring interacted with Glu510 in 1.49 Å (Figure 5c) and resulted in a slight shift of ring B of compound 4 from the sulfamoylphenyl motif of celecoxib (Figure 4a). It is worth noting that hydrogen bonding interactions within 2 Å between the hydroxyl group positioned at C3 of ring B and Arg499, as well as that at C1 of ring A and Tyr371, occurred not only in the case of compound 4, but also in oxyresveratrol (1) and acetoxylated compounds 3 and 6 ( Figure 5). In addition to key interactions with Arg499 and Tyr371, T-shaped pi-pi interactions of aryl ring B of these compounds with His75 and Tyr34 could be observed within 5 Å. For oxyresveratrol (1), aromatic pi interactions between its ring A with an alkyl chain of Leu338, Ala513 and Val335 were also observed as well as hydrogen bonding between its hydroxyl groups at C5 of ring B with Leu338 in 1.50 Å (Figure 5a). For acetoxylated analogue (3) of oxyresveratrol, the hydrogen bond of its acetoxy group at C3 of ring A with Ser516 could be observed in 2.80 Å (Figure 5b), whereas this occurred with Ser516 in 2.87 Å in the case of the acetoxylated analogue (6) of dihydrooxyresveratrol (Figure 5d). Hydroxy groups of dihydrooxyresveratrol 4 and its parent compound 1 can serve as either hydrogen-bond donors or hydrogen-bond acceptors, whereas methoxy groups of compounds 2 and 5 can function only as hydrogen-bond acceptors. The more negative the value of the ligand binding to COX-2, the greater the affinity of the receptor-ligand established at the binding site, indicating that the hydroxyl substituent considerably affects the ligand stability upon binding. This can be supported by other computational evidence that reported that the number of phenolic groups of ligand compounds and the ligand size influence their binding activity [65]. Considering the basis set superposition error upon noncovalent binding, we performed CP correction on the optimized complex geometries to obtain the corrected BE. From Table 1, the corrected BE of oxyresveratrol (1) was −68.47 kcal/mol in vacuo and −24.84 kcal/mol in implicit aqueous solvation, whereas that of its methoxylated stilbenoid (2), was −47.30 and −15.98 kcal/mol, respectively. Strong affinity, however, was observed for acetoxylated compounds 3 and 6, although their acetoxy groups act as hydrogen-bond acceptors, similar to the methoxy groups. It has been suggested that the acetoxylated analogue (6) of dihydrostilbenoid tightly binds to COX-2 with the corrected BE value of −70.34 kcal/mol in vacuo and −24.88 kcal/mol in implicit solvation, which is considerably stronger than that of −55.04 kcal/mol and −25.76 kcal/mol for celecoxib (Table 1). This arose from the fact that the acetoxy group is bulky, highly versatile and able to capture interactions at the binding site. A functional group with bulkiness that has rapidly transformed into a hydrophilic functional group through metabolism in circulation exhibits a good topical anti-inflammatory activity [66]. Therefore, binding affinities of polyphenolic compounds 1 and 4, as well as acetoxylated compounds 3 and 6, were found to be much stronger than that of the corresponding methoxylated compounds 2 and 5. The methoxylated analogues have also been reported in the literature as having poor inhibitory activity against COX-2 [67]. Hydroxy groups of dihydrooxyresveratrol 4 and its parent compound 1 can serve as either hydrogen-bond donors or hydrogen-bond acceptors, whereas methoxy groups of compounds 2 and 5 can function only as hydrogen-bond acceptors. The more negative the value of the ligand binding to COX-2, the greater the affinity of the receptor-ligand established at the binding site, indicating that the hydroxyl substituent considerably affects the ligand stability upon binding. This can be supported by other computational evidence that reported that the number of phenolic groups of ligand compounds and the ligand size influence their binding activity [65]. Considering the basis set superposition error upon noncovalent binding, we performed CP correction on the optimized complex geometries to obtain the corrected BE. From Table 1, the corrected BE of oxyresveratrol (1) was −68.47 kcal/mol in vacuo and −24.84 kcal/mol in implicit aqueous solvation, whereas that of its methoxylated stilbenoid (2), was −47.30 and −15.98 kcal/mol, respectively. Strong affinity, however, was observed for acetoxylated compounds 3 and 6, although their acetoxy groups act as hydrogen-bond acceptors, similar to the methoxy groups. It has been suggested that the acetoxylated analogue (6) of dihydrostilbenoid tightly binds to COX-2 with the corrected BE value of −70.34 kcal/mol in vacuo and −24.88 kcal/mol in implicit solvation, which is considerably stronger than that of −55.04 kcal/mol and −25.76 kcal/mol for celecoxib (Table 1). This arose from the fact that the acetoxy group is bulky, highly versatile and able to capture interactions at the binding site. A functional group with bulk- In brief, a series of dihydrogenated compound analogues (compound 4, 5 and 6) appears to be superior to the stilbenoid analogues (compound 1, 2 and 3) considering both of their computed binding energetics and biological activities (Table 1).

Fluorinated Analogues
Fluorine-containing compounds have been extensively investigated in drug development. Due to the unique chemical properties of fluorine, fluorine substitution in stilbenelike compounds has been reported to enhance the binding affinity of a compound to a protein target [68,69]. Our computational finding reveals that the incorporation of fluorine with the saturated bridge between both aromatic rings gives a strong affinity towards binding to COX-2. From Table 1, the corrected IEs of difluorooxyresveratrol (7) were −78.44 kcal/mol in vacuo and −37.90 kcal/mol in implicit aqueous solvation, the lowest values considering those values of other studied oxyresveratrol derivatives, including fluorinated compounds 8 and 9. Hydrogen bonding interactions between the hydroxyl group on ring B of compound 7 with Arg499 ( Figure 6) within 2 Å was still preserved in a similar fashion to celecoxib, a COX-2 specific NSAID. A negative electron density of difluorine atoms of compound 7 can facilitate its binding affinity without altering any steric bulk (Figure 7). a protein target [68,69]. Our computational finding reveals that the incorporation of fluorine with the saturated bridge between both aromatic rings gives a strong affinity towards binding to COX-2. From Table 1, the corrected IEs of difluorooxyresveratrol (7) were −78.44 kcal/mol in vacuo and −37.90 kcal/mol in implicit aqueous solvation, the lowest values considering those values of other studied oxyresveratrol derivatives, including fluorinated compounds 8 and 9. Hydrogen bonding interactions between the hydroxyl group on ring B of compound 7 with Arg499 ( Figure 6) within 2 Å was still preserved in a similar fashion to celecoxib, a COX-2 specific NSAID. A negative electron density of difluorine atoms of compound 7 can facilitate its binding affinity without altering any steric bulk (Figure 7).  With a similar binding mode of compound 7, acetoxylated analogue (9) of difluorooxyresveratrol has the corrected BE of −68.65 kcal/mol in vacuo and −24.04 kcal/mol implicit aqueous solvation, which lies in between that of difluorostilbenoid compound 7 and celecoxib (Table 1). This also demonstrates the significance of hydroxyl groups and the fluorine in a chemical structure of our ligands. Fluorine-containing compounds have been reported to have the ability to improve several pharmacokinetic and physicochemical properties and enhance their biological activities [70]. However, our studied fluorinated analogues appear to be cytotoxic to the normal cell lines. We have been investigating the fluorine substitution into oxyresveratrol derivatives to lower their in vitro cytotoxicity with the improved synthetic yield.
Among our studied compounds, the strongest binding affinity with COX-2 was found to be difluorooxyresveratrol (7). Regarding the cell cytotoxicity, we suggest dihydrooxyresveratrol (4) is the most potent and the second-best bound ligand for further lead optimization. been reported to have the ability to improve several pharmacokinetic and physicochem cal properties and enhance their biological activities [70]. However, our studied fluor nated analogues appear to be cytotoxic to the normal cell lines. We have been investiga ing the fluorine substitution into oxyresveratrol derivatives to lower their in vitro cyt toxicity with the improved synthetic yield. Among our studied compounds, the strongest binding affinity with COX-2 w found to be difluorooxyresveratrol (7). Regarding the cell cytotoxicity, we suggest dih drooxyresveratrol (4) is the most potent and the second-best bound ligand for further lea optimization.

Atomic Charges of Celecoxib and Oxyresveratrol Analogues
The electronic effect of ligand structures upon binding to COX-2 was further exam ined by considering a distribution of Mulliken atomic charges computed using M06-2 density functional and 6-31G+(d,p) basis set. We found the overall negative charge of ca bon atoms on a sulfamoylphenyl ring of celecoxib. This confirmed the fact that the lon pair electrons of pyrazolic nitrogen of celecoxib can be delocalized to the sulfamoylphen motif, which is helpful for its binding stability and anti-inflammatory activity. The pos tive charge of the sulfur was found to be a value of +1.514e for the unbound state celecoxib and +0.306e for its bound pose ( Table 2), indicating that the sulfonamide grou of celecoxib is an electron acceptor, which specifically contributes to a strong interactio with Arg499.

Atomic Charges of Celecoxib and Oxyresveratrol Analogues
The electronic effect of ligand structures upon binding to COX-2 was further examined by considering a distribution of Mulliken atomic charges computed using M06-2X density functional and 6-31G+(d,p) basis set. We found the overall negative charge of carbon atoms on a sulfamoylphenyl ring of celecoxib. This confirmed the fact that the lone-pair electrons of pyrazolic nitrogen of celecoxib can be delocalized to the sulfamoylphenyl motif, which is helpful for its binding stability and anti-inflammatory activity. The positive charge of the sulfur was found to be a value of +1.514e for the unbound state of celecoxib and +0.306e for its bound pose ( Table 2), indicating that the sulfonamide group of celecoxib is an electron acceptor, which specifically contributes to a strong interaction with Arg499. In comparison with celecoxib, a highly positive charge was also observed at C4 on ring B of the unbound pose of oxyresveratrol (1), dihydrooxyresveratrol (4) and difluorooxyresveratrol (7) as shown in Figure S4. It is important to note that only compound 4 preserved a highly positive charge when binding to COX-2, i.e., +0.642e (Table 2 and Figure S5). Our finding suggests that the characteristic of being an electron acceptor upon binding appears to be essential in design selective leads for COX-2 receptor.

Statistical Analysis
To investigate the effect of lipophilic substituents: methoxy and acetoxy groups of our oxyresveratrol analogues that have moved towards those with the ability of ligand binding, the absolute values of docking scores obtained from all three scoring functions in Table 1 were averaged out for each compound and compared with that of celecoxib. For any compounds that could remain in the COX-2 pocket as firmly as celecoxib, either their docking scores or estimated binding energetics were expected to be at least comparable to those of celecoxib. Therefore, the insignificance of the mean difference between each compound and celecoxib was statistically analyzed by using the Dunnett t-test, showing the mean square error and standard error of 59.086 and 6.276, respectively (Table S2). At 95% confidence, only acetoxylated analogues, i.e., compound 3, 6 and 9 showed no significant difference from celecoxib. This highlights the importance of acetoxy groups in ligand structure upon binding.
To confirm the significance of the mean difference of our computed binding affinity values from that of celecoxib (Table 1), we performed a Dunnett t-test statistical analysis of the average value from all four binding energetics of each compound against that of celecoxib at alpha level of 0.05. With the mean square error of 15.730 and a standard error of 2.804, the BE of the bound compounds 4, 6 and 7 were found to be significantly different from that of celecoxib at 95% confidence (Table S3). This agreed well with the top rank of dihydrooxyresveratrol (4) and difluorostilbenoid (7) suggested by the QM re-score.

Conclusions
Oxyresveratrol analogues can be exploited in the exploration of structurally modified phytochemicals for a novel design of COX-2 inhibitors. We extracted oxyresveratrol from the inner wood of Artocarpus lacucha Buch.-Ham. and modified its chemical structure to obtain its methoxylated and acetoxylated analogues, including those in a series of hydrogenated forms with excellent yields and fluorinated analogues. Our compounds exhibit anti-inflammatory potency against human COX-2 with low cytotoxicity to normal cell lines. Docking protocols displayed hydrogen bonds and hydrophobic pi interactions of oxyresveratrol analogues with key residues of COX-2 and some solvated water molecules. With the statistical significance at 95% confidence level of the mean difference, GoldScore and AutoDock scoring suggested acetoxylated analogues as the best leads whereas ChemPLP fitness function ranked dihydrooxyresveratrol (4) as the best one, agreed with its potency of IC 50 against COX-2. Other potential targets that regulate immune responses apart from COX could be further explored and identified by the inverse docking of this compound to members of a protein library and analysis of the similarity of binding patterns that occurred across the protein targets [71].
To predict a more accurate binding affinity, the binding energy of each COX-2/ligand complex was eventually computed at a quantum mechanical level, both in vacuo and with implicit aqueous solvation. The computed results demonstrated that dihydrooxyresveratrol (4), the hydrogenated compound, has the strongest binding affinity among non-fluorinated analogues, including its parent oxyresveratrol (1). This agreed with the screening results of the biological assay, showing the highest anti-inflammatory activity for dihydrooxyresveratrol (4) with very low cytotoxicity to the MRC-5 normal cell line. Analysis of atomic charges highlights the importance of the ligand characteristic of being an electron acceptor on the aromatic ring. Such electrostatic and van der Waals interactions have contributed to the enthalpic component of the binding affinity. However, entropic and solvation effects could also be captured and accounted for ligand-protein binding by free energy calculations such as a free energy perturbation, a linear response approximation (LRA) and a combination of molecular mechanical energies with the Poisson-Boltzmann continuum solvent approaches [72,73]. Incorporation of fluorine into hydrogenated analogues leads to a stronger binding affinity; however, its selective substitution with a low cytotoxicity has been challenging.

Chemistry
Oxyresveratrol (1) was chemically modified via semisynthetic pathways to produce stilbenoid compounds 1-3 and reductive compounds, dihydrostilbenoid 4-6 with high synthetic yields. Firstly, oxyresveratrol (1), a naturally occurring compound, was hydrogenated using H 2 on Pd/C in ethanol at room temperature to afford dihydrooxyresveratrol (2,3 ,4,5 -tetrahydroxybibenzyl, 4). Full O-acylations of 1 and 4 were carried out with acetic anhydride in pyridine, achieving compound 3 and 6, respectively. Full O-methyl analog 2 was produced from oxyresveratrol (1) with methyl iodide and potassium carbonate in acetone. The hydrogenation of compound 2 was subsequently performed under H 2 on Pd/C in ethanol at room temperature, giving compound 5.
The progress of all reactions was monitored on thin-layer chromatography using coated with silica gel 60 (Merck F254). TLC plates were examined under UV light (254 nm and 366 nm) and/or coloring with potassium permanganate stain. Column chromatography was performed on Merck silica gel 60 (0.063-0.200 mm, 70-230 mesh size). The chromatographic solvents were mentioned as v/v ratios. Air and moisture-sensitive reactions were conducted under an inert atmosphere of argon or nitrogen.

Isolation of Oxyresveratrol from A. lacucha
The heartwood of Artocarpus lacucha Buch.-Ham. (Synonym A. lakoocha Roxb.) was collected after botanical identification. The dried powder of A. lacucha heartwood was soaked overnight in ethanol for 72 h prior to each extraction. Three extracts were combined, filtered, and concentrated using a rotary evaporator. The ethanolic crude extract was then purified by column chromatography with 10-15% ethanol in dichloromethane as the eluent to give oxyresveratrol [74]. Oxyresveratrol

General Procedure for the Synthesis of Compound 2
A mixture of oxyresveratrol (1, 1.7 g, 7.0 mmol) and potassium carbonate (9.67 g, 70 mmol) in 60 mL of extra dry acetone was added with methyl iodide (35 mL, 70.0 mmol) at room temperature under an atmosphere of argon gas. The mixture was then refluxed for 72 h until the completion of the reaction. The solid in the reaction mixture was filtered off at room temperature. The filtrate was then evaporated under reduced pressure to obtain a yellow oil of crude product. A quick filtration of the crude product was carried out by using a short column with ethyl acetate as eluent. After removal of the solvent and hexane-recrystallization, the desired product of 2,3 ,4,5 -tetramethoxystilbene (2) was obtained as a white solid [30] with 75% yield, mp 81.2-81.9 • C. 1 1570, 1504, 1454, 1421, 1344, 1327, 1313, 1286, 1236, 1193, 1151, 1118, 1105, 1061, 1026 Oxyresveratrol (1, 1.00 mmol) or compound 2 (1 mmol) was dissolved in 5 mL of ethanol over 10 mol% Pd/C (10.64 mg, 0.10 mmol). The reaction mixture was stirred for 1 h under hydrogen gas in a balloon at room temperature. After the reaction had been complete, the catalyst was filtered off. The filtrate was evaporated under a reduced pressure to provide products 4 or 5 [31].  Oxyresveratrol (1, 1.00 mmol) or compound 4 (1.00 mmol) was dissolved in 15 mL of anhydrous pyridine. The mixture was stirred at ca. 0 • C in an ice-cold water bath before slowly adding excess acetic anhydride (10 mmol). The reaction mixture was stirred for 2 h. A complete reaction mixture was then quenched with 30 mL of ice-cold water and extracted with 30 mL of ethyl acetate. Three organic extracts were combined before washing with 30 mL of 10% HCl for pyridine removal and 30 mL of water three times. The organic layer was neutralized by 30 mL of sodium bicarbonate, then dried with anhydrous sodium sulfate and filtered afterwards. After the removal of the solvent, the desired product was obtained as compound 3 [30] or 6 [75].   Figure S16); 13

Characterization
Oxyresveratrol and its derivatives were characterized by 1 H and 13 C NMR recorded on Bruker NMR-systems (300 and 600 MHz). The deuterated NMR solvents were supplied from Merck for 99.8% chloroform-d 1 and Eurisotop for methanol-d 4 . Chemical shift (δ) is given in parts per million ppm referenced to TMS or specified solvent signals. The abbreviations of the signals are used to designate chemical shift multiplicities: s = singlet, d = doublet, t = triplet, q = quartet, m = multiplet, dd = doublet of doublets. Coupling constants are given in Hertz (Hz). 13 C NMR spectra were measured as ATP-or DEPTQexperiments. Infrared spectra were collected by using Perkin Elmer infrared spectrophotometer (spectrum GX) and Shimadzu IRAffinity-1S. The melting points were measured on a Kruss Optronic KSP1N apparatus. High resolution mass spectra (HRMS) were recorded on an orbitrap TM Exploris 120 Mass Spectrometer (Thermo Fisher Scientific, Bremen, Germany) using electrospray ionization (ESI) in positive ion mode.

COX-2 Inhibitory Activity
The in vitro inhibitory of tested compounds were evaluated using a human COX-2 human inhibitor screening assay kit (Item no. 701080) supplied by Cayman Chemicals USA. Celecoxib was used as a positive control. To determine the half-maximal inhibitory concentration (IC 50 ), the inhibitory activity of tested compounds was measured in triplicate at each concentration according to the procedure of COX-2 inhibitory immunoassay previously reported [76,77]. COX-2 enzyme and tested compounds at various concentrations were prepared in the reaction buffer (0.1 M Tris-HCl, pH 8.0, containing 5 mM EDTA and 2 mM phenol) containing heme. Prepared solutions were incubated at 37 • C for 10 min and added to arachidonic acid to initiate reactions. COX-derived prostaglandin H 2 was produced during a 2 min reaction incubated at 37 • C, then quenched with saturated stannous chloride solution. The quenched reactions were incubated at room temperature for 5 min to produce prostaglandin F 2α . Dilutions were made with the ELISA buffer and measured by using a PGE 2 ELISA kit according to the manufacturer's instruction. Absorbance was measured at the wavelength of 405-420 nm by a microplate reader.

Cytotoxicity Assay
The effect of oxyresveratrol and its derivatives on the viability of normal human fibroblasts (MRC-5 cells) were determined by MTT assay. MRC-5 cells were diluted to 2 × 10 5 cells/mL with a complete medium. The cells (100 µL) were seeded in the wells of the sterile 96-well cell culture plate. After incubation at 37 • C in CO 2 for 24 h, the cells were treated with different concentrations of sample compounds (100 µL) and doxorubicin as a positive control (100 µL). After incubation for 48 h, MTT solution (20 µL of 5 mg/mL) was added to each well and incubated for 3 h. After that, the complete medium was removed. Then, the cells were lysed with DMSO (100 µL). The yellow MTT dye was reduced by succinic dehydrogenase in the mitochondria of viable cells to purple formazan crystals. After that, absorbance (OD) was measured by a microplate reader at 570 nm. The percentage of cytotoxicity was determined with Equation (1) below [78].

Preparation of Enzyme Receptor
The x-ray crystallographic structure (PDB entry: 3LN1, chain A) of COX-2 with the resolution of 2.4 Å was retrived from the RCSB protein data bank. Cofactors and hemoglobin were removed from the protein-ligand complex. Missing atoms of amino acids were fixed by using Discovery Studio 4.1 [79]. The nitrogen-and carbon-terminal ends were, respectively, capped with acetoxy and methyl amino groups. The protonation state of amino acids was set at pH of 7.4 by using the PROPKA program [80]. To remove any high energy contacts and relax some of the distorted geometries, restrained molecular dynamics (MD) simulation [81] of chain A of COX-2 with bound celecoxib and 37 crystallized water molecules was carried out for one ns using GROMACS package version 2016.4 [82]. Atomic partial charges of bound celecoxib were detemined by using restrained electrostatic potential (RESP). Its topology for MD simulation was generated using a GAFF force field [83] and AnteChamber PYthon Parser interfacE (ACPYPE) [84], a wrapper script around the ANTECHAMBER package. By using the TIP3P model, the total number of 29,447 water molecules was generated and explicitly included in the system initially neutralized by 0.15 M NaCl.

Computational Models of Ligand Structures
Geometries of oxyresveratrol and its derivatives were built and optimized by using Gaussian 09 program [85] with M06-2X density functional [64] and 6-31G(d) basis set. For celecoxib, its input coordinate was extracted from the MD trajectories and added to the covalently bonded hydrogen atoms.

Molecular Docking Protocols
The automatic docking scheme was applied to the diverse conformations of the oxyresveratrol derivatives in the active site of the COX-2 receptor by using the AutoDock 4.2.5 [56] and GOLD Suite 2020.3 [58]. The size of a grid box was 40 × 40 × 40 Å with the dimensions of 0.375 Å grid spacing. The GOLD docking protocol [86], as previously described, was carried out by using the empirical GoldScore [57] and ChemPLP [59] scoring fitness functions. These two scores of the best-docked pose were examined and compared with the binding energy computed by using AutoDock. The docking performance of each scheme was validated against the crystallographic binding mode of COX-2/celecoxib prior to docking all oxyresveratrol analogues. Protein-ligand interactions were displayed using Discovery Studio 4.1 [79], Chimera 1.11.3 [87] and LigPlot+ v.2.2.4 [88]. The docked ligandprotein complexes predicted by GOLD Suite were then used for quantum mechanical (QM) optimization and calculations.

Computation of QM Energy and Binding Energy (BE)
Five water molecules and twenty amino acid residues: Val74, His75, Arg106, Phe191, Val330, Val335, Leu338, Ser339, Tyr341, Tyr371, Trp373, Arg499, Ala502, Phe504, Val509, Glu510, Ala513, Ser516, Leu517, Leu520 of the COX-2 active site with each of the bestdocked ligand poses obtained from the GOLD Suite was quantum-mechanically optimized using M06-2X [64,89] density functional theory (DFT). All dangling bonds were capped with the hydrogen atom. For the complex optimization, we applied the 6-31G(d) basis set to the binding region containing bound ligand and key amino acids: Arg106, Val335, Tyr341, Tyr371, Arg499, and Val509, as well as five water molecules. The lower basis set of 3-21G was applied for the rest. To confirm the binding stability, quantum mechanical (QM) energy of the optimized complexes of COX-2 with bound ligands was computed using M06-2X density functional and 6-31+G(d,p) basis set for both in vacuo and implicit aqueous solvent. The binding energy (BE) was calculated by Equation (2). The basis set superposition error on the calculated BE was also compensated by counterpoise (CP) correction on the optimized complex geometries. We performed all QM calculations using the Gaussian 09 package [85] with the ultrafine integration grid.

Statistical Analysis
A one-way analysis of variance was carried out for the docking scores and the estimated binding affinity as a dependent variable. The interaction of a dependent variable and different scoring functions were analyzed using a univariate general linear model with Dunnett's multiple comparison at the alpha level of 0.05. The dependence of QM energetics on different computational conditions was also investigated in the same fashion. We performed all statistical analysis by using IBM SPSS Statistics version 26 [90].
Supplementary Materials: The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/molecules27072346/s1, Figure S1: X-ray crystallographic pose of celecoxib bound to COX-2 receptor (PDB entry: 3LN1). Right panel shows the co-crystallized celecoxib and key residues of amino acids at the catalytic site of COX-2; Figure S2: The superposition of the one ns MD structure (grey) and the best-docked poses (yellow) of celecoxib from (a) AutoDock with the RMSD value of 0.68 Å, (b) ChemPLP with the RMSD value of 0.55 Å and (c) GoldScore with the RMSD value of 0.56 Å was within 6 Å of COX-2 binding pocket. Only the amino acid residues in proximity to the binding pocket are highlighted for clarity; Figure S3: Overlay of the lowest-QM-energy poses (stick representation) of celecoxib and oxyresveratrol derivatives in catalytic pocket of COX-2; Figure S4: Distribution of Mulliken atomic charges of unbound pose of (a) celecoxib, (b) oxyresveratrol (1), (c) dihydrooxyresveratrol (4) and (d) difluorooxyresveratrol (7); Figure S5: Distribution of Mulliken atomic charges of bound pose of (a) celecoxib, (b) oxyresveratrol (1), (c) dihydrooxyresveratrol (4) and (d) difluorooxyresveratrol (7) with COX-2; Table S1: Hydrophobic interactions of relaxed pose of celecoxib obtained from MD simulation in the catalytic pocket of COX-2; Table S2: Dunnett t-test of docking scores; Table S3: Dunnett t-test of QM interaction energy; Figure S6