Bioactivity and Molecular Docking Studies of Derivatives from Cinnamic and Benzoic Acids

Perez-Castillo, Yunierkis; Lima, Tamires C.; Ferreira, Alana R.; Silva, Cecília R.; Campos, Rosana S.; Neto, João B. A.; Magalhães, Hemerson I. F.; Cavalcanti, Bruno C.; Júnior, Hélio V. N.; de Sousa, Damião P.

doi:https://doi.org/10.1155/2020/6345429

BioMed Research International

On this page

Abstract Introduction Materials and Methods Results and Discussion Conclusions Data Availability Conflicts of Interest Acknowledgments Supplementary Materials References Copyright Related Articles

Research Article | Open Access

Volume 2020 | Article ID 6345429 | https://doi.org/10.1155/2020/6345429

Bioactivity and Molecular Docking Studies of Derivatives from Cinnamic and Benzoic Acids

Yunierkis Perez-Castillo,¹Tamires C. Lima,²Alana R. Ferreira,³Cecília R. Silva,⁴Rosana S. Campos,⁴João B. A. Neto,⁴Hemerson I. F. Magalhães,³Bruno C. Cavalcanti,⁵Hélio V. N. Júnior,⁴and Damião P. de Sousa³

Academic Editor: Mai S. Li

Received22 Feb 2020

Accepted23 Apr 2020

Published27 May 2020

Abstract

Over the last decade, there has been a dramatic increase in the prevalence and gravity of systemic fungal diseases. This study aimed therefore at evaluating the antifungal potential of ester derivatives of benzoic and cinnamic acids from three Candida species. The compounds were prepared via Fischer esterification, and the antifungal assay was performed by the microdilution method in 96-well microplates for determining the minimal inhibitory concentrations (MICs). The findings of the antifungal tests revealed that the analogue compound methyl ferulate, methyl o-coumarate, and methyl biphenyl-3-carboxylate displayed an interesting antifungal activity against all Candida strains tested, with MIC values of 31.25-62.5, 62.5-125, and 62.5 μg/ml, respectively. A preliminary Structure-Activity Relationship study of benzoic and cinnamic acid derivatives has led to the recognition of some important structural requirements for antifungal activity. The results of molecular docking indicate that the presence of the enoate moiety along with hydroxyl and one methoxy substitution in the phenyl ring has a positive effect on the bioactivity of compound 7 against Candida albicans. These observations further support the hypothesis that the antifungal activity of compound 7 could be due to its binding to multiple targets, specifically to QR, TS, and ST-PK. Additional experiments are required in the future to test this hypothesis and to propose novel compounds with improved antifungal activity.

1. Introduction

In the last decade, the number and severity of systemic fungal infections has increased intensely, and this has been linked specifically to an increased number of immunocompromised hosts, such as patients undergoing organ transplantation, receiving anticancer chemotherapy, suffering from AIDS, or in therapy with broad-spectrum antimicrobials [1, 2]. This factor, together with the limitations of available antifungal therapies and the crescent emergence of drug resistance, has created an urgent need for the development of new efficient broad-spectrum and safer antifungal agents, mainly with novel mechanisms for action [3, 4]. Candida species are broadly distributed in nature and constitute a part of the normal human microflora of the gastrointestinal tract, oral cavity, and vagina. However, a small percentage of the identified species is capable of causing opportunistic infections in humans, named as “candidiasis” [5–7]. Clinically, the most important pathogenic species of the Candida genus are C. albicans, C. parapsilosis, C. krusei, C. tropicalis, and C. glabrata [8].

Phenolic compounds comprise a broad class of secondary metabolites that possess a large spectrum of biological properties [9–12], and they have been successfully investigated for the treatment of systemic and superficial mycoses caused by different human pathogenic fungi [13–17]. Benzoic and cinnamic acids and their derivatives play an important role in antifungal activity, and several studies have focused on this bioactivity [18–25]. For example, Heleno and collaborators [26] reported the antifungal activity of the p-hydroxybenzoic and cinnamic acids as well as their protected glucuronide derivatives against different species of filamentous fungi, including species of the genus Aspergillus, Penicillium, and Trichoderma.

Thus, the objective of the present investigation was to investigate the antifungal potential of a set of 23 ester derivatives of benzoic and cinnamic acids from different Candida strains, including C. albicans, C. krusei, and C. parapsilosis. It also examined the Structure-Activity Relationship (SAR) of the screened compounds, thus, affording information for a better optimization of these classes of compounds as potential antimycotic agents.

2. Materials and Methods

2.1. General Procedure for the Synthesis of Compounds

A mixture of organic acid (0.5 g) and methanol (100 ml) was heated under reflux in the presence of sulfuric acid (0.8 ml) until the completion of the reaction which was checked by single spot TLC. Then, methanol was removed under reduced pressure and the solution was diluted with 20 ml of water. The product was extracted with ethyl acetate (30 ml). The organic phase was neutralized successively with NaHCO₃ 5% and water, dried over anhydrous Na₂SO₄, and filtered. The ethyl acetate phase was separated, which on evaporation yielded the ester derivatives [27]. The structural identification of the products was made by infrared, ¹H and ¹³C Nuclear Magnetic Resonance (NMR) analysis and comparison with the literature data [28–31].

2.2. In Vitro Antifungal Activity

Minimum inhibitory concentrations (MICs) for compounds benzoic and cinnamic derivatives were investigated by broth microdilution, according to the document M27-A3 [32], using RPMI broth (pH 7.0) buffered with 0.165 M MOPS (morpholinepropanesulfonic acid) (Sigma Chemical, St. Louis, MO). The compounds were dissolved in DMSO. To determine the susceptibility of the planktonic cells, drugs were tested at concentrations ranging from 0.98 to 1000 μg/ml for the compounds benzoic and cinnamic derivatives. Yeasts and compounds were incubated in 96-well culture plates at 35°C for 24 h, and the results were examined visually, as recommended by CLSI [32]. The minimum inhibitory concentration (MIC) of each compound was determined as the concentration that inhibited 50% of fungal growth.

2.3. Molecular Docking

2.3.1. Targets Selection

The selection of the potential targets of compound 7 was based on two criteria. First, there is evidence that coniferyl aldehyde derivatives can interfere with cell wall synthesis in Candida albicans. These previous results suggest that C-8 sterol isomerase (Uniprot ID A0A1D8PCB9, C-8-SI), Lanosterol 14-alpha demethylase (Uniprot ID P10613, L-14a-D), Thymidylate synthase (Uniprot ID P12461, TS), and Squalene monooxygenase (Uniprot ID Q92206, SQM) are potential targets of compounds bearing this scaffold. The second group of potential targets of compound 7 was selected based on a consensus target fishing approach. For this, the potential targets of the compound were predicted with five publicly available target fishing services: MolTarPred [33], SwissTargetPrediction [34], TargetNet [35], SuperPred [36], and Similarity Ensemble Approach (SEA) [37]. Any Candida albicans protein predicted as a potential target by any target fishing method was selected for modeling studies. In addition, the potential targets predicted by at least four methodologies were considered for homologs identification on Candida albicans. The latter consisted in the blast of the sequences of the predicted consensus targets against the Candida albicans proteome using the NCBI Blast server [38]. The predicted targets for which a homologous protein existed in Candida albicans were added to the list of potential targets of compound 7. This second group of included enzymes: Carbonic anhydrase (Uniprot ID Q5AJ71, CA), Polyamine oxidase (Uniprot ID Q5AMQ8, PO), Quinone reductase (Uniprot ID Q59Z38, QR), and Serine/threonine protein kinase (Uniprot ID A0A1D8PH81, ST-PK).

2.3.2. Molecular Docking

The X-ray structure of L-14a-D in complex with an inhibitor (PDB code 5tz1) and of CA (PDB code 6gwu) was obtained from the Protein Data Bank [39]. The rest of the potential targets had no structure determined, and homology models were obtained for them from the SWISS-MODEL server [40]. These models are searchable in the SWISS-MODEL repository using the above listed Uniprot IDs. The OMEGA software [41] was used to obtain the initial 3D structure of compound 7, and AM1-BCC charges were added to it with MOLCHARGE [42]. Molecular docking of compound 7 to its potential targets was performed with Gold [43] following the same procedure as in our previous research [44, 45]. Water molecules, cocrystalized ligands, and cofactors not involved in ligand binding were removed from the receptors. The receptor binding sites were defined based either on the cocrystalized ligands or functionally relevant residues. Primary docking was performed with the CHEMPLP scoring function, setting the Gold’s search efficiency to 200%. In total, 30 different docking solutions were explored and all of them were rescored using the ChemScore, GoldScore, and ASP scoring functions. The most probable binding modes of compound 7 to the selected targets were selected according to consensus scoring approach based on the scaling of each scoring function and their aggregation according to equation (1):

is the score of conformer according to scoring function and and are the mean and standard deviation of the scoring function across all conformers, respectively. Any ligand-receptor complex with an aggregated score higher than 1 was considered for further analysis.

2.3.3. Molecular Dynamics and Free Energy of Binding Calculations

All molecular dynamics (MD) simulations and free energy of binding calculations were performed with Amber 2018 [14] following the procedure of our previous publication [45]. All complexes were subject to the same modeling procedure which included two minimization steps, the heating of the systems, their equilibration, production runs, the extraction of representative snapshots, and the estimation of the free energies of binding from them through MM-PBSA calculations. All MD simulations were performed in explicit solvents using the ff14SB and gaff force fields for proteins and non-amino acidic residues, respectively.

In brief, the antechamber was used to obtain the topology and force field modifications of the complexes. The systems were enclosed in a truncated octahedron box and solvated with TIP3P water molecules. Either Na⁺ or Cl^- ions were added to the solvated systems to neutralize any excess charge. Long-range electrostatic interactions were treated with the PME method with cutoff distances of 10 Å and 12 Å during the first minimization step and the remaining steps listed above, respectively.

The systems formed by the complexes, water, and counterions were first minimized during 500 steps of the steepest descent method followed by 500 cycles of conjugate gradient at constant volume with all atoms except water molecules and ions constrained with a force constant of 500 kcal/mol.Å². During the second minimization step, no constraint was imposed on the system, and it took place for 1500 steps of the steepest descent method followed by 1000 cycles of conjugate gradient also at constant volume.

Next, the temperature of the systems was increased from 0 to 300 K at constant volume while keeping everything except water molecules and ions restrained with a force constant of 10 kcal/mol.Å². Heating took place during 10000 steps with a time step of 2 fs. Once the systems reached 300 K in temperature, they were equilibrated at constant pressure (1 bar) and temperature (300 K). Isotropic position scaling with a relaxation time of 2 ps was used to control the pressure during equilibration.

The last snapshot from the equilibration stage was used as input to 20 MD simulations lasting 2 ns each. To ensure a better exploration of the conformational space of the systems, each of these simulations was initialized with different random velocities. During the heating, equilibration, and production run stages, the bonds involving hydrogen atoms were constrained with the SHAKE algorithm and the temperature was controlled using a Langevin thermostat with a collision frequency of 1.0 ps^-1.

The MM-PBSA algorithm implemented in AmberTools 18 [46] was selected for the estimation of the free energies of binding of the ligand to the receptors. For these calculations, one MD snapshot every 200 ps was selected from every one of the previous production runs. This process was conducted by the selection of 200 diverse conformations of each system for MM-PBSA calculations. In addition, default implicit solvent parameters were used and the ionic strength was set to 100 mM.

3. Results and Discussion

3.1. Chemistry and Antifungal Activity of Compounds 1-23

In the present investigation, a series of 23 ester derivatives of benzoic and cinnamic acids (Figure 1) with different substitutions on the aromatic ring were assessed in vitro for their antifungal activity against three Candida species, namely, C. albicans, C. krusei, and C. parapsilosis. The antifungal assay was performed by the microdilution method in 96-well microplates for determining the minimal inhibitory concentrations (MICs). The results of the MIC testing, expressed in μg/ml, are reported in Table 1. Among the 23 evaluated ester derivatives, the methyl ferulate (7) was found to exhibit the best antifungal activity, with MIC values of 62.5, 31.25, and 62.5 μg/ml against C. albicans, C. parapsilosis, and C. krusei, respectively, followed by methyl biphenyl-3-carboxylate (22) (equal ) and methyl o-coumarate (2) (, 62.5, and 125 μg/ml, respectively). Methyl 3,5-di-tert-butyl-4-hydroxybenzoate (16), methyl biphenyl-4-carboxylate (23), and methyl benzoate (12) were the less toxic compounds, with MIC values of 500 μg/ml for the three Candida species. Furthermore, C. parapsilosis was the most susceptible species.

To obtain a better understanding of the relationship between the structural features of the screened ester derivatives and their antifungal effect, some structural and molecular characteristics which may contribute to the antifungal activity were recognized and will be discussed. In general, ester derivatives of the cinnamic acid were more efficient than their corresponding benzoic counterparts. For example, methyl ferulate (7) () had stronger antifungal activity than methyl vanillate (18) (), methyl cinnamate (1) (equal ) was more bioactive than methyl benzoate (12) (equal ), methyl p-coumarate (4) () was more potent than methyl 4-hydroxybenzoate (13) (equal ), and methyl sinapate (8) () was more effective than methyl syringate (15) ().

The addition of a hydroxyl group to the aromatic ring, in general, resulted in higher antifungal action. For example, when the bioactivity of the compounds methyl cinnamate (1), methyl o-coumarate (2), methyl m-coumarate (3), and methyl p-coumarate (4) was compared, it was found an increase in potency in the following order: (, 125-250, 250, and 250 μg/ml, respectively). A similar comparative trend can be seen when the compounds methyl benzoate (12) and methyl 4-hydroxybenzoate (13) were compared, once that 13 (equal ) was more bioactive than 12 (equal ). Furthermore, the position of the hydroxyl group on the benzene ring changes the activity, because compound 2 (ortho-hydroxyl; ) was more active than 3 (meta-hydroxyl; equal ) and 4 (para-hydroxyl; ). A study performed by Kim et al. (2011) [47] demonstrated that various benzaldehydes, like 2-hydroxy-3-methoxybenzaldehyde and 2-hydroxy-5-methoxybenzaldehyde, have been recognized to have a potent antifungal activity against different filamentous fungi, and a structure-activity relationship analysis showed that the antifungal effect was improved by the presence of an ortho-hydroxyl group on the aromatic ring.

However, the introduction of a second or third hydroxyl group into the benzene ring decreased the antifungal effect. For example, the compound methyl caffeate (5) (equal ) was less bioactive than methyl p-coumarate (4) (), and methyl gallate (14) () was less potent than methyl 4-hydroxybenzoate (13) (equal ). Ambobaré and coworkers (2002) [19] reported that the introduction of a supplementary hydroxyl group in different positions on the aromatic ring of the benzoic acid derivatives strongly decreased the antifungal activity against Eutypa lata compared to salicylic acid, while the addition of two hydroxyl groups did not noticeably enhance the effectiveness. Furthermore, 2-hydroxyl group was the best replacement for the benzene ring. In another study, Sánchez-Maldonado and collaborators (2011) [48] evaluated the antimicrobial effect of twelve hydroxybenzoic and hydroxycinnamic acids against four gram-positive bacteria: Lactobacillus plantarum, L. hammesii, Escherichia coli, and Bacillus subtilis. These authors described how the activity of the hydroxybenzoic acids diminished significantly with an increasing number of hydroxyl groups. For example, the antibacterial activity of gallic acid (three hydroxyl groups on the aromatic ring) was about 2-10-fold lesser when compared with other hydroxybenzoic acids.

The presence a 5-hydroxyl group together with an additional methoxy group at the 3-position on trans-cinnamic acid methyl moiety appears to be important to enhance the antifungal activity, because the methyl ferulate (7; ) was more potent than the compounds 1 (methyl cinnamate, equal ), 4 (methyl p-coumarate, ), 5 (methyl caffeate, equal ), and 6 (methyl 4-methoxycinnamate, equal ). Furthermore, the increase of methoxy groups on the benzene ring remarkably decreased the antifungal effect, because compound 7 (one methoxy group; ) was more potent than 8 (two methoxy groups; ), and this was more active than 9 (three methoxy groups; ).With regard to the aromatic substitution within class benzoic, it was found that the increase of methoxy groups on the benzene do not affect the antifungal effect, once the compounds methyl gallate (14), methyl syringate (15), and methyl vanillate (18) exhibited the same activity (). In contrast to the results reported herein, Sánchez-Maldonado and coworkers (2011) [48] observed that the antibacterial activity of hydroxybenzoic acids was enhanced through the replacement of a hydroxyl group by a methoxy group. Moreover, the presence of methoxy groups did not significantly influence the antimicrobial effect of cinnamic acids.

The introduction of bulky alkyl substituents into the aromatic ring resulted in a decrease in the antifungal effect. For example, the compound methyl 3,5-di-tert-butyl-4-hydroxybenzoate (16; two tert-butyl groups) was less potent (equal ) than compounds 14 and 15 (equal ). This diminution in the antifungal effect probably occurs due to steric hindrance caused by the presence of two tert-butyl groups. Similar results were also observed in others studies [49, 50].

The addition of a NO₂ group at the 2-position of the methyl cinnamate, resulting in compound 10 (methyl 2-nitrocinnamate), did not increase antifungal activity, as compounds 1 and 10 presented the same antifungal potency (equal ). Furthermore, the compound methyl 3-methyl-4-nitrobenzoate (20; ) was less bioactive than the compound methyl p-toluate (17; ), and the compound methyl 4-chlorocinnamate (11; ) was less potent than the compound methyl cinnamate (1; equal ). These observations suggest that the decrease in electron diversity caused by electron-withdrawing groups (such as nitro or chlorine) resulted in a reduction of the antifungal effect. Likewise, compounds possessing electron donating groups (like methyl) exhibited better antifungal activity, as methyl p-toluate (17; ) was more active than its corresponding methyl cinnamate (1; equal ). Similar results were also obtained by Fu and collaborators (2010) [51]. These authors found that caffeic acid anilides with electron donating groups (such as methyl or methoxy) in p-position of the benzene ring displayed better antibacterial activity than those with electron-withdrawing groups (like bromide, nitro, or fluoride).

Lipophilicity has been extensively studied and is considered as one of the most important parameters influencing the antifungal activity [21, 23, 34]; however, no correlation was observed in this set of tested compounds, as most lipophilic compounds, like 3,5-di-tert-butyl-4-hydroxybenzoate (16) and methyl biphenyl-4-carboxylate (23), exhibited weak antifungal effect.

3.2. Molecular Docking

The molecular docking was performed with the most potent compound of the series, methyl ferulate (7). Otherwise noted in the Materials and Methods section, default parameters were used for all calculations described from here on. Compound 7 was docked at the binding sites of C-8-SI, L-14a-D, TS, SQM, CA, PO, QR, and ST-PK following the procedure described above. The results of the docking studies are presented in Table 2. For SQM, Gold could not find any valid docking solution and hence this target was excluded from any further analysis. The docking to the rest of the potential targets reveals that for most of them more than one possible binding mode of compound 7 is predicted. Overall, the potential complexes summarized in Table 2 show favorable interactions between the ligand and the receptor. These interactions include contacts with functionally important residues in the receptors. In addition, the best docking scores are predicted for L-14a-D, PO, and ST-PK.

Despite the meaningful complexes predicted by molecular docking, this type of modeling tool is primarily optimized for the processing of large databases of chemical compounds in a short time. For this aim, molecular docking uses simplified models that take into account a limited number of factors associated with molecular recognition. In consequence, the docking predicted complexes are not necessarily stable. Previous studies have shown that many predicted ligand-receptor complexes, although presenting favorable interactions might not be feasible [44, 45]. To further refine the possible mechanism of action of compound 7, MD simulations and MM-PBSA calculations were performed as described in the Materials and Methods section. The results of the estimation of the free energy of binding of compound 7 in the complexes listed in Table 2 are summarized in Table 3.

As seen in Table 3, compound 7 is predicted to bind stably to ST-PK, TS, and QR. According to our results, the lowest free energy of binding is obtained for ST-PK, followed by TS and QR. Supporting the previous observation, despite having the highest docking scores, the complexes predicted with L-14a-D are the less stable ones among all studied systems. The most probable binding mode of compound 7 to ST-PK, TS, and QR as well as the network of interactions that it forms with the receptors are presented in Figure 2. The most probable binding mode of the compound to each receptor was selected as that having the lowest predicted ΔG of binding.

(a)

(b)

(c)

Figure 2

Predicted binding modes of compound 7 to ST-PK (a), TS (b), and QR (c) of Candida albicans (left). Compound 7 is depicted cyan and the receptors in gray, with noncarbon atoms following the scheme: blue for N, red for O, yellow for S, white for H, and orange for P. On the right side of the figure are represented the predicted interaction frequencies with the residues at the receptors binding sites. Darker lines indicate the higher frequencies of interaction. Only residues interacting with the ligand in more than 50% of the analyzed MD snapshots are labelled in the complexes’ structures.

The detailed description of the interactions observed between compound 7 and the ST-PK, TS, and QR receptors is discussed next. Figures of the predicted complexes of compound 7 with the receptors discarded due to its predicted poor ΔG of binding (L-14a-D, C-8-SI, PO, and CA) as well as their description are provided as Supporting Information. It must be considered that biomolecular complexes are dynamic systems, thus, the following analyses focus on the ligand-receptor interactions observed in more than 50% of the analyzed MD snapshots. In all cases, the complexes used for depiction are the centroids of the most populated clusters obtained from the 200 MD snapshots employed for the MM-PBSA calculations in each system.

Compound 7 is predicted to interact with QR through direct contact with the NADPH cofactor in all MD snapshots extracted for MM-PBSA calculations. This predicted binding mode would prevent the access of the enzyme’s substrate to its binding site and conduct of the blockage of the enzyme functioning. Specifically, compound 7 is predicted to hydrogen bond, one of the phosphate groups of NADPH in 84% of the selected snapshots through its hydroxyl group. Less frequent hydrogen bonds are also predicted with the backbones of I100 and Y325 in 28% and 12.5% of the snapshots, respectively. In addition, the complex is stabilized by the stacking of the ligand’s phenyl ring in front of the cofactor’s nicotinamide group. Compound 7 also contacts N98, F99, I100, L182, and Q183 in most of the selected MD snapshots. All of these residues are located close to the cofactor and stabilize the orientation of the aromatic moiety of the compound in the receptor-binding pocket. On the other hand, the enoate moiety of the ligand shows low frequency of interaction with the receptor and it mainly interacts with the solvent. These low-frequency interactions take place primarily with L324, L144, Y325, and Q323.

The results of the docking of compound 7 to the deoxyuridine monophosphate binding site of TS in the presence of tetrahydrofolate shows a large network of interactions between the ligand and the receptor. The most frequent ones along the selected MD snapshots are with the cofactor, R29, E68, W90, Y116, L173, C176, H177, D216, N224, H254, and Y256. Furthermore, the carbonyl and hydroxyl groups of the ligand can hydrogen bond the receptor in more than 75% of the selected complex conformations. Despite the enoate moiety of the ligand makes more contact with TS than to QR, it is also located at the mouth of the binding site and exposed to the solvent. In contrast, the phenyl ring locates close to Y116 in a position that favors the π-π stacking of both rings.

Finally, compound 7 is predicted to bind to ST-PK in its ATP binding grove, making contact mainly with L436, V444, A457, K459, E474, I491, F507, I509, L510, L561, I573, and D574. Similarly to the predicted binding mode to TS, the complex shows that both the hydroxyl and carbonyl groups of the ligand are able to form hydrogen bonds with the receptor in most of the studied MD snapshots. Specifically, the carbonyl moiety accepts a hydrogen bond from the backbone of L510, while the hydroxyl substituent is able to hydrogen bond the side chains of K459 and E474. This predicted mechanism of action is in agreement with that previously proposed [52], in which the ligand interacts with K459, E474, and D574, a feature common to potent kinase inhibitors [53].

Congeneric series of compounds such as those presented in this research should follow the similarity principle, that is, similar compounds should have similar bioactivity profiles. To further assess the potential targets of the studied series of compounds we performed the same calculations carried out with compound 7 for chemicals 2 and 22. These are the ones with the higher potency after compound 7. The full statistics of their modeling processes are provided as Supporting Information. Interestingly, compounds 2 and 22 follow the same pattern of predicted targets as compound 7. Specifically, compound 2 is predicted to bind stably to the TS and ST-PK proteins while its ΔG of binding predicts an unstable complex with TS. On the other hand, compound 22 is predicted to form stable complexes with the same three potential targets predicted for compound 7: TS, QR, and ST-PK. In addition, unfeasible complexes are predicted between compounds 2 and 22 and the L-14a-D, C-8-SI, PO, and CA receptors as for compound 7.

The most stable complexes formed by compounds 7, 2, and 22 (see Table 3 and Supporting Information) with TS, QR, and ST-PK are given in Figure 3. One common feature among the complexes predicted for compound 7 with TS, QR, and ST-PK is the presence of hydrogen bonds to the receptor mediated by the p-hydroxyl substituent at its phenyl ring. The importance of these interactions is highlighted by the analysis of the predicted binding modes of compounds 2 and 22 to the same targets. The lack of this group in the later makes these interactions impossible, thus changing in different degrees the geometry of the predicted ligand-receptor complexes. The importance of these hydrogen bonds can be illustrated by the predicted binding modes of compounds 2 and 22 to ST-PK. These compounds are rotated 180° relative to the predicted orientation of compound 7 on ST-PK. This rotation restores the hydrogen bond interactions of these ligands with K459, E474, and D574 which are predicted to stabilize the complex with compound 7. Interestingly, such an orientation is impossible for compound 7 since its bulky p-methoxy substituent would clash with ST-PK residues I491 and E508.

Another interesting observation is that the removal of the p-hydroxyl and m-methoxy groups in compounds 2 and 22 reduces the volume of the cavity occupied by these ligands at the bottom of the binding site of TS. In consequence, W90 is predicted to rotate its side chain toward the binding pocket to occupy the space filled by these substituents in compound 7. This large rotation of W90 reduces the volume of the TS binding cavity and still allows for favorable ligand-receptor interactions with compounds 2 and 22.

QR is a flavodoxin-like protein that is involved in the control of oxidative stress in Candida albicans [54]. In consequence, the inhibition of this enzyme could induce oxidative stress in this fungus, ultimately leading to cell death. In a similar way, TS is a validated target against Candida albicans. There is interference in DNA synthesis and repair if it is inhibited [55]. In addition, ST-PK has been shown to be involved in the hyphal growth in Candida albicans [56], and it has been proposed as an effective target for Candida albicans inhibition.

The proposed mechanisms of action are in agreement with the observed structure-activity relationship for the assayed series of compounds. For example, the lack of the enoate moiety present in compound 7 (see compound 18) reduces the number of favorable contacts with all receptors. More important, the removal of this moiety abolishes or reduces the number or predicted hydrogen bonds between the most active compounds and the receptors. In the same vein, the presence of the hydroxyl substitution in the phenyl ring allows for the formation of hydrogen bonds with functionally important residues in all cases. This is exemplified by the loss of activity observed for compound 9. Besides the positive influence of a hydroxyl group attached to the phenyl ring, the addition of a second methoxy substituent in the phenyl ring leads to a reduction in the bioactivity of the studied compounds. The latter can be explained by the steric constraints that this substitution imposes as the addition of a second methoxy group contributes with steric clashes with the receptors (see compounds 7, 8, and 9).

4. Conclusions

In summary, we have prepared a set of derivatives of benzoic and cinnamic acids and investigated for their antifungal activity against three Candida strains. Additionally, a preliminary Structure-Activity Relationship study also was performed, permitting the recognition of some important structural requirements for antifungal activity and providing information for a better optimization of these classes of compounds as potential antimycotic agents. Molecular docking and the estimation of the free energies of binding were also attributed to the most potent compounds against Candida strains: 7, 2, and 22. The results show that the antifungal activity of these compounds could be due to its binding on multiple targets and have suggested the contribution of substituent groups to bioactivity.

Data Availability

Part of the article data used to support the findings of this study has been deposited in the Federal University of Paraíba repository at https://repositorio.ufpb.br/jspui/handle/123456789/1356.

Conflicts of Interest

No potential conflict of interest was reported by the authors.

Acknowledgments

This research received financial supported from the Brazilian government agencies Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) (Process no. 421285/2016-8) and the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES).

Supplementary Materials

Predicted compound 7–L-14a-D complex. Predicted compound 7–C-8-SI complex. Predicted compound 7–PO complex. Predicted compound 7–CA complex. Docking results for compound 2. Docking results for compound 22 MM-PBSA results for compound 2 MM-PBSA results for compound 22. (Supplementary Materials)

References

W. Wang, G. Dong, J. Gu et al., “Structure–activity relationships of tetrahydrocarbazole derivatives as antifungal lead compounds,” Medicinal Chemistry Communications, vol. 4, no. 2, pp. 353–362, 2013.
View at: Publisher Site | Google Scholar
E. Rusu, M. Radu-Popescu, D. Pelinescu, and T. Vassu, “Treatment with some anti-inflammatory drugs reduces germ tube formation in Candida albicans strains,” Brazilian Journal of Microbiology, vol. 45, no. 4, pp. 1379–1383, 2014.
View at: Publisher Site | Google Scholar
K. S. Babu, X. C. Li, M. R. Jacob et al., “Synthesis, antifungal activity, and structure-activity relationships of Coruscanone A analogues,” Journal of Medicinal Chemistry, vol. 49, no. 26, pp. 7877–7886, 2006.
View at: Publisher Site | Google Scholar
R. Sharanappa and G. M. Vidyasagar, “Anti-Candida activity of medicinal plants: a review,” International Journal Pharmacy and Pharmaceutial Sciences, vol. 5, no. 4, pp. 9–16, 2013.
View at: Google Scholar
N. Papon, V. Courdavault, M. Clastre, and R. J. Bennett, “Emerging and emerged pathogenic Candida species: beyond the Candida albicans paradigm,” PLOS Pathogens, vol. 9, no. 9, article e1003550, 2013.
View at: Publisher Site | Google Scholar
J. C. O. Sardi, L. Scorzoni, T. Bernardi, A. M. Fusco-Almeida, and M. J. S. Mendes Giannini, “Candida species: current epidemiology, pathogenicity, biofilm formation, natural antifungal products and new therapeutic options,” Journal of Medical Microbiology, vol. 62, no. 1, pp. 10–24, 2013.
View at: Publisher Site | Google Scholar
E. Paramythiotou, F. Frantzeskaki, A. Flevari, A. Armaganidis, and G. Dimopoulos, “Invasive fungal infections in the ICU: how to approach, how to treat,” Molecules, vol. 19, no. 1, pp. 1085–1119, 2014.
View at: Publisher Site | Google Scholar
M. A. Pfaller, D. J. Diekema, G. W. Procop, and M. G. Rinaldi, “Multicenter comparison of the VITEK 2 antifungal susceptibility test with the CLSI broth microdilution reference method for testing amphotericin B, flucytosine, and voriconazole against Candida spp,” Journal of Clinical Microbiology, vol. 45, no. 11, pp. 3522–3528, 2007.
View at: Publisher Site | Google Scholar
K. Rashed, X. J. Zhang, M. T. Luo, and Y. T. Zheng, “Anti-HIV-1 activity of phenolic compounds isolated from Diospyros lotus fruits,” Phytopharmacology, vol. 3, no. 2, pp. 199–207, 2012.
View at: Google Scholar
C. A. L. Cardoso, R. G. Coelho, N. K. Honda, A. Pott, F. R. Pavan, and C. Q. F. Leite, “Phenolic compounds and antioxidant, antimicrobial and antimycobacterial activities of Serjania erecta Radlk. (Sapindaceae),” Brazilian Journal of Pharmaceutical Sciences, vol. 49, no. 4, pp. 775–782, 2013.
View at: Publisher Site | Google Scholar
M. Sova, Z. Zizak, J. Stankovic et al., “Cinnamic acid derivatives induce cell cycle arrest in carcinoma cell lines,” Medicinal Chemistry, vol. 9, no. 5, pp. 633–641, 2013.
View at: Publisher Site | Google Scholar
T. C. Lima, R. J. Souza, A. D. C. Santos et al., “Evaluation of leishmanicidal and trypanocidal activities of phenolic compounds from Calea uniflora less,” Natural Product Research, vol. 30, no. 5, pp. 551–557, 2015.
View at: Publisher Site | Google Scholar
A. Rao, Y. Zhang, S. Muend, and R. Rao, “Mechanism of antifungal activity of terpenoid phenols resembles calcium stress and inhibition of the TOR pathway,” Antimicrobial Agents and Chemotherapy, vol. 54, no. 12, pp. 5062–5069, 2010.
View at: Publisher Site | Google Scholar
N. C. G. Faria, J. H. Kim, L. A. P. Gonçalves, M. . L. Martins, K. L. Chan, and B. C. Campbell, “Enhanced activity of antifungal drugs using natural phenolics against yeast strains of Candida and Cryptococcus,” Letters in Applied Microbiology, vol. 52, no. 5, pp. 506–513, 2011.
View at: Publisher Site | Google Scholar
M. A. Ansari, A. Anurag, Z. Fatima, and S. Hameed, “Natural phenolic compound: a potential antifungal target,” in Microbial Pathogens and Strategies for Combating Them: Science, Technology and Education, A. Méndez-Vilas, Ed., vol. 2, pp. 1189–1195, Formatex Research Center, Spain, 2013 edition, 2013.
View at: Google Scholar
K. Rashed, A. Ciric, J. Glamoclija, and M. Sokovic, “Antibacterial and antifungal activities of methanol extract and phenolic compounds from Diospyros virginiana L.,” Industrial Crops and Products, vol. 59, pp. 210–215, 2014.
View at: Publisher Site | Google Scholar
J. Rodríguez-Morató, L. Xicota, M. Fitó, M. Farré, M. Dierssen, and R. de la Torre, “Potential role of olive oil phenolic compounds in the prevention of neurodegenerative diseases,” Molecules, vol. 20, no. 3, pp. 4655–4680, 2015.
View at: Publisher Site | Google Scholar
S. Tawata, S. Taira, N. Kobamoto, J. Zhu, M. Ishihara, and S. Toyama, “Synthesis and antifungal activity of cinnamic acid esters,” Bioscience, Biotechnology, and Biochemistry, vol. 60, no. 5, pp. 909-910, 2014.
View at: Publisher Site | Google Scholar
B. E. Amborabé, P. Fleurat-Lessard, J. F. Chollet, and G. Roblin, “Antifungal effects of salicylic acid and other benzoic acid derivatives towards Eutypa lata: structure–activity relationship,” Plant Physiology and Biochemistry, vol. 40, no. 12, pp. 1051–1060, 2002.
View at: Publisher Site | Google Scholar
A. López, D. S. Ming, and G. H. N. Towers, “Antifungal activity of benzoic acid derivatives from Piper lanceaefolium,” Journal of Natural Products, vol. 65, no. 1, pp. 62–64, 2002.
View at: Publisher Site | Google Scholar
B. Narasimhan, D. Belsare, D. Pharande, V. Mourya, and A. Dhake, “Esters, amides and substituted derivatives of cinnamic acid: synthesis, antimicrobial activity and QSAR investigations,” European Journal of Medicinal Chemistry, vol. 39, no. 10, pp. 827–834, 2004.
View at: Publisher Site | Google Scholar
S. S. Cheng, J. Y. Liu, E. H. Chang, and S. T. Chang, “Antifungal activity of cinnamaldehyde and eugenol congeners against wood-rot fungi,” Bioresource Technology, vol. 99, no. 11, pp. 5145–5149, 2008.
View at: Publisher Site | Google Scholar
M. Kratky and J. Vinsova, “Antifungal activity of salicylanilides and their esters with 4-(Trifluoromethyl)benzoic acid,” Molecules, vol. 17, no. 8, pp. 9426–9442, 2012.
View at: Publisher Site | Google Scholar
M. Kratky, J. Vinsova, and V. Buchta, “In vitro antibacterial and antifungal activity of salicylanilide benzoates,” The Scientific World Journal, vol. 2012, Article ID 290628, 7 pages, 2012.
View at: Publisher Site | Google Scholar
J. Guzman, “Natural cinnamic acids, synthetic derivatives and hybrids with antimicrobial activity,” Molecules, vol. 19, no. 12, pp. 19292–19349, 2014.
View at: Publisher Site | Google Scholar
S. A. Heleno, I. C. F. R. Ferreira, A. P. Esteves et al., “Antimicrobial and demelanizing activity of Ganoderma lucidum extract, p-hydroxybenzoic and cinnamic acids and their synthetic acetylated glucuronide methyl esters,” Food and Chemical Toxicology, vol. 58, pp. 95–100, 2013.
View at: Publisher Site | Google Scholar
M. O. Farias, T. C. Lima, A. L. A. L. Pérez et al., “Antifungal activity of ester derivatives from caffeic acid against Candida species,” International Journal of Pharmacy and Pharmaceutal Research, vol. 7, no. 1, pp. 151–159, 2016.
View at: Google Scholar
R. P. Mahajan, S. L. Patil, and R. S. Mali, “Convenient microwave assisted synthesis of naturally occurring methyl (E)-cinnamates,” Organic Preparations and Procedures International, vol. 37, no. 3, pp. 286–290, 2005.
View at: Publisher Site | Google Scholar
W. J. Huang, C. C. Chen, S. W. Chao et al., “Synthesis of N-hydroxycinnamides capped with a naturally occurring moiety as inhibitors of histone deacetylase,” ChemMedChem, vol. 5, no. 4, pp. 598–607, 2010.
View at: Publisher Site | Google Scholar
Y. Li, J. Wang, Z. Wang et al., “Palladacycle-catalyzed Suzuki–Miyaura reaction of aryl/heteroaryl halides with MIDA boronates in EtOH/H₂O or H₂O,” RSC Advances, vol. 4, no. 68, pp. 36262–36266, 2014.
View at: Publisher Site | Google Scholar
Z. Zhang and Z. Wang, “Diatomite-supported Pd nanoparticles: an efficient catalyst for Heck and Suzuki reactions,” Journal of Organic Chemistry, vol. 71, no. 19, pp. 7485–7487, 2006.
View at: Publisher Site | Google Scholar
Clinical and Laboratory Standards Institute (CLSI), Reference method for broth dilution antifungal susceptibility testing of yeasts, CLSI Standard M27-A3, 2008.
A. Peón, H. Li, G. Ghislat et al., “MolTarPred: a web tool for comprehensive target prediction with reliability estimation,” Chemical Biology & Drug Design, vol. 94, no. 1, pp. 1390–1401, 2019.
View at: Publisher Site | Google Scholar
D. Gfeller, A. Grosdidier, M. Wirth, A. Daina, O. Michielin, and V. Zoete, “SwissTargetPrediction: a web server for target prediction of bioactive small molecules,” Nucleic Acids Research, vol. 42, no. W1, pp. w32–w38, 2014.
View at: Publisher Site | Google Scholar
Z.-J. Yao, J. Dong, Y. J. Che et al., “TargetNet: a web service for predicting potential drug-target interaction profiling via multi-target SAR models,” Journal of Computer-Aided Molecular Design, vol. 30, no. 5, pp. 413–424, 2016.
View at: Publisher Site | Google Scholar
J. Nickel, B. O. Gohlke, J. Erehman et al., “SuperPred: update on drug classification and target prediction,” Nucleic Acids Research, vol. 42, no. W1, pp. W26–W31, 2014.
View at: Publisher Site | Google Scholar
M. J. Keiser, B. L. Roth, B. N. Armbruster, P. Ernsberger, J. J. Irwin, and B. K. Shoichet, “Relating protein pharmacology by ligand chemistry,” Nature Biotechnology, vol. 25, no. 2, pp. 197–206, 2007.
View at: Publisher Site | Google Scholar
S. F. Altschul, T. L. Madden, A. A. Schäffer et al., “Gapped BLAST and PSI-BLAST: a new generation of protein database search programs,” Nucleic Acids Research, vol. 25, no. 17, pp. 3389–3402, 1997.
View at: Publisher Site | Google Scholar
H. M. Berman, J. Westbrook, Z. Feng et al., “The Protein Data Bank,” Nucleic Acids Research, vol. 28, no. 1, pp. 235–242, 2000.
View at: Publisher Site | Google Scholar
S. Bienert, A. Waterhouse, T. A. P. de Beer et al., “The SWISS-MODEL Repository—new features and functionality,” Nucleic Acids Research, vol. 45, no. D1, pp. D313–D319, 2017.
View at: Publisher Site | Google Scholar
P. C. D. Hawkins, A. G. Skillman, G. L. Warren, B. A. Ellingson, and M. T. Stahl, “OMEGA, Santa Fe, NM: OpenEye Scientific Software,” http://www.eyesopen.com.
View at: Google Scholar
OpenEye Scientific Software, “QUACPAC, Santa Fe, NM: OpenEye Scientific Software,” http://www.eyesopen.com.
View at: Google Scholar
G. Jones, P. Willett, R. C. Glen, A. R. Leach, and R. Taylor, “Development and validation of a genetic algorithm for flexible docking,” Journal of Molecular Biology, vol. 267, no. 3, pp. 727–748, 1997.
View at: Publisher Site | Google Scholar
H. Turkez, F. R. . Nóbrega, O. Ozdemir et al., “NFBTA: a potent cytotoxic agent against glioblastoma,” Molecules, vol. 24, no. 13, article 2411, 2019.
View at: Publisher Site | Google Scholar
S. P. Lopes, Y. P. Castillo, M. L. Monteiro et al., “Trypanocidal mechanism of action and in silico studies of p-coumaric acid derivatives,” International Journal of Molecular Sciences, vol. 20, no. 23, article 5916, 2019.
View at: Publisher Site | Google Scholar
D. A. Case, I. Y. Ben-Shalom, S. R. Brozell et al., AMBER 2018, University of California, San Francisco, 2018.
J. H. Kim, K. L. Chan, N. Mahoney, and B. C. Campbell, “Antifungal activity of redox-active benzaldehydes that target cellular antioxidation,” Annals of Clinical Microbiology and Antimicrobials, vol. 10, no. 1, p. 23, 2011.
View at: Publisher Site | Google Scholar
A. F. Sánchez-Maldonado, A. Schieber, and M. G. Ganzle, “Structure–function relationships of the antibacterial activity of phenolic acids and their metabolism by lactic acid bacteria,” Journal of Applied Microbiology, vol. 111, no. 5, pp. 1176–1184, 2011.
View at: Publisher Site | Google Scholar
T. C. Lima, S. R. L. Santos, M. P. Uliana et al., “Oxime derivatives with larvicidal activity against Aedes aegypti L,” Parasitology Research, vol. 114, no. 8, pp. 2883–2891, 2015.
View at: Publisher Site | Google Scholar
K. Petnapapun, W. Chavasiri, and P. Sompornpisut, “Structure-Activity Relationships of 3,3-Phenylmethylene- bis-4-hydroxycoumarins: Selective and Potent Inhibitors of Gram-Positive Bacteria,” The Scientific World Journal, vol. 2013, Article ID 178649, 11 pages, 2013.
View at: Publisher Site | Google Scholar
J. Fu, K. Cheng, Z. M. Zhang, R. Q. Fang, and H. L. Zhu, “Synthesis, structure and structure–activity relationship analysis of caffeic acid amides as potential antimicrobials,” European Journal of Medicinal Chemistry, vol. 45, no. 6, pp. 2638–2643, 2010.
View at: Publisher Site | Google Scholar
M. Masłyk, M. Janeczko, A. Martyna et al., “The anti-Candida albicans agent 4-AN inhibits multiple protein kinases,” Molecules, vol. 24, no. 1, p. 153, 2019.
View at: Publisher Site | Google Scholar
A. D. Ferguson, P. R. Sheth, A. D. Basso et al., “Structural basis of CX-4945 binding to human protein kinase CK2,” FEBS Letters, vol. 585, no. 1, pp. 104–110, 2011.
View at: Publisher Site | Google Scholar
L. Li, S. Naseem, S. Sharma, and J. B. Konopka, “Flavodoxin-like proteins protect Candida albicans from oxidative stress and promote virulence,” PLoS Pathogens, vol. 11, no. 9, article e1005147, 2015.
View at: Publisher Site | Google Scholar
M. K. Kathiravan, A. B. Salake, A. S. Chothe et al., “The biology and chemistry of antifungal agents: a review,” Bioorganic & Medicinal Chemistry, vol. 20, no. 19, pp. 5678–5698, 2012.
View at: Publisher Site | Google Scholar
S. Goyard, P. Knechtle, M. Chauvel et al., “The Yak1 kinase is involved in the initiation and maintenance of Hyphal growth inCandida albicans,” Molecular Biology of the Cell, vol. 19, no. 5, pp. 2251–2266, 2008.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2020 Yunierkis Perez-Castillo et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

3781

Downloads

1640

Citations