Solvatochromism and Thermochromism of Fluconazole: An Experimental and Theoretical Study

This research aims to study fluconazole (FNZ)-solvent interactions in pure solvents through theoretical and experimental solvatochromism. This study showed a bathochromic shift as the polarity of the solvent increases, where the excited state of FNZ is of higher polarity and lower energy than that of the ground state. The multiparametric statistical analysis highlighted solvent polarizability, dispersion, and electronic induction as solvent parameters of great importance, with the basicity of solvent and hydrogen bond acceptor capacity as having certain relevance. The thermochromic effect was also determined by exposing the three pKa values of FNZ. In conclusion, this study shows the importance of the π-π stacking interaction of the FNZ dihalogenated phenyl ring, the solvent basicity for the hydroxyl group, and its acidity for the N4’ of one of the triazoles stand out, as well as the acid-base equilibria involving the -OH group and the two N4’.


Introduction
Invasive fungal infections represent a serious and ever-present threat to human health and are associated with at least 1.5 million deaths worldwide each year.Due to their broad spectrum of activity, azoles (imidazoles and triazoles) are the most common antifungal drugs in clinical use.][3] According to their date of discovery and subsequent structural modifications, the azoles are currently divided into four classes or generations.Second-generation triazoles, including fluconazole (FNZ) and itraconazole, exhibit a broader spectrum of antifungal activity compared to imidazoles.5][6] Because it is one of the most widely used antifungal drugs and a better understanding of its physicochemical properties could lead to new therapeutic targets, we selected fluconazole for this study.
Also, there is a new successful approach to discovering effective drugs to treat several existing diseases through drug repositioning, which means the use of a drug in an indication other than the one for which it was initially marketed. 7Several drugs have been successfully repositioned to a new indication, with the most prominent being sildenafil and thalidomide, which have generated historically high revenues. 8More and more companies are scanning the existing pharmacopeia for repositioning candidates, and the number of repositioning success stories is increasing. 9FNZ, although an old triazole, could be a great candidate to be repositioned as a new therapeutic drug, as it has been shown the increasing interest in FNZ, voriconazole and itraconazole due to the pandemic disease caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus. 10igure 1 shows the chemical structure of FNZ, 2-(2,4-difluorofenil)-1,3-bis(1,2,4-triazole-1-il)propan-2-ol.
Active pharmaceutical ingredients (APIs) interact with multiple environments throughout the entire manufacturing process and subsequent assimilation by the body. 11Therefore, it is essential to know the properties of APIs in different microenvironments to predict their pharmacokinetic and pharmacodynamic properties and to provide valuable information for all pharmaceutical steps. 12,13t is known that the microenvironment of the solution that surrounds a solute molecule, or cybotactic region, exerts a great influence on its physicochemical properties.The modifications that occur in the absorption spectra for certain organic molecules, depending on the polarity characteristics of the environment, are known as the solvatochromic effect. 12,14Its analysis makes it possible to measure the difference between the ground state and the excited state of a molecule, where the solvent acts as a constant electric field applied to the solute in question, with an intensity directly proportional to the magnitude of its dielectric constant.The extent and direction of solvatochromism (bathochromic effect, positive or red solvatochromism: shift to longer wavelengths; hypsochromic effect, negative or blue solvatochromism: shift to shorter wavelengths) depends on how stable the structure of the solute is in the ground state and in the excited state. 15,16It is known that the conformation with the highest dipole moment is the most stable in solvents with a high dielectric constant. 17][17][18][19] Thermochromism is the reversible transformation of a molecular structure or system, induced by changes in the temperature of a solution, liquid or solid, and which produces a spectral variation.The mechanism responsible for thermochromism varies with the molecular structure.This phenomenon can be caused by an equilibrium between two molecular species (acid-base, keto-enol, lactim-lactam), a ring opening, a thermally accessible triplet state or the formation of free radicals, among others.Thermochromic studies are carried out by determining the absorbance of a solution at different temperatures.The generation of isosbestic points indicates the presence of an intermolecular equilibrium.][22] The stability of the FNZ molecule has been extensively investigated through conformational studies, electronic delocalization, electrostatic potentials, and thermodynamic parameters, as well as the effect of some solvents on its UV-Vis absorption spectrum. 19,23Similarly, various polymorphs of FNZ have been analyzed in the crystalline state, proving its great conformational flexibility. 24owever, there is no evidence in the literature consulted of an analytical study that analyzes the interaction of solvents of different polarities with FNZ, nor its thermochromism at physiological pH.
In another way, theoretical approaches are common in the study of electronic structure to predict and understand the properties of molecular systems.A frequently used computational method is the density functional theory (DFT) which is based on the determination of the energy of an electronic state from the electronic density. 25The free energy of solvation calculates the energy change in the transfer of a molecule from the gas to the solvent and is used to calculate a variety of properties such as activity coefficients and solubilities. 268][29][30] In addition, the polarizable continuum model (PCM) is a method commonly used in computational chemistry to simulate the effects of solvation.While Integral Equation Formalism PCM (IEFPCM) is the current version of PCM that is applied in the usual quantum chemistry packages, the conductive polarizable continuum model (CPCM) variant of PCM is often considered one of the most popular successful solvation models. 29,31FT has introduced descriptors of global and local chemical reactivity.In the case of global descriptors, these provide information about the reactivity of the molecule and are calculated by using the Koopmans theorem, considering the energies of the frontier orbitals highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO). 26,32Some usual calculations to estimate the chemical reactivity of a molecule are energy gap, hardness (η), chemical potential (μ), electrophilicity (ω), electrodonating, and electroaccepting power. 33,34his research aims to study the solute-solvent interactions of the antimycotic azole FNZ in pure solvents through theoretical and experimental solvatochromism with multiparametric statistical analysis.The thermochromic effect on the UV-Vis spectrum in a buffer solution compatible with physiological pH is also determined.These two methodologies were not performed in the current published literature.The analysis of the physicochemical properties of FNZ may lead to a better understanding of its structural attributes and impact on biological action. 29,30,35

Quality control and preparation of experimental samples
The API and solvents used were of the highest analytical purity.FNZ was purchased from Parafarm (Buenos Aires, Argentina), and all the reagents and solvents used for spectroscopic studies were of the commercially available highest grade (Parafarm, Buenos Aires, Argentina).The absorption spectra were obtained by using a Varian Cary 50 Conc UV-Vis spectrophotometer equipped with a Peltier thermostatic and a fiber optic system (Varian, Mulgrave, Victoria, Australia) The radiation source was a xenon lamp, and a 1 cm optical path quartz cuvette was used.For FNZ quality controls, the melting point was determined in triplicate with an Opti Melt MPA100 automatic melt meter applying a heating ramp of 1 ºC min -1 and a temperature range of 130 to 160 ºC (Stanford Research System, Sunnyvale, California, USA).Additionally, analytical purity was determined by thin layer chromatography (TLC) and high-performance liquid chromatography (HPLC).For the TLC, a Merck silica gel 60 F254 aluminum plate (Darmstadt, Germany, Merck S. A.) and a mobile phase of chloroform-ethyl acetate-glacial acetic acid-water in a 4:4:4:1 ratio was used. 36The HPLC was performed using a Shimadzu LC-20AT equipment with an SPD-M20A diode array UV-Vis detector (Kyoto, Japan, Shimadzu), a Phenomenex Gemini C18 column 150 mm long stationary phase of 4.6 mm internal diameter and a particle diameter of 5 μm (Phenomenex, Torrence, California, USA).The mobile phase was a 19:81 acetonitrile-water solution at a flow rate of 0.5 mL min -1 .The preparation of FNZ samples was carried out by Argentinian Pharmacopoeia (FA) VII edition. 37The quality of the solvents was verified by measuring the refractive index with a Reichert AR200 pocket digital refractometer (Leica, Deerfield, Illinois, USA).
The buffer solution was prepared according to the FA 38 using 0.2 M monobasic potassium phosphate solution and NaOH until a solution of pH 7.4 was achieved.It was validated with a JP Selecta pH-2005 digital benchtop pH meter (Mettler Toledo, Leicester, United Kingdom).The solvatochromic study of FNZ was carried out with 1.0 × 10 -3 M solutions that comply with the Lambert-Beer law. 19Seven pure polar solvents, three protic (2-propanol (2-PrOH), ethanol (EtOH), methanol (MeOH)), and four aprotic (diethyl ether (DE), dichloromethane (DCM), acetonitrile (AcN), dimethyl sulfoxide (DMSO)) and a pH 7.4 buffer solution were used.The thermochromic study was carried out with a 1.0 × 10 -3 M FNZ solution in a pH 7.4 buffer.All spectra were determined at 20.0 ± 0.1 ºC.The experiments were performed in triplicate and average results were used for data analysis.

Solvatochromic and thermochromic studies
The distinction between specific and nonspecific solutesolvent interactions in the interpretation of experimental determinations of absorption spectra is difficult.Therefore, it is necessary to consider quantitative measures of polarity to differentiate between these two effects.Thus, several multiparametric equations have been developed that combine two or more solvent parameters to explain the particular and distinctive aspects of the solvent effect on a given compound, 27,[39][40][41] among them, those of Kamlet and Taft, 12,13 Catalán 40 and Laurence 41 are the most employed.
][44][45][46] The Kamlet and Taft 44 method rationalize the effect of the solvent in terms of a linear combination of three variables: the parameter p* refers to the shared effect of the polarity-polarizability of the solvent, which measures the ability of the environment to stabilize the charges of a dipole due to its dielectric effects.This parameter relates to nonspecific interactions.The parameter α describes the acidic capacity of the solvent to donate a proton through hydrogen bonding (HBD) and the parameter β measures the basic ability of the solvent to accept hydrogen bonds (HBA).The last two parameters imply specific interactions. 46hese solvatochromic parameters are related with the following equation: where A ( = 1/λ / cm -1 ) is the property of the solute to be correlated; A 0 is the value of this property for the same solute in the gas phase or in a hypothetical solvent for which α = β = p* = 0; 13,44 a, b and s are the coefficients that calculate the relative susceptibility of A to the parameters already defined: α, β and p*, respectively.
The Catalán 47 equation is an improvement in the multiparametric analysis concerning the Kamlet and Taft 45 equation. 13,39,47While the latter was made from data belonging to an average of around 40 determinations made by different authors in various laboratories, the former was created from data obtained in the same laboratory for pairs of homologous structures. 40oth the Catalán SPP descriptor and the Kamlet and Taft π* descriptor explain the polarity of the solvent, although the π* parameter was obtained from π-π* transitions of some azo dyes.The SPP descriptor was calculated from different transitions of many compounds.Even though multiparametric procedures have provided excellent results in the analysis of the effect of the solvent, there is also experimental evidence of some relevant deficiencies such as their inability to describe the solvatochromic behavior of nonpolar solutes that do not exhibit any specific interaction with the solvent. 47Therefore, Catalán 47 separated the SPP polarity scale into two parameters: the solvent polarizability (SP) and dipolarity (SdP) scales, resulting in a four-parameter model. 39According to this model, the specific interactions are described by the parameters SA, which shows the solvent acidity, and SB, which provides the measurement of the solvent basicity.Nonspecific interactions are now characterized by two independent parameters: SP and SdP.The equation is given by: where A ( = 1/λ / cm -1 ) is the solvent-dependent physicochemical property in a given solvent; A 0 is the statistical quantity of the value of the property in the gas phase; SA, SB, SP and SdP, represent the independent parameters, although complementary, that account for solute-solvent interactions; a SA , b SB , c SP and d SdP are the regression coefficients that measure the relative susceptibility of to the different interaction forces. 47n the other hand, the Laurence equation allows the estimation of the dispersion and electrostatic solute-solvent interactions of a given system with greater precision. 39,41It is based on both experimental measurements and theoretical calculations. 39According to this model, the nonspecific intermolecular forces are determined by the DI parameter, which describes the dispersion and induction interactions, and the ES parameter, which describes the electrostatic interactions between the permanent multipoles of the solute and the solvent.While the specific interactions are described by the parameter α 1 , which quantifies the solute HBA/solvent HBD interactions (the hydrogen-bond acidity of solvent); and the parameter β 1 , that accounts for the solute HBD/solvent HBA interactions (the hydrogen-bond basicity of solvents). 41his parametric analysis is determined by the equation: where A ( = 1/λ / cm -1 ) is the property of the solute to be correlated, A 0 is the statistical quantity of the value of the property in the gas phase.DI, ES, α 1 and β 1 represent the independent, although complementary, parameters that account for solute-solvent interactions; di, e, a 1 and b 1 are the regression coefficients that evaluate the relative susceptibility of A to the different interaction relationships. 27,41n solvatochromic and thermochromic studies, the UV-Vis absorption spectrum was measured from 225 to 325 nm, at a scanning speed of 1 nm s -1 .The solvatochromic samples were kept at a constant temperature of 20 ºC for 5 min before the spectroscopic determination.The thermochromic samples were thermostated for 5 min before each determination, taking the data with a gradient of 5 ºC in the temperature range between 5 and 75 ºC.The / cm -1 of each spectrum was measured by taking the midpoint between two positions of the spectrum where the absorbance is equal to 0.9 . 22,45he results were analyzed by the Kamlet and Taft method, the Catalán method, and the Laurence method.The wavenumber of maximum absorption, , was related to the empirical solvent parameters of Kamlet-Taft equation: α, β, p*, Catalan equation: SA, SB, SP, SdP, and Laurence equation: DI, ES, α 1 , β 1 .The correlation study was carried out using multiparameter linear regression data analysis, taking as the dependent variable and the former parameters as independent variables (Table 1).

Geometric optimization, conformational analysis, and molecular descriptors
The FNZ structure was traced using Gaussview 5.0 48 software and optimized with geometric calculations using Gaussian 9.0 software. 49Initially, the geometric optimization was carried out employing Molecular Mechanics universal force field (UFF) to obtain the most stable structure using the semi-empirical methods Austin Model 1 (AM1) and Parametric Method 3 (PM3).
The conformational analysis of FNZ reported by Chandrasekaran and Thilak Kumar 19 evidenced the relevance of the dihedral angle C5-C10-O3-H32 (Figure 2).Rotation around the C5-C10 single bond produces different conformers of the molecule since C10 is attached to a hydroxyl group, a phenyl ring, and methylene groups.Hence, the conformational scan of the C6-C5-C10-O31 dihedral angle was performed, rotating the eclipsed phenyl ring with the C10-O31 bond every 15º in the range of 0-360º and calculating the potential energy of every conformation.Subsequently, the minimum energy of this conformational scan was selected and the dihedral angle C5-C10-O31-H32 was analyzed, varying the position of H32 every 15º in the range of 0-360º until a new minimum was obtained.The FNZ structure obtained was optimized by ab-initio methods HF/3-21+G*, HF/6-31+G(d), DFT-B3LYP/3-21+G*, DFT-B3LYP/6-31+G(d) and DFT-CAMB3LYP/6-31+G(d).Since the geometric structure found by the aforementioned authors 19 was a conformer with the halogenated phenyl group in the opposite position to the hydroxyl group, this geometry was emulated by rotating it exactly 180º in relation to the initial molecule and optimized by semi-empirical methods. 19Then, optimizations were obtained by using the ab-initio methods already mentioned.Two slightly different FNZ conformers were obtained (Figure 2).
Frequency calculations were performed on the structures obtained by DFT-B3LYP and DFT-CAMB3LYP of FNZ conformers I and II to verify whether these conformations corresponded to energy minima according to their vibrational states.Next, the potential energy, the dipole moment, the wavelength of maximum absorbance, and the energies corresponding to the HOMO and LUMO frontier orbitals were calculated by Time-dependent selfconsistent field (TD-SCF) methods both in a vacuum and in solvation sphere with the Conductor-like Polarizable Continuum Model (CPCM), choosing ten solvents whose dipole moments are included in Gaussian. 49nalysis of the geometric structures of FNZ conformers I and II was performed by measuring specific bond lengths, angles, and dihedral angles.They were compared with those described by Chandrasekaran and Thilak Kumar 19 and some crystallized polymorphs of FNZ. 19,24 The global molecular descriptors hardness (η), chemical potential (μ), electronegativity (ω), electrondonating power (ω -) and electronaccepting power (ω + ) were calculated using the Koopmans procedure in DFT, in which EA (electronic affinity) = -E LUMO and IE (ionization energy) = -E HOMO , resulting in the following parameters: Theoretical solvatochromic study Solvatochromic analysis of FNZ was performed at the theoretical maximum absorbance wavelengths obtained by DFT-B3LYP and DFT-CAMB3LYP of conformers I and II.For this, the multiparametric equations of Kamlet and Taft, Catalán, and Laurence were used.First, the same solvents used in the experimental study were chosen, whose dipole moments were included by default in Gaussian (it was excluded 2-PrOH and buffer pH 7.4 and was added water instead).Second, cyclohexane (Cy), tetrachloromethane (Tet), and benzene (Ben) were incorporated into the database to include the FNZ solvation effect in non-polar solvents, constituting a set of ten solvents.Finally, the relative contributions of the specific and non-specific interactions were compared in each case.

Geometric optimization, conformational analysis, and molecular descriptors
Conformational scanning of the C6-C5-C10-O31 (θ) and C5-C10-O31-H32 (j) dihedral angles were performed, as well as the calculation of the potential energy of each FNZ conformer before molecular optimization by ab-initio methods.The energy minima were obtained with θ = 195º and j = 135º for conformer I and with θ = 0º and j = 15º for conformer II.
After geometric optimization, the potential energies in vacuum calculated at the DFT-B3LYP/6-31+G(d) level of the FNZ conformers were E I = -2,903,855.09kJ mol -1 a n d E II = -2 , 9 0 3 , 8 7 7 .4 6 k J m o l -1 , w h i l e a t the DFT-CAMB3LYP/6-31+G(d) level they were E I = -2,902,607.37kJ mol -1 and E II = -2,902,629.36kJ mol -1 .It can be noted that the lowest potential energy minima were obtained with the B3LYP calculation method, instead of with CAMB3LYP.
The frequency analysis of the vibrational modes belonging to the optimized molecular geometries for each FNZ conformer was calculated at the same theoretical level described.All observed frequencies were positive (real numbers), verifying that they correspond to true minima on the potential energy surface of the system.
The results of the structural analysis of the FNZ conformers are shown in Table 2.They are similar to those reported by Chandrasekaran and Thilak Kumar, 19 structurally optimized at the B3LYP/6-311++G(d,p) level and compared with X-ray diffraction analysis (XRD) crystallographic data. 19hortening of the N-N bonds of the triazoles would indicate conjugation of the semicarbazoles.Also, the elongation of the C-N bonds would suggest intense electronic delocalization throughout the molecule.It can be seen from the dihedral angles C10-C5-C6-C1 and C10-C5-C4-C3 that the hydroxyl group attached to the difluorophenyl ring is slightly rotated, with the B3LYP method, 2.224º for conformer I and 3.138º for conformer II, and with the CAMB3LYP method, 2.217º for conformer I and 3.340º for conformer II, respectively.The methylene group attached to the triazole rings is even more rotated due to electronic conjugation, reflected in the dihedral angles C11-N22-N27-C24, C11-N22-C23-N28, C14-N17-N29-C19 and C14-N17-C18-N30.The triazole ring of FNZ conformer II that contains the N17-N20 bond showed a 180º rotation.These results are in line with those reported by the cited authors. 19Intramolecular charge transfer (ICT) from the electron-donor substituents to the phenyl ring could also be corroborated due to the shortening of the C-N bond distances, the lengthening of all C-C distances, and the contraction of the C4-C5-C6 internal angle.This is correlated with the increase of the two angles C5-C6-C1 and C5-C4-C3 in the union of the difluorophenyl ring and the triazole groups through the methylene group.
Dihedral angles molecules corresponding to seven conformers of four polymorphs crystallized by Karanam et al. 24 were compared.Similarly, the theoretical calculations of two new conformers were contrasted with their corresponding sets of torsion angles, obtained by the same research group. 24None of our FNZ conformational structures fully agree in their torsion angles with those described in the literature. 24However, all the torsion angles vary in a very small range (2-50º) which is equivalent to differences of about 2 kcal mol -1 .This indicates that the FNZ has a high tendency to form numerous conformers in the gas phase.When these conformers crystallize in the form of crystal lattices, they can lead to the formation of different polymorphs by small changes in conformation. 24herefore, in theory, it would be possible to obtain crystalline polymorphs of conformations I and II of FNZ calculated as energetic minima.
Tables S1, S2, S3 and S4 (see Supplementary Information section) show the results of some molecular descriptors calculated based on the analysis of the frontier orbitals of the conformational structures I and II of FNZ, after being optimized with the bases B3LYP and CAMB3LYP.Data are ordered from vacuum to solvated states, by increasing solvent polarity, first nonpolar, then polar nonprotic, and finally polar protic.The effect of the various solvents on the FNZ molecule describes very well its energetic behavior, when passing from the gaseous phase to the solvent.A small HOMO-LUMO gap means that the molecule is more reactive (soft molecule), while a large gap implies less molecular reactivity (hard molecule).Similarly, a small gap implies that the electron density of the molecule will change more easily than that of a hard molecule.In concordance with the previous parameters, the electrophilicity index is an indicator of energy reduction due to the electronic jump between the donor HOMO and the acceptor LUMO.On the other hand, electronegativity measures the power of attraction of electrons, classifying molecules as strong, moderate, or marginal.
Similarly, a high value of the electroaccepting power descriptor is indicative of a high capacity for electronic acceptance, while a low value of electrondonating power involves a better molecular system for electronic donation.Finally, the chemical potential is the amount of resistance to electron density loss. 26,32It has been reported that, with increasing solvent polarity, the total molecular energy of the solvated molecule should decrease and thus increase the stability of its molecular structure, which is evidenced by the rise in solute polarity. 19t is observed that the energy gap, hardness, and softness of FNZ I and II increase slightly as the polarity of the solvent increases, regardless of the optimization method used, which is indicative of the hardness of the API's structure.Therefore, FNZ has high chemical stability.However, the chemical potential decreases when the API goes through the gas phase to a solvated one.This is caused when the energetic transformation between the HOMO and the LUMO takes place.Considering FNZ's electronegativity, it varies slightly from vacuum to more polar solvents and its reduction is so small that it is not indicative of the ionization of covalent bonds in the solution, classifying the molecule as marginal.
Both the electroaccepting and the electrodonating power decline with solvation and when the polarity of the solvents rises.This reveals a certain drop in the electronic acceptance capacity and an improvement in the tendency of the system to donate electrons in solution.Finally, the dipolar moment of the solute increases both, in the solvated phase and the increase of solvent polarity, presenting the highest maximum in DMSO and water, with important polarity within the selected non-protic and protic polar solvents, respectively.The evidence that FNZ has such a high dipolar moment in DMSO, a non-protic solvent, is indicative of the high polarizability inductive effect of the polar solvent affecting the API.These results are in agreement with the reference literature published by Chandrasekaran and Thilak Kumar. 190][41][42] The low solubility of FNZ in non-polar solvents prevents obtaining experimental data, therefore, they are used in the theoretical calculations of solvatochromism to have relevant information on the behavior of the API in non-polar solvents such as cyclohexane, with dipolar moment (µD) = 0 and parameters α = β = p* = 0, according to equation 1.The UV-Vis spectra of FNZ in solution are observed in Figure 3, with two absorption peaks at 261 and 266 nm.
Chandrasekaran and Thilak Kumar 21 carried out an exhaustive theoretical study describing the effects of the solvent on FNZ.However, they only experimentally determined absorbance in three polar solvents: DMSO, methanol, and water.The authors relate the absorbance at λ 266 nm due to π-π* transition and the one at λ 261 nm to a n-π* transition.We disagree in this matter and state that both correspond to π-π* transitions.This is because both energy absorbances at 261 and 266 nm are very similar, as can be seen by the molar absorption coefficient in Table 3.This means that the nature of the electron jump is equivalent, and it is in the range of allowed or partially allowed electronic transitions, that correspond to π-π* nature, as can be seen by log ε 266 ca.2.8 and log ε 261 ca.2.8.A n-π* electronic transition is mostly forbidden in the UV-Vis range analyzed by symmetry selection rules, as a result of the different spatial regions of each orbital that leads to poor overlapping between them.Thus, this kind of electronic excitation appears at a higher wavelength and usually with little ε max , showed as a smaller absorption peak in the UV-Vis spectrum than a π-π*. 50The typical values of molar absorption coefficients (ε) for some electronic transitions in organic molecules are around: 10 4 -10 5 cm -1 M -1 for π-π* in spin allowed of highly conjugated organic molecules, 200-500 cm -1 M -1 for π-π* in spin allowed with symmetry forbidden for small aromatic compounds, and 10-50 cm -1 M -1 for spin allowed n-π* for aromatic compounds with heteroatoms.Due to the high symmetry of FNZ, some π-π* transitions may be symmetry forbidden and show significantly reduced molar absorption coefficients. 51cording to the experimental data, a slightly bathochromic shift is observed as the polarity of the solvent increases.The change is greater when FNZ is surrounded by a polar solvent rather than non-polar one.It is caused by the polarization forces of attraction between the solvent and the API, which decreases the potential energy of the ground and excited states.This energy drop is higher for the excited state than for the ground state in π-π* electron excitation.Therefore, the difference between the two levels decreases, resulting in a bathochromic shift. 19owever, the solvent also affects the n-π* electronic transitions, causing a hypsochromic shift due to unpaired electron pairs.For this type of transition, the ground state is much more stable due to the effect of polar solvents on FNZ (hydrogen bonds forces).Consequently, the energies of the transition states are increased for polar solvents, with the subsequent hypsochromic shift of the spectrum.
In conclusion, a mixed effect results from the predominant interaction forces caused by a change in the solvent polarity and its ability to form hydrogen bonds.The signs and the relative contributions of the different coefficients of the solvation parameters in equations 1, 2, and 3, allow for discriminating opposite effects in the position of the molecular structure absorption maximum. 46onsidering the experimental solvatochromic shifts derived from the electronic transitions π-π* illustrated in Table 3 and taking acetonitrile solvent (λ max = 266.7 nm) as a reference, we observe that when FNZ goes from aprotic solvents of lower polarity such as diethyl ether and dichloromethane to acetonitrile, hypsochromic shifts of 0.1 nm are produced in both cases, while towards DMSO the bathochromic shift is 0.7 nm.
From acetonitrile to protic polar solvents (2-propanol, ethanol, and methanol) the bathochromic shifts for the first and the second are 0.7 and 0.1 nm, respectively, while the hypsochromic shift for the third is 0.2 nm.Finally, when FNZ moves to the diprotic solvent water (in buffer pH 7.4 solution) a hypsochromic shift of 0.4 nm can be observed.
Solvents may be classified by different physical constants.Using viscosity (η) as a criterion, solvents are of low viscosity when their dynamic viscosity is < 2 mPa s at 20 ºC, of medium viscosity between 2 and 10 mP s, and of high viscosity above 10 mPa s. 52 Figure 4 shows a linear relationship between the maximum wavelength of FNZ and the viscosity of the solvents, where the decreases as the η of the solvent 3. UV-Vis spectra of FNZ in different organic solvents.increases, although there are some exceptions as acetonitrile, dichloromethane, and 2-propanol.This data may indicate that the increasing viscosity of the solvent would be related to a hypsochromic behavior of the API, even though a mixed solvathochromism effect is caused by the change in the polarity of the solvent.Usually, no simple relationship with solvent viscosity was found when there has been an attempt to predict the solvent effect in the efficiency of a reaction due to an effective rate constant. 52able 4 shows the coefficients and errors associated with all the parameters of the multiple linear regressions in the FNZ solvatochromism analysis performed using the Kamlet and Taft, Catalán and Laurence equations.Also, the coefficient of determination, the relative contribution of each parameter, and the p-value of the Fisher statistic for the linear regression are reported.The results obtained both at λ 261 nm and λ 266 nm are quantitatively similar; however, all the parameters have statistical significance at 266 nm.
The triparametric equation 1 has a single parameter that accounts for the polarity-polarizability of the solvent, π*, while the parameters α and β describe the acidity and basicity capacities of the solvent to donate or accept a proton via a hydrogen bond, respectively.
Our results explain 87% of the variability of in different polar solvents.The greatest influence of the solvatochromism observed is due to the specific interactions, around 85% of the whole relative contributions.The β parameter accounts for around 50% of the total effect, while the α parameter explains the remaining 35% effect.Thus, the importance of the hydroxyl group attached to the methylene of FNZ, and the unpaired electron pairs of O and N can be deduced.The influence of non-specific interactions (given by the dihalogenated phenyl ring) is relatively low, although its influence on the global contribution to the solvatochromic effect is not neglected.The fact that the coefficients s and b have negative signs implies that the excited state is more stable as the polarity of the solvent increases, corresponding to a bathochromic shift.Contrariwise, the positive sign of the coefficient a shows the opposite influence. 53Consequently, the net effect observed is a bathochromic behavior, given by the magnitude and sign of the most relevant coefficients, that is, the increase in the polarity of the polar solvents stabilizes the excited state of the electronic transitions of FNZ more than the basal state.
Equations 2 and 3 are tetraparametric, therefore, both have the advantage of differentiating the contributions of the non-specific interactions of the solvent on the  API in two different parameters, improving the previous analysis.The outcomes obtained from the analysis of the Catalán equation at 266 nm allow us to explain 97% of the variability of in the solvation medium.It is remarkable that the coefficient of the parameter π* of equation 1 was the least relevant of all, while in this case the magnitude of the parameter SP, the solvent polarizability, represents the largest relative contribution to the position of the maximum absorbance band.Regarding the order of magnitude, the basicity (SB), the acidity (SA), and the dipolarity (SdP) of the solvent, respectively showed statistical significance.The relative contribution of 65% of the -c SP coefficient, 21% of -b SB and 6% of -d SdP contribute to the bathochromic effect, contrasted to 8% of the a SA coefficient, with a positive sign.Similar findings were achieved with Laurence's equation 3.In this case, the parameter π* is divided into two parameters, DI, representing the dispersion and electronic induction interactions, and ES, the electrostatic interactions between the permanent dipoles of the solute and the solvent.
The specific interactions are represented by α 1 and β 1 , the acidity and basicity of the solvent, respectively.The analysis yielded results like those obtained using the Catalán equation and allows us to explain 98% of the variability of in the solvation medium, with a 67% relative contribution of the DI parameter.The contributions of solvent basicity (β 1 ), acidity (α 1 ) and permanent dipole interactions (ES), showed better statistical significance of the linear regression.Similarly, the magnitude and sign of the multiparameter coefficients determine a net bathochromic shift of the absorbance peak maximum, with a relative contribution of 67% from the coefficient -di, 17% from -b 1 , 8% from -e and 8% from a 1 .
Thus, non-specific interactions are evidenced as the most influential between the API and its solvation medium, given by polarizability and dispersion interactions and electronic induction of the FNZ dihalogenated phenyl ring and the solvent basicity (HBA capacity) given by the hydroxyl group of the molecule.The solvatochromic analysis confirms that when the polarity of the solvent increases, the excited state stabilizes more than the basal state of FNZ, being the first one of higher polarity and lower energy, typical of the bathochromism observed. 13,27,41,44,52he International Union of Pure and Applied Chemistry (IUPAC) has defined the concept of pharmacophore or pharmacophore pattern as "the ensemble of steric and electronic features that is necessary to ensure the optimal supramolecular interactions with a specific biological target structure and to trigger (or to block) its biological response". 54 pharmacophore does not represent a real molecule or a real association of functional groups, but rather a purely abstract concept that accounts for the common molecular interaction capacities of a group of compounds towards their target structure.The pharmacophore can be considered as the largest common denominator shared by a set of active molecules.This definition discards a common misuse in medicinal chemistry that consists of naming as simple chemical functionalities such as guanidines, sulfonamides or dihydroimidazoles (formerly imidazolines), or typical structural skeletons such as flavones, phenothiazines, prostaglandins or steroids as pharmacophores). 55Ji et al. 56 have studied the pharmacophoric conformations of 14 azole antifungals, which have a similar docking mode at the active site of the enzyme lanosterol 14-α-demethylase from Candida albicans.
The halogenated phenyl group of azole inhibitors is deep in the same hydrophobic binding cleft as the 17-alkyl chain of the substrate.Interactions of π-π stacking could exist between the halogenated phenyl ring and the aromatic ring of residue Y132 in the hydrophobic cavity of the active site, which increases the fungicidal activity of the API. 56This information is consistent with the large relative contribution of the SP and DI parameters of the Catalán and Laurence equations and is in agreement with the experimental bathochromic behavior of FNZ in solution.
It has also been suggested that the oxygen atom attached to C10 is favorable to antifungal activity for two reasons.First, the chirality of that carbon atom gives it greater antifungal activity than other stereoisomers, for which FNZ led a new generation of more powerful, better tolerated, and metabolically more stable antifungals.Second, the oxygen atom interacts with the water molecules conserved in the active site forming hydrogen bonds to bind it to the H310 residue of the protein. 56ecently, crystallographic data of the lanosterol 14-α -demethylase-fluconazole complex from Saccharomyces cerevisiae provided information on a watermediated hydrogen bond network between the API and the enzyme.Two water molecules form networks of hydrogen bonds with FNZ.The first water molecule mediates the hydrogen bonds between the hydroxyl groups of FNZ and Y140, as well as a propionate of the heme cofactor.The second forms hydrogen bonds with the carbonyl oxygen of S382, the hydroxyl of Y126, and the atom N4' of the second triazole ring furthest from heme. 57The previous evidence agrees with the relative contribution of the SB and parameters of the Catalán and Laurence equations, with the same bathochromic effect on the solvation of FNZ.On the other hand, the relative contribution of the SA and α 1 parameters of the Catalán and Laurence equations in the inhibition function of the N4' in the second triazole ring is evident, thus indicating the HBD solvent capacity, although small and hypsochromic.Vol.34, No. 10, 2023   Theoretical solvatochromic study Table 5 summarizes the spectral behavior of both optimized conformers of FNZ in all the solvents studied.They are ordered by increasing polarity.The shifts of the absorption band in the UV-Vis in each solvent are similar to the experimental and theoretical data reported in the literature. 19All the theoretical λ in the UV-Vis region are lower than the experimental ones.The λ obtained with conformer II are lower than those of conformer I, and the λ obtained by CAMB3LYP are lower than those calculated by B3LYP.Every one of the values achieved corresponds to π-π* transitions, due to HOMO-LUMO electronic excitations.Taking the solvent acetonitrile as a reference, we observe that FNZ conformer I (λ B3LYP = 235.7 nm) shows on average hypsochromic shifts of 0.7 nm from nonpolar solvents such as cyclohexane, tetrachloromethane, and benzene.When FNZ moves from aprotic solvents of lower polarity, such as diethyl ether and dichloromethane to acetonitrile, hypsochromic shifts of 0.3 and 0.2 nm are produced, respectively, while towards DMSO no shift was observed.Similarly, a bathochromic shift of 0.1 nm is obtained when moving to the protic polar solvent ethanol, while there is no change for the protic solvents methanol or water.Similar behaviors are displayed in all the theoretical calculations described.
Tables 6 and 7 show the coefficients and errors associated with the parameters of the solvatochromism analysis of both FNZ conformers calculated by equations 1, 2, and 3. Also, the coefficient of determination, the relative contribution of each parameter, and the p-value of the Fisher statistic for the linear regression are reported.At first, the analysis is performed only by including polar solvents to compare them with the experimental results already obtained.Next, non-polar solvents data are included to consider the zero or nearly zero contributions of their solvatochromic parameters.
With respect to the theoretical solvatochromic analysis obtained by both the B3LYP and CAMB3LYP methods and evaluated by Kamlet and Taft with n = 7, it is observed that only the parameters π* and α are significant, with a similar percentage for both conformers.However, considering n = 10 for both methods, the significant parameters are π* and β, with a relative contribution of 44% for s and 42% for b, with 78% of the variability for conformer I by CAMB3LYP, displaying greater statistical significance.All coefficients have a positive sign, corresponding to a hypsochromic shift.
The theoretical multiparametric analysis of Catalán with n = 7 by the B3LYP method showed significance in all parameters in both conformers, with similar results.The relative contribution of the parameters corresponding to non-specific interactions changed drastically with respect to the experimental analysis, being almost the same in the theoretical study: 50% of d SdP and 39% of -c SP .The relative contribution of specific interactions is similar, around 6% in both a SA and b SB .By raising the number of solvents to n = 10, the SA contribution loses significance, increasing the statistical significance of the whole linear regression by two orders of magnitude.With the CAMB3LYP method, the results are like those obtained for both conformers, reducing the representativeness of the specific interactions in conformer II with n = 10 and minimizing the acid character of the solvent in conformer I.Although this  method better explains the variability ofv in the solvation environment, with 98%, the highest statistical significance is obtained with the B3LYP method.Finally, the theoretical analysis of the results obtained with the B3LYP method and carried out by applying the Laurence equation with n = 7 resulted in more than 99% of the variability for both conformers.The relative contribution of the non-specific interactions varied slightly, compared with Catalán equation.The electronic dispersion and induction were above the interactions between the permanent dipoles of the API and the solvent.As regards specific interactions, the solvent acidity was almost irrelevant.
With n = 10, the difference between the proportion of DI and ES interactions increased, where -di is greater than 60% and e is around 30%, without relevance of the basicity of the solvent and a greater statistical significance in the linear regression analysis.Using the CAMB3LYP method, similar results were obtained as with B3LYP, with a representativeness of more than 99% of the response variable for both conformers in the case of n = 7, and the highest statistical significance with n = 10.In the latter case, both conformers have a high relative contribution of the DI and ES parameters, as well as the ability to accept hydrogen bonds, α 1 , and without relevance of the capacity to give hydrogen bonds, β 1 .The net spectroscopic shift is bathochromic due to the magnitudes and signs of the coefficients: -di, e, -a 1 , and b 1 .It should be noted that although ES has a smaller magnitude than DI with respect to the type of interaction it represents, its statistical error is much lower, and its influence on the API solvation is more reliable.In short, the interactions of electronic polarizability, dispersion, and induction are equally important in theoretical and experimental solvatochromic behavior, and according to which equation is analyzed, Catalán or Laurence, and which method, B3LYP or CAMB3LYP, it takes on a certain relevance the acidity or basicity of the solvent.Both, B3LYP and CAMB3LYP methods, present the same percentage of explanation of the dependent variable with similar statistical significance.However, the theoretical results obtained with the B3LYP method are closer to the experimental values of λ max than with the CAMB3LYP method.
Figures S1 and S2 (see Supplementary Information section) show theoretical wavelength values for both FNZ conformers.The main difference between the two structures is that π-π* electronic transitions in polar solvent for FNZ conformer II are stronger than the ones in conformer I for the similar kind of solvents.The theoretical and experimental UV-Vis spectra were not depicted together in a single figure due to the differences in wavelength between the obtained data with the basis sets used for the analysis.makes evident the existence of three clearly identifiable isosbestic points at 262.4, 265.0, and 267.6 nm, without crossing the spectra in other positions.This correlates with intermolecular equilibria of the FNZ in the buffer, but not with intramolecular equilibria.
Three acid-base constants have been described for FNZ, according to its ability to give up H + , given by the -OH group of C10 or to accept H + , by any of the -N of the triazole groups 1,2,4-triazole. 49They are congruent with the three isosbestic points found in the API thermochromic study, where three intermolecular acid-base equilibria are proposed: pK a HA = 11.93,pK ab1 BH + = 5.23 and pK ab2 BH 2+ = 1.76 (Figure 6).
The stability of the molecule and charge delocalization of FNZ were studied by Chandrasekaran and Thilak Kumar, 19 who analyzed the number of electrons or natural population, the natural bonding orbitals, and the natural local molecular orbitals by molecular modeling calculations with the B3LYP/6-311++G(d,p) method.These results showed that the charge of the lone pairs of the N4' of each 1,2,4-triazole group is slightly higher than that of the other nitrogen atoms in the molecule.The N4' atom becomes more electronegative, since the quantity of charges transferred from the n (N4') orbital is very low, which makes the triazole ring a highly chelating agent, susceptible to protonation or hydrogen bonding.Therefore, the acid-base equilibria that are related to the isosbestic points probably would correspond to the loss of H + from the C10 hydroxyl group, and to the gain of one H + in each of the N4' of the azole rings, as illustrated in Figure 6 (boxed intermolecular equilibria).
The importance of the N4' of the azole ring as a chelating agent is highlighted by considering it as part of the azole pharmacophore, rather than N1' or N2', since the N4' of the azole ring is covalently bound to the heme of the target substrate.The compound would be deprived of activity if N4' is replaced by a carbon atom. 19More recent crystallographic studies of the binding of the enzyme to FNZ confirm the importance of the N4' of the triazoles by demonstrating coordination bonding of the nitrogen atom N4' in the triazole ring with the heme iron, replacing the water ligand and stabilizing the low spin bond.The triazole ring of FNZ coordinates with iron within the heme cofactor (Fe---N distance 2.13 Å), corroborated by spectrophotometric data. 55As already mentioned, the N4' of the other triazole participates in the network of hydrogen bonds with the residues of the active site outside the heme group.

Conclusions
The theoretical study of FNZ using the bases DFT-B3LYP/6-31+G(d) and DFT-CAMB3LYP/6-31+G(d) allows the geometric optimization of two structural conformers of the API corresponding to two true minima of energy potential.Both conformers show high electronic delocalization in the semicarbazoles, intramolecular charge transfer from the electron-donor substituents to the phenyl ring, the possibility of forming numerous intramolecular hydrogen bonds, and the same probability of crystallizing in different polymorphs when forming crystal lattices.The analysis of the molecular descriptors allows us to appreciate that both structures present great chemical stability and can be classified as marginal according to their electronegativity.Their dipolar moments rise in the solvated  phase as the polarity of the solvent increases, which indicates the high inductive effect of the polarizability of the polar solvent that affects the FNZ molecule.
The experimental solvatochromic study of FNZ results in a slightly bathochromic shift as the solvent polarity increases, given by π-π* electronic transitions.The multiparametric equations of Catalán and Laurence were more efficient than Kamlet and Taft's equation.The magnitudes and signs of the coefficients of the solvatochromic parameters of Catalán and Laurence agree in their relative relevance, around 65% of the SP parameter, solvent polarizability, and DI, electronic dispersion and induction, both related to π-π stacking type interactions between the dihalogenated phenyl ring of FNZ and the aromatic ring of residue Y132 in the hydrophobic cavity of the active site.Similarly, with a negative sign and a relative relevance of 20%, the parameters SB and β 1 become important, both indicative of the solvent basicity (HBA ability), consistent with the network of hydrogen bonds mediated by water between the API and the enzyme, specifically the hydroxyl group of FNZ and the Y140 residue.To a lesser extent, and with a positive sign, 8% of the SA and α 1 parameters are observed, related to the solvent acidity (HBD capacity).Even though it is a low value, it shows the importance of the N4' in the furthest triazole beyond the heme group within the water-mediated hydrogen bonding network.
The theoretical study of the solvatochromism of FNZ includes non-polar solvents in the multiparameter equations and confirms the bathochromic shift of the maximum of the spectral line as the polarity of the solvent increases.However, the relative relevance of the solvatochromic parameters changes.On average, 90% of the relevant interactions correspond to non-specific interactions, with polarizability/dispersion and electronic induction prevailing over dipolarity/interaction of permanent dipoles between solute and solvent.As the predominant interaction forces have opposite signs, a mixed solvatochromic resulting effect is observed, combined with the negative coefficients of the specific interaction parameters.Regarding the TD-DFT methods used, both B3LYP and CAMB3LYP result in shorter wavelengths than experimentally, matching the available bibliography.In our case, the theoretical results obtained with the B3LYP method are closer to the experimental values of λ max than with the CAMB3LYP method.
The thermochromic study shows intermolecular acidbase balances consistent with the three pK a of the FNZ, showing the deprotonation of the R-OH group and the protonation of each of the N4' of the azole rings, the most electronegative of each triazole.Similarly, the importance of N4' is evident by the coordination bond of the N4' atom in the triazole ring with heme iron, which is essential for the manifestation of antifungal activity.Both nonspecific interactions and hydrogen bonding of FNZ correlate with substrate binding and antifungal activity, which is inherent to the proposed pharmacophore for azoles.The suggested ionization equilibria at physiological pH justify a very good oral and intravenous bioavailability and effective tissue distribution of the API.
Besides the former pharmacological action fully analyzed against several kinds of fungi, the physicochemical properties of FNZ studied in this research work complement the available scientific data.In the future, this information could lead to reposit this drug as an API for new therapeutic targets.

Figure 1 .
Figure 1.Chemical structure of the antifungal FNZ.

Figure 2 .
Figure 2. Chemical structure of the conformer I (a) and II (b) of the antifungal FNZ.

Figure 4 .
Figure 4. Maximum wavelength of FNZ in function of increasing solvent viscosity.

R 2 :
coefficient of determination; F: p-value of Fisher's statistic.A 0 : property of the solute to be correlated.Correlation coefficients of polarity-polarizability of the solvent (s), acidity (a) and basicity (b) of the solvent in Kamlet and Taft analysis.Correlation coefficients of solvent polarizability (c SP ), dipolarity (d SdP ), solvent acidity (a SA ), solvent basicity (b SB ) in Catalán analysis.Regression coefficients of dispersion and electronic induction interactions (di), electrostatic interactions (e), acidity (a 1 ) and basicity (b 1 ) of the solvent in Laurence analysis.

aB3LYP:
functional Becke, 3-parameter, Lee-Yang-Parr; b CAM B3LYP: hybrid exchange correlation Coulomb Attenuated Method -Becke, 3-parameter, Lee-Yang-Parr; c n = 7; d n = 10; e coefficient of determination; f p-value of Fisher's statistic.A 0 : property of the solute to be correlated.Correlation coefficients of polarity-polarizability of the solvent (s), acidity (a) and basicity (b) of the solvent in Kamlet and Taft analysis.Correlation coefficients of solvent polarizability (c SP ), dipolarity (d SdP ), solvent acidity (a SA ), solvent basicity (b SB ) in Catalán analysis.Regression coefficients of dispersion and electronic induction interactions (di), electrostatic interactions (e), acidity (a 1 ) and basicity (b 1 ) of the solvent in Laurence analysis.

Figure 5
Figure 5 reveals the variation of the absorbance spectrum of FNZ in buffer solution at pH 7.4 due to the temperature gradient to which the solution was exposed.The thermal gradient plotted increases each line in 20 ºC from 5 to 65 ºC, although the measurements were performed every 5 ºC.The analysis of the thermochromic study

aB3LYP:
functional Becke, 3-parameter, Lee-Yang-Parr; b CAM B3LYP: hybrid exchange correlation Coulomb Attenuated Method -Becke, 3-parameter, Lee-Yang-Parr; c n = 7; d n = 10; e coefficient of determination; f p-value of Fisher's statistic.A 0 : property of the solute to be correlated.Correlation coefficients of polarity-polarizability of the solvent (s), acidity (a) and basicity (b) of the solvent in Kamlet and Taft analysis.Correlation coefficients of solvent polarizability (c SP ), dipolarity (d SdP ), solvent acidity (a SA ), solvent basicity (b SB ) in Catalán analysis.Regression coefficients of dispersion and electronic induction interactions (di), electrostatic interactions (e), acidity (a 1 ) and basicity (b 1 ) of the solvent in Laurence analysis.

Figure 6 .
Figure 6.The acid-base equilibria of the antifungal FNZ.

Table 1 .
Physical constants and empirical solvatochromic parameters of solvents

Table 4 .
Estimated coefficients and standard errors, correlations, significance, and the relative contributions percentage of the parameter's coefficients for multiple linear regression analysis of FNZ experimental data

Table 5 .
FNZ theoretical values of wavelength and wavenumber with f at the TD-DFT/6-31+G(d) level of theory

Table 6 .
Estimated coefficients and standard errors, correlations, significance, and the relative contributions percentage of the parameter's coefficients for multiple linear regression analysis of FNZ conformer I theoretical data

Table 7 .
Estimated coefficients and standard errors, correlations, significance, and the relative contributions percentage of the parameter's coefficients for multiple linear regression analysis of FNZ conformer II theoretical data