In silico and in vitro Studies Evaluating the Promising Antiureolytic Activity of Schiff’s Base 4-(3-Hydroxybenzylideneamino)phenol and Its Amine Derivative

Urea is the most widely used nitrogen fertilizer worldwide. However, ammonia volatilization, resulting from applying urea to the soil surface, causes economic and environmental losses; thus, urease inhibitors have been developed to mitigate these losses. In this work, the anti-ureolytic activity of Schiff’s base 4-(3-hydroxybenzylideneamino) phenol ( 3B4 ) and its amine-derived ( 3B4a ) was evaluated. The most promising urease inhibitor in soil was 3B4 (55.0 ± 3.9% inhibition), with comparable results to N -(butyl) thiophosphoric triamide ( p = 0.659). In the in vitro analysis ( Canavalia ensiformis ), the results of anti-ureolytic activity were similar, 22.6 ± 6.9% for 3B4 and 24.2 ± 9.6% for 3B4a. Biophysical interaction studies were also carried out through molecular docking studies and molecular fluorescence spectroscopy. These studies showed that both substances are preferentially competitive inhibitors, with the interaction between 3B4a and urease forming a more stable complex. In the analysis by Fourier transform infrared spectroscopy, no interaction was observed when 3B4 or 3B4a was mixed with urea (1:1) for 48 h, providing evidence of compatibility. Thus, the Schiff base 3B4 and its corresponding amine 3B4a may represent potential additives for urea fertilization aiming to assist in the urease inhibition process.


Introduction
Nitrogen (N) is essential to produce macromolecules and small key molecules by plants, such as chlorophyll.Although it is abundant in the atmosphere, representing about 80% of its volume, N is biologically unavailable due to its highly stable gas form. 1,2But, through the Haber-Bosch process, it is possible to obtain N fertilizers, contributing to high productivity in current crops, 1,2 as 30 to 50% of the agricultural production is attributed to introducing mineral fertilizers, including synthetic N. Thus, feeding about 48% of the world's population is due to using N fertilizers, 1,3 highlighting the importance of this element in agriculture to maintain the food supply.Urea is the most widely used synthetic nitrogen fertilizer, accounting for more than 70% of fertilizer use worldwide. 46][7] However, urea can undergo hydrolysis when applied to the soil, and this process is associated with the urease enzyme. 8,9he urease enzymes occur widely in nature, as in soils, playing a decisive role in N metabolism. 6This enzyme was crystallized for the first time by James B. Sumner, who received, in 1946, the Nobel Prize in Chemistry for showing that enzymes can be isolated in crystal form. 10,11reases were the first example of a nickel metalloenzyme cited in the literature. 6,12This enzyme catalyzes the hydrolysis of urea in ammonia and carbamate, which spontaneously decomposes into carbonic acid and more ammonia. 13,146][17][18][19] The use of urea can also contribute directly (nitrification and denitrification) and indirectly (NH 3 volatilization) to N 2 O emissions. 20,21This gas contributes to the greenhouse effect and destruction of stratospheric ozone. 1,9In addition, higher ammonia emissions can favor aerosol formations, decreasing visibility, 22 and favoring the formation of particles smaller than 2.5 µm in diameter (PM 2.5 ), 23 which negatively impacts human health. 24he degradation of urea by urease also generates economic losses for farmers since fewer nutrients are available to plants, affecting productivity. 5Pinto et al. 25 estimated that approximately 5.26% of the resource invested in applying 400 kg urea per ha per year in an area of 100 ha was lost due to ammonia volatilization and denitrification.Thus, the development of alternatives to improve the efficiency of N fertilizers is desirable from both an environmental and an economic point of view. 268][29] The N-(butyl) thiophosphoric triamide (NBPT) is the most commercially used urease inhibitor in agriculture, being able to delay urea hydrolysis and decrease ammonia volatilization. 5,21,29,30However, studies have shown that the stability of NBPT is a concern since its degradation occurs over time when applied to urea and is dependent on temperature and pH. 5,8,31,32Thus, it is necessary to increase the supply of urease inhibitors that are effective in different types of soils, compatible with urea, efficient in low concentrations, and safe for the environment and human health. 26mong the new substances synthesized by our research group, Schiff bases have stood out for their ease of synthesis, low cost, and higher thermal stability than NBPT, as reported in a patent document. 33Among 71 synthesized substances, eight stood out (Figure 1), and substance 3-hydroxybenzaldehyde (Scheme 1, 3B4) was chosen because it presented promising results in urease inhibition experiments in vitro and soil, in addition to not harming germination and plant development. 33,34n this work, an amine analog, 4-(3-hydroxybenzylideneamino) phenol (Scheme 1, 3B4a), derived from Schiff base 3B4, was compared regarding its anti-ureolytic activity in soil and in vitro.Molecular docking and molecular fluorescence spectroscopy studies were performed to understand better the structural features responsible for the anti-ureolytic activity and the interaction of the assets with the urease active site.Finally, the interaction of 3B4 and 3B4a with urea was also evaluated, aiming at the application of these substances as additives to increase the efficiency of fertilizers.

Chemicals
All chemicals were obtained from commercial suppliers and used without further purification.Urease (from jack beans, nominal activity 40.1 units mg -1 solid) was obtained from Sigma-Aldrich (St. Louis, USA).Uncorrected melting points were determined in a melting point meter PF 1500 FARMA (Gehaka, São Paulo, Brazil).Reactions under microwave irradiation (MW) were carried out in a CEM-Microwave-Enhanced Life reactor (Matthews, USA).The 1 H nuclear magnetic resonance (NMR) and 13 C NMR spectra were recorded on a Bruker Avance DRX/400 (Bruker, Billerica, USA).Chemical shift values (d) were given in parts per million (ppm).Fourier transform infrared spectroscopy (FTIR) spectra were recorded in a Spectrometer One PerkinElmer (PerkinElmer, Waltham, USA).
A microplate spectrophotometer (Multiskan FC, Thermo Scientific, Waltham, USA), a centrifuged model Heraeus Fresco 17 Centrifuge (Thermo Scientific, Waltham, USA), and a Shaker MaxQ4000 (Thermo Scientific, Waltham, USA) were used to study the enzymatic activity of urease inhibition.

Synthesis of 3B4 and 3B4a
Synthesis of 3B4 (Scheme 1) was performed using the 3-hydroxybenzaldehyde (10 mmol) (Sigma-Aldrich, St. Louis, USA) and 4-hydroxyaniline (10 mmol) (Sigma-Aldrich, St. Louis, USA) in ethanol.The mixture was placed in a microwave reactor and irradiated at reflux temperature, with a ramp time of 4 min and reaction time of 10 to 20 min under maximum agitation.Purification was performed by recrystallization using hexane:ethyl acetate (1:1).The product was characterized by melting point analysis, NMR ( 1 H and 13 C), and FTIR spectroscopy.
The 3B4 amine derivative (Scheme 1), compound 3B4a, was prepared according to Byung and Sang, 35 with modifications.Briefly, 5 mmol of 3B4 were solubilized in ethanol, followed by a slow addition of sodium borohydride (NaBH 4 ; in a 1:2 ratio).The mixture was placed under a microwave irradiation reactor for 10 min and then subjected to extraction with ethyl acetate and distilled water.The organic phase was collected, dried with sodium sulfate (Na 2 SO 4 ), and evaporated after filtration.The product was characterized by melting point analysis, NMR ( 1 H and 13 C), and IR spectroscopy.
The soil was sieved using 2 mm mesh, and samples of 0.5 g, previously activated with water for 24 h, were placed in 15 mL tubes.To evaluate urease inhibition in soil, 500 μM of the test substances were solubilized in ethanol and 80 mM of urea (Merck, Darmstadt, Germany).Samples were incubated in the Shaker (Thermo Scientific, MaxQ4000, Waltham, USA) at 37 °C for 1 h.Reactions were stopped by adding 5 mL of 1 M potassium chloride (KCl) in 1 M hydrochloric acid (HCl).Each system was incubated for a further 30 min and then centrifuged (Thermo Scientific, Heraeus Fresco 17 Centrifuge, Waltham, USA) for 3 min at 10,000 × g.After centrifugation, 100 μL of a solution containing 3.4% (m/v) sodium salicylate, 2.5% (m/v) sodium citrate, 2.5% (m/v) sodium tartrate, and 0.012% (m/v) sodium nitroprusside (SNP) were added to the 20 μL of the supernatant, which was incubated for 15 min at 25 °C and 600 rpm.Each mixture received 0.1 vol of a solution containing 3.0% (m/v) sodium hydroxide (NaOH) in 1% (v/v) sodium hypochlorite, followed by incubation for 1 h at 25 °C in the dark under 600 rpm.The ammonium formed was quantified (λ max = 660 nm) (Multiskan FC, Thermo Scientific, Waltham, USA) with the aid of a calibration curve using ammonium chloride as a standard.All experiments were performed in triplicate (n = 3) with five technical replicates for each substance.NBPT was used as a reference for urease inhibition.

In vitro inhibitory activity against urease (Cavanalia ensiformis)
In vitro urease inhibition studies were conducted to understand the interactions of synthesized substances (3B4 and 3B4a) with urease presented in theoretical molecular analyses.Reactions were performed in 20 mM phosphate buffer (pH 7.0) containing 1.0 mM EDTA (ethylenediamine tetraacetic acid), 12.5 mU urease of C. ensiformis enzyme (Sigma-Aldrich, St. Louis, USA), 10 mM urea in the presence or absence of each compound-test at 50 μM (solubilized in ethanol). 39The reactions were incubated at 25 °C for 10 min and stopped by adding 0.5 vol of 1% m/v phenol in 5 mg L -1 SNP and 0.7 vol of 0.5% m/v NaOH in 0.1% v/v sodium hypochlorite.The systems were incubated at 50 °C for more 5 min, followed by quantification of NH 4 + (λ max = 630 nm) (Multiskan FC, Thermo Scientific, Waltham, USA) with the aid of an analytical curve using ammonium chloride as a standard.Hydroxyurea (HU) (Sigma-Aldrich, St. Louis, USA) was used as the inhibitor reference.The experiments were performed in quadruplicate (n = 4) with four repetitions for each substance.In experiments evaluating anti-ureolytic activity, the online software Chemicalize, developed by ChemAxon, 40 was used to predict the pK a and logP of substances 3B4 and 3B4a.

Molecular modeling
The docking study was performed with compounds 3B4 and 3B4a to predict binding conformation, interactions with essential residues, and coordination with the bi-nickel center in the urease's active site.The lowest energy binding poses for derivatives were chosen as the initial conformations for the molecular dynamics (MD) simulations of the proteinligand complexes for 100 ns to evaluate the complexes' stabilities.Root-mean-square deviation (RMSD) of ligands and root-mean-square fluctuation (RMSF) of protein were calculated based on the analysis of the MD trajectories.Finally, MM-PBSA calculations were performed using the trajectory files obtained after the MD simulations to determine the protein-ligand complexes' binding energy.All simulations were performed according to our research group's previous procedures. 41

Interaction between inhibitors-urease by molecular fluorescence spectroscopy
The spectrofluorimetric measurements were performed using the Varian Cary Eclipse spectrofluorimeter (Agilent, Santa Clara, USA) with a thermostat provided by a single cell peltier, with 1 cm optical path quartz cuvettes.Substances 3B4, 3B4a, NBPT, HU, and urea were used in interaction studies with the enzyme urease (C.ensiformis) (Sigma, St. Louis, USA).The stock solutions of the substances were prepared in ethanol, and the urease solutions alone or containing the inhibitors were prepared in 20 mM phosphate buffer (pH = 7.40) and ultrapure water (Milli-Q ® purification system, Merck Millipore, Burlington, USA).
Spectrophotometric measurements were performed on a Hitachi U-2010 UV-Vis spectrophotometer (Hitachi, Tokyo, Japan).Based on the absorbance signal (λ max = 280 nm) and the molar extinction coefficient for urease enzyme (54780 M -1 cm -1 ), 42 the urease concentration in the stock solution was calculated.
The spectrofluorimetric interaction studies were performed using 1.0 µM urease, λ ex = 280 nm, λ em = 290 to 420 nm, and the slit was 5/10 nm (λ ex /λ em ).Thermodynamic parameters and intermolecular forces involved in the interaction of enzyme and inhibitors were determined based on spectrofluorimetric titrations.The urease concentration was fixed at 1.0 μM, and increasing inhibitor concentrations (5-80 μM) at 23, 30, and 38 °C.The reference signal was established considering the system without the inhibitors under the same conditions.The quenching mechanism was assessed using the Stern-Volmer equation 1: (1) where F 0 and F are the fluorescence intensities in the absence and presence of the ligand, respectively; [Q] is the concentration of the ligand; K SV is the Stern-Volmer quenching constant; k q is the maximum rate constant for diffusional quenching bimolecular constant in biopolymer systems (2.0 × 10 10 M -1 s -1 ); and τ 0 is the average lifetime (10 -8 s).To evaluate the degree of interaction between the compounds through the binding constant (K b ) and to determine the stoichiometry of the complex formed between urease and the inhibitors (n) were used equation 2: (2) where, F 0 and F are the fluorescence intensities in the absence and presence of the ligand, respectively; [Q] is the concentration of the ligand; n is the stoichiometric proportion of the supramolecular complex formed between enzyme and inhibitor (number of sites occupied by the ligand in the urease structure); and K b is the binding constant (expresses the strength of the interaction between the ligand and the macromolecule).The thermodynamic parameters involved in binding were also determined.The values of enthalpy (ΔH) and entropy (ΔS) were calculated by linearizing the Van't Hoff equation 3: (3) R is the gas constant (8.3144621J mol -1 K -1 ); T is the experimental temperature (Kelvin); and is the binding constant at a given temperature.To evaluate the process spontaneity, Gibbs free energy (ΔG) was calculated by equation 4: The 3D fluorescence was used to assess conformational changes in the enzyme structure in the presence of the inhibitors.The spectrofluorimetric measurements were performed under the following conditions: 1.0 μM enzyme urease and 3B4 or 3B4a at 70 μM.The instrumental conditions were the same as in the previous experiments.The interaction of the substrate (urea), the standard inhibitors (HU and NBPT), and the proposed inhibitors (3B4 and 3B4a) with the enzyme was evaluated through a competition assay.Spectrofluorimetric titrations of the enzyme (1 μM) with 3B4 and 3B4a (5-80 μM) were performed in the presence and absence of inhibitors or substrate (10 μM).

Interaction between 3B4 and 3B4a with urea
The FTIR technique was used to verify the interaction between the proposed inhibitors (3B4 and 3B4a) and urea.FTIR spectra were recorded using Spectrum One PerkinElmer equipment with an attenuated total reflectance (ATR) system (PerkinElmer, Waltham, USA), and absorption was measured in wavenumber (cm -1 ).The analysis was performed with the mixture (1:1) of urea and 3B4 or 3B4a, crushed with a pestle in an agate mortar.The samples (protected from light) were prepared and analyzed after 0, 24, and 48 h.

Statistical analysis
The results were presented as means ± standard deviations.The Kolmogorov-Smirnov test assessed normality, while homoscedasticity was by Brown-Forsythe.If the data followed a normal distribution and homoscedastic, parametric statistics were used (one-way analysis of variance (ANOVA) test followed by Tukey post-test).Outlier analysis was performed with a 10% Q ROUT test.The s statistical difference was considered significant when p-values < 0.05.Statistical analyzes were performed using GraphPad Prism version 6.01. 43

Synthesis
The imine was obtained from a condensation reaction between 3-hydroxybenzaldehyde and 4-hydroxyaniline with a high yield (72%) and easy preparation.The respective amine is the product of reductive amination of 3B4, thus guaranteeing a high yield (76%).Both urease inhibitors were characterized by melting point, FTIR, and 1 H and ¹³C NMR (spectra of both substances can be accessed in the Supplementary Information (SI) section), and the values were consistent with the literature.

Evaluation of urease inhibition in soil and in vitro
Among the urease inhibitors tested, 3B4 was as efficient as NBPT and inhibited soil ureases by 55.0% on average (Figure 2a, p = 0.659).Notably, the imine 3B4 was roughly twice as much potent as the corresponding amine 3B4a (Figure 2a, p = 0.006).The most potent substance in soil, 3B4, and the corresponding amine were submitted to in vitro test with urease from C. ensiformis.The antiureolytic activity of the substances 3B4 and 3B4a at 50 µM (22.6 ± 6.9 and 24.2 ± 9.6%, respectively) were similar to that of the reference inhibitor HU (Figure 2b, p = 0.988 and p = 0.991, respectively).
The difference observed between the activity of 3B4 and 3B4a, when comparing the in vitro and soil results, suggests that the 3B4a substance may be protonated at the soil pH, becoming more retained in soil colloids, and less available to act as a urease inhibitor.Colloids are mainly responsible for the chemical activity of the soil and have a net negative charge that attracts and retains particles/substances with a positive charge. 44Using the online software Chemicalize 40 it was possible to obtain the pK a values of the substances 3B4 and 3B4a and, from the pK a values, calculate the distribution of micro species.As can be seen in Figure 3, at the pH of the soil used in the test (5.8),around 98.7% of the substance 3B4 is in its neutral form (Figure 3a, H 2 A), compared to 55.8% of the 3B4a (Figure 3b, H 2 B).Thus, the protonation of 3B4a (43%, Figure 3b, H 3 B + ) favors its retention in the soil, which may justify the lower activity (2 times lower) in relation to the 3B4 imine (Figure 2).
In the in vitro assay, the influence of the soil matrix is inexistent, showing only the activity of 3B4a in contact with the enzyme urease.Thus, molecular modeling and spectroscopy biophysical studies were proposed to understand which groups are the most important for 3B4 and 3B4a to act as urease inhibitors and how the interaction between these inhibitors and the enzyme works.

Molecular modeling
The docking simulations revealed that the position of the hydroxyl group in the aromatic ring (3-and 4-position for 3B4 and 3B4a, respectively) is a critical factor for enzyme inhibition. 45For amine derivative 3B4a, the binding conformation showed 4-hydroxyl coordination with the urease's bi-nickel center (Figure 4b).Further stability was due to the formation of some interactions, including conventional hydrogen bonds between the -OH group with carbamylated lysine KCX490 and His409 residues.The analysis of the imine derivative 3B4 revealed traditional hydrogen bonds between the -OH group with His519 and Ala440 residues, and coordination with only one Ni II center, which suggests the formation of a less stable complex, and slightly lower activity.Finally, both compounds performed π-alkyl interactions with Leu589 and π-cation with His585, as well as hydrogen bonds between CME592 and imine or amine groups of 3B4 and 3B4a, respectively.
These results are in line with the results obtained in in vitro tests.When evaluated in vitro, the substance  3B4a showed slightly greater antiureolytic activity than 3B4, since it was completely bioavailable, favoring the formation of a stable complex with urease.It is important to highlight that, although both substances have a hydrophobic profile (solubility studies, Table S1, SI section), in the study of the antiureolytic effect the substances were completely solubilized in a hydroalcoholic solution.Thus, in all studies, substances 3B4 and 3B4a were bioavailable.However, the substance 3B4a interacts with the colloids present in the soil, as it is protonated at the soil pH, which reduces its bioavailability and activity.
To evaluate the 3B4-and 3B4a-urease supramolecular complexes' stabilities at the urease catalytic site, MD simulations were performed using the docking poses as initial conformation.For both derivatives, the system acquires stability after 20 ns of simulation and remains until the end of the simulations.The RMSD values (Figure 5a) from their docking poses showed minimal deviations (between 0.1 and 0.3 nm), indicating that the molecular docking results are reliable.In addition, the RMSF plots (Figure 5b) show low fluctuations (between 0.1 and 0.5 nm), with emphasis on smaller changes for the residues that coordinate with nickel (His 407 , His 409 , KCX 490 , His 492 , His 519 , and His 545 ).Finally, MM-PBSA calculations were applied to determine the binding energy of the complexes with the urease.The results showed that the binding energy value for derivative 3B4a (ΔG bind = -46.65 ± 0.43 kJ mol -1 ) is about two times lower than the binding energy value for 3B4 (ΔG bind = -20.38 ± 0.83 kJ mol -1 ), which agrees with binding constants (K b ) which were experimentally calculated for both compounds.
Binding and thermodynamic parameters of macromoleculeligand   Figure S7, SI section).In Figure 6 it is also possible to observe the Stern-Volmer quenching plot and double logarithmic curves for calculating the binding constant of substances (3B4 and 3B4a).As can be seen, urease showed an intense fluorescence band between 334 nm when excited at 280 nm (Figures 6a and 6d).It is also possible to observe that increasing concentrations of the inhibitors (3B4 and 3B4a) led to a reduction in the fluorescence intensity (quenching process) associated with a bathochromic displacement (Figures 6a and 6d).This spectral profile may be related to conformational alterations of the enzyme, resulting from the interaction with inhibitors, and changes in the microenvironment of Trp and Tyr residues, causing an increase in polarity (water presence) around these residues, 46 and will be covered in the 3D fluorescence study.Different interactions at the molecular level can result in quenching, such as excited-state reactions, molecular rearrangement, transfer of energy, ground-state complex formation, and collisional quenching.Consequently, changes or perturbations in fluorescence properties like quantum yield, intensity, and/or lifetime are observed. 47t is possible to determine the quenching process by evaluating the constant k q , whose values are illustrated in Table 1.When k q is less than 2.0 × 10 10 M −1 s −1 , quenching tends to be predominantly dynamic, while higher values of k q indicate that quenching will be mostly static. 48From these results, it is possible to infer that the suppression of the intrinsic fluorescence of the enzyme by the two inhibitors (3B4 and 3B4a) occurs through static quenching.Static quenching occurs in the ground state, resulting in the formation of a non-fluorescence-emitting complex between the fluorophore and the quenching agent. 47The k q range for substance 3B4 was between 0.51 and 0.85 × 10 12 M −1 s −1 and for substance 3B4a from 0.84 to 1.16 × 10 12 M −1 s −1 , exceeding the limiting diffusional constant, showing the predominance of static quenching.
The degree of interaction between the compounds was evaluated using the K b (binding constant), and the stoichiometry of the complex urease:inhibitors (n) were determined (Table 1).The increase in K b and temperature between the enzyme and 3B4 or 3B4a suggests that the rise in temperature favors the formation of the complex.The 3B4-urease complex presented a stoichiometry that varied between 0.73 to 0.95; while the 3B4a-urease complex varied from 1.53 to 1.63 in relation to increasing temperature.In general, 3B4a showed higher K b values (9.55-41.69× 10 5 M −1 ) than 3B4 (0.004-0.04 × 10 5 M -1 ), suggesting a stronger interaction with the enzyme.This result agrees with the molecular modeling study, where a more stable 3B4a-urease complex was observed.Literature data also point out to interaction involving urease and other organic compounds, such as K b around 6.6-84.9× 10 5 M -1 for Biginelli adducts, 45 around 5.71-7.95× 10 3 M -1 for N-carbamothioylbenzamide, 41 and around 18.6-52.5× 10 5 M -1 for (+)-usnic acid. 46Thus, the substance 3B4a appears more similar to Biginelli adducts and (+)-usnic acid, while substance 3B4 is more similar to N-carbamothioylbenzamide.
The thermodynamic parameters of 3B4 and 3B4a (Table 1) were calculated using equation 3, described in the Experimental section.Figure S8 (SI section) shows a graphical representation of the linearized Van't Hoff equation for 3B4 and 3B4a.From these parameters, it was possible to determine the predominant types of intermolecular bonds in the interaction process.When ΔH > 0 and ΔS > 0, hydrophobic forces are predominant; for ΔH < 0 and ΔS < 0, there is a predominance of van der Waals forces and hydrogen bonds; and when ΔH < 0 and ΔS > 0, there is a predominance of electrostatic interactions. 46Thus, there is evidence that hydrophobic interactions are the preferential forces in the formation of the urease-3B4 and urease-3B4a complexes, which could be expected given the high hydrophobicity of these substances, observed in solubility studies (Table S1, SI section).The interaction process is entropically favorable, and the negative values observed for ΔG indicate that the interaction between urease and 3B4 or 3B4a is thermodynamically favored.Additionally, the ΔG values obtained in fluorescence studies are in agreement with the ΔG bind obtained in theoretical molecular docking studies, showing higher values for 3B4a than for substance 3B4.
Substances 3B4 and 3B4a have functional groups that favor hydrogen bonds with the active site of the urease enzyme, as demonstrated in the molecular docking study.However, the presence of aromatic groups results in preferential interactions of hydrophobic forces, observed in the calculated thermodynamic parameters.A qualitative study to evaluate the solubility of 3B4 and 3B4a showed that both substances are insoluble in water and aqueous solutions (Table S1, SI section).Furthermore, calculations performed on the Chemicalize 40 platform indicated a logP of 3.24 for 3B4 and 2.56 for 3B4a (Figure S9, SI section), providing further evidence of the lipophilicity of these substances.
Nonpolar molecules tend to self-associate in aqueous solution.The association of two nonpolar molecules in water decreases the nonpolar surface area, which reduces the amount of structured water and thus generates a favorable entropy of association.The enthalpy contribution to hydrophobic interactions is due to the presence of water molecules that occupy hydrophobic binding sites, interfering with the formation of hydrogen bonds with the receptor. 49hus, it is suggested that hydrophobic forces are also relevant in substance-enzyme and substance-solvent interactions.
Evaluation of conformational changes in the urease: 3D fluorescence studies 3D fluorescence studies were performed to investigate the bathochromic displacement observed from the addition of increasing concentrations of 3B4 or 3B4a, which may be evidence of possible alterations in the urease's secondary structure and the fluorophore residue's microenvironment.Figures 7 and S10 (SI section) illustrate the 3D fluorescence spectra for free urease and the complexes with 3B4 (Figure 7b) or 3B4a (Figure S10b).These spectra have three emission peaks: (i) scattering Rayleigh that occurs at λ ex = λ em ; (ii) fluorescence from the amino acid residues Trp and Tyr; and (iii) fluorescence from the highly excited electronic states of the aromatic residues composing the protein structure. 41As can be seen, there is a significant decrease in the fluorescence of peaks 2 and 3.The main parameters obtained for 3D fluorescence spectra and the Stokes displacement are summarized in Table 2, with the position and fluorescence intensity of the maximum emission peaks.From the results presented above, it was possible to observe that 3B4 and 3B4a led to similar spectral profiles (Figures 7 and S10), with a reduction of fluorescence intensity of peaks 2 and 3 by 78 and 38%, respectively, in the presence of 3B4 and by 65 and 44% for 3B4a, indicating conformational changes in the region of Tyr and Trp residues of urease enzyme.The decrease in emission at peak 3 could be linked to alterations in molecular structure that reduce the visibility of aromatic amino acid residues.This alteration makes energy absorption more intricate, subsequently reducing transitions to excited states that typically cause strong fluorescence emission, resembling the behavior observed in the free enzyme. 41Stokes displacement (peak 2), observed for longer wavelengths in the presence of 3B4 and for shorter wavelengths for 3B4a, represents conformational changes in the native protein structure, making the microenvironment of interaction more polar and apolar, respectively.

Competition assay
The inhibition mechanism of the 3B4 and 3B4a was evaluated in the presence of urea (urease enzyme substrate) and classical/commercial inhibitors (HU and NBPT).The ratio between binding constants in the presence and absence of competitors and substrate (K b '/K b ) was used to assess the influence on urease-3B4 and urease-3B4a interactions.K b ' and K b are the binding constants in the presence and absence of a competitor/substrate, respectively.The formation of the urease-inhibitor complex is favored in the presence of the competitor/substrate when K b '/K b > 1; and disadvantaged when K b '/K b < 1 (Table 3).
Analyzing the ratios between K b ' and K b , it can be observed that the binding constant decreased notably in the presence of the enzyme-substrate (urea).Therefore, complex formation decreased in the presence of urea, which strongly indicates that 3B4 and 3B4a interact in the active site of the enzymes, competing with urea for it.NBPT is a mixed-type commercial inhibitor. 50Thus, the ligand can bind to the enzyme's active site or other regions, modifying its conformation and inhibiting the enzymatic activity.HU is a competitive type inhibitor, interacting directly at the enzyme's active site. 51In the presence of these competitors, the formation of the urease-3B4 or 3B4a complex is  favored, which suggests that both substances do not interact only in the active site of the enzyme.However, in NBPT presence, the constants ratio (K b '/K b ) is much larger for 3B4a than for 3B4.As NBPT is a mixed inhibitor, there is less competition, which might suggest that the preferred binding site for 3B4a might be the active site.There is also the possibility that the interaction of NBPT with other parts of the enzyme is causing a change in the conformation of the active site, favoring the binding with 3B4 and 3B4a, or these ligands, in particular 3B4a, are interacting with the active site of the enzyme, as well as with a different allosteric site other than NBPT.
Interaction between urea and 3B4 or 3B4a FTIR is a technique commonly used to verify interactions between actives and excipients, where the reduction in the intensity or suppression of a peak and the appearance of new peaks are evidence of an interaction. 52Figure 8a illustrates the spectra obtained for urea (U), 3B4 (I-1) and the urea + 3B4 (1:1) mixture at times 0 (UI1-0 h),   24 (UI1-24 h) and 48 h (UI1-48 h). Figure 8b illustrates the corresponding spectra for 3B4a.For urea, there is a characteristic stretching band of the N-H bond at 3428-3429 and 3328-3329 cm -1 , a bending band of this bond at 1588-1589 cm -1 , and stretching bands for the C=O and C-N bond at 1673-1674 and 1457 cm -1 , respectively.The spectrum obtained for 3B4 (Figure 8a) shows a characteristic stretching band of the O-H bond at 3321 cm -1 and stretching bands from the C=N and C=C bonds at 1622, 1583, 1505, and 1456 cm -1 .For 3B4a (Figure 8b), there is a typical stretching band of the N-H bond at 3256 cm -1 , a bending band of the N-H aromatic secondary amine bond at 1516 cm -1 , and stretching bands from C=C bonds at 1472-1451 cm -1 .The N-H stretching band probably covers the stretching band for the O-H bond.
No significative changes in the spectrum profile were recorded when equal amounts of 3B4 or 3B4a were added to urea, regardless of the incubation time of the systems at room temperature.However, some bands were masked by others of greater intensity, such as the stretching of the O-H bond in the imine 3B4 and the stretching bands from the C=C bond in the amine 3B4a.Thus, no interaction between urea and 3B4 or 3B4a was observed in the experimental conditions tested, which may be evidence of these substances' compatibility for future agriculture applications.

Conclusions
Investment in the development of new urease inhibitors is a worldwide demand.In this work, we evaluated the anti-urease activity of substances 3B4 and 3B4a, showing that both are promising inhibitors.The results obtained in the molecular modeling showed that the functional groups 3-OH in 3B4 and 4-OH in 3B4a are essential for the inhibition of the enzyme urease; and that 3B4a forms a more stable supramolecular complex with the enzyme due to the formation of a hydrogen bond with the amine group.The molecular fluorescence study also proved that the enzyme's interaction with 3B4a is stronger.Hydrophobic forces are also preferably formed in the inhibitor-urease complex studied and both inhibitors are preferably competitive.Finally, FTIR analyses did not show interactions with urea, providing evidence of compatibility, which is essential for developing a new product for agronomic use.Finally, as perspectives, we propose to evaluate the efficiency and safety of these inhibitors in field studies.

Figure 2 .
Figure 2. Effect of 3B4 and 3B4a on soil (a) and C. ensiformis (b) ureases.Hydroxyurea (HU) and N-(butyl)thiophosphoric triamide (NBPT) were used as reference urease inhibitors.Data were analyzed by one-way ANOVA followed by Tukey post-test at a significance level of 5%.Distinct letters represent statistical differences.

Figure 6
Figure 6 illustrates the spectral emission profile of urease (1 μM) at different concentrations of 3B4 and 3B4a (5-80 μM) at pH 7.4 and 30 °C (the spectral profiles for the other temperatures evaluated can be accessed in

Figure 4 .
Figure 4. Interactions of 3B4 (a) and 3B4a (b) with Ni metals and residues from urease.The binding conformation of the ligands is shown in stick representation, and orange spheres represent the two Ni II ions.Coordination with Ni II and hydrogen bonding are dotted lines.

Figure 5 .
Figure 5. RMSD plots for the ligands (a) and RMSF plots for the backbone Cα atoms of the urease (b) with sequential numbering from the PDB.