Next Article in Journal
Multiple Targets Directed Multiple Ligands: An In Silico and In Vitro Approach to Evaluating the Effect of Triphala on Angiogenesis
Next Article in Special Issue
Phytochemicals in Garlic Extract Inhibit Therapeutic Enzyme DPP-4 and Induce Skeletal Muscle Cell Proliferation: A Possible Mechanism of Action to Benefit the Treatment of Diabetes Mellitus
Previous Article in Journal
The Guanidine Pseudoalkaloids 10-Methoxy-Leonurine and Leonurine Act as Competitive Inhibitors of Tyrosinase
Previous Article in Special Issue
Alpha-Amylase and Alpha-Glucosidase Enzyme Inhibition and Antioxidant Potential of 3-Oxolupenal and Katononic Acid Isolated from Nuxia oppositifolia
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Discrimination of Naturally-Occurring 2-Arylbenzofurans as Cyclooxygenase-2 Inhibitors: Insights into the Binding Mode and Enzymatic Inhibitory Activity

by
Ericsson Coy-Barrera
Bioorganic Chemistry Laboratory, Facultad de Ciencias Básicas y Aplicadas, Universidad Militar Nueva Granada, Cajicá 250247, Colombia
Biomolecules 2020, 10(2), 176; https://doi.org/10.3390/biom10020176
Submission received: 9 December 2019 / Revised: 12 January 2020 / Accepted: 20 January 2020 / Published: 23 January 2020

Abstract

:
2-arylbenzofuran-containing compounds are chemical entities that can be naturally produced by several organisms. A wide-range of activities is described for several compounds of this kind and they are, therefore, valuable moieties for a lead finding from nature. Although there are in-vitro data about the activity of 2-arylbenzofuran-related compounds against cyclooxygenase (COX) enzymes, the molecular level of these COX-inhibiting constituents had not been deeply explored. Thus, 58 2-arylbenzofurans were initially screened through molecular docking within the active site of nine COX-2 crystal structures. The resulting docking scores were statistically analyzed and good reproducibility and convergence were found to discriminate the best-docked compounds. Discriminated compounds exhibited the best performance in molecular dynamics simulations as well as the most-favorable binding energies and the lowest in-vitro IC50 values for COX-2 inhibition. A three-dimensional quantitative activity-structure relationship (3D-QSAR) was also demonstrated, which showed some crucial structural requirements for enhanced enzyme inhibition. Therefore, four hits are proposed as lead structures for the development of COX-2 inhibitors based on 2-arylbenzofurans in further studies.

Graphical Abstract

1. Introduction

Cyclooxygenase, COX, (or Prostaglandin H2 synthase), is a well-known enzyme that catalyzes the rate-limiting step in the formation of prostanoids such as prostaglandins (PGs) and thromboxane A2 (TxA2). These products are considered as important mediators of fever, pain, and inflammation [1]. Three COX isozymes are known such as the constitutive form including COX-1 (which exists in several tissues), the inducible form, COX-2 (which is expressed throughout inflammation), and the COX-1 splice variant, COX-3 [1]. In recent decades, the fact for understanding the distinct functions of COX isozymes has encouraged many studies, including those dedicated to the search for anti-inflammatory agents. In this sense, it is well-recognized that selective COX-2 inhibitors have been developed as securer substitutes for the nonsteroidal anti-inflammatory drugs (NSAIDs), mainly due to the adverse effects caused by the inhibition of the COX-1 enzyme. Therefore, a pronounced interest for finding selective anti-COX-2 agents still remains [2].
Several types of compounds (natural or synthetic) with anti-COX-2 properties have been discovered in recent years. Several chemical entities have been studied as selective COX-2 inhibitors with interesting and important potential [3]. However, the most important ones are those named coxibs (such as celecoxib, rofecoxib, and valdecoxib), which have been approved due to the lower adverse effects in comparison to the traditional NSAIDS (non-steroidal anti-inflammatory drugs) [3]. However, it has recently been reported that some coxibs and other selective COX-2 inhibitors cause non-favorable gastrointestinal and cardiovascular events [3,4]. Because of that, novel agents for COX-2 inhibition are still required.
A huge interest on benzofuran-containing compounds has been originated several years ago because of their important biological and pharmacological properties and their extensive natural occurrence [5]. Benzofuran moiety is considered to be an indole isostere (a potent pharmacodynamic moiety exhibiting anti-inflammatory activity) and, consequently, several anti-inflammatory benzofuran-based compounds have been studied [6,7]. However, despite the plethora of reported benzofuran derivatives, 2-arylbenzofurans have not been sufficiently explored at a molecular level for COX-2 inhibition in spite of some of them exhibiting selective in-vitro anti-COX-2 activity [8,9]. Thus, as part of the research on anti-inflammatory compounds of natural origin, a combined in-silico study comprising different computational approaches/tools were used on the basis of three motivations: (1) discriminating putative COX-2 binders through a molecular docking-based screening of 2-arylbenzofuran-related compounds, within the active site of nine reported crystal structures of COX-2, (2) describing some insights into the binding mode of the best-docked compounds using molecular dynamics and binding energy calculations, and (3) recognizing some structural requirements of the best-docked compounds toward COX-2 inhibition by means of comparative molecular field analysis (CoMFA). This last purpose required the expansion of the in-silico predictions to experimental results through the in-vitro COX-2 inhibitory assay performed on an in-house collection of 2-arylbenzofurans. Results revealed four very important hits to be used in future studies for developing novel selective COX-2 inhibitors based on 2-arylbenzofurans.

2. Materials and Methods

2.1. Protein Preparation

The X-ray crystallographic structure of nine COX-2 proteins (from Mus musculus) were obtained from the Protein Data Bank (PDB) (http://www.rcsb.org) at a resolution < 3.00 Å (Table 1). Water molecules, ligands, and other heteroatoms were removed from the protein molecules. The addition of hydrogen atoms to the proteins was then performed. These crystallographic structures were retained without any processing for molecular docking. The co-crystallized inhibitors/substrates were employed to define the corresponding active site, the flexible residues within this active site, and a validation criterion of docking calculations (re-docking). By means of a literature survey, a maximum of five common residues involved in the binding site were defined as flexible residues: Arg106, Val335, Ser339, Tyr341, and Ser530. PDB files of proteins were effectively prepared using the AutoDock/Vina (Molecular Graphics Lab, The Scripps Research Institute, La Jolla, CA, USA) plugin under PyMOL 2.3 (Schrödinger LLC, New York, NY, USA), and saved as pdbqt files. COX-2 structures were also superimposed by sequence alignment for calculating the root-mean-square deviation (RMSD) and sequence similarity using Discovery Studio Client 16.1 (Biovia, San Diego, CA, USA).

2.2. Ligand Preparation

A targeted, custom-made library comprising only 2-arylbenzofurans (n = 58) was compiled from the information of an in-house collection of isolated compounds as well as different published works reporting anti-inflammatory activity of structurally-related 2-arylbenzofurans of natural origin. Compiled compounds were three-dimensionally (3D)-sketched using the molecular builder and visualization tool Avogadro 1.2.0 (open-source, http://avogadro.cc). A Monte-Carlo randomized conformational search was then performed, without any geometrical restrictions, using the semi-empirical AM1 parametrization implemented in the Spartan’14 software (Wavefunction, Irvine, CA, USA) with a limit of 500 conformers. Energetically-lowest stable conformers, within a 6 kcal/mol energy range, were geometrically optimized at a density functional theory (DFT) level using B3LYP hybrid functional and 6–311++G (2df, 2p) basis set that is also implemented in Spartan’14. After structural optimization, harmonic vibrational frequencies were calculated at the same level of theory to verify the reliability of the stationary point at a minimum and all structures converged successfully. During geometry optimization, the highest energy occupied molecular orbital (HOMO) and lowest energy unoccupied molecular orbital (LUMO) were also computed for further analysis. The verified, optimized ligand structures were saved as PDB-files for further molecular docking and molecular dynamics analysis. Additionally, structural similarity analysis of 2-arylbenzofurans was performed using the substructure fragment dictionary-based binary fingerprint descriptor (FragFp) implemented in the open-source data visualization and analysis program DataWarrior [17]. The absorption, distribution, metabolism, and excretion–toxicity (ADMET) properties of test compounds were computationally predicted using the Free ADMET Filtering tool (FAF-Drugs4, http://fafdrugs4.mti.univ-paris-diderot.fr).

2.3. Molecular Docking

Molecular docking calculations were performed by means of Autodock/Vina (Molecular Graphics Lab, The Scripps Research Institute, La Jolla, CA, USA) using PyMOL 2.3 (Schrödinger LLC, New York, NY, USA) and running on a Ubuntu 12.04 server equipped with a dual Intel Xeon® processor CPU @ 2.6 GHz (32 CPU), 64 GB DDR3 RAM. Prepared proteins and ligands were used for simulations using flexible residues. The selection of such flexible residues within the active site was based on a position at 4.0 Å from co-crystallized ligand (i.e., substrate or inhibitor, as presented in Table 1 and Figure A1). Thus, the active site was placed into a 24 × 24 × 24 Å cubic box at the geometrical center of the set of selected flexible residues, which has 0.375 Å as grid point spacing. Ten different pdbqt-files of each ligand were used as replicates to be docked into the flexible active site of prepared enzymes to evaluate convergence. Co-crystallized ligands were also docked in order to assess the docking performance through a re-docking strategy. Additionally, the specificity and sensitivity of the docking protocol, 15 compounds of diverse chemical nature, which have IC50 < 4 μM, against COX-2 (3LN1) were compiled from ChEMBL [18]. Furthermore, 50 decoys per each active compound were compiled from the directory of useful decoys—enhanced (DUD-E) [19]. Thus, decoys and bioactive compounds were processed using the same docking protocol. The resulting data was then assessed through a receiver operating characteristic (ROC) and scores enrichment using the screening explorer webserver [20]. Once the docking parameters were validated, the docking simulations of test 2-arylbenzofurans were performed. Docking results were initially examined in PyMOL 2.3 and the resulting poses were top-ranked regarding the Vina scores (kcal/mol). Molegro Virtual Docker 6.0 (CLC Bio Company, Aarhus, Denmark) was secondly used to evaluate the best-poses of best-docked compounds as a rescoring strategy. RMSD of atomic positions of the best poses between replicates was calculated through an in-house bash script. Discovery Studio Client 16.1 (Biovia, San Diego, CA, USA) and Maestro 11.8 (Schrödinger, Cambridge, MA, USA) were used to visualize the respective 3D models and 2D residual interactions diagrams, respectively, using the pdbqt outputs for the highest-scored ligands.

2.4. Statistical Analysis

Descriptive and inferential statistics tests were carried out using R project software version 3.0.2 (R Foundation, Vienna, Austria). Principal component analysis (PCA) was also carried out in SIMCA 14.0 software (Umetrics, Umeå, Sweden) using the docking scores in order to observe plausible relationships with structural variations among test 2-arylbenzofurans.

2.5. Molecular Dynamics Simulations

The ligand’s best poses from selected docked (ligands) and crystal (COX-2, 3LN1) structures were chosen as inputs for molecular dynamics (MD) simulations of ligand-enzyme complexes and the apoenzyme. Such simulations were carried out on the systems comprised of COX-2+celecoxib, COX-2 + 11, COX-2 + 29, COX-2 + 42, and COX-2 + 49 (obtained after molecular docking). These systems were individually subjected to an additional 150-ns conventional MD simulation for exploring the performance over the time of these ligands within the active site of the enzyme. These MD simulations were run in Gromacs 5.0.5 (open source, http://www.gromacs.org) on an Ubuntu 12.04 server, using NPT (constant pressure and temperature) and periodic boundary conditions, as previously reported [21]. Hence, docked ligands were prepared by adding hydrogen atoms in UCSF Chimera (UCSF, CA, USA) and the resulting pdb-file was uploaded to the atb server (http://compbio.biosci.uq.edu.au/atb/) to add the respective Gromo53a6 force field and get the itp-type topology file. COX-2 topologies were obtained in Gromacs using the Gromos53a6 force fields, due to the presence of the Heme group. The simple point charge (SPC) water model was then implemented for solvation in a triclinic box using a 1.0-nm margin distance. 0.10 M NaCl was added to the simulation systems and water molecules were randomly replaced until neutrality. An energy minimization through a 2000-steps steepest descent method was then used. NVT (constant volume and temperature) equilibration at 310 K during 500 ps, followed by NPT equilibration during 1000 ps using the Parrinello-Rahman method at 1 bar as a reference were done on the systems using position restraints. Lastly, solute position restraints were released and a production run for 90 ns was performed. Temperature and pressure were kept constant at 310 K and 1 bar, respectively. Coordinates were recorded in a 1 fs time step. Electrostatic forces were calculated using the particle-mesh Ewald (PME) method. Periodic boundary conditions were used in all simulations and covalent bond lengths were constrained by the LINear Constraint Solver (LINCS) algorithm.

2.6. Binding Free Energy Analyses

Binding-free energy was calculated using the g_mmpbsa tool [22] (open source drug discovery consortium (OSDD), New Dehli, India). This tool calculated components of the free energy of the protein–substrate binding (ΔGBind) using the molecular mechanic/poisson-boltzmann surface area (MM/PBSA) method [23,24]. In this method, ΔGbind calculation between a protein and a ligand is carried out by ∆Gbind = ∆H − T∆S ≈ ∆EMM + ∆Gsol − T∆S, ∆EMM = ∆Einternal + ∆Eelectrostatic − ∆EvdW, ∆Gsolv = ∆Gelecsolv + DGvdWsolv, where the total gas phase energy on the binding of MM energy is ΔEMM, the free energy of solvation is ΔGsolv, and the entropy contribution is TΔS. A Poisson-Boltzmann model was used to compute the electrostatic solvation energy in a continuum solvent. The derivation of a non-polar solvation energy term was computed as a solvent-accessible surface area (SASA). ΔEMM were calculated using the Lennard-Jones and Coulomb potential [23,24]. ΔGBind was used to analyze the binding associations between COX-2 and selected ligands (i.e., celecoxib, 11, 29, 42, and 49) by decomposing the total binding-free energy into each residue. The binding energy calculations of the selected ligands were performed for 100 snapshots taken at an interval of 500 ps during the last stable 40-ns MD simulations.

2.7. COX-2 Inhibition Assay

The ability of selected compounds to inhibit bovine COX-2 was determined using an enzyme immunoassay (EIA) kit (Cayman Chemical Company, Ann Arbor, MI, USA), according to the reported method [8]. The activity was expressed as half-maximal inhibitory concentrations, IC50 (in μM), obtained after non-linear regression, using the program GraphPad Prism version 5.00 (GraphPad Software, San Diego, CA, USA).

2.8. Comparative Molecular Field Analysis (CoMFA)

The best-docked poses of those test 2-arylbenzofurans experimentally evaluated against COX-2 were merged and aligned by means of particular tethers placed on benzofuran and aromatic rings moieties, using the molecular overlay tool included in the software Discovery Studio Client 16.1 (Biovia, San Diego, CA, USA). Compound 11 was selected as a template to select tethers. The aligned molecules set (in an sdf-file) was randomly divided into two subsets (training and test sets, corresponding to 70% and 30%, respectively). Comparative molecular field analysis (CoMFA) analysis was then performed using Open3DQSAR, using the standard protocol [25]. The experimental in-vitro COX-2 inhibition (expressed as IC50 in M) was converted into a negative logarithmic form (pIC50) and then used as an independent variable. The models were validated by leave-one-out (LOO) and leave-many-out (LMO) methods. The quality of the models after validation was evaluated on predicting the independent variable for the test set.

3. Results and Discussion

3.1. 2-Arylbenzofuran-Related Chemical Space

The selection of test 2-arylbenzofuran-type compounds was based on three parameters: (1) natural origin, (2) enlarging the chemical space of an in-house collection (n = 26) of 2-arylbenzofurans, and (3) having an anti-inflammatory effect (preferably, inhibition of enzymatic activity or expression of COX-2). Hence, a targeted, custom-made library, comprising 58 2-arylbenzofurans, was then compiled. Common names of these compounds (158) are found in Table A1 and their structures are presented in Figure A2. The 2-arylbenzofuran-related chemical space was examined to characterize it and facilitate the data analysis according to structural fragments. Therefore, a structural similarity analysis between test compounds was performed using the FragFp descriptor available in DataWarrior software [17]. Compound 11 was selected as a reference to compute the similarity index due to its good COX-2 inhibitory activity, as previously reported [8]. The resulting similarity chart is presented in Figure 1.
After such an analysis, some clusters were shown, which involved similarity FragFp indexes between 0.93 and 0.28, according to the heatmap based on the structure similarity index between 0 (red) and 1 (light blue) (Figure 1). Clusters were related to different classes of substituted benzofurans, which were then subdivided into moracin-type (including isopentyl-substituted and dihydrooxepine-containing) (FragFp 0.93–0.69), 2H-pyran- and benzofuro[6.5-b]furan-containing (FragFp 0.67–0.38), propyl(en)-susbtituted (FragFp 0.90–0.71), dihydropyran-containing (FragFp 0.66–0.62), furocyclohexadienone-type (FragFp 0.28), licarin-type (i.e., 2-aryldihydrobenzofuran) (FragFp 0.50–0.46), and other four ungrouped compounds (FragFp 0.67–0.59). This structure similarity analysis clustered quantitatively the test compounds into six classes, which indicated these compounds have particular moieties/fragments to have further analysis.

3.2. Molecular Docking Simulations

Molecular docking was chosen as a first-line discrimination of bioactives from the above-mentioned chemical space related to 2-arylbenzofuran compounds due to the capacity to simulate the binding of small compounds within the active site of enzymes, as a structure-based discrimination procedure for the in-silico prediction of putative competitors [26,27,28]. In this regard, some concerns should be considered using accurate parameters for suitable selection of sampling and scoring procedures during structure-based screening [29,30]. Therefore, the molecular docking of the co-crystallized ligands for each enzyme structure (ES1–ES9, Table 1) indicated a good performance of the present docking protocol since these re-docking calculations resulted into low conformational RMSD values (<0.5 Å), such as the docked/co-crystallized superposition of celecoxib presented in Figure 2a. Furthermore, the docking performance was also evaluated through a benchmarking approach from docking results of 15 known selective COX-2 inhibitors in comparison to the binding scores of a group of decoys (n = 750). Calculation of the area under the curve (AUC) from ROC curves was used to estimate the sensitivity and specificity of the docking protocol [31]. AUC and the Boltzmann-enhanced discrimination of the receiver operating characteristic (BEDROC) resulted in 0.908 and 0.827, respectively. Therefore, the validation of the docking protocol was considered successful (Figure 2b). This evaluation reached true positive identification over 92% of the active compounds involving a recognition within 10% of the test compounds (Figure 2c).
Once the parameters were adequately validated, each geometrically-optimized structure 158 was docked into the active site of each COX-2 structures (ES1–ES9, Table 1) using AutoDock/Vina. These enzyme structures ES1–ES9 were retrieved from the PDB and their respective data are shown in Table 1. These M. musculus-derived COX-2 structures were selected to compare the performance of selective (S) and non-selective (NS) COX-2 inhibitors and 2-arylbenzofurans, as a discriminating initiative using the re-docking strategy, to assess the predictive potential of the simulation protocol. It was, therefore, required to retrieve COX-2 structures bound to common anti-inflammatory drugs (i.e., S/NS COX-2 inhibitors) and the natural ligand. In addition, in spite of the importance to include human COX-2, they were not used in the present docking study. Other MmCOX-2 PBD structures were discarded due to some issues related to sequence, residues, and a bound ligand. Such facts were taken as criterion to select the test COX-2 structures. Furthermore, the active site was adopted according to the information described in the referenced studies. After a detailed scrutiny of such information, five common residues involved in the binding site were initially defined (i.e., Arg106, Val335, Ser339, Tyr341, and Ser530). They were then used in the present docking study as flexible residues as well as those located at 4.0 Å from the co-crystallized ligand. These structures used for docking consist of very similar COX-2 protein sequences in a different solution when the enzyme was co-crystallized with a ligand or unbounded. These structures ES1–ES9 presented small alterations on the tertiary structure depending on the crystallization protocol. In addition, the presence of the inhibitor bound to the active site can also induce changes in the 3D-conformation of the binding site and modify the accessibility of some residues for ligand interaction. These facts were examined on superimposing the 9 chain-A 3D-structures of COX-2 after sequence alignment (Figure 3), which exhibited slight differences in a primary and tertiary structure (i.e., sequence similarity 99.2%–99.8% and RMSD 0.35–0.50 Å (Figure 3a), respectively). The differences between these COX-2 sequences were found to be related to particular residues such as Asn34, Ala33, His90, His208, Gln310, and Lys333 (Figure 3b), but the active site remained identical. Therefore, the same ligands docked into the same binding site of different crystal structures of the same enzyme is very important to statistically support the docking results.

3.2.1. Performance and Trends of the Docking Results

Docking results were examined by descriptive statistics to explore the trends and performance of the docking protocol, on one hand, and visualize and rationally select the best-docked compounds to discriminate putative binders of COX-2 on the other. After 10 replicates per calculation of benzofuran-related compounds 158 with nine crystal structures of COX-2 (ES1–ES9), the docking results were expressed as Vina Scores (in kcal/mol) and the best pose into the active site of each enzyme structure was also achieved. The mean docking scores, the standard deviation (SD) among replicates, and the RMSD (in Å) of atomic positions between the best poses after 10 replicates for each test benzofuran docked with each COX-2 structure are listed in Table A2. In general, the SD values were found to be within the 0.02–0.55 range, which indicates good reproducibility for docking calculations. However, the RMSD values exhibit high variability (i.e., 0.049–1.582 Å range). However, the most RMSD values were below 0.8 Å, which demonstrated good convergence in the calculated complexes formation. As expected, arachidonic acid (AA) exhibited the highest variability due to the long chain-derived conformations (14 rotatable bonds), but planar (more rigid) structures exhibited the lowest SD and RMSD values. 2-arylbenzofurans possessing an additional ring (e.g., pyran, furan, or oxepine) exhibited overall stronger docking scores (<−10 kj/mol). The previously mentioned facts were confirmed with the global mean Vina scores (GMVS) for the docking results with all COX-2 structures (Table 2).
Additionally, docking scores were found to be adequately distributed and reproducible. The relative standard deviation (RSD) percentages resulted in the 2.7%–33.2% range, but only 16.3% of the RSD values were above 10% and 3.6% of the RSD values above 20%. These facts are a good indication of the reproducibility. It is well-recognized that the docking procedure has a limitation related to distinguishing the correct pose among the generated scoring positions. However, the difference in the best poses among replicates were considered to be convergent since 76.3% of the RMSD values were below 0.6 Å. A similar trend was observed with inhibitors and the respective co-crystallized structure after re-docking calculations (RMSD < 0.5 Å).
Regarding the behavior of Vina scores among COX-2 structures (ES1–ES9), a similar distribution was found as shown in the respective boxplot (Figure 4a), which was built from the resulting mean docking scores. Mean values per COX-2 structure were similar for ES1–ES5, as a first group, and ES6, ES8, and ES9 as a second one. ES7 exhibited the lowest mean Vina score among enzymes and the lowest docking score (−12.7 ± 0.2 kcal/mol) corresponding to celecoxib (cel). These facts mean that the particular behavior of the ligands into each binding pocket was slightly different, which can be understood as a consequence of structural adaptations of the active site by the presence of the co-crystallized ligands. This issue should be considered a very important factor in the moment to get more reliable results during docking studies.
Trends within docking results were also inspected by multivariate statistics in order to identify the plausible relationship between docking scores and structural variations of 2-arylbenzofurans. A PCA was performed using the Vina scores dataset. The first two components explained the 82.3% of the total variance with good model fitting (R2 = 0.823 and Q2 = 0.778 for the principal component 1 (PC1) and PC2). The resulting score plot (PC1 vs. PC2, Figure 4b) demonstrated good differentiation of the compounds by variance according to the Vina scores along COX-2 structures. Hence, this unsupervised analysis clustered the compounds in four groups, which is shown by discriminating across PC1 depending on mean Vina scores (66.2% explained variance). Group 1 (green) consisted of the best-docking 2-arylbenzofurans (such as 11, 28, 29, 30, 40, 42, 49, 51, and 54) with a comparable profile to the selective COX-2 inhibitors celecoxib (cel) and SC558. 2-arylbenzofurans possessing an additional cycle (such as pyran, furan, and oxepine) were the predominant structural feature in group 1. In contrast, group 2 (dark blue) clustered the weakest-docking 2-arylbenzofurans with profiles related to those of arachidonic acid (AA), indomethacin, and naproxen. This group was characterized by unsubstituted moracin-type compounds as well as some dihydrobenzofurans. Groups 3–4 (yellow and red) exhibited their separation by dual incidence of PC1 and PC2, whose separation along PC2 was mainly dependent by the scoring behavior with ES6 and ES8. These results indicated that the PCA led to discriminated bioactive compounds, according to the molecular docking performance.

3.2.2. Selection of the Best-Docked Compounds

An examination of the previously mentioned docking results for compounds 158 (Table 2) was performed to select the best-docked compounds as putative COX-2 binders. Hence, from the global mean Vina scores (GMVS < −10.5 kj/mol)), best-pose convergence (RMSD < 0.5) and scoring reproducibility (%RSD <6%) along nine COX-2 structures, the compounds 9′-nor-7, 8-dehydro-isolicarin-B (11), moracin H (29) psoralidin (42), and maximol (49) were established as the best-docked compounds among the test set. These selected compounds were, therefore, clustered in the green group of the PCA-derived score plot (Figure 4b), which exhibits the lowest-scoring behavior. Other compounds in this group were not selected because they did not meet completely the previously mentioned features (i.e., GMVS, RMSD, and %RSD). In addition, Molegro Virtual Docker (MVD) was also used to simulate the intermolecular interaction as a rescoring strategy. An identical trend was obtained, so the best-docked compounds used in Vina exhibited the best MolDock scores (i.e., cel = −163.2, 49 = −156.9, 42 = −156.6, 11 = −139.1, 29 = −136.2), which is another indication of the good performance of the docking protocol and the filtering by descriptive and multivariate statistics. These four compounds were used for the first time to simulate their interaction within the active site of COX-2 using a validated docking protocol. However, compound 11 was previously evaluated against COX-2 enzyme (affording an IC50 of 3.32 μM) [8] and compound 42 was reported to exhibit good ability to reduce the COX-2 expression at 50 µM [32].

3.2.3. Binding Mode and Residual Interactions of the Best-Docked Compounds

Binding modes of the best-docked 2-arylbenzofurans with the COX-2 structure ES7 (PDB code: 3LN1) was examined looking for important/crucial interactions between active site residues and ligand moieties, according to the individual 2D-residual interaction diagrams (Figure 5).
Analysis of these diagrams was useful to delineate some insights into the binding mode of the simulated complexes, which showed some important hydrophobic and polar interactions. These interactions were taken as key contacts and an indication of the importance of the presence of the hydroxiaryl moiety for interacting with the respective active site residues to stabilize the simulated enzyme-ligand complexes. Particularly, celecoxib exhibited H-bonds between Phe504, Arg499 (as donors), and Gln178 (as acceptor) and sulfonamide as well as Arg106 (as donor) with 1H-pyrazol-2-ium moiety. In contrast, other chemical interactions were found on analyzing the binding modes of the lowest-scored poses for the four strongest-docking compounds (i.e., 11, 29, 42, and 49; Figure 5). These compounds resulted in the lowest mean docking scores (−10.48 to −11.26 kcal/mol range) and were found to be well-positioned into the COX-2 active site. However, they involve different orientations depending on the type of interactions. In this regard, residues Tyr341, Gln178, and Ser339 were found to be key polar contacts to stabilize the enzyme ligand complex for those compounds possessing a resorcinyl substituent, through a π-π-stacking for 29 and acting as an H-donor for 49. Furthermore, the p-hydroxyphenyl group in 49 interacted with Trp373 via π-π-stacking. Psoralidin (42) exhibited another binding mode involving His75 as an H-donor to its phenolic OH at benzofuran moiety. Ala513 also interacted with phenolic OH at aryl moiety in 42. In addition, the methylenedioxy moiety in 11 was important to orient the molecule toward Arg499, acting as an H-donor acceptor. The combination of these structural features of such benzofurans could serve as an important starting point to design a novel series of COX-2 inhibitors.

3.3. Stability and Chemical Reactivity of the Best-Docked 2-Arylbenzofurans

Frontier molecular orbitals HOMO and LUMO can reflect the successful biological interaction of ligands within the binding pocket of proteins [33]. In this sense, the DFT level-derived calculations of HOMO and LUMO energies were a suggestion of kinetic stability and chemical reactivity of these best docked compounds. The resulting B3LYP/6–311++G (2df, 2p) models are presented in Figure 6. These compounds exhibited lower HOMO-LUMO energy gaps to that of celecoxib and other 2-arylbenzofurans, which infer high stability but also high reactivity. Pearson’s correlation between COX-2 docking scores and HOMO-LUMO energy gaps for 2-arylbenzofurans 158 was also evaluated. A positive correlation among scores and energy gaps was found (0.812, p < 0.001), which indicated that the self-reactivity of compounds might influence the binding mode. Hence, among the best-docked compounds, 11 and 29 exhibited higher Vina scores and HOMO-LUMO energy gaps, whereas 42 and 49 showed an opposite behavior (Figure 6). This relationship can be rationalized since these energy gaps and eigenvalues defined the biological activity of ligands via high reactivity and good stability. Therefore, enzyme inhibition by these best-docked 2-arylbenzofurans might be afforded through an electron donating ability of test compounds with the enzyme-binding site pocket of COX-2.

3.4. Molecular Dynamics Studies of the Best-Docked 2-Arylbenzofurans

In order to extend the information on the binding mode of the best docked compounds, 150-ns MD simulations were separately accomplished using the COX-2 (ES7) alone (apoenzyme) and docked distinctly with compounds 11, 29, 42, 49, and celecoxib. Ligand-enzyme trajectories for resulting complexes were examined through the variation of geometric properties over the time. Thus, RMSD of the COX-2 backbone reflected the receptor frame stability by computing the time-dependent distance (Å) among different positions of the atom set (Figure 7a). Apoenzyme exhibited a normal evolution during the simulation but revealed a slight perturbation at 30 and 50 ns (RMSD 0.50–0.65 nm). A normal stabilization was then reached throughout the rest of the MD simulation (RMSD 0.63 nm). The COX-2·celecoxib complex evolved typically but showed a late perturbation at 80 ns (until RMSD 0.99 nm). Although compounds 11 and 49 exhibited the most-perturbed profile, both complexes (with celecoxib and 11) attained stability after 90 ns and 49 after 110 ns. The steady progress was maintained over the remaining time for these three ligands. On the other hand, compounds 29 and 42 exhibited the least-perturbed MD performance and the evolution was found to be similar to that of apoenzyme, which indicates good geometric properties on interacting with COX-2. The computed variations in RMSD values for these complexes showed that the COX-2 structure is affected differentially by the interaction with each ligand, which achieved a stable condition at the end of each MD simulation.
Fluctuations of the Cα atomic positions for each residue of the enzyme was explored through the root mean square fluctuation (RMSF). It was used to scrutinize the flexibility and secondary structure of the COX-2 enzyme under binding with the best-docked 2-arylbenzofurans. Thus, lower RMSF values implied constrained regions whereas higher RMSF values denoted more flexibility. All MD-simulated complexes showed similar performance (Figure 7b) involving fluctuations within the range of 0.2 and 1.6 nm. Such alterations were found to be near the formerly recognized crucial interactions (e.g., Arg106, Phe186, Leu292, Val330, Leu377, and Ser516). However, in comparison to the apoCOX-2, the binding with test ligands did not substantially modify the flexibility of these regions to maintain an overall steady state (RMSF differences > 0.2 nm), excepting the most-fluctuating residues around Arg106. Therefore, the COX-2 inhibition mode of compounds 11, 29, 42, and 49 might be interpreted through a complex stabilization, which is in agreement with the previously reported MD study for other known COX-2 inhibitors [34].

3.5. Binding-Free Energy Calculations of the Best-Docked 2-Arylbenzofurans

The binding-free energies (ΔGbind) for the best-docked compounds (i.e., 11, 29, 42, 49, and celecoxib) during the interaction of COX-2 for the last 40 ns of MD trajectory was estimated by the MM/PBSA approach to evaluate the global stability of the resulting ligand-enzyme complexes. The calculated binding energies are presented in Table 3. All five ligands exhibited negative binding energies. However, 42 showed comparable binding energy to that of celecoxib (−195.7 ± 4.5 kJ/mol vs. −197.8 ± 4.5 kJ/mol, respectively) and significant differences to that of other ligands, which justify the attained docking performance. The main contribution to the binding energy was due to van der Waal (ΔEvdW) energies (< −205 kJ/mol), which resulted in similar values among test compounds, while the contribution of electrostatic energy exhibited higher differences between them. Celecoxib and 49 exhibited the lowest electrostatic energy whereas 11 exhibited the highest one. In addition, the contribution of the polar solvation energy was found to be more unfavorable for 49 and celecoxib. The solvent accessible surface area (SASA) energy was similar for test ligands, even though compound 49 exhibited the highest ΔGsasa. Therefore, according to these results, electrostatic and polar solvation energy contributions could rationalize the difference in the binding mode of test ligands, which plays hydrophobic interactions as an important role in stabilization and even binding of test 2-arylbenzofurans within the active site of COX-2.
As a consequence, non-polar electrostatic interactions could be implied as the main driving force for the molecular recognition of COX-2 by test ligands. This fact was then confirmed after per-residue decomposition of the total binding energy of the simulated COX-2 ligand complexes. In this context, most binding energy contributing residues were found to be different depending on the ligand, as presented in Table 4. Regarding ligands, 42 exhibited similar binding energy to that of celecoxib (<−80 kJ/mol), which has the highest values among test ligands.
Regarding residues, Arg106, Leu338, and Val509 were found to be those residues that contribute the most (<−8 kJ/mol) to the binding energy for the COX-2·celecoxib complex (as described previously [34]), whereas Arg106 and Arg499 were observed for COX-2·11, Val335 for COX-2·29, Val335, Val509, and Ala513 for COX-2·42, and Glu510 and Ala513 for the COX-2·49 complex. Ser339 (<−2 kJ/mol) and Val509 (<−6 kJ/mol) were the common residues that contributed to the binding energy for all test ligands. Such residues with higher binding energy contributions represented those contact points for electrostatic interactions, as described in molecular docking and Cα atom fluctuations. Similar information was described for the binding mode of stilbene analogs within the active site of COX-2 [35].

3.6. In-Vitro COX-2 Activity of an in-House Collection of 2-Arylbenzofurans

The information on the potential as COX-2 inhibitors for the best-docked compounds (i.e., 11, 29, 42, 49) was expanded by measuring the IC50 (in μM) against COX-2 of an in-house library of 26 2-arylbenzofurans, which several of them were evaluated in a previous study [8]. The resulting IC50 values are presented in Table 5. Test compounds exhibited IC50 values in the range of 438 and 0.752 µM. In general, furocyclohexadienone and dihydrobenzofuran-containing compounds exhibited the highest IC50 values, except compound 49, whose IC50 was 1.25 µM, which can be attributed to the resorcinyl and styryl substitutions at C3 and C5, respectively. These moieties influenced substantially to the binding of 49 since some interactions were favored with hydrophobic (i.e., styryl) and polar (i.e., resorcinyl) regions of the active site of COX-2. No significant differences were found for the IC50 values of compounds 11 and 29, whose performances in molecular docking and molecular dynamics were very similar. The most active compound among the compound set was 42 (IC50 = 0.752 µM), which exhibited an interesting docking behavior (mean Vina score −11.25 ± 0.59 kcal/mol), the highest stability and reactivity (ΔEHOMO-LUMO = 3.77 eV), and an excellent performance over the time for stabilizing the COX-2·42 complex. This predicted compound was previously reported to exhibit anti-inflammatory properties as a dual inhibitor of the expression of COX-2 and 5-LOX [32]. These results can be understood as an adequate validation of the molecular docking-derived discrimination to find COX-2 inhibitors based on 2-arylbenzofurans. Previously mentioned most-active compounds also exhibited favorable ADMET properties (calculated using the FAF-Drugs4 web server [36]), which is afforded as accepted since the overall result with no Lipinski’s rule violation, good solubility, and oral bioavailability (Table A3).

Comparative Molecular Field Analysis (CoMFA)

A CoMFA was applied to correlate molecular interaction fields (MIF) and COX-2 inhibition in order to consistently rationalize the inhibitory activity of 2-arylbenzofurans regarding electrostatic and steric properties [38] as well as describing and predicting some structural requirements for COX-2 inhibition by 2-arylbenzofurans. For this purpose, the best-docked pose of compound 11 was used as a template to outline tethers and overlay the test compounds (Table 5), due to the structural features. Aligned structures and experimental data of COX-2 inhibition (as pIC50 values) were segmented into a training set (70%) to create the respective CoMFA model and a test set (30%) for the external validation [38]. MIFs were computed by means of electrostatic and steric probes. Partial least squares (PLS) regression was performed (employing up to five PLS components). Linear relationships between MIF fluctuations according to variations of the experimental pIC50 were accomplished. Thus, the best model achieved good correlation between MIF values and experimental pIC50 of test compounds using four PLS components, which comprises a correlation coefficient r2 = 0.816, a F-test = 162.321, cross-validated LOO coefficient q2LOO = 0.766, and a cross-validated LMO coefficient q2LMO = 0.701. These parameters were found to be enough for a statistically robust, predictive model [39]. The adequate correlation was examined through a Y-scrambling procedure [40]. Hence, 30 scramblings and 10 runs were employed. After that, no correlation appeared since the models dropped significantly (Rscr2 and Qscr2 < 0.3). This fact suggested the model was not reached as a consequence of a coincidental correlation. Therefore, this CoMFA-based model predicted the activity for the entire compound set, expressed as pIC50 (pred) (Table 5), which demonstrated a reasonable correlation for both training and test datasets.
Steric and electrostatic field outputs (stdev * coeff) reached values of 31.1% and 37.4%, respectively. Figure 8 showed the corresponding translated contour surfaces from the field contributions. The electrostatic field contour map (Figure 8a) displayed the positively and negatively-charged regions to favor COX-2 inhibition as blue and red contours, respectively. Accordingly, phenolic hydroxyl groups at C4′ and C7, resorcinyl moiety, and electron-withdrawing groups at C6 enhanced the COX-2 activity to interact with some residues located in the polar region of the active site of this enzyme (e.g., Hiz75, Arg499, and Gln178), as observed for compounds 29, 42, and 49. Alkyl (such as allyl, styryl, and prenyl) and methoxyl substitutions on 2-aryl moiety also favored the COX-2 inhibition (e.g., compounds 11 and 42) to interact positively or negatively with residues in the hydrophobic zone of the COX-2 active site (i.e., Val102, Leu103, and Leu345).
On the other side, the enzymatic inhibitory activity might also be enhanced/decreased by the presence of bulky groups at benzofuran and 2-aryl moieties, as observed in the steric field contour maps (Figure 8b). Unfavorable and favorable zones by steric effects were then depicted as yellow and green contours, respectively. COX-2 inhibitory activity could be improved if a bulky steric group is located on C3′, C5, and C7, such as substituted prenyl groups and dihydropyran ring, as observed in most active compounds 29 and 42. In contrast, a bulky group in C2′ (e.g., hydroxypropenyl) and C3 (e.g., phenyl group) negatively influenced the enzymatic inhibition. Consequently, the CoMFA model indicated that the particular presence of both bulky and electronegative substituents (oriented sterically and electrostatically toward polar and hydrophobic zones within the active site) could be considered as important structural features to improve the COX-2 inhibitory activity by 2-arylbenzofurans.

4. Conclusions

An in-silico study was primarily performed on a custom-made library comprising 58 benzofuran-containing compounds of natural origin in order to discriminate binders and non-binders of COX-2. Thus, those 2-arylbenzofurans, having an additional ring (such as oxepine) or bulky substitutions in benzofuran or aryl moieties, can be considered as potential hits for COX-2 inhibition due to the favorable characteristics for an adequate binding within the active site of this enzyme. In this regard, the best-docked compounds exhibited the best performance in MD simulations, the most-favorable binding energies, and the lowest in-vitro IC50 values for COX-2 inhibition. Consequently, the molecular docking behavior of 2-arylbenzofurans within the active site of this enzyme might be used as an appropriate strategy to predict/analyze their binding mode and potentiality as COX-2 inhibitors. Lastly, the findings through the CoMFA model demonstrated a structure-activity relationship, which shows some crucial structural requirements for enhanced enzyme inhibition. These structural features must be oriented sterically and electrostatically toward both polar and hydrophobic zones within the active site to promote crucial interactions. These oriented interactions are the main driving force for the molecular recognition of COX-2 by test ligands. In this context, four hits (i.e., 11, 29, 42, and 49) are then proposed as lead structures for the development of COX-2 inhibitors. Therefore, such information should be conserved during further studies by focusing on the development of COX-2 inhibitors based on 2-arylbenzofurans.

Funding

This study was financially supported by The Vicerrectoría de Investigaciones at Universidad Militar Nueva Granada – Validity 2018, through the Project IMP-CIAS-2924 financially supported this study.

Acknowledgments

Authors thank MU Nueva Granada for the financial support.

Conflicts of Interest

The author declares no conflict of interest. The funders had no role in the design of the study, in the collection, analyses, or interpretation of data, in the writing of the manuscript, or in the decision to publish the manuscript.

Appendix A

Table A1. Common names of test 2-arylbenzofuran-related compounds.
Table A1. Common names of test 2-arylbenzofuran-related compounds.
No.Common NameReferenceNo.Common NameReference
1moracin Y[7]30moracin L[7]
2ocophyllal A[8]31moracin U[7]
3ocophyllal B[8]32moracin P[7]
4moracin A[7]33moracin Q[7]
5moracin B[7]34moracin V[7]
6moracin F[7]35moracin O[7]
7moracin J[7]36moracin X[7]
8moracin M[7]37substituted benzofuran 2 b[41]
9(+)-obtusafuran[42]38isoparvifuran[43]
10substituted benzofuran 1 a[44]39moracin D[7]
119′-nor-7,8-dehydro-isolicarin-B[8]40moracin K[7]
12egonol[37]41kanzonol U[45]
13acetylegonol[37]42psoralidin[46]
14homoegonol[37]43liliflol B[8]
15zanthocapensol[47]44(−)-licarin B[8]
16zanthocapensate[47]45(+)-licarin B[8]
17acetylhomoegonol[37]46(−)-licarin A[8]
18moracin C[7]47(+)-licarin A[8]
19moracin I[7]48(+)-acuminatin[8]
20moracin N[7]49maximol [48]
21moracin S[7]50picrasmalignan A[49]
22moracin T[7]51salvianolic acid B[50]
23dinklagein A[9]52substituted pterocarpan 1 c[51]
24dinklagein B[9]53substituted pterocarpan 2 d[51]
25moracin R[7]54(+)-glyceollin I[52]
26moracin Z[7]55sulfuretin[53]
27moracin W[7]56(+)-mirandin A[8]
28moracin G[7]57kadsurenone[8]
29moracin H[7]58(+)-denudatin B[8]
a 6-hydroxy-2-(2-hydroxy-4-methoxyphenyl)-benzofuran, b ethyl 2-methyl-5-phenoxybenzofuran-3-carboxylate, c 2S,3S-3,8-dihydroxy-2,9-dimethoxy-pterocarpan, d 2R,3R-3,8-dihydroxy-2,9-dimethoxy-pterocarpan.
Table A2. Molecular docking results (Vina scores) of 2-arylbenzofuran-related compounds with COX-2.
Table A2. Molecular docking results (Vina scores) of 2-arylbenzofuran-related compounds with COX-2.
No.ES1 (1CVU) aES2 (1CX2) aES3 (4COX) aES4 (5COX) aES5 (1PXX) aES6 (3NT1) aES7 (3LN1) aES8 (4RS0) aES9 (4PH9) a
VS ± SD bRMS cVS ± SD bRMS cVS ± SD bRMS cVS ± SD bRMS cVS ± SD bRMS cVS± SD bRMS cVS ± SD bRMS cVS ± SD bRMS cVS ± SD bRMS c
1−8.8 ± 0.10.435−9.1 ± 0.10.337−8.6 ± 0.10.808−8.4 ± 0.20.566−8.5 ± 0.20.361−8.3 ± 0.10.879−9.7 ± 0.20.699−8.1 ± 0.10.299−8.2 ± 0.10.303
2−9.6 ± 0.20.308−8.9 ± 0.20.225−10.0 ± 0.20.319−9.7 ± 0.10.269−9.7 ± 0.10.240−8.3 ± 0.20.282−10.2 ± 0.10.674−8.5 ± 0.10.183−7.8 ± 0.10.264
3−7.3 ± 0.20.239−9.0 ± 0.20.272−8.9 ± 0.20.255−8.6 ± 0.10.848−8.7 ± 0.10.441−7.6 ± 0.10.264−8.9 ± 0.10.178−7.6 ± 0.10.165−6.3 ± 0.20.210
4−8.7 ± 0.10.933−8.7 ± 0.10.877−8.5 ± 0.10.936−8.1 ± 0.40.492−8.3 ± 0.20.960−8.2 ± 0.20.358−9.5 ± 0.10.447−8.1 ± 0.20.334−7.8 ± 0.11.075
5−8.6 ± 0.20.973−9.1 ± 0.10.262−8.1 ± 0.10.699−8.4 ± 0.30.311−8.4 ± 0.10.355−8.2 ± 0.20.238−9.6 ± 0.10.128−8.1 ± 0.10.325−8.0 ± 0.10.344
6−8.5 ± 0.10.342−8.3 ± 0.10.344−8.2 ± 0.20.587−8.1 ± 0.10.288−8.0 ± 0.10.885−7.3 ± 0.10.927−9.8 ± 0.20.038−7.4 ± 0.10.351−7.7 ± 0.11.182
7−8.5 ± 0.10.346−9 ± 0.20.352−7.9 ± 0.20.279−8.4 ± 0.10.880−8.2 ± 0.10.334−8.4 ± 0.10.374−9.5 ± 0.10.673−8.0 ± 0.10.344−8.0 ± 0.11.201
8−8.5 ± 0.10.432−8.8 ± 0.10.854−8.4 ± 0.10.924−8.5 ± 0.10.869−8.4 ± 0.20.385−8.4 ± 0.10.400−9.3 ± 0.20.370−8.2 ± 0.10.650−7.9 ± 0.10.298
9−8.7 ± 0.20.869−8.8 ± 0.20.503−8.9 ± 0.20.986−9.8 ± 0.20.304−9.0 ± 0.10.308−8.6 ± 0.10.455−9.2 ± 0.20.734−8.6 ± 0.10.664−8.7 ± 0.10.619
10−8.4 ± 0.10.346−8.8 ± 0.10.314−8.7 ± 0.20.358−8.6 ± 0.10.707−8.8 ± 0.10.329−8.1 ± 0.10.490−9.3 ± 0.20.941−8.1 ± 0.10.362−7.8 ± 0.10.274
11−10.6 ± 0.20.276−10.5 ± 0.40.802−10.8 ± 0.20.299−10.9 ± 0.10.924−11.3 ± 0.10.178−10.8 ± 0.20.661−10.6 ± 0.10.935−10.4 ± 0.20.118−10.9 ± 0.20.220
12−9.2 ± 0.10.395−9.9 ± 0.40.797−9.7 ± 0.10.499−9.2 ± 0.10.789−9.2 ± 0.20.324−6.9 ± 0.10.252−9.9 ± 0.20.860−8.0 ± 0.10.213−8.9 ± 0.10.313
13−9.8 ± 0.10.652−10.6 ± 0.20.662−10.6 ± 0.10.287−9.5 ± 0.10.422−9.6 ± 0.10.273−7.6 ± 0.30.255−9.8 ± 0.20.594−7.4 ± 0.31.110−8.3 ± 0.20.947
14−8.7 ± 0.31.054−8.7 ± 0.20.577−9.4 ± 0.20.362−8.2 ± 0.10.716−8.8 ± 0.10.285−7.7 ± 0.10.307−9.0 ± 0.20.918−7.8 ± 0.10.275−8.9 ± 0.10.267
15−9.6 ± 0.11.134−10.3 ± 0.20.605−9.0 ± 0.10.327−9.3 ± 0.20.318−9.4 ± 0.10.295−6.9 ± 0.30.277−9.8 ± 0.10.829−7.1 ± 0.11.393−9.0 ± 0.20.244
16−9.6 ± 0.10.276−10.6 ± 0.20.767−9.6 ± 0.10.433−9.1 ± 0.30.258−9.6 ± 0.10.477−7.4 ± 0.20.252−9.6 ± 0.20.477−7.5 ± 0.20.107−8.8 ± 0.10.260
17−9.2 ± 0.10.587−9.2 ± 0.30.572−10.2 ± 0.20.424−8.9 ± 0.20.414−9.4 ± 0.10.253−7.6 ± 0.10.293−9.4 ± 0.10.573−7.7 ± 0.10.555−8.2 ± 0.10.974
18−9.0 ± 0.10.485−8.7 ± 0.30.350−9.5 ± 0.10.316−9.4 ± 0.10.373−9.2 ± 0.10.314−7.7 ± 0.10.351−11.0 ± 0.10.827−7.0 ± 0.50.391−7.4 ± 0.20.418
19−9.3 ± 0.20.385−8.5 ± 0.10.353−9.2 ± 0.20.768−8.9 ± 0.10.872−8.6 ± 0.10.386−8.5 ± 0.20.860−10.7 ± 0.10.144−9.1 ± 0.10.311−8.4 ± 0.10.360
20−9.5 ± 0.11.199−9.6 ± 0.20.900−9.6 ± 0.20.317−9.1 ± 0.20.333−9.8 ± 0.20.325−7.4 ± 0.10.924−11.1 ± 0.10.586−7.3 ± 0.10.393−8.6 ± 0.10.388
21−9.9 ± 0.10.361−9.6 ± 0.40.913−9.2 ± 0.10.673−10.3 ± 0.10.789−9.0 ± 0.10.344−8.1 ± 0.10.845−10.8 ± 0.20.532−8.6 ± 0.10.305−8.5 ± 0.10.403
22−9.1 ± 0.10.274−10.0 ± 0.20.489−8.8 ± 0.10.520−8.9 ± 0.10.394−9.1 ± 0.10.343−8.4 ± 0.10.331−10.7 ± 0.30.958−8.1 ± 0.11.139−8.3 ± 0.10.359
23−9.4 ± 0.20.345−9.6 ± 0.30.563−9 ± 0.10.569−10.4 ± 0.20.338−9.2 ± 0.10.784−8.5 ± 0.20.632−10.9 ± 0.10.358−8.9 ± 0.20.478−8.4 ± 0.20.412
24−9.3 ± 0.20.756−9.2 ± 0.30.478−9.2 ± 0.20.656−10.6 ± 0.20.589−9.4 ± 0.10.779−8.2 ± 0.10.456−10.5 ± 0.30.665−8.7 ± 0.20.851−8 ± 0.20.458
25−8.4 ± 0.10.372−8.7 ± 0.30.391−9.3 ± 0.10.401−8.7 ± 0.20.885−9.6 ± 0.10.953−6.7 ± 0.30.382−10.8 ± 0.10.249−7.6 ± 0.20.391−6.7 ± 0.10.412
26−9.0 ± 0.10.950−9.5 ± 0.20.919−9.4 ± 0.10.473−8.8 ± 0.30.776−9.3 ± 0.10.322−6.8 ± 0.70.395−10.7 ± 0.20.719−7.5 ± 0.10.554−8.0 ± 0.10.392
27−10.7 ± 0.10.535−9.7 ± 0.30.561−9.9 ± 0.10.977−10.0 ± 0.30.571−9.7 ± 0.50.823−10.0 ± 0.10.507−11.1 ± 0.30.413−9.8 ± 0.20.373−10.5 ± 0.20.403
28−10.7 ± 0.10.859−9.8 ± 0.10.344−10.4 ± 0.10.275−10.3 ± 0.10.336−9.9 ± 0.10.329−10.4 ± 0.50.303−11.3 ± 0.10.660−10.6 ± 0.10.290−10.8 ± 0.10.303
29−10.3 ± 0.20.351−10.9 ± 0.50.692−10.3 ± 0.10.328−10.0 ± 0.10.350−10.1 ± 0.10.336−10.5 ± 0.20.288−11.4 ± 0.10.883−10.6 ± 0.10.257−10.3 ± 0.53.621
30−10.2 ± 0.10.864−10.5 ± 0.20.329−10.2 ± 0.10.472−10.4 ± 0.10.366−9.3 ± 0.20.363−10.4 ± 0.40.346−11.3 ± 0.10.604−10.7 ± 0.10.354−10.8 ± 0.10.355
31−9.1 ± 0.10.340−9.9 ± 0.20.372−10.0 ± 0.40.396−9.7 ± 0.10.404−10.7 ± 0.10.932−7.2 ± 0.30.353−11.3 ± 0.10.622−8.7 ± 0.10.404−8.8 ± 0.10.419
32−8.8 ± 0.10.915−9.3 ± 0.10.316−10.3 ± 0.10.353−9.7 ± 0.10.918−9.5 ± 0.10.335−6.6 ± 1.80.547−11.3 ± 0.20.999−6.3 ± 0.20.300−7.8 ± 0.10.390
33−7.0 ± 0.10.311−9.2 ± 0.20.301−10.2 ± 0.20.302−9.0 ± 0.10.297−9.4 ± 0.10.263−7.1 ± 0.10.213−9.6 ± 0.20.662−6.4 ± 0.10.206−7.0 ± 0.10.288
34−8.6 ± 0.20.340−9.1 ± 0.10.359−10.3 ± 0.10.409−9.8 ± 0.30.931−9.0 ± 0.10.672−6.4 ± 1.70.380−11.2 ± 0.10.918−6.3 ± 0.10.364−7.0 ± 0.10.358
35−8.2 ± 0.20.481−9.6 ± 0.10.314−10.0 ± 0.10.353−9.2 ± 0.10.358−9.8 ± 0.10.893−7.2 ± 0.10.381−9.8 ± 0.10.683−6.9 ± 0.20.352−7.6 ± 0.11.113
36−9.1 ± 0.10.792−9.7 ± 0.30.321−9.0 ± 0.10.507−8.3 ± 0.30.898−9.1 ± 0.10.735−8.2 ± 0.10.296−10.7 ± 0.20.815−8.6 ± 0.10.302−8.5 ± 0.11.144
37−9.4 ± 0.10.592−9.4 ± 0.10.474−9.3 ± 0.10.630−9.1 ± 0.20.887−9.5 ± 0.10.773−7.9 ± 0.10.225−9.6 ± 0.20.865−8.4 ± 0.10.140−8.7 ± 0.20.216
38−9.0 ± 0.10.463−8.6 ± 0.20.839−8.8 ± 0.20.502−9.4 ± 0.10.301−8.5 ± 0.10.486−8.2 ± 0.10.327−9.8 ± 0.10.049−8.2 ± 0.10.538−8.1 ± 0.10.328
39−10.1 ± 0.10.306−10.3 ± 0.50.334−9.7 ± 0.10.922−9.2 ± 0.10.299−10.2 ± 0.10.312−10.2 ± 0.30.339−11.8 ± 0.20.826−10.2 ± 0.10.316−10.1 ± 0.20.331
40−10.0 ± 0.10.562−10.5 ± 0.10.283−10.1 ± 0.10.954−10.1 ± 0.10.492−9.8 ± 0.20.335−9.8 ± 0.50.886−11.7 ± 0.10.712−10.6 ± 0.10.305−10.1 ± 0.20.398
41−10.3 ± 0.10.322−10.3 ± 0.10.258−10.3 ± 0.20.307−9.6 ± 0.10.793−9.3 ± 0.10.780−8.7 ± 0.10.242−10.7 ± 0.10.747−8.9 ± 0.20.222−9.2 ± 0.10.352
42−10.5 ± 0.10.302−11.6 ± 0.10.328−10.9 ± 0.10.324−10.9 ± 0.10.315−10.6 ± 0.10.303−11.0 ± 0.10.258−12.0 ± 0.20.683−11.9 ± 0.10.264−11.9 ± 0.10.241
43−9.2 ± 0.10.221−9.6 ± 0.20.324−10.2 ± 0.10.301−10.0 ± 0.10.310−9.9 ± 0.20.283−8.0 ± 0.10.296−9.7 ± 0.10.665−8.6 ± 0.10.226−8.7 ± 0.10.275
44−9.1 ± 0.30.255−9.7 ± 0.20.249−9.8 ± 0.20.289−10.2 ± 0.20.239−9.6 ± 0.10.243−7.1 ± 0.10.278−9.8 ± 0.20.915−8.3 ± 0.10.172−8.6 ± 0.10.239
45−7.4 ± 0.20.246−9.9 ± 0.40.236−10.1 ± 0.30.234−9.3 ± 0.20.303−9.9 ± 0.20.287−8.5 ± 0.20.411−10.2 ± 0.20.904−8.7 ± 0.10.142−8.9 ± 0.10.253
46−9.0 ± 0.20.286−9.2 ± 0.10.908−9.7 ± 0.10.264−9.5 ± 0.20.295−9.7 ± 0.20.290−6.8 ± 0.10.738−9.8 ± 0.20.943−7.9 ± 0.10.292−7.5 ± 0.20.262
47−7.3 ± 0.10.267−9.4 ± 0.20.300−9.6 ± 0.10.295−9.4 ± 0.10.282−9.5 ± 0.20.225−7.3 ± 0.10.261−10.1 ± 0.20.891−7.9 ± 0.10.254−7.5 ± 0.10.330
48−7.1 ± 0.20.274−9.8 ± 0.30.298−9.5 ± 0.10.405−9.5 ± 0.20.245−9.7 ± 0.10.229−7.9 ± 0.10.264−10.1 ± 0.10.912−8.1 ± 0.10.187−7.4 ± 0.10.179
49−10.2 ± 0.10.582−11.4 ± 0.10.329−11.6 ± 0.10.513−12.2 ± 0.20.319−10.8 ± 0.20.671−11.1 ± 0.60.512−12.7 ± 0.20.878−10.6 ± 0.10.352−10.7 ± 0.10.308
50−8.3 ± 0.10.281−9.1 ± 0.10.515−11.3 ± 0.20.712−9.2 ± 0.20.846−6.5 ± 0.20.428−7.1 ± 0.30.258−10.1 ± 0.10.910−7.1 ± 0.10.240−7.6 ± 0.10.488
51−9.2 ± 0.10.510−10.5 ± 0.50.459−10.6 ± 0.20.808−9.9 ± 0.20.643−10.2 ± 0.10.372−10.8 ± 0.10.726−11.2 ± 0.20.978−10.4 ± 0.20.568−10.5 ± 0.20.525
52−8.4 ± 0.10.494−8.6 ± 0.20.334−9.0 ± 0.10.352−8.7 ± 0.20.336−8.8 ± 0.20.793−6.9 ± 0.20.297−9.1 ± 0.30.916−7.4 ± 0.10.256−6.7 ± 0.10.293
53−9.5 ± 0.10.358−9.4 ± 0.10.375−9.6 ± 0.10.556−9.1 ± 0.20.337−9.8 ± 0.10.653−8.4 ± 0.10.389−9.4 ± 0.10.889−8.9 ± 0.10.281−7.5 ± 0.20.368
54−10.5 ± 0.10.365−10.5 ± 0.10.312−10.7 ± 0.10.348−10.7 ± 0.10.354−8.9 ± 0.10.307−10.7 ± 0.60.324−10.3 ± 0.20.453−10.4 ± 0.10.242−10.8 ± 0.10.372
55−8.6 ± 0.10.376−9.1 ± 0.10.361−8.7 ± 0.20.297−8.9 ± 0.10.374−8.4 ± 0.20.338−8.3 ± 0.10.364−9.8 ± 0.10.998−8.4 ± 0.10.374−8.3 ± 0.20.412
56−7.8 ± 0.10.256−8.8 ± 0.20.252−9.0 ± 0.10.307−9.0 ± 0.10.300−8.4 ± 0.20.641−8.2 ± 0.10.257−8.8 ± 0.30.856−8.5 ± 0.10.193−7.9 ± 0.20.250
57−6.4 ± 0.10.275−8.5 ± 0.40.436−8.8 ± 0.20.281−8.2 ± 0.40.232−9.0 ± 0.10.260−6.7 ± 0.10.240−9.8 ± 0.20.643−5.9 ± 0.10.172−6.6 ± 0.20.228
58−7.6 ± 0.10.247−8.4 ± 0.10.471−9.6 ± 0.10.230−9.3 ± 0.10.246−7.8 ± 0.20.508−8.3 ± 0.10.173−8.8 ± 0.20.995−8.4 ± 0.10.186−7.9 ± 0.20.249
Natural Ligand d and Known Inhibitors e
AA−8.6 ± 0.21.381−8.4 ± 0.31.259−7.8 ± 0.31.331−7.5 ± 0.61.236−8.1 ± 0.41.562−7.7 ± 0.20.417−8.2 ± 0.30.533−7.8 ± 0.41.252−7.7 ± 0.32.993
CEL−10.7 ± 0.10.493−11.6 ± 0.10.421−10.9 ± 0.10.551−10.8 ± 0.10.313−10.9 ± 0.10.434−11.0 ± 0.10.411−12.7 ± 0.20.885−11.1 ± 0.10.500−11.4 ± 0.10.384
DCL−8.7 ± 0.30.289−8.2 ± 0.20.270−8.3 ± 0.40.256−8.1 ± 0.10.370−8.7 ± 0.10.174−8.5 ± 0.20.264−8.6 ± 0.10.838−8.0 ± 0.10.152−7.9 ± 0.10.248
IMD−8.2 ± 0.10.530−10.3 ± 0.30.446−10.7 ± 0.30.825−10.5 ± 0.10.298−9.2 ± 0.10.457−9.5 ± 0.10.289−9.9 ± 0.40.821−9.6 ± 0.10.169−9.1 ± 0.10.205
NPX−8.6 ± 0.10.854−8.5 ± 0.20.339−8.7 ± 0.10.275−8.4 ± 0.20.266−8.3 ± 0.10.309−9.3 ± 0.10.276−8.6 ± 0.10.765−8.7 ± 0.20.075−8.8 ± 0.10.221
SC58−10.7 ± 0.10.310−12 ± 0.20.313−10.9 ± 0.10.613−10.6 ± 0.20.349−11.2 ± 0.10.2631−10.6 ± 0.10.449−12.0 ± 0.20.786−10.6 ± 0.10.268−11.2 ± 0.10.356
S−IBP−7.5 ± 0.10.222−7.6 ± 0.10.650−8.2 ± 0.10.199−8.2 ± 0.20.235−8.5 ± 0.10.592−8.2 ± 0.10.180−7.5 ± 0.10.528−7.7 ± 0.10.589−7.8 ± 0.22.549
R−IBP−6.5 ± 0.10.377−6.0 ± 0.11.065−6.7 ± 0.10.550−6.1 ± 0.11.376−6.5 ± 0.20.267−6.9 ± 0.11.434−6.7 ± 0.10.665−5.5 ± 0.10.574−5.8 ± 0.10.751
aES = Enzyme Structure (PDB Code). b VS = Vina Scores (docking energies or affinities in kcal/mol) ± Standard Deviation (SD) of the affinities (best poses) after ten replicates per molecular docking calculation. c Root-mean-square deviation (RMSD) of atomic positions between the best poses after 10 replicates. d Arachidonic acid (AA). e Reported inhibitors: CEL = celecoxib (selective (S)). DCL = diclofenac (selective (S)). IMD = indomethacin (S). NPX = naproxen (N-S). SC558 = 1-phenylsulfonamide-3-trifluoromethyl-5-parabromophenylpyrazole (S). S-IBP = (S)-ibuprofen (N-S). R-IBP = (R)-ibuprofen (non-inhibitor).
Table A3. Absorption, distribution, metabolism, and excretion–toxicity (ADMET) properties of best-docked compounds within the active site of COX-2.
Table A3. Absorption, distribution, metabolism, and excretion–toxicity (ADMET) properties of best-docked compounds within the active site of COX-2.
Descriptors11294249
MW316.39348.43348.43431.54
logP2.04041.01720.95672.0612
logD3.21.051.473.72
logSw−2.82−2.51−2.41−3.48
tPSA36.9268.1579.1569.92
RotatableB4235
RigidB20222028
Flexibility0.170.080.130.15
HBD0233
HBA4554
HBD_HBA4787
Rings2214
MaxSizeRing914179
NumCharges0000
TotalCharge0000
HeavyAtoms23252532
CarbonAtoms19202028
HeteroAtoms4554
ratioH/C0.210.250.250.14
Lipinski_Violation0000
Sol (mg/L)18,800.8328,341.4231,451.2713,240.01
Sol Forecast IndexGood SolubilityGood SolubilityGood SolubilityGood Solubility
OB (VEBER)GoodGoodGoodGood
OB (EGAN)GoodGoodGoodGood
TrafficLights0001
PhospholipidosisNon InducerNon InducerNon InducerNon Inducer
Overall ResultAcceptedAcceptedAcceptedAccepted
MW: Molecular weight in Da. log P = Octanol-water partition coefficient. tPSA = topological polar surface area. HBA = Hydrogen bond atoms. HBD = Hydrogen bond donors. Sol = solubility. OB = Oral bioavailability.
Figure A1. Structures of natural ligand and selective and non-selective COX-2 inhibitors.
Figure A1. Structures of natural ligand and selective and non-selective COX-2 inhibitors.
Biomolecules 10 00176 g0a1
Figure A2. Structures of the test 2-arylbenzofuran-related compounds.
Figure A2. Structures of the test 2-arylbenzofuran-related compounds.
Biomolecules 10 00176 g0a2

References

  1. Clària, J. Cyclooxygenase-2 biology. Curr. Pharm. Des. 2003, 9, 2177–2190. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Praticò, D.; Dogné, J.-M. Selective cyclooxygenase-2 inhibitors development in cardiovascular medicine. Circulation 2005, 112, 1073–1079. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Rayar, A.M.; Lagarde, N.; Ferroud, C.; Zagury, J.F.; Montes, M.; Sylla-Iyarreta Veitia, M. Update on COX-2 selective inhibitors: Chemical classification, side effects and their use in cancers and neuronal diseases. Curr. Top. Med. Chem. 2017, 17, 2935–2956. [Google Scholar] [CrossRef] [PubMed]
  4. White, W.B.; Faich, G.; Borer, J.S.; Makuch, R.W. Cardiovascular thrombotic events in arthritis trials of the cyclooxygenase-2 inhibitor celecoxib. ACC Curr. J. Rev. 2003, 12, 14–15. [Google Scholar] [CrossRef]
  5. Khanam, H. Shamsuzzaman Bioactive Benzofuran derivatives: A review. Eur. J. Med. Chem. 2015, 97, 483–504. [Google Scholar] [CrossRef] [PubMed]
  6. Ragab, F.A.E.-F.; Eid, N.M.; Hassan, G.S.; Nissan, Y.M. Synthesis and anti-inflammatory activity of some benzofuran and benzopyran-4-one derivatives. Chem. Pharm. Bull. (Tokyo) 2012, 60, 110–120. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Naik, R.; Harmalkar, D.S.; Xu, X.; Jang, K.; Lee, K. Bioactive benzofuran derivatives: Moracins A–Z in medicinal chemistry. Eur. J. Med. Chem. 2015, 90, 379–393. [Google Scholar] [CrossRef]
  8. Coy, E.D.; Cuca, L.E.; Sefkow, M. COX, LOX and platelet aggregation inhibitory properties of Lauraceae neolignans. Bioorg. Med. Chem. Lett. 2009, 19, 6922–6925. [Google Scholar] [CrossRef]
  9. Kapche, D.W.F.G.; Lekane, N.M.; Kulabas, S.S.; Ipek, H.; Tok, T.T.; Ngadjui, B.T.; Demirtas, I.; Tumer, T.B. Aryl benzofuran derivatives from the stem bark of Calpocalyx dinklagei attenuate inflammation. Phytochemistry 2017, 141, 70–79. [Google Scholar] [CrossRef]
  10. Kiefer, J.R.; Pawlitz, J.L.; Moreland, K.T.; Stegeman, R.A.; Hood, W.F.; Gierse, J.K.; Stevens, A.M.; Goodwin, D.C.; Rowlinson, S.W.; Marnett, L.J.; et al. Structural insights into the stereochemistry of the cyclooxygenase reaction. Nature 2000, 405, 97–101. [Google Scholar] [CrossRef]
  11. Kurumbail, R.G.; Stevens, A.M.; Gierse, J.K.; McDonald, J.J.; Stegeman, R.A.; Pak, J.Y.; Gildehaus, D.; iyashiro, J.M.; Penning, T.D.; Seibert, K.; et al. Structural basis for selective inhibition of cyclooxygenase-2 by anti-inflammatory agents. Nature 1996, 384, 644–648. [Google Scholar] [CrossRef] [PubMed]
  12. Rowlinson, S.W.; Kiefer, J.R.; Prusakiewicz, J.J.; Pawlitz, J.L.; Kozak, K.R.; Kalgutkar, A.S.; Stallings, W.C.; Kurumbail, R.G.; Marnett, L.J. A novel mechanism of cyclooxygenase-2 inhibition involving interactions with Ser-530 and Tyr-385. J. Biol. Chem. 2003, 278, 45763–45769. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Duggan, K.C.; Walters, M.J.; Musee, J.; Harp, J.M.; Kiefer, J.R.; Oates, J.A.; Marnett, L.J. Molecular basis for cyclooxygenase inhibition by the non-steroidal anti-inflammatory drug naproxen. J. Biol. Chem. 2010, 285, 34950–34959. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Wang, J.L.; Limburg, D.; Graneto, M.J.; Springer, J.; Hamper, J.R.B.; Liao, S.; Pawlitz, J.L.; Kurumbail, R.G.; Maziasz, T.; Talley, J.J.; et al. The novel benzopyran class of selective cyclooxygenase-2 inhibitors. Part 2: The second clinical candidate having a shorter and favorable human half-life. Bioorg. Med. Chem. Lett. 2010, 20, 7159–7163. [Google Scholar] [CrossRef]
  15. Blobaum, A.L.; Xu, S.; Rowlinson, S.W.; Duggan, K.C.; Banerjee, S.; Kudalkar, S.N.; Birmingham, W.R.; Ghebreselasie, K.; Marnett, L.J. Action at a distance: Mutations of peripheral residues transform rapid reversible inhibitors to slow, tight binders of cyclooxygenase-2. J. Biol. Chem. 2015, 290, 12793–12803. [Google Scholar] [CrossRef] [Green Version]
  16. Orlando, B.J.; Lucido, M.J.; Malkowski, M.G. The structure of ibuprofen bound to cyclooxygenase-2. J. Struct. Biol. 2015, 189, 62–66. [Google Scholar] [CrossRef] [Green Version]
  17. Sander, T.; Freyss, J.; von Korff, M.; Rufener, C. DataWarrior: An open-source program for chemistry aware data visualization and analysis. J. Chem. Inf. Model. 2015, 55, 460–473. [Google Scholar] [CrossRef]
  18. Gaulton, A.; Hersey, A.; Nowotka, M.; Bento, A.P.; Chambers, J.; Mendez, D.; Mutowo, P.; Atkinson, F.; Bellis, L.J.; Cibrián-Uhalte, E.; et al. The ChEMBL database in 2017. Nucleic Acids Res. 2017, 45, D945–D954. [Google Scholar] [CrossRef]
  19. Mysinger, M.M.; Carchia, M.; Irwin, J.J.; Shoichet, B.K. Directory of Useful Decoys, Enhanced (DUD-E): Better ligands and decoys for better benchmarking. J. Med. Chem. 2012, 55, 6582–6594. [Google Scholar] [CrossRef]
  20. Empereur-Mot, C.; Zagury, J.-F.; Montes, M. Screening explorer–An interactive tool for the analysis of screening results. J. Chem. Inf. Model. 2016, 56, 2281–2286. [Google Scholar] [CrossRef]
  21. Guerrero-Perilla, C.; Bernal, F.A.; Coy-Barrera, E. Insights into the interaction and binding mode of a set of antifungal azoles as inhibitors of potential fungal enzyme-based targets. Mol. Divers. 2018, 22, 929–942. [Google Scholar] [CrossRef]
  22. Kumari, R.; Kumar, R.; Lynn, A. g_mmpbsa—A GROMACS tool for high-throughput MM-PBSA calculations. J. Chem. Inf. Model. 2014, 54, 1951–1962. [Google Scholar] [CrossRef] [PubMed]
  23. Sun, H.; Li, Y.; Tian, S.; Xu, L.; Hou, T. Assessing the performance of MM/PBSA and MM/GBSA methods. 4. Accuracies of MM/PBSA and MM/GBSA methodologies evaluated by various simulation protocols using PDBbind data set. Phys. Chem. Chem. Phys. 2014, 16, 16719–16729. [Google Scholar] [CrossRef] [PubMed]
  24. Sun, H.; Li, Y.; Shen, M.; Tian, S.; Xu, L.; Pan, P.; Guan, Y.; Hou, T. Assessing the performance of MM/PBSA and MM/GBSA methods. 5. Improved docking performance using high solute dielectric constant MM/GBSA and MM/PBSA rescoring. Phys. Chem. Chem. Phys. 2014, 16, 22035–22045. [Google Scholar] [CrossRef]
  25. Tosco, P.; Balle, T. Open3DQSAR: A new open-source software aimed at high-throughput chemometric analysis of molecular interaction fields. J. Mol. Model. 2010, 17, 201–208. [Google Scholar] [CrossRef]
  26. Morris, G.M.; Lim-Wilby, M. Molecular docking BT. In Molecular Modeling of Proteins; Kukol, A., Ed.; Humana Press: Totowa, NJ, USA, 2008; pp. 365–382. ISBN 978-1-59745-177-2. [Google Scholar]
  27. Sliwoski, G.; Kothiwale, S.; Meiler, J.; Lowe, E.W., Jr. Computational methods in drug discovery. Pharmacol. Rev. 2013, 66, 334–395. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Li, Q.; Shah, S. Structure-based virtual screening BT. In Protein Bioinformatics: From Protein Modifications and Networks to Proteomics; Wu, C.H., Arighi, C.N., Ross, K.E., Eds.; Springer: New York, NY, USA, 2017; pp. 111–124. ISBN 978-1-4939-6783-4. [Google Scholar]
  29. Westermaier, Y.; Barril, X.; Scapozza, L. Virtual screening: An in silico tool for interlacing the chemical universe with the proteome. Methods 2015, 71, 44–57. [Google Scholar] [CrossRef]
  30. Chen, Y.-C. Beware of docking! Trends Pharmacol. Sci. 2015, 36, 78–95. [Google Scholar] [CrossRef]
  31. Ma, H.; Bandos, A.I.; Gur, D. On the use of partial area under the ROC curve for comparison of two diagnostic tests. Biom. J. 2015, 57, 304–320. [Google Scholar] [CrossRef]
  32. Yang, H.J.; Youn, H.; Seong, K.M.; Yun, Y.J.; Kim, W.; Kim, Y.H.; Lee, J.Y.; Kim, C.S.; Jin, Y.-W.; Youn, B. Psoralidin, a dual inhibitor of COX-2 and 5-LOX, regulates ionizing radiation (IR)-induced pulmonary inflammation. Biochem. Pharmacol. 2011, 82, 524–534. [Google Scholar] [CrossRef]
  33. Pang, X.; Zhou, L.; Zhang, L.; Xu, L.; Zhang, X. Two rules on the protein-ligand interaction. Nat. Preced. 2008. [Google Scholar] [CrossRef]
  34. Chaudhary, N.; Aparoy, P. Deciphering the mechanism behind the varied binding activities of COXIBs through Molecular Dynamic Simulations, MM-PBSA binding energy calculations and per-residue energy decomposition studies. J. Biomol. Struct. Dyn. 2017, 35, 868–882. [Google Scholar] [CrossRef] [PubMed]
  35. Bouaziz-Terrachet, S.; Toumi-Maouche, A.; Maouche, B.; Taïri-Kellou, S. Modeling the binding modes of stilbene analogs to cyclooxygenase-2: A molecular docking study. J. Mol. Model. 2010, 16, 1919–1929. [Google Scholar] [CrossRef] [PubMed]
  36. Lagorce, D.; Bouslama, L.; Becot, J.; Miteva, M.A.; Villoutreix, B.O. FAF-Drugs4: Free ADME-tox filtering computations for chemical biology and early stages drug discovery. Bioinformatics 2017, 33, 3658–3660. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Bertanha, C.S.; Braguine, C.G.; Moraes, A.C.G.; Gimenez, V.M.M.; Groppo, M.; Silva, M.L.A.; Cunha, W.R.; Januário, A.H.; Pauletti, P.M. Cyclooxygenase inhibitory properties of nor-neolignans from Styrax pohlii. Nat. Prod. Res. 2012, 26, 2323–2329. [Google Scholar] [CrossRef] [PubMed]
  38. Roy, K.; Kar, S.; Das, R.N. Introduction to 3D-QSAR. In Understanding the Basics of QSAR for Applications in Pharmaceutical Sciences and Risk Assessment; Elsevier: London, UK, 2015; pp. 291–317. ISBN 9780128015056. [Google Scholar]
  39. Doweyko, A.M. Three-dimensional quantitative structure–activity relationship: The state of the art. In Comprehensive Medicinal Chemistry II; Elsevier: London, UK, 2007; pp. 575–595. ISBN 9780080450445. [Google Scholar]
  40. Clark, R.D.; Fox, P.C. Statistical variation in progressive scrambling. J. Comput. Aided Mol. Des. 2004, 18, 563–576. [Google Scholar] [CrossRef]
  41. Yadav, P.; Singh, P.; Tewari, A.K. Design, synthesis, docking and anti-inflammatory evaluation of novel series of benzofuran based prodrugs. Bioorg. Med. Chem. Lett. 2014, 24, 2251–2255. [Google Scholar] [CrossRef]
  42. Gregson, M.; Ollis, W.D.; Redman, B.T.; Sutherland, I.O.; Dietrichs, H.H.; Gottlieb, O.R. Obtusastyrene and obtustyrene, cinnamylphenols from Dalbergia retusa. Phytochemistry 1978, 17, 1395–1400. [Google Scholar] [CrossRef]
  43. Muangnoicharoen, N.; Frahm, A.W. Arylbenzofurans from Dalbergia parviflora. Phytochemistry 1981, 20, 291–293. [Google Scholar] [CrossRef]
  44. Jurd, L.; Wong, R.Y. Diarylheptanoid and other phenolic constituents of Centrolobium species. Aust. J. Chem. 1984, 37, 1127. [Google Scholar] [CrossRef]
  45. Fukai, T.; Sheng, C.-B.; Horikoshi, T.; Nomura, T. Isoprenylated flavonoids from underground parts of Glycyrrhiza glabra. Phytochemistry 1996, 43, 1119–1124. [Google Scholar] [CrossRef]
  46. O’Neill, M.J.; Adesanya, S.A.; Roberts, M.F.; Pantry, I.R. Inducible isoflavonoids from the lima bean, Phaseolus lunatus. Phytochemistry 1986, 25, 1315–1322. [Google Scholar] [CrossRef]
  47. Luo, X.; Pedro, L.; Milic, V.; Mulhovo, S.; Duarte, A.; Duarte, N.; Ferreira, M.-J. Antibacterial Benzofuran Neolignans and Benzophenanthridine Alkaloids from the Roots ofZanthoxylum capense. Planta Med. 2011, 78, 148–153. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Zawawi, N.K.N.A.; Ahmat, N.; Mazatulikhma, M.Z.; Shafiq, R.M.; Wahid, N.H.; Sufian, A.S. Bioactive oligostilbenoids from Shorea maxwelliana King and their chemotaxonomic significance. Nat. Prod. Res. 2012, 27, 1589–1593. [Google Scholar] [CrossRef]
  49. Zhao, F.; Chen, L.; Bi, C.; Zhang, M.; Jiao, W.; Yao, X. In vitro anti-inflammatory effect of picrasmalignan A by the inhibition of iNOS and COX-2 expression in LPS-activated macrophage RAW 264.7 cells. Mol. Med. Rep. 2013, 8, 1575–1579. [Google Scholar] [CrossRef] [Green Version]
  50. Zhao, L.-M.; Liang, X.-T.; Li, L.-N. Prionitisides A and B, two phenolic glycosides from Salvia prionitis. Phytochemistry 1996, 42, 899–901. [Google Scholar] [CrossRef]
  51. Máximo, P.; Lourenço, A. A pterocarpan from Ulex parviflorus. Phytochemistry 1998, 48, 359–362. [Google Scholar] [CrossRef]
  52. Burden, R.S.; Bailey, J.A. Structure of the phytoalexin from soybean. Phytochemistry 1975, 14, 1389–1390. [Google Scholar] [CrossRef]
  53. Shu, L.; Shen, Y.; Yang, L.; Gao, X.; Hu, Q.-F. Flavonoids derivatives from arundina graminifolia and their cytotoxicity. Asian J. Chem. 2013, 25, 8358–8360. [Google Scholar] [CrossRef]
Figure 1. Similarity chart of 2-arylbenzofurans 1–58 after structure similarity analysis using the substructure fragment dictionary-based binary fingerprint descriptor (FragFp) ([17]. Colors according to the heatmap based on the structure similarity index between 0 (red) and 1 (light blue). Clusters according to: (a) moracin-type, (b) 2H-pyran- and benzofuro[6.5-b]furan-containing, (c) propyl(en)-susbtituted, (d) dihydropyran-containing, (e) furocyclohexadienone-type, and (f) licarin-type.
Figure 1. Similarity chart of 2-arylbenzofurans 1–58 after structure similarity analysis using the substructure fragment dictionary-based binary fingerprint descriptor (FragFp) ([17]. Colors according to the heatmap based on the structure similarity index between 0 (red) and 1 (light blue). Clusters according to: (a) moracin-type, (b) 2H-pyran- and benzofuro[6.5-b]furan-containing, (c) propyl(en)-susbtituted, (d) dihydropyran-containing, (e) furocyclohexadienone-type, and (f) licarin-type.
Biomolecules 10 00176 g001
Figure 2. (a) Docked (blue sticks) and co-crystallized (red sticks) superposition of celecoxib after re-docking within the active site of COX-2. (b) Receiver operating characteristic (ROC) and (c) enrichment curve for the benchmarking of the docking protocol using a set of active compounds and decoys against COX-2 (3LN1).
Figure 2. (a) Docked (blue sticks) and co-crystallized (red sticks) superposition of celecoxib after re-docking within the active site of COX-2. (b) Receiver operating characteristic (ROC) and (c) enrichment curve for the benchmarking of the docking protocol using a set of active compounds and decoys against COX-2 (3LN1).
Biomolecules 10 00176 g002
Figure 3. (a) Superposition of chains-A from nine COX-2 structures (represented as solid ribbons) after sequence alignment. The active site is highlighted by the H-bond donor/acceptor surface, based on the heatmap (orange = donor, blue = acceptor). The HEME group as grey sticks. (b) Sequence alignment chains-A from 9 COX-2 crystal enzymes. The figure was produced using ESPript (http://espript.ibcp.fr).
Figure 3. (a) Superposition of chains-A from nine COX-2 structures (represented as solid ribbons) after sequence alignment. The active site is highlighted by the H-bond donor/acceptor surface, based on the heatmap (orange = donor, blue = acceptor). The HEME group as grey sticks. (b) Sequence alignment chains-A from 9 COX-2 crystal enzymes. The figure was produced using ESPript (http://espript.ibcp.fr).
Biomolecules 10 00176 g003
Figure 4. (a) Boxplot from Vina scores after docking between 2-arylbenzofuran-related compounds and nine crystal structures of COX-2. (b) Score plot from Principal Components Analysis (PCA) on Vina scores after docking of test compounds with COX-2 structures. Clustering: green (group 1), dark blue (group 2), Yellow (group 3), and red (Group 4). <VS = lowest Vina Scores. >VS = Highest Vina Scores.
Figure 4. (a) Boxplot from Vina scores after docking between 2-arylbenzofuran-related compounds and nine crystal structures of COX-2. (b) Score plot from Principal Components Analysis (PCA) on Vina scores after docking of test compounds with COX-2 structures. Clustering: green (group 1), dark blue (group 2), Yellow (group 3), and red (Group 4). <VS = lowest Vina Scores. >VS = Highest Vina Scores.
Biomolecules 10 00176 g004
Figure 5. Three-dimensional (3D) models and ligand interaction diagrams (LID) of COX-2·ligand complexes of the lowest-scored pose of the best docked compounds 11, 29, 42, and 49, as well as celecoxib (cel). In LID, light green, aquamarine, and purple lines depict the active site surface. Dark green arrows indicate the π-π stacking between residues and aromatic moieties of compounds. Magenta arrows indicate the H-bonding. Residues are differentiated by colors according to the interaction type as indicated in the LID legend: hydrophobic (green), polar (aquamarine), and a positive charge (purple). 3D-models showed docked compounds in bold sticks and enzymes as a light green cartoon.
Figure 5. Three-dimensional (3D) models and ligand interaction diagrams (LID) of COX-2·ligand complexes of the lowest-scored pose of the best docked compounds 11, 29, 42, and 49, as well as celecoxib (cel). In LID, light green, aquamarine, and purple lines depict the active site surface. Dark green arrows indicate the π-π stacking between residues and aromatic moieties of compounds. Magenta arrows indicate the H-bonding. Residues are differentiated by colors according to the interaction type as indicated in the LID legend: hydrophobic (green), polar (aquamarine), and a positive charge (purple). 3D-models showed docked compounds in bold sticks and enzymes as a light green cartoon.
Biomolecules 10 00176 g005
Figure 6. HOMO (highest occupied molecular orbital) and LUMO (lowest unoccupied molecular orbital) energies (in eV), HOMO-LUMO gaps (ΔE), and corresponding orbitals obtained by density functional theory (DFT) level for the best docked compounds (11, 29, 42, 49, and celecoxib (cel)).
Figure 6. HOMO (highest occupied molecular orbital) and LUMO (lowest unoccupied molecular orbital) energies (in eV), HOMO-LUMO gaps (ΔE), and corresponding orbitals obtained by density functional theory (DFT) level for the best docked compounds (11, 29, 42, 49, and celecoxib (cel)).
Biomolecules 10 00176 g006
Figure 7. Molecular dynamics (MD) simulations during 150 ns for COX-2 apoenzyme (black line) and docked separately with 11 (purple line), 29 (green line), 42 (blue line), 49 (red line), and celecoxib (orange line). (a) Root mean square deviations (RMSD) along MD-simulated trajectories. (b) Root mean quadratic fluctuations (RMSF) of COX-2 residues (chain A) along MD-simulated trajectories.
Figure 7. Molecular dynamics (MD) simulations during 150 ns for COX-2 apoenzyme (black line) and docked separately with 11 (purple line), 29 (green line), 42 (blue line), 49 (red line), and celecoxib (orange line). (a) Root mean square deviations (RMSD) along MD-simulated trajectories. (b) Root mean quadratic fluctuations (RMSF) of COX-2 residues (chain A) along MD-simulated trajectories.
Biomolecules 10 00176 g007
Figure 8. Contour maps from the field contributions from the comparative molecular field analysis (CoMFA). (a) Electrostatic contour maps. Blue and red contours depict zones where positive and negative charges, respectively, will favor the COX-2 inhibitory activity. (b) Steric contour maps. Yellow and green contours indicate regions where a bulky group is unfavorable and favorable, respectively, for COX-2 inhibition.
Figure 8. Contour maps from the field contributions from the comparative molecular field analysis (CoMFA). (a) Electrostatic contour maps. Blue and red contours depict zones where positive and negative charges, respectively, will favor the COX-2 inhibitory activity. (b) Steric contour maps. Yellow and green contours indicate regions where a bulky group is unfavorable and favorable, respectively, for COX-2 inhibition.
Biomolecules 10 00176 g008
Table 1. Data of the crystal structures of COX-2.
Table 1. Data of the crystal structures of COX-2.
Code aPDB ID bResolution [Å]BoundSourceReference c
ES11CVU2.40arachidonic acidMus musculus[10]
ES25COX3.00unboundedM. musculus[11]
ES34COX2.90indomethacinM. musculus[11]
ES41CX23.00SC-558M. musculus[11]
ES51PXX2.90diclofenacM. musculus[12]
ES63NT11.73naproxenM. musculus[13]
ES73LN12.40celecoxibM. musculus[14]
ES84RS02.81S-ibuprofenM. musculus[15]
ES94PH91.81rac-ibuprofenM. musculus[16]
aES = Enzyme Structure. b Protein Data Bank (PDB) code. c Reference including the active site.
Table 2. Global mean values of Vina scores, relative standard deviation percentage (%RSD), and root mean square deviation (RMSD) (in Å) of 2-arylbenzofurans docked with cyclooxygenase-2 (COX-2) structures.
Table 2. Global mean values of Vina scores, relative standard deviation percentage (%RSD), and root mean square deviation (RMSD) (in Å) of 2-arylbenzofurans docked with cyclooxygenase-2 (COX-2) structures.
No.GMVS a%RSD bRMSD cNo.GMVS a%RSD bRMSD cNo.GMVS a%RSD bRMSD c
1−8.63 ± 0.505.850.52123−9.37 ± 0.838.880.51345−9.21 ± 0.9310.060.335
2−9.19 ± 0.849.140.30724−9.23 ± 0.899.670.62646−8.79 ± 1.1112.600.475
3−8.10 ± 0.9411.610.31925−8.50 ± 1.3515.850.49347−8.67 ± 1.1413.140.345
4−8.43 ± 0.505.900.71226−8.78 ± 1.184.880.57448−8.79 ± 1.1513.110.333
5−8.50 ± 0.536.250.40427−10.16 ± 0.5013.400.61149−11.26 ± 0.815.480.496
6−8.14 ± 0.749.090.54928−10.47 ± 0.469.380.71150−8.48 ± 1.5818.660.520
7−8.43 ± 0.526.190.53129−10.49 ± 0.434.140.45051−10.37 ± 0.577.160.621
8−8.49 ± 0.394.580.57630−10.42 ± 0.5412.750.79052−8.18 ± 0.9211.300.452
9−8.92 ± 0.384.300.60531−9.49 ± 1.215.210.47153−9.07 ± 0.725.600.467
10−8.51 ± 0.469.030.49032−8.84 ± 1.6618.800.56454−10.39 ± 0.587.950.842
11−10.76 ± 0.275.400.45833−8.32 ± 1.4217.120.31655−8.72 ± 0.495.610.433
12−8.99 ± 0.9810.880.49434−8.63 ± 1.7420.140.52656−8.49 ± 0.455.320.368
13−9.24 ± 1.2012.960.57835−8.70 ± 1.2314.150.54857−7.77 ± 1.3817.810.307
14−8.58 ± 0.576.590.52936−9.02 ± 0.798.710.64658−8.46 ± 0.677.980.367
15−8.93 ± 1.1713.080.60237−9.03 ± 0.586.410.534AA−7.98 ± 0.374.601.329
16−9.09 ± 1.0511.530.36738−8.73 ± 0.586.650.426CEL−11.23 ± 0.625.520.488
17−8.87 ± 0.869.750.51639−10.2 ± 0.696.790.443DCL−8.33 ± 0.303.650.318
18−8.77 ± 1.2414.140.42540−10.3 ± 0.597.420.547IMD−9.67 ± 0.798.130.449
19−9.02 ± 0.717.890.49341−9.70 ± 0.725.740.447NPX−8.66 ± 0.293.320.376
20−9.11 ± 1.2013.180.59642−11.25 ± 0.595.280.335SC558−11.09 ± 0.575.130.412
21−9.33 ± 0.899.550.57443−9.32 ± 0.758.010.322S−IBP−7.71 ± 0.334.220.638
22−9.04 ± 0.849.250.53444−9.13 ± 0.9810.730.320R−IBP−7.71 ± 0.273.520.784
a Global mean values of Vina scores (GMVS) (kcal/mol) ± standard deviation (SD) for each test 2-arylbenzofuran docked with nine different COX-2 structures. b Relative Standard Deviation (RSD) percentages for each test 2-arylbenzofuran. c Mean RMSD values for each test 2-arylbenzofuran.
Table 3. Binding-free and related energies calculated using molecular mechanic/Poisson Boltzmann surface area (MM/PBSA) method for test compounds.
Table 3. Binding-free and related energies calculated using molecular mechanic/Poisson Boltzmann surface area (MM/PBSA) method for test compounds.
LigandsΔEvdW aΔEele bΔGsol cΔGsasa dΔGbind e
11−211.7 ± 2.5−96.3 ± 1.6185.6 ± 7.5−19.0 ± 0.2−141.5 ± 5.8
29−221.2 ± 5.4−109.2 ± 3.9175.6 ± 2.7−20.2 ± 0.3−175.0 ± 4.3
42−205.6 ± 3.5−102.5 ± 3.6132.2 ± 4.8−19.8 ± 0.2−195.7 ± 4.5
49−226.9 ± 4.5−120.8 ± 4.6189.2 ± 3.7−24.4 ± 0.4−182.9 ± 2.0
celecoxib−225.5 ± 8.9−138.2 ± 5.7184.5 ± 8.5−18.6 ± 0.3−197.8 ± 7.5
a Van der Waal (vdW) energy. b Electrostatic energy. c Polar solvation energy. d Solvent accessible surface area (SASA) energy. e Binding-free energy. Energies are expressed as mean values (in kJ/mol) ± standard deviation.
Table 4. Molecular mechanic/Poisson Boltzmann surface area (MM/PBSA) per-residue decomposition energies for selected residues of those complexes between COX-2 and 11, 29, 42, 49, and celecoxib using the last 40-ns molecular dynamics (MD) simulations.
Table 4. Molecular mechanic/Poisson Boltzmann surface area (MM/PBSA) per-residue decomposition energies for selected residues of those complexes between COX-2 and 11, 29, 42, 49, and celecoxib using the last 40-ns molecular dynamics (MD) simulations.
Residue11294249Cel
Val74−0.22 ± 0.42−0.15 ± 0.07−0.08 ± 0.05−4.12 ± 0.76−0.37 ± 0.08
His75−0.57 ± 0.40−2.12 ± 0.58−0.24 ± 0.82−2.26 ± 0.73−1.93 ± 0.37
Leu103−1.07 ± 0.15−1.02 ± 0.26−3.63 ± 0.83−0.77 ± 0.05−0.54 ± 0.04
Arg106−9.64 ± 3.76−0.35 ± 0.51−6.67 ± 0.99−0.36 ± 2.84−8.94 ± 1.33
Gln178−1.77 ± 0.30−1.68 ± 0.32−1.48 ± 0.21−0.53 ± 0.29−4.97 ± 0.38
Val330−2.08 ± 0.06−0.08 ± 0.08−2.13 ± 0.38−0.02 ± 0.02−0.28 ± 0.15
Val335−5.43 ± 0.46−9.27 ± 0.34−11.49 ± 0.72−2.17 ± 0.61−4.65 ± 0.81
Leu338−1.41 ± 0.90−3.29 ± 1.09−4.13 ± 0.74−5.39 ± 0.77−8.87 ± 1.19
Ser339−2.59 ± 0.81−4.39 ± 0.56−4.86 ± 0.39−4.98 ± 0.91−6.37 ± 1.01
Tyr341−0.07 ± 0.63−4.58 ± 0.55−1.56 ± 0.54−1.18 ± 0.75−3.01 ± 0.52
Leu345−2.29 ± 0.41−1.59 ± 0.26−1.42 ± 0.43−0.85 ± 0.15−1.80 ± 0.08
Leu370−0.16 ± 0.04−0.13 ± 0.03−0.17 ± 0.03−1.91 ± 0.15−1.32 ± 0.38
Tyr371−0.08 ± 0.11−0.55 ± 0.16−0.58 ± 0.32−0.46 ± 0.26−1.49 ± 0.92
Arg499−5.32 ± 1.33−3.83 ± 1.39−1.35 ± 1.41−6.24 ± 1.26−7.90 ± 1.17
Phe504−3.12 ± 0.44−2.93 ± 0.91−2.37 ± 0.51−2.41 ± 0.48−6.23 ± 1.03
Val509−6.17 ± 0.58−6.41 ± 1.03−8.11 ± 0.33−7.36 ± 0.91−9.91 ± 0.43
Glu510−3.01 ± 1.03−0.61 ± 0.99−1.27 ± 0.55−13.3 ± 1.76−6.89 ± 0.58
Gly512−0.61 ± 0.17−1.32 ± 0.25−1.07 ± 0.07−3.68 ± 0.56−3.29 ± 0.66
Ala513−3.84 ± 1.03−5.55 ± 0.71−9.86 ± 0.24−7.57 ± 0.37−5.09 ± 0.83
Ser516−2.17 ± 0.91−0.52 ± 0.61−2.48 ± 0.58−2.17 ± 0.71−2.78 ± 0.49
Leu517−1.50 ± 1.24−1.36 ± 0.58−5.96 ± 0.59−1.54 ± 0.28−1.89 ± 0.32
Ligand−60.9 ± 1.5−61.5 ± 1.97−83.4 ± 2.44−69.6 ± 2.24−92.6 ± 3.70
Table 5. In-vitro COX-2 inhibition of an in-house collection (n = 26) of 2-arylbenzofurans.
Table 5. In-vitro COX-2 inhibition of an in-house collection (n = 26) of 2-arylbenzofurans.
NoCOX-2 IC50pIC50 apIC50 bNoCOX-2 IC50pIC50 apIC50 bNoCOX-2 IC50pIC50 apIC50 b
(µM)(exp)(pred)(µM)(exp)(pred)(µM)(exp)(pred)
216.7 c3.784.04293.754.433.94491.254.904.71
313.2 c3.493.35381582.802.865045.63.343.21
828.73.543.72420.7525.125.185110.53.984.14
91192.923.3243 c1022.993.455410.93.963.84
113.27 c4.494.0944 c26.23.583.125512.63.904.17
124.33 c4.363.9845 c52.73.282.9556 c1252.902.51
138.56 c4.074.4646 c32.93.483.5957 c4382.362.21
147.37 c4.134.5547 c71.53.153.5558 c2402.622.72
1712.3 c3.914.2248 c29.33.533.78cel0.102--
apIC50(exp) = −log(IC50(exp) in M), b predicted from CoMFA model. c Compounds previously evaluated for COX-2 inhibitory activity [8,37].

Share and Cite

MDPI and ACS Style

Coy-Barrera, E. Discrimination of Naturally-Occurring 2-Arylbenzofurans as Cyclooxygenase-2 Inhibitors: Insights into the Binding Mode and Enzymatic Inhibitory Activity. Biomolecules 2020, 10, 176. https://doi.org/10.3390/biom10020176

AMA Style

Coy-Barrera E. Discrimination of Naturally-Occurring 2-Arylbenzofurans as Cyclooxygenase-2 Inhibitors: Insights into the Binding Mode and Enzymatic Inhibitory Activity. Biomolecules. 2020; 10(2):176. https://doi.org/10.3390/biom10020176

Chicago/Turabian Style

Coy-Barrera, Ericsson. 2020. "Discrimination of Naturally-Occurring 2-Arylbenzofurans as Cyclooxygenase-2 Inhibitors: Insights into the Binding Mode and Enzymatic Inhibitory Activity" Biomolecules 10, no. 2: 176. https://doi.org/10.3390/biom10020176

APA Style

Coy-Barrera, E. (2020). Discrimination of Naturally-Occurring 2-Arylbenzofurans as Cyclooxygenase-2 Inhibitors: Insights into the Binding Mode and Enzymatic Inhibitory Activity. Biomolecules, 10(2), 176. https://doi.org/10.3390/biom10020176

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop