Next Article in Journal
Zeaxanthin, a Molecule for Photoprotection in Many Different Environments
Next Article in Special Issue
Selenoxide Elimination Triggers Enamine Hydrolysis to Primary and Secondary Amines: A Combined Experimental and Theoretical Investigation
Previous Article in Journal
Strawberry Puree Functionalized with Natural Hydroxytyrosol: Effects on Vitamin C and Antioxidant Activity
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Molecules
Volume 25
Issue 24
10.3390/molecules25245828
molecules-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Molecular Docking Study on Several Benzoic Acid Derivatives against SARS-CoV-2

by
Amalia Stefaniu
*,
Lucia Pirvu
*,
Bujor Albu
and
Lucia Pintilie
National Institute for Chemical-Pharmaceutical Research and Development, 112 Vitan Av., 031299 Bucharest, Romania
*
Authors to whom correspondence should be addressed.
Molecules 2020, 25(24), 5828; https://doi.org/10.3390/molecules25245828
Submission received: 15 November 2020 / Revised: 30 November 2020 / Accepted: 1 December 2020 / Published: 10 December 2020
(This article belongs to the Special Issue From a Molecule to a Drug: Chemical Features Enhancing Pharmacological Potential)
Browse Figures
Review Reports Versions Notes

Abstract

:
Several derivatives of benzoic acid and semisynthetic alkyl gallates were investigated by an in silico approach to evaluate their potential antiviral activity against SARS-CoV-2 main protease. Molecular docking studies were used to predict their binding affinity and interactions with amino acids residues from the active binding site of SARS-CoV-2 main protease, compared to boceprevir. Deep structural insights and quantum chemical reactivity analysis according to Koopmans’ theorem, as a result of density functional theory (DFT) computations, are reported. Additionally, drug-likeness assessment in terms of Lipinski’s and Weber’s rules for pharmaceutical candidates, is provided. The outcomes of docking and key molecular descriptors and properties were forward analyzed by the statistical approach of principal component analysis (PCA) to identify the degree of their correlation. The obtained results suggest two promising candidates for future drug development to fight against the coronavirus infection.

1. Introduction

Severe acute respiratory syndrome coronavirus 2 is an international health matter. Previously unheard research efforts to discover specific treatments are in progress worldwide. Virtual screening of existing compound databases against three protein targets (main protease, RNA dependent RNA polymerase and spike protein) to inhibit coronavirus replication is one of the actual approaches that allows researchers to quickly select best drug candidates for further in vitro assays [1,2,3,4,5].
Boceprevir is a direct-acting antiviral agent (DAA), acting as an inhibitor of NS3/4A, a serine protease enzyme encoded by hepatitis C virus (HCV) [6]. The serine protease enzyme plays a vital role in the replication and cleavage of viral proteins. We found boceprevir as a cocrystallized ligand in a complex with a protein named 3C-like proteinase, the main protease found in coronaviruses, characterized by X-ray diffraction, introduced in the protein data bank with the entry ID: 6WNP [7], at 1.45 Å resolution. The main protease in coronaviruses contains a cysteine-histidine dyad able to achieve catalytic cleavage of the coronavirus polyprotein [8]. Cysteine acts as a nucleophile due to its free electron pair on the sulfur atom, donated to form intramolecular bonds and histidine, respectively, that acts as a general base by its imidazole heterocycle [9]. The main protease of SARS-CoV-2 consists of three domains and a characteristic CYS145-HIS41 dyad in the active site [10].
Due to their low toxicity and antioxidant activity, phenolic acids and flavonoids appear as the most feasible and secure natural antiviral compounds. The potential antiviral activity of vegetal polyphenols is generally based on their capacity to alter virus replication and functional protein synthesis [11]. Another aspect to be mentioned is the capacity of volatile oils, saponins and triterpenic acids, to act as effective solvents and detergents, therefore to solubilize and destroy the lipid layer of the enveloped viruses. Regarding the general chemical aspects associated with antiviral activity, the presence of phenyl ring(s), vinyl and carboxyl moieties and ester, hydroxyl and methoxy groups appear to be the basis of plant phenolics antiviral efficacy [11,12,13]. For example, it is considered that phenolics with free hydroxyl groups interfere with viral adsorption and further cell penetration [11], which is sustained by the fact that high polar phenolics create a protective coating on the cell’s surface. On the other hand, the plant compounds’ bioavailability allows them to reach the circulation and be able to manifest antiviral activity. Accordingly, clinical studies have revealed that gallic acid, catechins, flavones and quercetin glucosides have the best bioavailability in humans [14]; also, data suggests that anthocyanins are fully absorbed in humans [15]. The most potent antiviral compounds proved to interfere with virus replication and/or viral essential protein synthesis are some of the most common vegetal polyphenols, namely quercetin, kaempferol and apigenin, ellagic, rosmarinic and gallic acids, catechin and epicatechin, and various alkyl gallates [11]. Similary, chrysin, acacetin and apigenin inhibited viral transcription of the human rhinovirus 14 [16], while proanthocyanidins from Myrothamnusf labellifolia have proved antiviral activity against herpes simplex virus type-1 by viral adsorption and cell penetration inhibition [17].
Furthermore, molecular docking studies in the last two years, made on dozens of natural and synthesized antiviral compounds, associate their antiviral activity with the capacity of their active groups (phenyl groups and phenyl moieties such as hydroxyl, carbonyl, amino, azo, nitrile, sulfonyl) to bind the active groups of several amino acids found in the active site of SARS-CoV-2 main protease. GLU166, HIS41, ASN142, GLY143, HIS163 and THR 190, are on the top of the most frequent amino acids bound [10,18,19,20].
Recent virtual screening of 22 U.S.-FDA approved antiviral drugs in the parallel with 24 natural plant-based molecules, indicated theaflavin digallate (from green tea) with the best docking score (−10.574), a result obtained using the Glide (Schrödinger) module on the COVID-19 main protease (structure 6LU7) [20]. It must be noted that excepting HIS41, which bound Osp2 from the carboxyl moiety of gallic acid, all other interactions occurred at hydroxyl groups from the flavan ring. Therefore, it is confirmed that the high degree of hydroxylation is associated with a good antiviral effect.
Another aspect to be mentioned is the conclusion of a computational survey to a drug repurposing study claiming that five neutral antiviral drugs have inhibitory activities against SARS-CoV-2 main protease [21]. As it is well known, the logP value describes lipophilicity for neutral compounds. Therefore, the compounds with an octanol-water partition coefficient (logP) value between one to three have good passive absorption across lipid membranes (bioavailability), while those with a logP value greater than three or less than one usually have lesser bioavailability [22]. In this context, further calculations can be made in connection with the formulation of antiviral products to obtain targeted bioavailability in the human intestine.
Gallic acid, one of the most abundant phenolic acids in plants, has various health benefits including anti-inflammatory, antioxidant and antimicrobial activities [23] and is viewed as the lead compound with promising pharmacological properties to design and develop new drugs [24]. Regarding semisynthetic derivatives of gallic acid, the alkyl gallates with antimicrobial activity were reported too [25,26]. Their biological activity is associated with the length with their alkyl side chain, which affects membrane binding capability [25]. Generally, membrane binding of the alkyl gallates increases with increasing alkyl chain length and is correlated with their antiviral activity.
Starting from such intriguing premises, in molecular docking approach we conducted in silico screening on several benzoic acid derivatives and on a homologue series of alkyl gallates, starting from the lead compound gallic acid, against SARS-CoV-2 main protease.

2. Results and Discussion

We selected the following compounds: benzoic acid, 4-aminobenzoic acid, 4-hydroxybenzoic acid, 3,4,5-trimethoxybenzoic acid (eudesmic acid), 3,4-dihydroxybenzoic acid (protocatehuic acid), 2,5-dihydroxybenzoic acid (gentisic acid), 4-hydroxy-3-methoxybenzaldehyde (vanillin), 4-hydroxy-3,5-dimethoxybenzoic acid (syringic acid), 4,5-dihydroxy-3-oxocyclohex-1-ene-1-carboxylic acid, epicatechin, 3,4,5-trihydroxybenzoic acid (gallic acid), methyl 3,4,5-trihydroxybenzoate, ethyl 3,4,5-trihydroxybenzoate, propyl 3,4,5-trihydroxybenzoate, isopropyl 3,4,5-trihydroxybenzoate, butyl 3,4,5-trihydroxybenzoate, isobutyl 3,4,5-trihydroxybenzoate, pentyl 3,4,5-trihydroxybenzoate, isopentyl 3,4,5-trihydroxybenzoate, octyl 3,4,5-trihydroxybenzoate. Their 3D optimized structures obtained with the Spartan program, along with the atomic numbering scheme, arbitrary chosen by the software, are given in Figure 1.

2.1. Results of Molecular Docking Simulations

Intermolecular interactions occurring in boceprevir, benzoic acid and alkyl gallate derivatives in complex with the 6WNP protein fragment were identified. The lengths of hydrogen-bonding interactions were measured. The results are given in terms of docking score function and RMSD (root mean square deviation).
To validate the molecular docking protocol, boceprevir was initially docked into the crystal structure of the main protease fragment and its interactions with the target 3C-like proteinase were analyzed. As can be seen Figure 2a, the native ligand forms eight hydrogen bonding interactions with residues CYS145, SER144, GLY143, HIS41, HIS164, GLU166, GLU166 and GLU166. In Figure 2b, the superposition of the binding pose of boceprevir, obtained by redocking, is displayed. As illustrated in Figure 2c, all docked ligands were found to have similar binding poses to the native ligand, thus validating the chosen docking approach.
Figure 3 reveals the obtained docking scores for the cocrystallized and selected ligands. Boceprevir exhibits the greatest score (−63.95), due to its numerous interactions with the amino acid residues from the protein’s active binding site, i.e., eight hydrogen bonding interactions with nitrogen or oxygen atoms of amino acids residues N sp2 CYS145 (2.900 Å), N sp2 SER144 (3.053 Å), N sp2 GLY143 (2.783 Å), N sp2 HIS41 (2.604 Å), O sp2 HIS164 (3.103 Å), N sp2 GLU166 (3.118 Å), O sp2 GLU166 (2.908 Å and O sp2 GLU166 (3.229 Å). Remarkable is the docking score for octyl-gallate (−60.22, RMSD 1.12), close to boceprevir, and for epicatechin (−49.57). Benzoic acid derivatives, in the first half of the graph exhibit moderate scores ranging between −29.59 (benzoic acid) and −37.25 (syringic acid). We observed an increasing trend of score with the increasing number of hydroxyl groups. For instance, gallic acid docking simulations resulted in a −38.31 score. Regarding the gallates, the increasing length of alkyl chain led to docking scores in the order gallic acid < methyl gallate < ethyl gallate < isopropyl gallate < propyl gallate < butyl gallate < isobutyl gallate < isopentyl gallate < pentyl gallate < octyl gallate. The ramification of lateral side chain led to a slight decrease of binding affinity, noticeable with an increasing number of -CH2 (e.g., isopentyl: −46.17 versus pentyl gallate: −48.77).
In Table S1 from the Supplementary files available online are listed the obtained values for docking score and RMSD (root mean square deviation), the amino acids group interactions and type and length of hydrogen-bonding interactions formed by each ligand in complex with SARS-CoV-2 main protease.
Figure 4 illustrates the intramolecular interactions by H-bonding (a), and amino acid group interactions occurring in the complex formed by octyl gallate and the 6WNP main protease fragment.
Amino acid residues CYS145 and SER144 interacted by hydrogen-bonding both with boceprevir and with octyl-gallate (see Figure 4a); similar interactions resulted in similar docking results. Propyl gallate and pentyl gallate revealed the same ten intramolecular interactions with N sp2 GLU166, N sp2 HIS163, O sp3 SER144, O sp3 SER144, O sp2 LEU141, N sp2 SER144, N sp2 GLY143, N sp2 CYS145, N sp2 GLY143 and O sp2 ASN142. Their obtained docking scores were −42.13 and −48.77, respectively; lower than for octyl-gallate but greater than the results for gallic acid (−38.31). These observations are in good agreement with experimental findings of Takai E. et al., 2019 [25], anticipating increasing antiviral effects of alkyl gallates with increasing alkyl chain length, except for cetyl and stearyl gallate (which are not included in this study), a fact experimentally validated by fluorescence analysis of the binding of alkyl gallates to phospholipid membranes. Beyond a certain alkyl chain length (8–11), a reduction of antibacterial and antiviral activities of the alkyl gallates was observed, probably due to a self-association process [25,26,27]. This was the reason for choosing to break off the gallates screening at octyl. Increasing docking results were observed with increasing length of the n-alkyl side chain. The good docking result for octyl gallate recommends it as good alternative for developing new therapeutic antiviral agents.

2.2. Results of Oral Bioavailability Evaluation

In Table 1 are listed key molecular descriptors and properties to evaluate the oral bioavailability [28,29] and Veber’s [30] rules, where: MW is the molecular weight that should be less than 500 Daltons, HBD is the number of hydrogen bond donors (recommended to be lower than 5), HBD is the number of hydrogen bond acceptors with acceptable values less than 10 and log P is the water-octanol partition coefficient, that should be less than 5. Veber D.F. et al., 2002 [30] imposed additional restrictions related with the molecular descriptor PSA (polar surface area), namely, no larger than 140 Å2 and with a maximum 10 rotatable bonds for good oral bioavailability.
Therefore, the boceprevir antiviral exhibited two violations of Lipinski’s criteria, namely molecular mass (521.69 g mol−1) and six hydrogen bond donors. The structure of boceprevir presents the maximum allowed number of flexible bonds (10) and maximum hydrogen bond acceptors (10). Although it had these exceptions, the docking score was the highest among the investigated ligands, suggesting strong interactions and stability of the complex formed with the SARS-CoV-2 main protease. The results indicated that all proposed ligands met the restrictive criteria for good oral bioavailability. Increased hydrophilicity was observed for all compounds due to their hydroxyl, carboxyl and/or methoxy groups on their skeleton. These functional groups offer good absorption and permeation properties. Thus, by means of NH/OH/N/O groups, hydrophilic interactions were favored and further reflected in good and high docking scores. Concerning logP values, there were observed positive values for benzoic acid (0.79), 2,5-dihydroxybenzoic acid (0.81) and 0.72 for octyl 3,4,5-trihydroxy benzoate. A combination of molecular factors and properties, mainly due to the increased hydrophobicity of the lateral n-octyl chain of octyl gallate, was also reflected by its best docking score function, indicating this compound as the best antiviral candidate among all screened compounds in the study.

2.3. Results of Quantum Reactivty Analysis

Frontier molecular orbitals, the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) localization and energy levels for octyl-gallate, are illustrated in Figure 5.
The resulting band gap (ΔE) can provide useful information on the chemical reactivity and kinetic stability of each ligand. The same values for energy gap were given for alkyl gallates, starting at ethyl to octyl-gallate, suggesting the same stability. Slight differences were found for the values of ionization potential (I = −EHOMO) and electron affinity (A = −ELUMO), according Koopmans’ theorem [31]. The theorem allows estimation of quantum global reactivity parameters, starting from calculated energies of frontier molecular orbitals, and describes the molecules in terms of chemical hardness (η), global softness (σ), ionization (I), electron affinity (A), electronegativity (χ) and electrophyliciy index (ω) [32,33]. Obtained quantum reactivity parameters for all investigated ligands are given in Table S2 of the Supplementary Materials.
The global reactivity parameters analysis provides deep structural insights, a holistic characterization for revealing the properties of interest leading to strong binding affinity related to the protease target. The quantum reactivity parameters are related to relative nucleophilicity and electrophilicity. Ionization potential (I), refers to the energy needed to remove an electron from a molecule, and electron affinity (A) measures the ability of a molecule to accept electrons and form anions species [34,35]. Such parameters are useful to estimate further chemical reactivity behavior. Some of the investigated molecules can also be seen as lead compounds for a new series of (semi)synthetic molecules. Therefore the data on their reactivity are useful.

2.4. Results of Principal Component Analysis (PCA)

Principal component analysis (PCA) is a statistical tool for the identification of linear combinations of the variables which account for certain proportions of the variance of a set of variables. The selection is based on the eigenvalues of the dispersion matrix of the variables. The principal components are associated with decreasing eigenvalues and, therefore, share the amount of the variance. Typically, the first few principal components account for virtually all the variance. PCA also represents the pattern of similarity of the observations and the variables by displaying them as points in maps [36,37,38,39]. PCA analysis of all properties was calculated with Spartan software, along with docking scores, and data are listed in Table S4 of the Supplementary Materials.
PCA analysis was employed to find the degree of correlation between molecular descriptors and properties and their involvement in the resulting docking score.
The PCA correlation matrix (Table S4 in Supplementary Materials) revealed fairly good correlations between area and mass (r = 0.95), area and ovality (r = 0.96), docking score and polarizability (r = 0.97), volume and area (r = 0.97), volume and ovality (r = 0.98) and a moderate correlation (r = 0.55–0.66) between the dipole moment and docking score, mass, area, volume, and ovality, and between PSA and polarizability and mass, respectively.
Table 2 and Figure 6 are related to the eigenvalues which reflect the quality of the projection from the N–dimensional initial (n = 6) to a lower number of dimensions.
From Table 2 it can be seen that the first eigenvalue equals 6.529 and represents 65.29 % of the total variability. Each eigenvalue corresponds to a factor, and each factor to a one dimension. A factor is a linear combination of the initial variables and all the factors are uncorrelated (r = 0). The eigenvalues and the corresponding factors are sorted by descending order of how much of the initial variability they represent (converted to %). Therefore, the first two factors allow us to represent 70.69% of the initial variability of the data.

3. Methods

3.1. Methods for Molecular Docking Simulations

The docking simulation was carried out using CLC Drug, Discovery Work Bench (QIAGEN, Aarhus, Denmark,). SARS-CoV-2 main protease bound to boceprevir at 1.45 Å (PDB ID: 6WNP) [7], which was imported from the Protein Data Bank. Ligands were constructed using the Spartan’16 program [40,41] and their geometries were optimized to obtain the lowest energy conformers. The viral protein fragment contains three binding pockets: 48.13 Å3, 40.45 Å3, 36.62 Å3. Redocking of the cocrystallized ligand (boceprevir) was realized to validate the docking protocol. The amino acid residues forming the binding site were protonated and water molecules were removed.

3.2. Methods for Energy Minimization and Computation of Molecular Properties

The examined ligands were first generated in 3D by importing their corresponding files from the Pubchem database (National Center for Biotechnology Information, U.S. National Library of Medicine Rockville Pike, Bethesda, MD, USA, https://pubchem.ncbi.nlm.nih.gov) or directly constructed (long chain alkyl gallates) in the Spartan’18 program [40,41]. Their geometry was optimized in a multistep procedure by molecular mechanics force fields, developed at Merck Pharmaceuticals (Merck Research Labs, Kenilworth, NJ, USA), to realize energy minimization [42]. Molecular descriptor and properties were calculated using the DFT, B3LYP algorithm [43] and the 6-31G (d, p) polarization basis set [40]. Calculations were performed for equilibrium geometry at ground state in gas for neutral molecules. No solvent corrections were done.

3.3. Methods for Principal Component Analysis (PCA)

The predicted ligand data were processed to compute principal component analysis using the free Excel add-in Real Statistics Resource Pack software (Release 7.2, www.real-statistics.com), Copyright (2013–2020) Charles Zaiontz (www.real-statistics.com).

4. Conclusions

In this work we analyzed and predicted the antiprotease activity of natural derivatives of benzoic and semisynthetic alkyl gallate acids on SARS-CoV-2 main protease. The investigation was corroborated with drugability and quantum reactivity evaluations. The docking results of the two studied series (benzoic acid series versus gallic acid series) suggested 2,5-dihydroxybenzoic acid (gentisic acid) and octyl gallate as the best potential candidates among the investigated structures. The two compounds had similar logP values (0.81 and 0.72). Octyl gallate had the best docking score (−60.22), but decreased dipole moment (1.31). 2,5-dihydroxybenzoic acid (gentisic acid) had a lower score (−33.84) but an increased dipole moment (4.62) which means higher polarity and also higher reactivity. As is known, the dipole moment represents a measure of net molecular polarity and, therefore, the larger the difference in electronegativities of bonded atoms the larger the reactivity of the molecule. Accordingly, a combination of the two compounds can be considered. On the other hand, if the benzoic acid series indicated that the position of hydroxyl groups on the ring was more important than the ester, hydroxyl and methoxy groups’ number or type, the galllic acid series clearly indicated similar results for all studied compounds, apart from octyl gallate which was proved to be the best potential candidate among investigated structures, exhibiting antiviral activity against the SARS-CoV-2 protease and, therefore, could be studied further for drug development.

Supplementary Materials

The following are available online. Table S1: The list of intermolecular interactions between the compounds docked with 6WNP, Table S2: Quantum chemical reactivity parameters calculated with Koopman’s relationships, Table S3: Molecular properties for the investigated ligands, calculated with Spartan Software, Table S4: Correlation matrix of PCA.

Author Contributions

Conceptualization, A.S. and L.P. (Lucia Pirvu); methodology, A.S.; software, L.P. (Lucia Pintilie) and B.A.; validation, A.S., B.A. and L.P. (Lucia Pirvu); formal analysis, L.P. (Lucia Pintilie); investigation, L.P. (Lucia Pintilie) and B.A.; resources, L.P. (Lucia Pintilie) and B.A.; data curation, L.P. (Lucia Pirvu); writing—original draft preparation, A.S.; writing—review and editing, L.P. (Lucia Pirvu); visualization, A.S.; supervision, L.P. (Lucia Pirvu); project administration, L.P. (Lucia Pintilie). All authors have read and agreed to the published version of the manuscript.

Funding

This research work was funded by the Romanian National Authority for Scientific Research ANCS, Competitiveness Operational Programme COP-A1-A1.2.3-G-2015, Project title “Innovative technologies for new, natural health products”, ID P_40_406, SMIS 105542.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

AA-electron affinity
B3LYPBecke, 3-parameter, Lee-Yang-Parr
DAAdirect-acting antiviral agent
DFTDensity Functional Theory
EHOMOenergy of the highest occupied molecular orbital
ELUMOenergy of the lowest unoccupied molecular orbital
FDAFood and Drug Administration (U.S.)
HBAhydrogen bond acceptor
HBDhydrogen bond donor
HCVhepatitis C virus
HOMOthe highest occupied molecular orbital
Iionization potential
logPoctanol-water partition coefficient
LUMOthe highest occupied molecular orbital
LVLipinski’s violations
MWMolecular weight
PCAPrincipal component analysis
PSApolar surface area
rPearson correlation coefficient
rbrotatable bond
RMSDRoot mean square deviation
RNAribonucleic acid
SARSCoV-2-Severe acute respiratory syndrome coronavirus 2
ΔEgapenergy gap between frontier molecular orbitals
µchemical potential
ηglobal chemical hardness
σglobal softness
χelectronegativity
ωglobal electrophilicity index

References

  1. Vijayakumar, B.G.; Ramesh, D.; Joji, A.; Prakasan, J.J.; Kannan, T. In silico pharmacokinetic and molecular docking studies of natural flavonoids and synthetic indole chalcones against essential proteins of SARS-CoV-2 proteins. Eur. J. Pharmacol. 2020, 886, 173448. [Google Scholar] [CrossRef] [PubMed]
  2. Cherrak, S.A.; Merzouk, H.; Mokhtari-Soulimane, N. Potential bioactive glycosylated flavonoids as SARS-CoV-2 main protease inhibitors: A molecular docking and simulation studies. PLoS ONE 2020, 15, e0240653. [Google Scholar] [CrossRef] [PubMed]
  3. Kumar, C.S.; Ali, D.; Alarifi, S.; Radhakrishnan, S.; Akbar, I. In silico molecular docking: Evaluation of coumarin based derivatives against SARS-CoV-2. J. Infect. Public Health 2020, 13, 1671–1677. [Google Scholar]
  4. El-hoshoudy, A.N. Investigating the potential antiviral activity drugs against SARS-CoV-2 by molecular docking simulation. J. Mol. Liq. 2020, 318, 113968. [Google Scholar] [CrossRef] [PubMed]
  5. Narkhede, R.R.; Pise, A.V.; Cheke, R.S.; Shinde, S.D. Recognition of Natural Products as Potential Inhibitors of COVID-19 Main Protease (Mpro): In-Silico Evidences. Nat. Prod. Bioprospect. 2020, 10, 297–306. [Google Scholar] [CrossRef] [PubMed]
  6. Kiser, J.J.; Flexner, C. Direct-acting antiviral agents for hepatitis C virus infection. Annu. Rev. Pharmacol. Toxicol. 2013, 53, 427–449. [Google Scholar] [CrossRef] [Green Version]
  7. Anson, B.; Mesecar, A. 6WNP X-ray Structure of SARS-CoV-2 Main Protease Bound to Boceprevir at 1.45 A. Available online: https://www.rcsb.org/structure/6WNP (accessed on 3 December 2020).
  8. Fan, K.; Wei, P.; Feng, Q.; Chen, S.; Huang, C.; Ma, L.; Lai, B.; Pei, J.; Liu, Y.; Chen, J.; et al. Biosynthesis, purification, and substrate specificity of severe acute respiratory syndrome coronavirus 3C-like proteinase. J. Biol. Chem. 2004, 279, 1637–1642. [Google Scholar] [CrossRef] [Green Version]
  9. Ryu, Y.B.; Park, S.J.; Kim, Y.M.; Lee, J.Y.; Seo, W.D.; Chang, J.S.; Park, K.H.; Rho, M.C.; Lee, W.S. SARS-CoV 3CLpro inhibitory effects of quinone-methide triterpenes from Tripterygium regelii. Bioorg. Med. Chem. Lett. 2010, 20, 1873–1876. [Google Scholar] [CrossRef]
  10. Kumar, Y.; Singh, H.; Patel, C.N. In silico prediction of potential inhibitors for the Main protease of SARS-CoV-2 using molecular docking and dynamics simulation-based drug-repurposing. J. Infect. Public Health 2020. [Google Scholar] [CrossRef]
  11. Kamboj, A.; Saluja, A.K.; Kumar, M.; Atri, P. Antiviral activity of plant polyphenols. J. Pharm. Res. 2012, 5, 2402–2412. [Google Scholar]
  12. Park, E.S.; Moon, W.S.; Song, M.J.; Kim, M.N.; Chung, K.H.; Yoon, J.S. Antimicrobial activity of phenol and benzoic acid derivatives. Int. Biodeterior. Biodegrad. 2001, 47, 209–214. [Google Scholar] [CrossRef]
  13. Kratky, M.; Konecna, K.; Janousek, J.; Brablikova, M.; Jand’ourek, O.; Trejtnar, F.; Stolarikova, J.; Vinsova, J. 4-Aminobenzoic Acid Derivatives: Converting Folate Precursor to Antimicrobial and Cytotoxic Agents. Biomolecules 2019, 10, 9. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Manach, C.; Scalbert, A.; Morand, C.; Rémésy, C.; Jimenez, L. Polyphenols: Food sources and bioavailability. Am. J. Clin. Nutr. 2004, 79, 727–747. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Passamonti, S.; Terdoslavich, M.; Franca, R.; Vanzo, A.; Tramer, F.; Braidot, E.; Petrussa, E.; Vianello, A. Bioavailability of flavonoids: A review of their membrane transport and the function of bilitranslocase in animal and plant organisms. Curr. Drug Metab. 2009, 10, 369–394. [Google Scholar] [CrossRef] [Green Version]
  16. Smith, T.J.; Kremer, M.J.; Luo, M.; Vriend, G.; Arnold, E.; Kamer, G.; Rossman, M.G.; McKinlay, M.A.; Diana, G.D.; Otto, M.J. The site of attachment in human rhinovirus 14 for antiviral agents that inhibit uncoating. Science 1986, 233, 1286–1293. [Google Scholar] [CrossRef]
  17. Kirsten, G.; Joachim, K.; Eva, L.; Wali, H.; Derksena, A.; Detersa, A.; Andreas, H. Proanthocyanidin-enriched extract from Myrothamnus flabellifolia Welw. exerts antiviral activity against herpes simplex virus type-1 by inhibition of viral adsorption and penetration. J. Ethnopharmacol. 2011, 134, 468–474. [Google Scholar]
  18. Shah, B.; Modi, P.; Sagar, S.R. In silico studies on therapeutic agents for COVID-19: Drug repurposing approach. Life Sci. 2020, 252, 117652. [Google Scholar] [CrossRef]
  19. Cardoso, W.B.; Mendanha, S.A. Molecular dynamics simulation of docking structures of SARS-CoV-2 main protease and HIV protease inhibitors. J. Mol. Struct. 2021, 1225, 129143. [Google Scholar] [CrossRef]
  20. Peele, K.A.; Durthi, C.P.; Srihansa, T.; Krupanidhi, S.; Sai, A.V.; Babu, D.J.; Indira, M.; Reddy, A.R.; Venkateswarulu, T.C. Molecular docking and dynamic simulations for antiviral compounds against SARS-CoV-2: A computational study. Inform. Med. Unlocked. 2020, 19, 100345. [Google Scholar] [CrossRef]
  21. Wang, J. Fast identification of possible drug treatment of coronavirus disease-19 (COVID-19) through computational drug repurposing study. J. Chem. Inf. Model. 2020, 60, 3277–3286. [Google Scholar] [CrossRef]
  22. Hetal, T.; Bindesh, P.; Sneha, T. A review on techniques for oral bioavailability enhancement of drugs. Int. J. Pharm. Sci. Rev. Res. 2010, 4, 203–223. [Google Scholar]
  23. Kahkeshani, N.; Farzaei, F.; Fotouhi, M.; Alavi, S.S.; Bahramsoltani, R.; Naseri, R.; Momtaz, S.; Abbasabadi, Z.; Rahimi, R.; Farzaei, M.H.; et al. Pharmacological effects of gallic acid in health and diseases: A mechanistic review. Iran J. Basic Med. Sci. 2019, 22, 225–237. [Google Scholar] [CrossRef] [PubMed]
  24. Nayeem, N.; Asdaq, S.M.B.; Salem, H.; Alfqy, S.A. Gallic Acid: A Promising Lead Molecule for Drug Development. J. Appl. Pharm. 2016, 8, 1000213. [Google Scholar] [CrossRef] [Green Version]
  25. Takai, E.; Hirano, A.; Shirak, K. Effects of alkyl chain length of gallate on self-association and membrane binding. J. Biochem. 2011, 150, 165–171. [Google Scholar] [CrossRef] [PubMed]
  26. Król, E.; de Sousa Borges, A.; da Silva, I.; Polaquini, C.R.; Regasini, L.O.; Ferreira, H.; Scheffers, D.J. Antibacterial activity of alkyl gallates is a combination of direct targeting of FtsZ and permeabilization of bacterial membranes. Front. Microbiol. 2015. [Google Scholar] [CrossRef]
  27. Fujita, K.; Kubo, I. Plasma membrane injury induced by nonyl gallate in Saccharomyces cerevisiae. J. Appl. Microbiol. 2002, 92, 1035–1042. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Lipinski, C.A. Lead-and drug-like compounds: The rule-of-five revolution. Drug Discov. Today Technol. 2004, 1, 337–341. [Google Scholar] [CrossRef] [PubMed]
  29. Lipinski, C.A.; Lombardo, F.; Dominy, B.W.; Feeney, P.J. Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings. Adv. Drug Deliv. Rev. 2001, 46, 3–26. [Google Scholar] [CrossRef]
  30. Veber, D.F.; Johnson, S.R.; Cheng, H.Y.; Smith, B.R.; Ward, K.W.; Kopple, K.D. Molecular properties that influence the oral bioavailability of drug candidates. J. Med. Chem. 2002, 45, 2615–2623. [Google Scholar] [CrossRef]
  31. Koopmans, T. Über die Zuordnung von Wellenfunktionen und Eigenwerten zu den Einzelnen Elektronen Eines Atoms. Physica 1934, 1, 104–113. [Google Scholar] [CrossRef]
  32. Sastri, V.S.; Perumareddi, J.R. Molecular orbital theoretical studies of some organic corrosion inhibitors. Corrosion 1997, 53, 617–622. [Google Scholar] [CrossRef]
  33. Yankova, R.; Genieva, S.; Halachev, N.; Dimitrova, G. Molecular structure, vibrational spectra, MEP, HOMO-LUMO and NBO analysis of Hf(SeO3)(SeO4)(H2O)4. J. Mol. Struct. 2016, 1106, 82–88. [Google Scholar] [CrossRef]
  34. Landeros-Martinez, L.L.; Orrantia-Borunda, E.; Flores-Holguin, N. DFT Chemical Reactivity analysis of biological molecules in the presence of silver ion. Org. Chem. Curr. Res. 2015, 4, 153. [Google Scholar] [CrossRef]
  35. Pirvu, L.; Neagu, G.; Terchescu, I.; Albu, B.; Stefaniu, A. Comparative studies of two vegetal extracts from Stokesia laevis and Geranium pratense: Polyphenol profile, cytotoxic effect and antiproliferative activity. Open Chem. 2020, 18, 488–502. [Google Scholar] [CrossRef]
  36. Jolliffe, I.T. Principal Component Analysis; Springer: New York, NY, USA, 2002; ISBN 978-0-387-22440-4. [Google Scholar]
  37. Jackson, J.E. A Use’s Guide to Principal Components; John Wiley & Sons: New York, NY, USA, 1991; ISBN 9780471725336. [Google Scholar] [CrossRef]
  38. Saporta, G.; Niang, N. Chapter 1: Principal component analysis: Application to statistical process control. In Data Analysis; Govaert, G., Ed.; John Wiley & Sons: London, UK, 2009; pp. 1–23. [Google Scholar]
  39. Abdi, H.; Williams, L.J. Principal component analysis. WIREs Comput. Stat. 2010, 2, 433–459. [Google Scholar] [CrossRef]
  40. Shao, Y.; Molnar, L.F.; Jung, Y.; Kussmann, J.; Ochsenfeld, C.; Brown, S.T.; Gilbert, A.T.B.; Slipchenko, L.V.; Levchenko, S.V.; O’Neill, D.P.; et al. Advances in methods and algorithms in a modern quantum chemistry program package. Phys. Chem. Chem. Phys. 2006, 8, 3172–3191. [Google Scholar] [CrossRef]
  41. Hehre, W.J. A Guide to Molecular Mechanics and Quantum Chemical Calculations; Wavefunction, Inc.: Irvine, CA, USA, 2003. [Google Scholar]
  42. Halgren, A.T. Merck molecular force field. I. Basis, form, scope, parameterization, and performance of MMFF94. J. Comput. Chem. 1996, 17, 490–519. [Google Scholar] [CrossRef]
  43. Lee, C.; Yang, W.; Parr, R.G. Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. Phys. Rev. B 1988, 37, 785–789. [Google Scholar] [CrossRef] [Green Version]
Molecules 25 05828 g001a 550Molecules 25 05828 g001b 550
Figure 1. 3D optimized structures of investigated ligands and their numbering atomic labels.
Figure 1. 3D optimized structures of investigated ligands and their numbering atomic labels.
Molecules 25 05828 g001aMolecules 25 05828 g001b
Molecules 25 05828 g002 550
Figure 2. Hydrogen bonding interactions of the native ligand (boceprevir); (a) superposition of the native ligand; (b) superposition of all docked ligands in the active binding site of 3C-like proteinase (6WNP) (c).
Figure 2. Hydrogen bonding interactions of the native ligand (boceprevir); (a) superposition of the native ligand; (b) superposition of all docked ligands in the active binding site of 3C-like proteinase (6WNP) (c).
Molecules 25 05828 g002
Molecules 25 05828 g003 550
Figure 3. Docking scores for investigated ligands against SARS-CoV-2 main protease (6WNP).
Figure 3. Docking scores for investigated ligands against SARS-CoV-2 main protease (6WNP).
Molecules 25 05828 g003
Molecules 25 05828 g004 550
Figure 4. Hydrogen bonding interactions of octyl-gallate (a) and interactions with amino acid residues from the active binding site of 6WNP protein fragment (b).
Figure 4. Hydrogen bonding interactions of octyl-gallate (a) and interactions with amino acid residues from the active binding site of 6WNP protein fragment (b).
Molecules 25 05828 g004
Molecules 25 05828 g005 550
Figure 5. Highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) molecular frontier orbitals and their energy gap for octyl-gallate.4
Figure 5. Highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) molecular frontier orbitals and their energy gap for octyl-gallate.4
Molecules 25 05828 g005
Molecules 25 05828 g006 550
Figure 6. Results of PCA analysis: scree plot of the eigenvalues and cumulative variability versus the F1-F6 components (a) and PC1 PC2 PC3 3D Loading plot (b). The scree plot (Figure 6a) is a useful visual aid for determining the number of the principal components, which depends on the elbow point at which the remaining eigenvalues are relatively small and all about the same size. This point is not so evident in the scree plot, but we may say that our elbow point is the third point. In conclusion Plot 4 and the Table 2 indicate that the first three PCs are sufficient to explain most of the variance (more than 90.56%) of the data set without overfitting the model. Detailed information about the best principal component locations of the extracted eigenvectors PC1-PC2-PC3 3D loading plots are also presented (Figure 6b). It is worth mentioning that the dipole moment is located in the PC2 × PC3 space, PSA in the PC1 × PC2 space and the other variables in the PC1 × PC3 space.
Figure 6. Results of PCA analysis: scree plot of the eigenvalues and cumulative variability versus the F1-F6 components (a) and PC1 PC2 PC3 3D Loading plot (b). The scree plot (Figure 6a) is a useful visual aid for determining the number of the principal components, which depends on the elbow point at which the remaining eigenvalues are relatively small and all about the same size. This point is not so evident in the scree plot, but we may say that our elbow point is the third point. In conclusion Plot 4 and the Table 2 indicate that the first three PCs are sufficient to explain most of the variance (more than 90.56%) of the data set without overfitting the model. Detailed information about the best principal component locations of the extracted eigenvectors PC1-PC2-PC3 3D loading plots are also presented (Figure 6b). It is worth mentioning that the dipole moment is located in the PC2 × PC3 space, PSA in the PC1 × PC2 space and the other variables in the PC1 × PC3 space.
Molecules 25 05828 g006
Table
Table 1. Lipinski and Veber’s parameters for drugability assessment.
Table 1. Lipinski and Veber’s parameters for drugability assessment.
LigandMWPSAHBDHBALogPrbLV
Benzoic acid122.12333.690120.7910
4-Aminobenzoic acid137.13858.47133−0.9310
4-Hydroxybenzoic acid138.12253.44423−0.2910
3,4,5-Trimethoxybenzoic acid212.20153.22315−2.1440
3,4-Dihydroxybenzoic acid154.12171.21734−1.3710
2,5-Dihydroxybenzoic acid154.12171.262340.8110
4-Hydroxy-3-methoxybenzaldehyde152.14941.01213−1.5320
4-Hydroxy-3,5-dimethoxybenzoic acid198.17464.70625−2.2430
4,5-Dihydroxy-3-oxocyclohex-1-ene-1-carboxylic acid172.13683.67135−0.9210
Epicatechin290.271101.29456−3.7210
3,4,5-Trihydroxybenzoic acid170.1289.40845−2.4610
Methyl 3,4,5-trihydroxybenzoate184.14775.75235−2.1920
Ethyl 3,4,5-trihydroxybenzoate198.17475.42535−1.8630
Propyl 3,4,5-trihydroxybenzoate212.20175.43335−1.3740
i-Propyl 3,4,5-trihydroxybenzoate212.20175.06835−1.5430
Butyl 3,4,5-trihydroxybenzoate226.22875.43335−0.9550
i-Butyl 3,4,5-trihydroxybenzoate226.22875.14935−0.9740
Pentyl 3,4,5-trihydroxybenzoate240.25575.42635−0.5460
i-Pentyl 3,4,5-trihydroxybenzoate240.25575.41535−0.6250
Octyl 3,4,5-trihydroxybenzoate282.33675.390350.7290
MW—molecular weight (g mol−1); PSA—polar surface area (Å2) HBD—hydrogen bond donor; HBA—hydrogen bond acceptor; rb—rotatable bonds count; LV—Lipinski’s violations.
Table
Table 2. Eigenvalues from the Principal Component Analysis (PCA) analysis.
Table 2. Eigenvalues from the Principal Component Analysis (PCA) analysis.
F1F2F3F4F5F6
Eigenvalue6.5291.4401.0870.4800.2910.084
Variability, %65.29%14.40%10.87%4.80%2.91%0.84%
Cumulative, %65.29%79.69%90.56%95.36%98.27%99.11%
Sample Availability: Samples of the compounds are not available from the authors.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Stefaniu, A.; Pirvu, L.; Albu, B.; Pintilie, L. Molecular Docking Study on Several Benzoic Acid Derivatives against SARS-CoV-2. Molecules 2020, 25, 5828. https://doi.org/10.3390/molecules25245828

AMA Style

Stefaniu A, Pirvu L, Albu B, Pintilie L. Molecular Docking Study on Several Benzoic Acid Derivatives against SARS-CoV-2. Molecules. 2020; 25(24):5828. https://doi.org/10.3390/molecules25245828

Chicago/Turabian Style

Stefaniu, Amalia, Lucia Pirvu, Bujor Albu, and Lucia Pintilie. 2020. "Molecular Docking Study on Several Benzoic Acid Derivatives against SARS-CoV-2" Molecules 25, no. 24: 5828. https://doi.org/10.3390/molecules25245828

APA Style

Stefaniu, A., Pirvu, L., Albu, B., & Pintilie, L. (2020). Molecular Docking Study on Several Benzoic Acid Derivatives against SARS-CoV-2. Molecules, 25(24), 5828. https://doi.org/10.3390/molecules25245828

Article Metrics

Back to TopTop