Urease Inhibitors of Agricultural Interest Inspired by Structures of Plant Phenolic Aldehydes

Figure S1. Thermal stability of phenolic aldehyde derivatives 2A7 and 2B10 (PA derivatives), 2A9 (SA derivative), 2D2 (VN derivative) and the reference urease inhibitor N-(butyl)thiophosphorictriamide (NBPT) assessed by monitoring compounds mass changes as a function of increments in temperature. _______________ *e-mail: adefatima@qui.ufmg.br, lvmodolo@icb.ufmg.br 100 200 300 400 500 0 20 40 60 80 100 B


Introduction
Nitrogen (N) is a key nutrient absorbed by plants mostly as nitrate (NO 3 − ) and/or ammonium (NH 4 + ).Despite the great abundance of N in nature, less than 2% is bioavailable to plants.Biological N fixation, soil organic matter mineralization and lightening are natural processes known to increase the input of absorbing N in soil. 1,2However, these sources of bioavailable N are not enough to guarantee food supply for the growing world's population predicted to reach over 9.5 billion people within the next 35 years. 3Then, N fertilizers have been widely used to improve crop productivity, in particular urea, due to its high N content (46%), low price per N unit and easy management. 4rease (EC 3.5.1.5;urea amidohydrolase) catalyzes the hydrolysis of urea furnishing the gaseous products ammonia (NH 3 ) and carbon dioxide (CO 2 ). 5 It occurs in a variety of organisms including plants, fungi, bacteria and some vertebrates.Soil ureases play essential roles in N global cycle, also contributing to agriculture with respect to the availability of N as NH 4  + to plant growth in soils supplemented with urea-based fertilizers. 6,7On the other hand, soil ureases can also be detrimental for crop production specially when using the technique of covering fertilization.In this case, the urea applied on soil surface is rapidly hydrolyzed by soil ureases releasing NH 3 far away from rhizosphere allowing for N losses to the atmosphere due to the volatile nature of NH 3 .In fact, depending on climate and soil physicochemical properties, more than 50% of the N-urea applied to soil surface can be lost mainly by NH 3 volatilization. 4,8Besides negatively affecting plant N nutrition, excessive N losses to atmosphere as NH 3 remarkably impact natural ecosystems by contributing either directly or indirectly to acid rain, lakes and rivers eutrophication and formation of nitrous oxide, an atmospheric pollutant. 8ne of the strategies that have been adopted to minimize N losses as NH 3 is the use of urea-based fertilizers supplemented with urease inhibitors to slow down urea hydrolysis on soil surface and further increase the possibility of urea incorporation to soil by rain and irrigation. 9he potential of several classes of substances to inhibit the ureolytic activity of soil ureases have been investigated.][15] Indeed, the N-(butyl)thiophosphoric triamide (NBPT) was found to become a very effective inhibitor when transformed to its corresponding oxo-derivative (oxo-NBPT) by soil microorganisms. 16][19][20] Additionally, NBPT is more efficient in neutral soils with limited organic matter. 17,21Tropical soils exhibit organic matter and microbial biomass dynamics different from temperate soils. 22The amendment of soil with organic matter demanded from 2-to 4-fold NBPT to alleviate N volatilization by 20% in comparison to soils devoid of crop residues. 23Although other works investigated the NBPT effectiveness in tropical soils, [24][25][26] further research is needed for the development of novel and cost-effective urease inhibitors with improved efficiency in tropical soils and different environmental conditions.
Nature is undoubtedly a source of metabolites with potential to interfere with the activity of ureases as determined by in vitro assays with pure enzymes from Helicobacter pylori or Canavalia ensiformis (jack bean). 27][30] Although there is no report on the ability of natural phenolic aldehydes to inhibit ureases, it is likely that such secondary metabolites may work on ureases or be good prototypes for the design of urease inhibitors.Examples of natural phenolic aldehydes that have been explored as health promoters include protocatechuic aldehyde (PA), syringaldehyde (SA) and vanillin (VA).
The aim of the study herein presented was to use the natural products PA, SA and VA as building blocks for the development of four urease inhibitors of agricultural interest (Figure 1).Urea (urease substrate) or thiourea (urease inhibitor) core was also introduced to the structure of phenolic aldehyde derivatives synthesized (Figure 1).Then, in vitro assays were performed with pure jack bean urease to check the potential of synthesized compounds as inhibitors of ureolytic activity and disclose the mechanism of action of promising molecules.The effect of such phenolic aldehyde derivatives on soil ureases was addressed to confirm the potential of synthesized compounds for use as additives in urea-based fertilizers.

Preparation of phenolic aldehyde derivatives
An ethanolic mixture containing protocatechuic aldehyde (PA), syringaldehyde (SA) or vanillin (VA) individually (1 mmol; Sigma-Aldrich), ethyl acetoacetate (1.5 mmol) and urea or thiourea (1.5 mmol), here referred to as (thio)urea, and p-sulfonic acid calix [4]arene (0.5 mol%) was maintained under reflux and vigorous stirring for 8 h.After then, the mixture was filtered and the phenolic aldehyde derivative formed was recrystallized using ethanol.The phenolic aldehyde derivatives synthesized based on urea (2A7) or thiourea (2A9, 2B10 and 2D2) structure were characterized by 1 H and 13 C nuclear magnetic resonance (NMR), infrared, melting point and elemental analysis and the data compared to those reported elsewhere. 31,32The phenolic aldehyde derivatives were obtained in 49 to 80% yield.

In vitro urease activity assay
Initially, the phenolic aldehydes and derivatives 2A7, 2A9, 2B10 and 2D2 were screened for the ability to inhibit in vitro the ureolytic activity of purified Canavalia ensiformis (jack bean) type III urease (Sigma, St. Louis, Mo, USA).Each reaction medium containing 20, 1 and 10 mmol L -1 of phosphate buffer (pH 7.0), ethylene diamine tetra acetic acid (EDTA) and urea, respectively, 12.5 mU urease and compounds-test at 0 or 1.6 mM was incubated for 10 min at 25 °C.Reactions were stopped by adding 0.5 volume of 1% m/v phenol in 5 mg L -1 sodium nitroprusside (SNP) followed by the addition of 0.7 volume of 0.5% m/v NaOH in 0.1% v/v NaOCl solution.After samples incubation at 50 ºC for 5 min, the absorbance was measured at 630 nm to determine the amount of ammonium (NH 4 + ) formed. 33Hydroxyurea (HU) was used as a reference of urease inhibitor.Urease inhibition was determined in terms of percentage of NH 4 + formed in compounds-test reactions in relation to total urease activity in reactions without compounds.Three independent experiments were performed, each with four replicates.

Effect of phenolic aldehyde derivatives on the kinetic parameters of jack bean urease
The inhibition profile exhibited by the natural product derivatives 2A7, 2A9, 2B10 and 2D2, synthesized in this study, was determined by incubating inhibitors at concentrations necessary to inhibit jack bean urease activity between 30 and 40% (from 0.3 to 1.6 mM) in reaction medium containing 20 mmol L -1 phosphate buffer (pH 7.0), 1 mmol L -1 EDTA, urea (ranging from 1 to 32 mmol L -1 ) and 12.5 mU urease.The stoppage of reactions, NH 4 + quantification and urease inhibition calculation were done as described previously.Jack bean urease kinetic parameters such as initial velocity (V o ), K M (Michaelian constant) and maximum velocity (V max ) were obtained using Hyper32 software. 34The OriginPro8 (Origin Lab, Northamptom, MA) software was used to pursue Michaelis-Menten hyperbolas and Lineweaver-Burk plots.The equilibrium dissociation constants for ureaseinhibitor complex (K i ) and for urease-urea-inhibitor complex (K' i ) were determined from the α and the α' values. 35
Sieved soil samples (0.5 g; particles smaller than 2 mm) were incubated with 72 mmol L -1 urea in the absence or 3.2 mmol L -1 of 2A7, 2A9, 2B10, 2D2 or NBPT (used as a reference of soil ureases inhibitor) at 37 o C for 1 h.Ureases activity was stopped by incubating the systems with 5 mL of 1 mol L -1 KCl in 10 mmol L -1 HCl for 30 min at 25 °C.A supernatant aliquot was taken after soil decantation and added to a solution containing 3.4, 2.5 and 2.5% of sodium salicylate, sodium citrate and sodium tartrate, respectively, and 120 mg L -1 SNP.After 15 min incubation at 25 °C (under darkness), 0.1 volume of 3.0% NaOH in 1.0% sodium hypochlorite was added to each reaction system following incubation under darkness for 1 h at 25 °C and 600 rpm.The NH 4 + formed was detected by spectrophotometric measurements at 660 nm.
Then, assays using different concentrations (from 0.05 to 3.2 mM) of phenolic aldehyde derivatives or NBPT were performed to determine the concentration of compound-test that causes 50% inhibition of soil ureases (IC 50 ).Independent experiments were performed, each with at least five replicates.

Statistical analysis
Data were subjected to analysis of variance (ANOVA) by general linear model (GLM) procedure and contrast analysis at 5% significance level using the software R (Software Foundation, Boston, MA, USA).

Inhibition of ureolytic activity of jack bean urease
The in vitro assay with purified jack bean type III urease showed that, among the natural products tested, only protocatechuic aldehyde (PA; at 1.6 mM) effectively inhibited the enzyme activity (68% inhibition) while vanillin (VN) and syringaldehyde (SA) marginally reduced the production of NH 4 + (Figure 2).The derivatives 2A7 and 2B10, originated from PA, were the most potent urease inhibitors showing results (94% urease inhibition) comparable to that observed for the standard inhibitor hydroxyurea (HU; Figure 2).Compounds 2A9 and 2D2, derived from VN and SA, respectively, caused enzyme inhibition that averaged 58.6% (Figure 2).

Mechanism of action of phenolic aldehyde derivatives toward jack bean urease
Urease is categorized as a Michaelian enzyme since the graph of initial velocity (V o ) versus urea concentration exhibits a typical hyperbolic behavior (Figure 3).The average urea K M (Michaelian constant) and urease maximum velocity (V max ) in reactions free of urease inhibitor were, respectively, 3.4 ± 0.4 mM and 8.1 ± 0.4 µmol NH 4 + min -1 mg -1 protein.The addition of 2A7, 2A9, 2B10 or 2D2 to the reaction medium caused a concentration-dependent increment of urea K M (apparent Michaelis constant: K M (app) ) and decrease of urease V max (apparent maximum velocity: V max (app) ) (Figure 3).All the phenolic aldehyde derivatives tested behaved as mixed inhibitors, as attested by the lines intersection in the second quadrant of Lineweaver-Burk plots (Figure 3; right column).
The derivative 2A7 was the most potent mixed inhibitor since the K i and K' i values for complexes formed with this compound were the lowest in comparison with the others inhibitors (Table 1).In general, K i values were lower than the K' i values for complexes related to the same inhibitor by 2.4-to 15.5-fold.

Soil ureases activity assay
When tested at 3.2 mM, all phenolic aldehyde derivatives were able to inhibit soil ureases at different extents; 2A7 (PA derivative) and 2D2 (SA derivative) were found to be as efficient as NBPT (commercial inhibitor; 40% enzyme inhibition) while the VN-derived 2A9 and the PA-derived 2B10 inhibited soil ureases by up to 30% (Figure 4).
The concentration of 2A7 and 2D2 necessary to cause the inhibition of soil ureases by 50% (IC 50 ) were, in average, 3.25 mM.The derivative 2A9 exhibited a maximum inhibitory activity of 16% when used at 0.05 mM or higher concentrations.There was not a pattern in the behavior of results observed for the replicates of independent experiments performed with 2B10 and NBPT, which did not allowed for determining the IC 50 values for such inhibitors.

Thermal stability of natural phenolic-derived urease inhibitors
The thermal stability of the natural phenolic aldehydederived urease inhibitors, assessed by mass changes of compounds as a function of fast increments in temperature, revealed that the first event of mass loss (decrease by 5%) for derivatives 2A7, 2D2, 2A9 and 2B10 occurred at 254, 253, 244 and 226 ºC, respectively (Figure S1, Supplementary Information section).This same event was observed when NBPT was subjected to 151 ºC (Figure S1, Supplementary Information section).The second event, characterized by 20% mass loss, took place at 200.5 ºC for NBPT while similar percentage of mass loss for the phenolic aldehyde derivatives was registered at 241, 269.1, 271.2 and 276.4 ºC for 2B10, 2A7, 2A9 and for 2D2, respectively, (Figure S1, Supplementary Information section).

Discussion
The potential of a series of plant natural products as urease inhibitors of clinical and/or agricultural interest has been documented. 27Among them, methyl gallate (phenolic ester) and its glycosylated derivative isolated from Paeonia lactiflora roots were shown to be promising with respect to the inhibition of H. pylori urease. 37The promising effect of these phenolic esters prompted us to investigate the potential of phenolic aldehyde derivatives as urease inhibitors of agricultural interest.Thus, the plant natural products protocatechuic aldehyde (PA), vanillin (VN) and syringaldehyde (SA) were selected as prototypes for the design of new urease inhibitors based on urea or thiourea scaffolds, urease substrate and inhibitor respectively.
In vitro assays revealed that PA per se decreased the activity of jack bean type III urease, while VN and SA were found to be inactive as they caused less than 1% enzyme inhibition (Figure 2).Besides methyl gallate and related derivatives, 37 other plant phenolic compounds, such as (iso) quercitrin, avicularin, guaijaverin, flavonoid glucosides and shoreaphenol, were also reported to inhibit jack bean urease. 29,30,38,39Interestingly, structural modifications on PA, VN and SA dramatically improved their ability to inhibit the ureolytic activity of jack bean urease (Figure 2).Indeed, the conversion of phenolic aldehydes to derivatives bearing urea or thiourea core (Figure 1) yielded the VN derivative  2D2 that is 230-fold more potent than VN by itself, the SA derivative 2A9 of potency 66-fold higher than that of SA by itself and the PA derivatives 2A7 and 2B10 were about 40% more potent in comparison with PA.Under our experimental conditions, the novel urease inhibitors 2A7 and 2B10 were as potent as hydroxyurea (HU; known urease inhibitor), while 2D2 (also novel) and 2A9 was less effective than HU (Figure 2).The potential of the derivative 2A9 as an inhibitor of one of the jack bean ureases was recently reported, 40 although the experimental conditions were different from the one reported herein.Thus, the outstanding performance of these phenolic derivatives, compared to the natural products they originated from, might be attributed to the combination of a catechol skeleton with urea or thiourea core.We performed kinetic experiments varying urea concentration in urease-catalyzed reactions containing each inhibitor at fixed concentrations (Figure 3).The urea K M value obtained from reactions carried out in 20 mmol L -1 phosphate buffer (pH 7.0) was, in average, 3.4 mM and urease V max 8.1 µmol NH 4 + min -1 mg -1 prot.Other studies with jack bean urease, under experimental conditions distinct from the one used here, reported urea K M values ranging from 1.0 to 4.0 mM. 7,41he kinetic behavior of jack bean urease in the presence of the phenolic aldehyde derivatives is consistent with the one expected for an enzyme in the presence of a mixed inhibitor.Mixed inhibitors are known to be capable of binding both the free enzyme (forming an enzyme-inhibitor complex) and the enzyme-substrate complex (forming an enzyme-inhibitor-substrate complex). 42The values obtained for the dissociation constants for both urease-inhibitor and urease-urea-inhibitor complexes indicate that the phenolic aldehyde derivatives synthesized bind more efficiently to the urease active site in comparison to allosteric ones (Table 1).The potency of compounds with respect to the binding to urease active site is 2A7 > 2B10 > 2D2 > 2A9.As for the binding to allosteric site(s) the order of potency is 2A7 >> 2B10 = 2D2 = 2A9.
The derivatives 2A7 and 2D2 were the most efficient compounds to inhibit soil ureases, clustering together with the reference inhibitor NPBT (Figure 4).These synthesized compounds are able to inhibit soil urease activity by 50% when used at 3.25 mM.Notably, the maximum inhibition of soil ureases exhibited by the SA derivative 2A9 was 16% when employed at 0.05 mM, no matter higher concentrations would be applied on soil.In fact, in the case of soil ureases, (a)biotic conditions such as temperature, pH, moisture and the presence of different types of clays, organic matter and viable microorganisms particularly are known to affect ureases performance. 14Moreover, soil matrix comprises complex physicochemical features and biological processes that may culminate in the chemical modifications of xenobiotic substances, 43 as is the case of synthetic urease inhibitors.Such chemical transformations triggered by soil microbiota may result in loss or increment of the function of a certain compound. 43Alternatively, the complex nature of soil matrix may affect the bioavailability of the xenobiotic for interaction with the target enzymes.Taking these into account it is likely that the PA derivative 2B10 undergoes some structural transformation caused by soil microbiota as it is known to occur with NBPT. 16t is well known that NBPT is sensitive to heat. 44or the most active derivatives on soil (2A7 and 2D2), thermogravimetric analysis shows no decomposition of such compounds up to 170 o C (Figure S1; Supplementary Information section).In addition to the thermal stability, the derivatives 2A7 and 2D2 are obtained in a single synthesis step (78% average yield) after a simple purification procedure (recrystallization) that furnishes compounds as  solid materials.These are desirable features for obtaining urease inhibitors to be used as additive in urea-based fertilizers.

Conclusions
Overall, the hybridization of structures of the natural products PA, SA and VN with (thio)urea core furnished derivatives with inhibitory effect on ureases activity displaying mechanisms of action typical of mixed inhibitors.The interesting physicochemical features of phenolic aldehyde derivatives herein studied, together with their ability to inhibit soil ureases, make these compounds, especially 2A7 and 2D2, promising candidates for further studies as additive in urea-based fertilizers.

Figure 2 .
Figure 2. Inhibition of jack bean urease by phenolic aldehydes and its derivatives.The compounds hydroxyurea (HU), protocatechuic aldehyde (PA), syringaldehyde (SA), vanillin (VA), 2A7 and 2B10 (PA derivatives), 2A9 (SA derivative) and 2D2 (VN derivative) were employed at 1.6 mM in reactions containing 10 mmol L -1 urea and 12.5 mU urease.Results are representative of three independent experiments, each with four replicates.Different letters indicate significant difference (p < 0.05 by contrast analysis) among the compounds.

Figure 3 .
Figure 3. Representative Michaelis-Menten hyperbola and Lineweaver-Burk plots for jack bean urease in the presence of phenolic aldehyde derivatives.The compounds 2A7 and 2B10 (PA derivatives), 2A9 (SA derivative) and 2D2 (VN derivative) were employed at different concentrations (0.3 to 1.6 mM) in reactions containing urea ranging from 1 to 32 mmol L -1 and 12.5 mU urease.A V o versus urea concentration plot obtained from data of assays with 2A7 is shown to exemplify the Michaelian behavior of urea catalysis.

Figure 4 .
Figure 4. Effect of phenolic aldehyde derivatives and NBPT on the activity of soil ureases.The compounds N-(butyl) thiophosphorictriamide (NBPT), 2A7 and 2B10 (PA derivatives), 2A9 (SA derivative) and 2D2 (VN derivative) were applied to soil at 3.2 mM in the presence of 72 mmol L -1 urea.Results are representative of independent experiments, each with at least five replicates.Distinct letters indicate significant difference (p < 0.05 by contrast analysis) among the compounds.

Table 1 .
Inhibition constants for phenolic aldehyde derivatives towards jack bean type III urease : equilibrium dissociation constant for urease-phenolic aldehyde derivative complex; K' i : equilibrium dissociation constant for urease-ureaphenolic aldehyde derivative complex; values are the mean ± standard deviation of triplicate determinations from a representative experiment. i