A Novel Series of Acylhydrazones as Potential Anti-Candida Agents: Design, Synthesis, Biological Evaluation and In Silico Studies

In the context of an increased incidence of invasive fungal diseases, there is an imperative need of new antifungal drugs with improved activity and safety profiles. A novel series of acylhydrazones bearing a 1,4-phenylene-bisthiazole scaffold was designed based on an analysis of structures known to possess anti-Candida activity obtained from a literature review. Nine final compounds were synthesized and evaluated in vitro for their inhibitory activity against various strains of Candida spp. The anti-Candida activity assay revealed that some of the new compounds are as active as fluconazole against most of the tested strains. A molecular docking study was conducted in order to evaluate the binding poses towards lanosterol 14α-demethylase. An in silico ADMET analysis showed that the compounds possess drug-like properties and represent a biologically active framework that should be further optimized as potential hits.


Introduction
The continuous rise of fungal infections due to the increased number of immunocompromised patients is currently a major concern. Different Candida strains, including C. albicans, C. krusei and C. parapsilosis are known to be etiological agents of human infections. C. albicans is the most prevalent cause of invasive fungal diseases, being a serious public health problem due to high rates of mortality [1]. Although C. albicans accounted for more than half of the cases of candidemia in European countries, a progressive increase of non-albicans infections has been reported. The incidence rates for candidemia are 16.3% for C. parapsilosis and 2-5% for C. krusei [2,3]. The occurrence of C. parapsilosis as a nosocomial pathogen continues to rise, this non-albicans Candida strain being commonly associated with prosthetic surfaces, parenteral nutrition and central venous catheters [4].
The mechanism of action of azole antifungals involves the blockage of ergosterol biosynthesis through inhibition of lanosterol 14α-demethylase (CYP51), a member of the cytochrome The mechanism of action of azole antifungals involves the blockage of ergosterol biosynthesis through inhibition of lanosterol 14α-demethylase (CYP51), a member of the cytochrome P450-dependent superfamily. These antifungal agents bind in a strong and reversible manner to the heme iron atom situated in the catalytic site of the enzyme. As a result of CYP51 inhibition, the normal function of fungal membrane is affected due to reduced rigidity and increased permeability and fluidity [5]. Development of different mutations on CYP51 has led to a decreased binding affinity of the drug to the target enzyme and to an increased frequency of drug resistance [6]. The emerging fungal resistance constitutes a widespread problem for the treatment of invasive infections caused by Candida species.
The important role of CYP51 in fungal growth explains the actual interest in maintaining this enzyme as a target for the development of new anti-Candida drugs, but, in the same time, efforts are being made in order to discover compounds with a different inhibitory mechanism.
To accelerate the progress of development of more effective antifungal drugs, the identification of hit or lead compounds needs to be done. Recently published studies revealed that substituted thiazoles exhibit anti-Candida activity [7][8][9]. Furthermore, hydrazone and acylhydrazone derivatives, which possess an azomethine group −N=CH−, have been developed and evaluated for their inhibitory activity on different Candida spp. [10][11][12][13][14].
Gidaro et al. recently developed a ligand-based pharmacophore hypothesis, based on structurally similar molecules with confirmed in vitro anti-Candida activity. In their study, the importance of several chemical features of the theoretical pharmacophore hypothesis was observed: a central feature represented by a heterocycle (aromatic ring and hydrophobic feature-Ar-HF) containing a hydrogen bond acceptor-HBA (thiazole containing the nitrogen atom) linked to a hydrogen bond donor site-HBD and a second HBA site (hydrazone and carbonylic moiety) and to other aromatic or hydrophobic features (phenyl) [15].
Based on this ligand-based pharmacophore hypothesis and on several molecules with in vitro anti-Candida activity, herein we focused on designing a new scaffold that contains some of the above-mentioned chemical features, among others, that are found in different active compounds, as illustrated in Figure 1. Association of various structural moieties, which have proved to possess an important antifungal effect, on the same chemical structure represents an important strategy for the development of new drugs, leading to compounds with improved biological activity.  Following that, nine new 1,4-phenylene-bisthiazole acylhydrazones were synthesized and evaluated in vitro for their inhibitory activity against different Candida strains. Moreover, computational approaches in drug discovery domain, including molecular docking and prediction of ADMET properties, are of great importance because they could reduce the number of experiments required, which are money and time-consuming processes [16]. The pharmacological activity of compounds is closely related to their pharmacokinetics and pharmacotoxicology. Among the interaction of a drug with its molecular target, properties such as absorption, distribution, metabolism, excretion and toxicity (ADMET) influence its efficacy and safety. In order to predict their binding poses towards CYP51 and to evaluate the drug-like properties of the synthesized compounds, a molecular docking study and an in silico ADMET prediction were also carried out.
computational approaches in drug discovery domain, including molecular docking and prediction of ADMET properties, are of great importance because they could reduce the number of experiments required, which are money and time-consuming processes [16]. The pharmacological activity of compounds is closely related to their pharmacokinetics and pharmacotoxicology. Among the interaction of a drug with its molecular target, properties such as absorption, distribution, metabolism, excretion and toxicity (ADMET) influence its efficacy and safety. In order to predict their binding poses towards CYP51 and to evaluate the drug-like properties of the synthesized compounds, a molecular docking study and an in silico ADMET prediction were also carried out.
The target compounds, 1,4-phenylene-bisthiazole acyl-hydrazones 5a-i, were obtained according to the route illustrated in Scheme 2. Following that, nine new 1,4-phenylene-bisthiazole acylhydrazones were synthesized and evaluated in vitro for their inhibitory activity against different Candida strains. Moreover, computational approaches in drug discovery domain, including molecular docking and prediction of ADMET properties, are of great importance because they could reduce the number of experiments required, which are money and time-consuming processes [16]. The pharmacological activity of compounds is closely related to their pharmacokinetics and pharmacotoxicology. Among the interaction of a drug with its molecular target, properties such as absorption, distribution, metabolism, excretion and toxicity (ADMET) influence its efficacy and safety. In order to predict their binding poses towards CYP51 and to evaluate the drug-like properties of the synthesized compounds, a molecular docking study and an in silico ADMET prediction were also carried out.
The hydrazide derivative 4 was subsequently synthesized by refluxing ethyl 4-methyl-2-(4-(2-methylthiazol-4-yl)phenyl)thiazole-5-carboxylate (3) with hydrazine hydrate in absolute ethanol (yield 69%). The final acylhydrazone derivatives 5a-i were obtained in good Analyzing the IR spectra of the newly synthesized acylhydrazones, we could identify specific signals such as: a sharp signal with medium intensity at 3112-3094 cm −1 , corresponding of the stretching vibration of the C 5 -H from the thiazole ring, a strong signal of high intensity at 1682-1658 cm −1 due to the presence of the C=O bond, an absorption band at 3281-3241 cm −1 due to the stretching vibration of the N-H bond, as well as a characteristic band at 1611-1598 cm −1 , which was attributed to the stretching of the C=N double bond.
Signals that are specific for each final compound could also be observed: two strong signals found between 1265-1252 cm −1 and 1039-1025 cm −1 , corresponding to the asymmetric and symmetric stretching vibration of the C-O-C ether bond in compounds 5a-c, a wide absorption band characteristic for the stretching vibration of the O-H bond at 3419 cm −1 for compound 5d and at 3418 cm −1 for compound 5e, a characteristic signal at 1231 cm −1 due to the C-F bond in compound 5f and strong signals between 743-701 cm −1 due to the C-Cl vibrations in compounds 5g and 5h.
The mass spectra recorded for all newly synthesized compounds revealed the presence of the molecular ion peaks (M + 1), as expected from their molecular formulas.
In the 1 H-NMR spectra of the key intermediate 4, we can observe two broad signals due to the NH and NH 2 protons at 9.61 ppm and 4.63 ppm, respectively. The absence of NH 2 protons signal, the shift of the signal corresponding to the NH proton from 9.61 ppm, in compound 4, to 11.75-11.93 ppm, in compounds 5a-i, and the presence of the N=CH proton and additional aromatic protons, confirmed the condensation between the carbohydrazide 4 and the corresponding aromatic aldehydes. In the 13 C-NMR spectra of compounds 5a-i, the presence of the C=N peaks from the azomethine group at 139.3-149.0 ppm proved once again that the synthesis of the acylhydrazones has been successfully accomplished.

Anti-Candida Activity Assay
The newly synthesized compounds were evaluated for their in vitro anti-Candida inhibitory activity, using four human opportunistic pathogenic Candida strains. The obtained results were compared with those of fluconazole, used as reference drug, and are illustrated in Tables 1 and 2.  The results obtained for the most active compounds are marked in bold.
The best inhibitory activity on both Candida albicans strains was obtained when the Ar substituent was a substituted phenyl ring, in orto position with a hydroxyl group (compound 5d) or in para position with a fluoro atom (compound 5f). Both compounds presented the same anti-Candida activity as fluconazole on C. albicans strains (MIC = 15.62 µg/mL, MFC = 31.25 µg/mL). Compounds 5d, 5f, 5g and 5h were the most active in our series on Candida krusei, showing MIC and MFC values equal to those of fluconazole (MIC = 15.62 µg/mL, MFC = 31.25 µg/mL). The Candida parapsilosis strain seemed to be less susceptible to almost all our tested compounds, compound 5g being the only one with an inhibitory activity equal to the one of our reference drug (MIC = 7.81 µg/mL, MFC = 15.62 µg/mL).
Some structure-activity relationships (SAR) could be observed based on the obtained MIC and MFC values. Compound 5d (2-OH-phenyl derivative) showed a 2-fold greater anti-Candida activity, when compared with compound 5e (4-OH-phenyl derivative), suggesting that the presence of the hydroxyl group in orto position of the phenyl ring increases the anti-Candida activity.
When Ar substituent group was represented by 4-F-phenyl or 4-Cl-phenyl, the corresponding target compounds 5f and 5g showed superior anti-Candida activity compared to compounds substituted with 4-OH-phenyl or 4-OCH 3 -phenyl. Therefore, para-substitution of the phenyl ring with a halogen atom could enhance the anti-Candida activity, fact which was also observed by Secci et al. [18] and by our group in our previous studies [12,19,20]. When Ar substituent group was 2,4-di-Cl-phenylene (compound 5h), a decrease in the inhibitory activity on Candida parapsilosis was observed. No improvement in the anti-Candida activity was observed when the phenyl nucleus from the Ar substituent group was replaced by the thienyl moiety (compound 5i).
Antifungal agents are considered to have fungicidal activity if the MFC/MIC ratio is ≤4 and fungistatic effect if the MFC/MIC ratio is >4. All our tested compounds have a MFC/MIC ratio equal to 2, indicating that these compounds behave as fungicidal agents [21].

Molecular Docking
In the last few years, lanosterol 14α-demethylase has been intensively studied for the discovery and development of novel antifungal agents, which are capable of inhibiting it [22,23]. The present study is focused on analyzing the binding affinity of our tested compounds 5a-i towards both the active site and the access channel of the enzyme. The predicted best binding conformation of each compound to the catalytic site of lanosterol 14α-demethylase, in terms of binding energy, expressed as the variation of Gibbs free energy (∆G), as well as the calculated inhibition constant (Ki) and the statistical analysis of the cluster containing the top binding conformation are presented in Table 3. All tested compounds 5a-i have good to very good binding affinity to lanosterol 14α-demethylase of Candida albicans, their binding energy values being under −11 kcal/mol. Statistically speaking, for many compounds in our series and especially for those with the best binding affinity to the target enzyme, it could be observed that most of the resulted conformations are in the same 2Å cluster with the best binding conformation.
A detailed energy analysis of the cluster of the top binding conformation provided that the average value of the Gibbs free energy is high and close to the best binding conformation found, fact that is well correlated with a low standard deviation of the values around the average value. A low standard deviation of the root-mean-square deviation of the three-dimensional Cartesian coordinates of each conformation from the top binding cluster found, led us to conclude that each individual pose is relatively close to others found in the same cluster. Moreover, the presence of a reduced number of clusters with residual poses, that are not found in the top cluster, is correlated with the idea that the high number of conformations found for each compound, in the top cluster, relates to the constant binding of the compounds in the same zone of the protein, interacting with the same amino acid residues. Thus, we can conclude that this in silico assay has good reproducibility.
In order to evaluate how different substituents, located on the Ar group, influence the binding affinity of the tested compounds, we have visually analyzed and compared the binding mode of the compounds to the target lanosterol 14α-demethylase. It could be noticed that the 1,4-phenylenebisthiazole moiety enters the depth of the catalytic site along the access channel, oblique to the active center of the enzyme, parallel to it at some distance. It is worth highlighting the difference in polarity of the binding regions from the access channel to the catalytic site of the enzyme. Two distinct areas of different polarity can be identified: a lipophilic, non-polar zone in the depth of the binding pocket, which contains the 1,4-phenylene-bisthiazole moiety, and a more polar one, located at the exit of the access channel to the active site of the enzyme, which contains the acylhydrazone group and the diverse substituted phenyl ring or thienyl moiety.
The polar area of the access channel contains polar amino acid residues (Tyr132, Tyr118, His377), with which the ligands can interact if the phenyl ring (Ar) is optimally substituted. Substitution of the nucleus with polar groups (o-OH, p-F) will lead to an optimal interaction between the molecules and His377 through a hydrogen bond, as shown in Figures 2 and 3.  Molecules 2019, 24, x FOR PEER REVIEW 7 of 15  Substitution of the aryl group with non-polar or voluminous substituents will lead to repulsion due to His377 residue, and a mild but significant repositioning of the molecule in the binding pocket. For instance, compound 5h is located nearer to the exit of the access channel, compared to 5d or 5f. Due to the particular positioning of the molecule and the consequent conformational changes in its structure, the occurrence of hydrogen bonds between Tyr118 and C=O from the acyl-hydrazone group could be observed as a novelty (Figure 4).   Substitution of the aryl group with non-polar or voluminous substituents will lead to repulsion due to His377 residue, and a mild but significant repositioning of the molecule in the binding pocket. For instance, compound 5h is located nearer to the exit of the access channel, compared to 5d or 5f. Due to the particular positioning of the molecule and the consequent conformational changes in its structure, the occurrence of hydrogen bonds between Tyr118 and C=O from the acyl-hydrazone group could be observed as a novelty (Figure 4).  Substitution of the aryl group with non-polar or voluminous substituents will lead to repulsion due to His377 residue, and a mild but significant repositioning of the molecule in the binding pocket. For instance, compound 5h is located nearer to the exit of the access channel, compared to 5d or 5f. Due to the particular positioning of the molecule and the consequent conformational changes in its structure, the occurrence of hydrogen bonds between Tyr118 and C=O from the acyl-hydrazone group could be observed as a novelty (Figure 4).  Substitution of the aryl group with non-polar or voluminous substituents will lead to repulsion due to His377 residue, and a mild but significant repositioning of the molecule in the binding pocket. For instance, compound 5h is located nearer to the exit of the access channel, compared to 5d or 5f. Due to the particular positioning of the molecule and the consequent conformational changes in its structure, the occurrence of hydrogen bonds between Tyr118 and C=O from the acyl-hydrazone group could be observed as a novelty (Figure 4). When the phenyl ring (Ar) is substituted with voluminous groups, the molecule gathers a long length and can no longer fit between the depth of the pocket and His377. Thus, different interactions could be observed between the compounds and the amino acid residues of the binding pocket, depending on how the phenyl ring is substituted, which could explain the differences in binding affinity of the tested compounds.
The newly synthesized compounds 5a-i are predicted to bind to lanosterol 14α-demethylase in a different manner than clinically approved azoles, as long as there was no interaction between the tested compounds and the heme Fe 2+ from the active site of the enzyme. All interactions which took place between the tested compounds and lanosterol 14α-demethylase involve amino acid residues from the access channel of the enzyme.

In Silico ADMET Evaluation
Drug-likeness of a molecule is a qualitative concept in designing new drugs, based on the theoretical calculation of certain parameters, giving a better assumption and understanding of the pharmacological activity of the compound [24]. Moreover, many drug candidates fail in late-phase studies due to an inadequate pharmacokinetic and pharmacotoxicological profile. Therefore, ADMET properties prediction has become an important approach in the drug discovery domain [25].
The drug-like properties of the synthesized compounds 5a-i, according to the Lipinski Rule of Five, are summarized in Table 4. The predicted water solubility, gastrointestinal absorption, blood-brain barrier permeability, P-glycoprotein and CYP2D6 enzyme inhibition are presented in Table 5.  Abbreviations: 1 logarithm of water solubility; 2 gastrointestinal absorption; 3 blood-brain barrier permeability; 4 P-glycoprotein.
All newly synthesized acylhydrazone derivatives obey, without any violations, the Lipinski Rule of Five, which provides that, for a good absorption and permeation, compounds should have a MW less than 500 Da, no more than 5 HBDs, 10 HBAs and mLogP value under 4.15 [26]. The molar refractivity for all compounds is less than 130.
Topological polar surface area (TPSA) is a descriptor that is correlated with passive transport across the cell membranes [27]. Compounds 5a, 5b, 5c, 5f, 5g and 5h, that contain a phenyl ring substituted with a methoxy group or one or two halogen atoms, have TPSA values under 140 Å 2 , being associated with higher gastrointestinal absorption, compared to compounds 5d, 5e and 5i, containing a phenyl ring substituted with a hydroxyl group or a thienyl residue, which have TPSA values above 140 Å 2 .
Percentage of gastrointestinal absorption (%Abs) varies from 56 to 66%, which indicates a moderate to good oral bioavailability. The ability of the compounds to pass the blood-brain barrier is absent, this being correlated with a reduced risk of developing side effects on the central nervous system, upon administration. Compounds are assumed to have a poor solubility in water, with logS values between −7.76 and −9.01. P-glycoprotein (P-gp), a member of the ABC superfamily, plays an important role in intestinal absorption after oral administration, in biliary and renal excretion, as well as in avoiding brain penetration of drugs. Being a substrate to P-glycoprotein can generate a lot of drug-drug interactions and seriously modify the bioavailability and safety of a certain drug [28]. Our tested compounds did not seem to be substrates for P-gp.
Moreover, azole antifungal agents are known for their lack of specificity for the fungal enzyme. CYP51 is part of cytochrome P450 superfamily, found both in human and fungal cells. Therefore, it is important to know the possible pharmacokinetic interactions in which the newly synthesized compounds can be involved [29]. All tested compounds act as non-inhibitors of CYP2D6, so drug-drug interactions with compounds that are metabolized by this isoform of CYP450 superfamily are less likely to appear.

General Information
Reagents and solvents used in this study were of analytical grade purity and purchased from Sigma-Aldrich Chemicals (Buchs, Switzerland) and Alfa Aesar (Karlsruhe, Germany). The progress of all reactions was checked by thin layer chromatography (TLC) using precoated Silica Gel 60F 254 sheets (Merck, Darmstadt, Germany). The spots on plates were visualized under UV light (254 nm). Melting points (m.p.) were determined using an MPM-H1 melting point apparatus (Schorpp Gerätetechnik, Überlingen, Germany) and are uncorrected. Elemental analyses for C, H, N and S were carried out on a Vario El CHNS instrument (Hanau, Germany). FT-IR spectra were recorded on a JASCO FT-IR-6100 spectrophotometer (Cremella, Italy). 1 H-NMR and 13 C-NMR analyses were performed on an Avance NMR spectrometer (Bruker, Karlsruhe, Germany), operating at 500 MHz and 125 MHz respectively, in DMSO-d 6 , using tetramethylsilane (TMS) as internal standard. Chemical shifts (δ values) are expressed in parts per million (ppm). M + 1 peaks were determined on an Agilent 1100 series and an Agilent Ion Trap SL mass spectrometer (Santa Clara, CA, USA), operating at 70 eV.

Anti-Candida Activity Assay
Evaluation of the anti-Candida activity was done according to the guidelines of Clinical Laboratory Standards Institute (CLSI) [30]. In order to determine the minimum inhibitory concentration (MIC) and minimum fungicidal concentration (MFC) values of the tested compounds, Candida albicans ATCC 10231, Candida albicans ATCC 18804, Candida krusei ATCC 6258 and Candida parapsilosis ATCC 22019 were used. All the tested fungal strains were obtained from Food Biotechnology Laboratory, Life Sciences Institute, University of Agricultural Sciences and Veterinary Medicine Cluj-Napoca, Romania. The cultures were stored on potato dextrose agar (Sifin, Germany). Prior to antifungal susceptibility testing, each strain was inoculated on potato dextrose agar plates to ensure optical growth characteristics and purity. The medium used for susceptibility testing was Roswell Park Memorial Institute (RPMI) 1640 with L-glutamine, adjusted to pH 7.0 with 3-(n-morpholino)propanesulfonic acid. The initial density of Candida spp. was approximately 2 × 10 6 colony forming units/mL (CFU/mL). Inoculums (density of 0.5 in McFarland scale) were prepared in a 0.9% NaCl sterile solution. Then, tested strains were suspended in RPMI 1640 medium, to give a final density of 2 × 10 3 CFU/mL. Stock solutions (1 mg/mL) were obtained by dissolving the newly synthesized compounds and the reference antifungal, fluconazole, in sterile DMSO. These solutions were stored at 4 • C. A series of double diluting solutions of the above compounds were prepared in RPMI 1640 medium obtaining final concentrations in the range of 500-0.015 µg/mL. The broth microdilution method was employed for the minimum inhibitory concentration test. The medium was placed into each of the 96 wells of the microplates. Sample solutions at high concentration were added into the first rows of the microplates and two-fold dilutions of the compounds were made by dispensing the solutions into the remaining wells. 10 µL of the culture suspensions were inoculated into the wells. The growth control, sterility control and control of antifungal compound were used. Plates were incubated under normal atmospheric conditions at 30 • C for 48 h. Minimum inhibitory concentration (MIC) values have been determined by adding resazurin (20 µL, 0.02%) followed by incubation for 2 h. The MIC was defined as the lowest concentration required arresting the growth of the fungal strain. For determination of minimum fungicidal concentration (MFC), 0.01 mL of the medium withdrawn from the culture tubes, showing no macroscopic growth at the end of the 24 h, was subcultured on potato dextrose agar plates to determine the number of vital organisms and incubated further at 30 • C for 48 h. The MFC was defined as the lowest concentration of the agent at which no colonies are observed.

Molecular Docking
The molecular docking study targeted the compounds 5a-i into the active site of fungal lanosterol 14α-demethylase, using AutoDock 4.2.6 [31], in order to evaluate their binding affinity for the protein and to predict the binding poses of the new acylhydrazone derivatives. Furthermore, focusing on understanding the differences between compounds in terms of substituents' variability and how different groups might influence the protein-ligand interactions, we decided to keep the same target protein and docking protocol as reported in our previous work [17]. The macromolecular target was built by homology modeling using the three-dimensional structure 5EQB from Protein Data Bank (www.rcsb.org). The search space size (x, y, z) was set to 75 × 75 × 75 points with a 0.375 Å grid spacing, with the center Cartesian coordinates set to x = −23.134, y = 13.943, z = 19.959.
The unique change in our protocol was the increasing of the generated conformations to 100 for each ligand, in order to increase the precision and to obtain more realistic results. The inhibition constant (Ki) was calculated using the formula K i = e ∆G×1000 R×T , where ∆G represents the computed binding affinity energy, R = 1.98719 cal mol −1 K −1 (Regnault constant) and T = 298.15 K = 25 • C.
Clustering analysis of the resulted conformations was performed in order to analyze and to evaluate the homogeneity of binding of the ligands to the catalytic site of the enzyme [17,32]. Visualization of the docking results was performed using UCSF Chimera [33].

In Silico ADMET Evaluation
Theoretical in silico evaluation of ADMET properties of compounds 5a-i was performed using the Swiss-ADME online platform [34,35]. In this study, we assessed the following descriptors, taking into consideration the Lipinski Rule of Five: molecular weight (MW), number of rotatable bonds (RoB), hydrogen bond acceptors (HBA) and hydrogen bond donors (HBD), Moriguchi LogP (mLogP) and molar refractivity [26]. Water solubility, topological polar surface area and gastrointestinal absorption, which are predictors of the oral bioavailability of a certain compound, the capacity to inhibit CYP2D6 and P-glycoprotein, which are involved in many drug-drug interactions, and the ability of the compounds to cross the blood-brain barrier were also investigated [27]. Percentage of gastrointestinal absorption (%Abs) was calculated using the formula: %Abs = 109 − 0.345 × TPSA [36].

Conclusions
Nine new acylhydrazone derivatives 5a-i were designed, synthesized, fully characterized and evaluated for their inhibitory activity against different Candida spp. In the in vitro anti-Candida activity assay performed, compounds 5d and 5f exhibited the best inhibitory activity against both C. albicans strains and compound 5g displayed the highest inhibition against C. parapsilosis. The most susceptible fungal strain to the newly synthesized compounds was C. krusei, four out of nine compounds being as active as the reference drug used, fluconazole, against this Candida strain.
A molecular docking study was performed in order to investigate the binding poses of the new compounds 5a-i towards lanosterol 14α-demethylase. The obtained results from the docking simulations showed that our compounds interact with the amino acid residues from the access channel to the active site of the enzyme, with no interaction with amino acids from the catalytic site.
Analyzing the results of the in silico ADMET screening, it can be observed that the new synthesized acylhydrazone derivatives displayed favorable drug-like properties according to Lipinski Rule of Five, making them suitable for development as oral drug candidates. Moreover, the compounds are not predicted to pass the blood-brain barrier or to inhibit CYP2D6 and P-glycoprotein, so the risk for central nervous system side effects and drug-drug interactions is considerably low. All in all, the above-mentioned observations provide a promising basis for the development of new anti-Candida agents.