Electronic structure and molecular properties of nitisinone and mesotrione in water

Context Nitisinone is a medium-sized organic molecule that is used in treating hereditary tyrosinemia type 1 (HT-1). The structurally analogous mesotrione, however, is used as a pesticide/herbicide. What molecular properties are responsible for the similarity/dissimilarity of these molecules is investigated here. The solvent effect reduces the electron affinity to rather negative values and causes the negative electron affinity which manifests itself in a very high positive absolute reduction potential. Methods B3LYP method was utilized for a geometry optimization of nitisinone and mesotrione in their neural and ionized (L0, L+, L−) forms of 6 structures. The calculations were conducted in water as a solvent using conductor-like polarizable continuum model (CPCM), nitisinone also in vacuo. The complete vibrational analysis at the true energy minimum allows evaluating the thermodynamic functions with focus to the zero-point energy and overall entropic term. The change of the Gibbs energy on reductions and/or oxidation facilitates evaluating the absolute reduction and absolute oxidation potentials. Also, DLPNO-CCSD(T) method that involves the major part of the correlation energy has been applied to nitisinone and mesotrione and their molecular ions.

Nitisinone and mesotrione belong to the class of β-triketone herbicides; these are utilized instead of triazine herbicides such as atrazine.Triketones act as inhibitors of the 4-hydroxyphenylpyruvate dioxygenase (HPPD); they, however, increase tyrosine level in plasma.HPPD is an enzyme playing an important role in catabolism of tyrosine.
HPPD, expressed mainly in the liver, converts 4-hydroxyphenyl-pyruvate to homogentisic acid (HGA); it is used in treating hereditary tyrosinemia, type 1 [1].In 2020, the European Commission has issued the approval of the extended indication for nitisinone to include treatment of adult patients with rare disease alkaptonuria (AKU) based on successful results of the SONIA2 trial [2][3][4][5][6].Nitisinone is used as a medicament (Orfadin) along with restriction of phenylalanine and tyrosine on diet.
Mesotrione is used as a selective herbicide effective mainly to maize.HGA acts as a precursor for synthesis of plastoquinone and α-tocopherol in plants; the inhibition effect to HPPD causes bleaching and death in weed [7].Mesotrione can cause eye irritation and could contribute to obesity but it is non-toxic by oral consumption.A relationship was found between levels of tyrosine and accumulation of lipid droplets in the non-alcoholic fatty liver disease [8].
The present study focuses to the calculations of the electronic structure of nitisinone and mesotrione and to determine molecular descriptors that cause similarity and/or dissimilarity of these species in vacuo and in solvent (water) with focus to redox properties.
In the first step, the structures of nitisinone [13] and mesotrione [14] have been used for the geometry optimization.After the gradient criteria indicated the global energy minimum, a set of molecular properties was evaluated: the energies of HOMO (the highest occupied molecular orbital) and LUMO (the lowest unoccupied molecular orbital), the permanent dipole moment (p), isotropic value of the quadrupole moment (Q), and the isotropic value of the dipole polarizability (α).The vibrational analysis facilitates calculation of the partition function from which the zero-point vibration energy and standard thermodynamic functions (inner energy U ø , enthalpy H ø , entropy S ø , and Gibbs energy G ø ) were evaluated.
Using energies of neutral molecule, molecular cation, and anion, the vertical ionization energy E i = E + − E 0 and electron affinity E eg = E − − E 0 were calculated in the frozen geometry.These quantities are further processed to get the molecular electronegativity = (E i − E eg )∕2 ) according to Mulliken, the chemical hardness = (E i + E eg )∕2 accord- ing to Pearson, and the electrophilicity index = 2 ∕2 according to Parr [15][16][17].The adiabatic variants of the ionization energy and electron affinity were evaluated after the geometry optimization of the respective molecular ions L + and L − .The reaction Gibbs energy Δ r G ø served for the evaluation of the absolute oxidation potential and absolute reduction potential via thermodynamic equation for one-electron process, e.g., E abs ø (L 0 /L q ) [V] = -Δ r G ø [J mol −1 ]/F with the Faraday constant F. The ab initio calculations were conducted also by DLPNO-CCSD(T) method (domain-based local pair natural orbitals-coupled cluster singlets, doublets, and triplets) that account to the major part of the correlation energy.The augmented basis set aug-cc-pVTZ for S-{13,13,13,111/77111/111/11}, F-, O-, N-, C-{88,111/3111/111/11}, and H-{3111/111/11} was applied with 1426 and 1518 basis functions for nitisinone and mesotrione, respectively; the auxiliary basis set was augcc-pVTZ/C [18,19].The tedious numerical evaluation of the gradient, vibrational displacements, and the polarizability prevents fast geometry optimization and vibrational analysis.Thus, the calculations were conducted in the fixed geometry optimized by the B3LYP method.

Results and discussion
The set of molecular properties calculated for nitisinone and mesotrione is presented in Table 1.There are reports about three polymorphs of mesotrione confirmed by the X-ray powder diffraction and infrared spectra in the solid state [20][21][22].Nitisinone was considered in in vacuo as well as in water yielding the following comparisons:  Table 2 brings electronic energies and standard Gibbs energies of nitisinone and mesotrione after the geometry optimization followed by the complete vibrational analysis.These data serve for evaluating the adiabatic ionization   energy E i and adiabatic electron affinity E eg ; these allow determining the derived electronic properties such as electronegativity , hardness η, and electrophilicity ω.The reaction Gibbs energies serve for determining the absolute oxidation and reduction potentials at standard conditions.According to Table 3, the electronic redox properties (E i , E eg , χ, η, ω) and thermodynamic properties (E ox ø , E red ø ) of nitisinone are essentially the same as for mesotrione.The chemical hardness is a measure of the willingness of the molecule against the electron transfer; for nitisinone and mesotrione, it is η = 39 kcal mol −1 (rather low value).The electrophilicity index is extremely high: w = 200 kcal mol −1 .
Consequently, the absolute reduction potential is very high E red ø = + 3.80 V; for amino acids and small biogenic molecules, it is typically ~ + 1 V when the same methodology (B3LYP) is used.Noticeably, ω is defined on the basis of electronic energies of L 0 , L + , and L − species and it is common for oxidation and reduction processes.On the contrary, redox potential is derived from the standard reaction Gibbs energies and is different for oxidation and reduction processes.
An inspection to Table 3 confirms that the overall temperature-dependent correction to Gibbs energy G ø (L q ) − E el (L q ) is about 111 and 133 kcal mol −1 which in fact is only 0.01% of the total energy.Therefore, it might be possible to express the redox potential by an approximate formula E r * = − Δ r E el /F when the Gibbs energy is not available due to a very tedious numerical evaluation of the vibrational frequencies.The ionization energy and/or the electron affinity then serve for the reaction electronic energies.Some of the molecular properties of nitisinone and mesotrione are classified as bulky properties and they increase with increasing molar mass: absolute value of the quadrupole moment |Q|, isotropic dipole polarizability α, zero-point energy E zpe , and the total entropic term S•T ø belong to them.The dipole moment reflects a delicate separation of the barycenters of positive-negative charges: p = 6.3 and 8.9 D for nitisinone and mesotrione, respectively.
Application of a more exact DLPNO-CCSD(T) method (M2) to above-mentioned neutral and ionized species in geometries optimized by B3LYP method (M1) brings results presented in Table 3.The ionization energy is high (E i = 198 and 222 kcal mol −1 , higher than by M1), and the electron affinity is very negative (about E eg ~ − 83 kcal mol −1 , similar to M1). High ionization energy is reflected into increased molecular electronegativity χ and increased chemical hardness η.The resulting electrophilicity index is very high (ω ~ 170 kcal mol −1 ) which predetermines high reduction potential.
The electrostatic potential distributed along the surface of the molecule is taken against a unit charge and is known as the molecular electrostatic potential (MEP) [23,24].It is visualized as a three-dimensional contour map plotted on the isovalue surface of charge density.It identifies sites possessing positive or negative values that are suitable for nucleophilic and/or electrophilic interactions along the molecule.HyperChem software has been exploited for generating MEP [25] that is drawn in Fig. 2 along with the optimized molecular structures.With expectations, high negative potential occurs in the sites of electronegative atoms such as fluorine and oxygen (blue colored); positive potential is delocalized over the rings and localized at the sulfur atom and/or carbon od the -CF 3 moiety (red colored).
The calculated IR spectra over the whole wavenumber range are shown in Fig. 3 for nitisinone and in Fig. 4 for mesotrione.The DFT-B3LYP is not a perfect tool for reproducing the vibrational spectra quantitatively; however, after some scaling (scaling factor about 0.96) this method yields acceptable results [26].Both spectra can be divided into three domains: above 3000 (I), below 1700 (II), and below 250 (III) cm −1 .The experimental window is 500-4000 cm −1 .For nitisinone, the high-frequency domain is limited by medium-intensity bands at I, 2950-3250 cm −1 where a number of low-intensity transitions exist in between.This feature is well reproduced by calculations.There is a broad silent gap between domains I and II.Domain II starts with highintensity peaks at 1700 cm −1 followed by a huge number of resolved, weakly resolved, or oddly resolved transitions [20].For the far-IR domain III, experimental data are missing.
Last but not least, some notes regarding biological activity are given.Both studied compounds, nitisinone and mesotrione, are potent inhibitors of 4-hydroxyphenylpyruvate  It could be hypothesized that a direct oxidation of Fe II to Fe III by strong oxidizing agents (mesotrione and nitisinone) may destroy the functionality of HPPD.However, experiments show that the {HPPD•Fe II• NTBC} adduct under aerobic conditions maintains typical electronic transition at 450 nm for a long time, which confirms the presence of Fe II [29].There is an important difference in price of nitisinone and mesotrione (at the moment mesotrione is 15 times cheaper).

Conclusions
In conclusion, DFT-B3LYP calculations show that two derivatives of β-triketones, nitisinone and mesotrione, possess essentially molecular parameters in water such as energies of HOMO and LUMO, the ionization energy and electron affinity, and derived electronic properties such as molecular electronegativity, chemical hardness, and electrophilicity index.They differ in the dipole moment (high value p = 7.3 and 8.9 debye) and bulk properties such as isotropic polarizability (α = 270 and 309 a 0 3 ), solvated surface and volume, zero-point vibration energy (143 and 163 kcal mol −1 ), and the total entropic contribution (S•T ø = 42 and 44 kcal mol −1 ).The value of the absolute reduction potential is very high (E red ø = 3.80 and 3.82 V) so that these species are effective oxidizing agents.DLPNO-CCSD(T) method, that includes explicitly the correlation energy, applied to nitisinone and mesotrione confirms qualitative predictions of more approximate B3LYP in water: high adiabatic ionization energy, very negative electron affinity (E eg ~ − 83 kcal mol −1 ), high electronegativity, low chemical hardness, high electrophilicity index (ω = 170 kcal mol −1 ), and high dipole moment (p = 7.8 and 9.5 debye).

370
Page 2 of 7 (i) HOMO and LUMO energies are the same in vacuo and in water [items 1, 2]; (ii) as a solvent effect, the total energies of the neutral molecule and molecular ions are lowered (items 3 through 5); (iii) upon solvation, the ionization energy is lowered and electron affinity adopts more negative values (items 6, 7); (iv) the chemical hardness is reduced as well (item 8); (v) the polar solvent rises the dipole moment substantially (item 10); (vi) the isotropic value of the quadrupole moment is solvent-independent (item 11); (vii) the isotropic dipole polarizability is raised in water (item 12); (viii) the contributions to thermodynamic functions (translational, rotational, and vibrational) are almost insensitive to solvation (items 15-17 and 20-23); (ix) the electronic contributions

Table 1
Molecular properties of nitisinone and mesotrione by DFT-B3LYP a ). Abbr.str structure

Table 3
[27,28]nase (HPPD), 40-50 kDa.In active sites, this enzyme contains Fe II in hexacoordinate form[27,28].In plants, HPPD exists in its monomeric form but dimers or tetramers were identified in bacteria.In plants, HPPD transforms tyrosine into other products necessary for living: after several steps, the lack of tyrosine prevents the formation of chlorophyll, which in turn initiates plant death.HPPD inhibitors include not only nitisinone and mesotrione but