Hui
Xing‡
a,
Sarah K. M.
McGregor‡
a,
Bruna D.
Batista
ab,
Cassidy
Whitefield
a,
Isobella S. J.
Stone
a,
Cody-Ellen
Murray
c,
Rebecca M.
Hurst
c,
Yizhou
Liu
a,
Sharon
Chow
a,
Tyler
Fahrenhorst-Jones
a,
Qi
Zhao
a,
Sevan D.
Houston
a,
Shu-Hong
Hu
a,
Thierry
Lonhienne
a,
Amanda
Nouwens
a,
Jed M.
Burns
a,
G. Paul
Savage
d,
Gimme H.
Walter
c,
Luke W.
Guddat
a,
Michelle A.
Rafter
b and
Craig M.
Williams
*a
aSchool of Chemistry and Molecular Biosciences, University of Queensland, Brisbane, 4072, Queensland, Australia. E-mail: c.williams3@uq.edu.au
bHealth and Biosecurity, CSIRO, Ecosciences Precinct, Brisbane, 4102, Queensland, Australia
cSchool of the Environment, University of Queensland, Brisbane, 4072, Queensland, Australia
dCSIRO Manufacturing, Ian Wark Laboratory, Melbourne, 3168, Victoria, Australia
First published on 10th December 2024
Weed management is an essential intervention for maintaining food security and protecting biodiversity but is heavily reliant on chemical control measures (i.e., herbicides). Concerningly, only one herbicide has been developed with a new mode of action (MOA) since the 1980s. Therefore, alternative strategies for preventing weed growth need to be explored. The lesser-known concept of halting weed growth through herbistasis could be one strategy to alleviate the lack of success in obtaining new MOA leads, but this type of activity has rarely been investigated. Herein reported is a bioisosteric cyclooctatetraene (COT) for phenyl ring replacement tactic, using the commercial acetolactate synthase (ALS) inhibitor metsulfuron methyl, that has unearthed a rare agent displaying herbistatic properties against the weed, Cryptostegia grandiflora (rubber vine).
The term “herbistatic agent” was initially introduced in 1948 in the context of the herbicide, 2,4-D (1),12 associating it with its synthetic auxin herbicidal activity. However, this term has garnered only sporadic mention since that time, with only three additional examples of herbistatic agents identified: (1) glaucarubolone (2) in 1996 – a naturally occurring quassinoid which was found to inhibit NADH oxidation in soybean cells,13 (2) diflufenzopyr (3) in 1999 – assumed to block the polar transport of naturally occurring auxin indole-3-acetic acid (IAA) (4),14 and (3) thaxtomin (5) in 2013 – a bacterial phytotoxin presumed to inhibit cellulose biosynthesis15 (Fig. 1). However, parallels can be drawn to the related “plant dwarfism” effect, where application of brassinazole (6),16 or related systems such as fenarimol (7)17 (Fig. 1), can induce the inhibition of the brassinosteroid biosynthetic pathway. Notably, the effects of brassinosteroid inhibition are reported to be reversible. The only identifiable definition for “herbistatic” is reported as “inhibiting the growth of weeds, which can be reversible under certain conditions”.15 This seemingly implies that both inhibition of the target enzyme and stunted plant growth are key features of herbistasis. This definition is closely analogous to that used for fungistatics (e.g., azoles and echinocandins), which “halt the growth of, but does not kill the fungus”18 and for bacteriostatics (e.g., chloramphenicol, ethambutol, tetracycline).19
Conceptually, from a weed control perspective, treatment with a herbistatic also provides a beneficial growth advantage to agricultural or native species (e.g., increased availability of essential nutrients through stunted weed growth), but a herbistatic would likely require a different mode of resistance development as the weed “killing” mechanism is absent (i.e. possibly extending the utility of known MOAs). The latter might be beneficial for prolonging resistance or treating resistance species that have blocked “killing” mechanisms, but little is known about herbistatics as there are so few, and whether they have select advantages over herbicides [i.e. as do some bacteriostatics (e.g., clindamycin20) over bactericides]. To that end, the herbicide metsulfuron methyl (MM) (8) was chosen as a model template to explore herbistasis. MM is a highly potent broad-spectrum herbicide (Fig. 2), and a member of the sulfonylurea herbicide class.21 This herbicide is a powerful inhibitor of acetolactate synthase (ALS) [also known as acetohydroxyacid synthase (AHAS)], which is the first enzyme in the branched-chain amino acid biosynthesis pathway and critical to plant growth.22–25 Furthermore, anecdotal reports claim that MM limits plant growth at sub-lethal dosages without affecting the biomass through its ability to limit cell division.26 Consequently, it may already have structural features conducive to herbistatic influences.
Reported herein is the utilisation of cyclooctatetraene (COT, 9),27,28 as a more dynamic biological probe to pursue any underlying herbistatic properties of MM through the synthesis of COT-MM (10). The COT ring system has previously been deployed as a phenyl ring (11) bioisostere/biomotif.27,28 However, unlike other phenyl ring bioisosteres, such as bicyclo[1.1.1]pentane, bicyclo[2.2.2]octane, and cubane,29 the COT ring contains both steric and π character components. Furthermore, the COT ring system has shape shifting properties arising from the dynamic equilibrium that exists between ring inversion (e.g., 12), π-bond shifts, and to a much lesser extent valence tautomerization (i.e., bicyclo[4.2.0]octatriene 13, and corresponding 14) (Fig. 2).
Silyl deprotection of the sulfone (23), followed by conversion of the intermediate sulfinic acid to the sulfonamide 24 was achieve in 60% yield.35 To introduce the urea functionality and complete the synthesis, sulfonamide 24 was reacted with phenylcarbamate (26), derived from 1,3,5-triazin-2-amine (25), to provide COT-MM (10) in 37% yield (Scheme 1).
Scheme 2 (a) Reaction of singlet oxygen with COT via the valence tautomer. (b) Synthetic route employed to access COT-BT (28), with its parent herbicide BT (29) also shown. |
In the view that 1,4-disubstituted COT ring systems are most efficiently obtained by rhodium catalysed valence isomerisation of 1,4-disubstituted cubanes,42 this strategy was adopted for the synthesis of COT-BT. Half hydrolysis of cubane diester 30 was performed with sodium hydroxide, and the resulting acid converted to the corresponding acid chloride (31).43 Barton chlorination, using 32, gave 33 on irradiation in carbon tetrachloride. Subsequent hydrolysis, followed by borane reduction, gave the cubane alcohol 34.42 Conversion to the thioacetate 35 proceeded smoothly using a Mitsunobu reaction.44 The thioacetate was then hydrolysed to thiol 36, and subsequent exposure to sodium hydride followed by treatment with N,N-dimethylcarbamoyl chloride (37) afforded 38 in 66% yield. COT-BT (28) was then simply obtained through gentle heating of 38 in toluene at 60 °C using the commercially available rhodium(I) catalyst [Rh(nbd)Cl]2 (Scheme 2b).
Initially, an in vitro study was performed over an 18 day post treatment period to assess the plant growth inhibition of COT-MM and COT-BT in relation to their corresponding commercial herbicides MM and BT (i.e. phenyl ring containing). Compared to the vehicle control, the application of MM and COT-MM reduced plant vertical growth by 100% and 96%, respectively (Fig. 3a). By comparison, treatment using COT-BT resulted in no reduction in growth, and BT by only 21% compared to the negative vehicle control (Fig. 3a).
Fig. 3 The effect of COT-MM and COT-BT on vertical growth of rubber vine in vitro (a) and grown under glasshouse conditions (b–e). (a) The experiment was conducted over 18 days and the results (plant height 18 DAT–1 DAT) are compared to the commercial herbicides MM and BT and two vehicle only controls. Control 1 (1% DMSO, 1% Tween 80), BT (0.1 mM BT, in 1% DMSO, 1% Tween 80), COT-BT (0.1 mM COT-BT, in 1% DMSO, 1% Tween80), Control 2 (1% Tween 80), MM (0.1 mM MM K+ in 1% Tween 80), and COT-MM (0.1 mM COT-MM K+ in 1% Tween 80). n = 13 independent replicates. The diamond symbol represents the mean, the horizontal line represents the median, the box represents interquartile range (IQR), and the whiskers represent 1.5 × IQR. (b) Plant mortality over time (top) and the impact on rubber vine height (bottom), with reference lines at 35 mm and 155 mm, n = 16 independent replicates. Control (0.10 mM K2CO3 in 1% Tween 80), MM (0.10 mM MM K+ in 1% Tween 80), and COT-MM (0.10 mM COT-MM K+ in 1% Tween 80) measured periodically over 20 DAT. (c) The visual effects of treatments on rubber vine plants 20 DAT, and the relative impact on shoot dry biomass (SDB) and root dry biomass (RDB) for MM and COT-MM compared to the control. (d) Long-term impact of chemical treatment of COT-MM on rubber vine grown under glasshouse conditions. Rubber vine plants (n = 24) were treated with COT-MM (0.1 mM COT-MM K+ in 1% Tween 80) and the plant height and plant mortality were recorded weekly from 0 to 70 DAT. (e) The visual effects of the treatment on the rubber vine plants at 70 DAT grown in 5 cm square pots. Dead plants (n = 6) marked with red asterisk. For (a), (b), and (d) the different letters above columns indicate significant differences in average plant height over time (p < 0.05, Tukey's HSD test), error bars represent the standard deviation, and data points outside the whiskers represent outliers. [Note: due to the instability of sulfonylurea functional groups in dimethylsufoxide (DMSO), conversion of COT-MM (and MM) to its potassium salt was necessary to improve aqueous solubility for plant experiments and enzymatic assays (see ESI†)]. |
Application of both MM and COT-MM led to a statistically significant decrease of 15% in shoot dry biomass (SDB), while root dry biomass (RDB) remained unaffected by these treatments. COT-BT showed no significant effect on any plant growth parameter assessed, and BT did not significantly impact SDB or RDB (Fig. S1 and Table S1, ESI†). Both negative control treatments (1% DMSO with 1% Tween 80, and 1% Tween 80) also demonstrated indistinguishable effects on plant vertical growth, as well as on SDB and RDB (Table S1, ESI†).
In a subsequent glasshouse study, the effect of COT-MM was compared with that of MM over a 20 day post-treatment period, with both demonstrating a statistically significant negative impact on rubber vine growth (Fig. 3b) (Table S2 and Fig. S2, ESI†). Control-treated (i.e., K2CO3 added) plants exhibited a 4-fold increase in vertical growth from 0 to 20 days after treatment (DAT) (from roughly 35 mm to 150 mm), whereas plants treated with MM and COT-MM remained the same height as when first treated (Fig. 3b). A mortality rate of 6% for plants treated with MM was observed at 7 DAT, which increased to 19% at 11 DAT, and peaked at 63% by 15 DAT.
Comparably, treatment with COT-MM induced a 6% mortality rate from 11 DAT, but this remained constant until the end of the experiment at 20 DAT. Control-treated plants displayed no mortality through the duration of the experiment (Fig. 3b).
Shoot and root dry biomass (SDB, RDB) were equally and significantly affected by MM and COT-MM at 20 DAT compared to the control. MM and COT-MM reduced SDB by 82% and 71%, respectively, and RDB was decreased by 94% and 89%, respectively (Fig. 3c, Table S2 and Fig. S2, ESI†).
The experimental outcomes reported above support the view that COT-MM is a herbistatic agent, as the inhibition to plant growth was no longer observed to be comorbid with plant death, as is typical of most commercial herbicides. Applying MM and COT-MM at concentrations ranging from 0.01 mM to 0.40 mM showed comparable behaviour in respect to both plant vertical growth and RDB, with minor variation to SDB as a function of concentration (Table S3 and Fig. S3, ESI†). Minimal effects to rubber vine plants grown in vitro from BT and COT-BT indicates that the inclusion of the COT functional group did not have an influence on plant singlet oxygen or plant P450 oxidation pathways.
A longer-term glasshouse study was performed with COT-MM to determine the extent of herbistasis, and whether it was reversible. Rubber vine mortality was observed to peak at 25% at 28 DAT and remained constant until the end of the experiment at 70 DAT (Fig. 3d and e). Similarly to the 20 day glasshouse study, no plant death was observed within the first 7 DAT, and a comparable 8% mortality rate was observed at 14 DAT. While no plant death was observed after 28 DAT, the 18 surviving plants exhibited stunted vertical growth up to 42 DAT, with the first noticeable increase in height measured at 49 DAT. By 56 DAT plant height had increased by 57% compared to the height at the time of treatment (0 DAT), rising to 87% at 63 DAT, and reaching 119% at the end of the experiment at 70 DAT (Fig. 3d). These results provide additional support that COT-MM is a herbistatic agent, as rubber vine plants appear to recover after sustaining more than six weeks of inhibited growth (Fig. S4, ESI†).
To test whether COT-MM has the same MOA as its parent template MM, a detailed kinetic study was performed (Fig. 4). The ALS from Arabidopsis thaliana (AtALS) was chosen as a model enzyme on which to perform inhibition studies. Its sequence identity and sequence homology with rubber vine ALS are 82.1% and 96.4%, respectively. The amino acids in the herbicide binding site are also completely conserved (Fig. S5, ESI†). This therefore represents an excellent interpretative model for analysis.
Fig. 4 Inhibition plots of MM and COT-MM for AtALS. Ki measurement of (a) MM fitted with the tight binding equation, and (b) COT-MM fitted with the simple inhibition equation. The time-dependent accumulative inhibition curves at 25 °C for (c) MM and (d) COT-MM, under normal oxygen aerobic conditions (teal and blue respectively), and nitrogen anaerobic conditions (black). The individual replicates are fitted with grey dotted lines, and the averaged fitted curves (bolded lines) use the respective F values ([E]F/[EI]) and eqn (S3) (ESI†) to determine the average values of kiapp and k3. All measurements were performed in triplicate with error bars representing the standard deviations (SD) of the data. |
The inhibition assay was performed after nitrogen purging of the master solution to minimize any effects of the side-reaction involving molecular oxygen. The continuous assay method was applied to monitor the disappearance of pyruvate directly and to determine when the lag phase of the reaction had been completed. Initial testing showed that MM was an extremely tight binding inhibitor while COT-MM was a weak binding inhibitor, and therefore required a different range of inhibitor concentrations between the two experiments (Fig. 4). Applying the tight binding equation (eqn (S2), ESI†) a Ki value of 47 nM was obtained for MM (Fig. 4a), which is a similar value to that previously published (80 nM)46 taking into account a discontinuous assay was used for the latter value. A Ki value of 19 mM was obtained for COT-MM (Fig. 4b), but since tight binding was not observed, the equation for non-competitive inhibition was used (eqn (S1), ESI†).25 No inhibition of AtALS was observed for BT or COT-BT, as was expected since the target of BT is reported to be long-chain fatty acid biosynthesis.39–41
Next, the accumulative inhibition of AtALS by MM and COT-MM was measured under standard aerobic conditions as well as after nitrogen purging to investigate the influence of dissolved oxygen on the inhibition kinetics. The key values for the assessment of accumulative inhibition are the apparent rate of inactivation in the presence of the inhibitor (kiapp), and the enzyme recovery rate that occurs in the absence of inhibitor (k3) ((eqn (S3) and Fig. S6, ESI†).25 These values are dependent on the ratio of free enzyme to inhibitor–enzyme complex, so were kept similar in all experiments by determining the F value from the Ki values and eqn (S4) (see ESI† for details).
Under aerobic assay conditions, the kiapp, k3, and kiapp/k3 values for MM were 0.25 min−1, 0.005 min−1, and 50 respectively, and for COT-MM they were 0.10 min−1, 0.003 min−1, and 31, respectively (Fig. 4c, d and Table S4, ESI†). The higher kiapp/k3 for MM suggests a slower recovery rate, and therefore an increase in accumulative inhibition. The decreased kiapp value for COT-MM compared to MM was initially suspected to arise from the reactive oxygen quenching abilities of COT (Scheme 2a) protecting ThDP and FADH2 from oxidative inactivation. However, no oxygenated COT-MM by-products were detected by mass spectrometry (see below), hence it is likely that this may result from a modified inhibitor orientation within the active site (see below).
To further explore the importance of oxygen towards ALS enzyme inhibition, the master solution was purged with nitrogen and supplemented with the known oxygen quenching reagent, 2-mercaptoethanol, to simulate anaerobic conditions. No time-dependent inhibition could be measured for COT-MM or MM under deoxygenated conditions (k3 = 0), and both displayed a 5-fold decrease of kiapp to 0.06 min−1 and 0.02 min−1, respectively (Fig. 4c and d). The maximum enzyme velocity (Vmax) was constant across all samples (Table S4, ESI†). The apparent non-reversible inhibition under deoxygenated conditions is tentatively suggested to result from degradation of ThDP through the physical aspects of competitive/non-competitive inhibition when MM and COT-MM occupy the enzyme active site. This contribution to ALS inhibition is likely still present under oxygenated conditions, although only a minor component compared to the dominating accumulative pathway.
In a subsequent experiment the reaction temperature was increased from 25 °C to 30 °C to investigate the relative inhibition rates with elevated Vmax at temperatures closer to those of glasshouse plant experiments (Table S4 and Fig S7, ESI†). The five degree increase in temperature approximately doubled the Vmax, and the kiapp under aerobic conditions was also doubled to 0.57 min−1 for MM, and to 0.23 min−1 for COT-MM. The enzyme recovery rate in the absence of inhibitor again appeared equivalent under oxygenated conditions between COT-MM and MM (k3 = 0.02 min−1), and an order of magnitude higher than that at 25 °C. It is important to note that this recovery rate represents holoenzyme and [ES*] → [ES#] (after [ES*I] → [ES*]), and not the original enzyme without cofactors (Fig. S6, ESI†). As the experimental methods and consistent F values ensure that the observed rate constants were derived at equivalent [ESI] concentrations, inhibitors with the same binding inhibition mechanism should also demonstrate comparable k3 values. Therefore, COT-MM and MM are expected to act via the same mechanism, as would be expected due to their similar chemical composition and that ALS herbicides do not typically induce structural changes to the enzyme active site as a consequence of herbicide inhibition.
While the addition of 2-mercaptoethanol was not found to influence the binding affinity or Ki values of MM or COT-MM (data not shown), it was omitted under nitrogen sparged conditions at 30 °C. As was expected, COT-MM again did not display any accumulative inhibition (k3 = 0), with a kiapp of 0.02 min−1. Notably, this kiapp value is the same as the 25 °C-condition using 2-mercaptoethanol, supporting the theory that ThDP degradation via physical inhibition is dominating under anaerobic conditions, with minimal contribution from reactive oxygen pathways (Table S4 and Fig S7, ESI†). Any oxygen quenching from the COT ring in the absence of 2-mercaptoethanol is assumed to be minimal. A small amount of accumulative inhibition was observed for MM in the nitrogen sparged condition (kiapp/k3 = 5), which may be due to the fact that MM is an extremely tight binder, and that residual oxygen is likely remaining in the system. The AtALS binding affinity of MM (47 nM) is almost 500 times stronger than COT-MM (19 mM), so it is anticipated that this variation to binding affinity would also modify the ability to trap reactive oxygen species.
The lower Ki and lower kiapp values of COT-MM across the kinetic experiments indicates a decreased active site presence in AtALS compared to MM. Nevertheless, plant experiments using equal concentrations of COT-MM and MM (i.e. 0.05, 0.20 and 0.40 mM) showed comparable inhibition to plant growth (Fig. S3, ESI†). This highlights that oxidative inactivation of this enzyme is a potent mode of action, which could dominate over the herbicide binding affinity and the apparent rate of enzyme inactivation as the leading predictor of herbicide activity. Coincidentally, the reported rate of enzyme activation (kobs(act))47 for WT AtALS is comparable to the kiapp observed for COT-MM under 30 °C aerobic conditions. This serendipitous result corroborates the observed herbistatic behaviour under glasshouse conditions, with BCAA production at steady state due to the ALS inactivation rate from COT-MM being equal to its intrinsic reactivation rate. In plants, this could sustain the essential biological functions by constant BCAA production, but prevent further plant growth due to restricted output. It is assumed that this inhibition mechanism would persist until the inhibitor is metabolised or degraded, or until secondary processes result in plant death. The kiapp for MM at both 25 °C and 30 °C is higher than the enzyme reactivation rate (i.e., kiapp> kobs(act)), so the enzyme equilibrium moves towards the inactivated “off” state, yielding accumulative inhibition and increased plant death.
Fig. 6 Binding mode of MM and COT-MM to AtALS. (a) Crystal structure of MM (cyan, PDB:1YHY) in complex with AtALS (grey) and cofactors (green), overlayed with the docking pose of COT-MM (magenta). (b) Docking pose of COT-MM with ligand interactions to residues (pale blue), by hydrogen bonds (yellow dashes), π–π stacking (green dashes), salt bridge (orange dashes), and hydrophobic interactions (black dashes). Nitrogen atoms in dark blue, sulfur in yellow, and oxygen in red. |
However, the higher kiapp/k3 ratio for MM suggests increased accumulative inhibition, likely due to its decreased size and compact bent configuration (Fig. 6a). Therefore, it was considered that the enzyme—MM complex can more readily trap the reactive ThDP-carbanion intermediate without access to a second acceptor keto-acid that is required for product formation. Analogous to the previously proposed “waiting room” theory,51 where both donor and acceptor keto-acids are bound simultaneously in the initial stages of the reaction. Closure of the active site after MM binds prevents further exchange with solvent or substrate, leaving the ThDP-carbanion exposed and susceptible to degradation or side-reactions.22 Consequently, oxygenation within the enzyme—MM complex can cause the formation of ThDP-peracetate, which then oxidises FADH2 to the inactive FAD form, with re-reduction by pyruvate oxidase for enzyme activation being the rate limiting step (kobs(act) = 0.2 min−1 for AtALS).47 A comparable process is postulated to occur for the enzyme–COT-MM complex, however its extended structure may allow the active site to stay open, enabling the migration of oxygen or keto-acids and decreasing the possibility of ThDP-peracetate formation. As reflected by the increased Ki, the loss of stabilizing phenyl ring interactions may also prevent the ThDP-peracetate side-product from being constrained in the correct orientation for reaction with FADH2. While the lower kiapp values of COT-MM compared to MM demonstrates less efficient enzyme inactivation, the exact mechanism behind accumulative inhibition is yet to be confirmed. However, both the molecular docking and the time-dependent enzyme kinetic studies strongly suggest the target for the herbistasis activity of COT-MM is through the inhibition of ALS.
This observation is remarkable in that the ALS inhibition MOA is the same as MM as confirmed by AtALS kinetics, and supported by molecular docking, but the mechanism of action has switched from herbicidal to herbistatic for COT-MM. Mass spectroscopy, and the COT moiety control agent, COT-BT, suggests the herbistasis effect is not the result of oxidative modification of COT-MM due to plant or ALS reactive oxygen species. Molecular docking reveals that the COT-MM adopts a fully extended structure compared to that for the other sulfonylureas and MM, with the COT moiety in the “tub-shaped” conformation (D2d) (Fig. 2) expanding beyond the binding access channel. This supports why COT-MM does not bind as tightly to AtALS as MM, enabling it to move in and out of the binding channel more freely. This would result in a decrease in BCAA synthesis, rather than the potent decrease which occurs for tight binding inhibitors such as MM.
Inducing plant herbistasis is rare, and even more so for synthetic chemicals. As explored herein, the herbistatic features displayed by COT-MM were realised by targeting the well-established ALS MOA and applying a bioisostere/biomotif chemical modification protocol using a pre-existing commercial herbicide template. Deploying these methods to reinvigorate the concept of herbistasis could provide additional tools for tackling the weed crisis. For example, (1) weeds may need to evolve alternate mechanisms of resistance, which would potentially prolong resistance development; or (2) the low-mortality and herbistasis induced by COT-MM could be employed in combination with weed biocontrol agents at crucial times of the weed lifecycle to suppress plant growth and seed set while enabling the agent population to persist on the weed in the environment. Lastly, potential applications outside of weed control can also be envisaged, such as placing crop seedlings into stasis to avoid the impact of impending deleterious weather like drought.
Footnotes |
† Electronic supplementary information (ESI) available: Experimental procedures, copies of NMR spectra. See DOI: https://doi.org/10.1039/d4sc06923d |
‡ Joint first authors. |
This journal is © The Royal Society of Chemistry 2025 |