Development of a general quantum-chemical descriptor for steric effects: density functional theory based QSAR study of herbicidal sulfonylurea analogues

Design, synthesis and insecticidal activity of novel Isoxazoline Acylhydrazone compounds

BACKGROUND

Nowadays, the diamondback moth has ascended to become one of the most formidable pests plaguing cruciferous vegetables. Consequently, the exigency for the development of efficacious pesticide candidates for crop protection has never been more paramount. In response to this pressing need, this study presents a compendium of novel isoxazoline derivatives, incorporating acylhydrazone moieties, synthesized with the express purpose of serving as potential insecticides.

RESULTS

The structures of these derivatives were confirmed using Proton nuclear magnetic resonance (1 H NMR), Carbon-13 nuclear magnetic resonance (13 C NMR), and high-resolution mass spectrometry (HR-MS). Most of these derivatives demonstrated effective insecticidal activities against Plutella xylostella. Notably, compound E3 exhibited exceptional insecticidal activity against Plutella xylostella (LC50  = 0.19 mg L-1 ), surpassing the effectiveness of ethiprole (LC50  = 3.28 mg L-1 ), and comparable to that of fluxametamide (LC50  = 0.22 mg L-1 ). Interestingly, compound E3 also displayed potent insecticidal activity against Pyrausta nubilalis (LC50  = 0.182 mg L-1 ) and Chilo suppressalis (LC50  = 0.64 mg L-1 ), and the LC50 values of fluxametamide were 0.23 mg L-1 (P. nubilalis) and 2.26 mg L-1 (C. suppressalis), respectively. The molecular docking results revealed that the compound E3 can form a hydrogen bond and two Pi-Pi bonds with the active sites of GABA receptors. In addition, the DFT calculations were also performed to study the relationship between insecticidal activities. The structure-activity relationships suggested that the identity of the R substituent was crucial for their pesticidal activities.

CONCLUSION

The results of the present study suggest that isoxazoline acylhydrazone derivatives could be promising candidates against P. xylostella and other Lepidopteran pests. © 2023 Society of Chemical Industry.

Novel Sulfonylurea Derivatives as Potential Antimicrobial Agents: Chemical Synthesis, Biological Evaluation, and Computational Study

Methicillin-resistant Staphylococcus aureus (MRSA) is a worldwide health threat and has already tormented humanity during its long history, creating an urgent need for the development of new classes of antibacterial agents. In this study, twenty-one novel sulfonylurea derivatives containing phenyl-5-vinyl and pyrimidinyl-4-aryl moieties were designed and synthesized, among which, nine compounds exhibited inhibitory potencies against Gram-positive bacterial strains: MRSA (Chaoyang clinical isolates), S. aureus ATCC6538, vancomycin-resistant Enterococci-309 (VRE-309), and Bacillus subtilis ATCC 6633. Especially, 9i and 9q demonstrated inhibitory activities against the four bacterial strains with minimum inhibitory concentrations (MICs) of 0.78-1.56 μg/mL, and quite a few of other MRSA clinical strains with MICs of 0.78 μg/mL, superior to those of the positive controls vancomycin (MIC of 1 μg/mL) and methicillin (MIC of >200 μg/mL). This is the very first time that sulfonylurea derivatives have been identified as promising inhibitors against different MRSA clinical isolates. In addition, all the MIC values of the synthesized compounds against Candida albicans were greater than 100 μg/mL. Since the reported anti-Candida activities of sulfonylureas were due to acetohydroxyacid synthase (AHAS) inhibition, the molecular target against MRSA for the target sulfonylureas was thought to be a different mode of action. Density functional theory (DFT) calculations were finally performed to understand the structure-activity relationships, based on which, significant differences were observed between their HOMO maps for compounds with strong antibacterial activities and weak anti-MRSA effects. The present results hence provide valuable guidance for the discovery of novel agents to treat bacterial infections, especially against MRSA.

A Quantitative Structure-Activity Relationship Approach to Determine Biotoxicity of Amide Herbicides for Ecotoxicological Risk Assessment

Amide herbicides have been widely applied in agriculture and found to be widespread and affect nontarget organisms in the environment. To better understand the biotoxicity mechanisms and determine the toxicity to the nontarget organisms for the hazard and risk assessment, five QSAR models were developed for the biotoxicity prediction of amide herbicides toward five aquatic and terrestrial organisms (including algae, daphnia, fish, earthworm and avian species), based on toxicity concentration and quantitative molecular descriptors. The results showed that the developed models complied with OECD principles for QSAR validation and presented excellent performances in predictive ability. In combination, the investigated QSAR relationship led to the toxicity mechanisms that eleven electrical descriptors (E_HOMO, E_LUMO, α_xx, α_yy, α_zz, μ, qN^-, Q_xx, Q_yy, qH⁺, and q^-), four thermodynamic descriptors (Cv, S^θ, H^θ, and ZPVE), and one steric descriptor (V_m) were strongly associated with the biotoxicity of amide herbicides. Electrical descriptors showed the greatest impacts on the toxicity of amide herbicides, followed by thermodynamic and steric descriptors.

Chemical synthesis, biological activities, and molecular simulations of novel sulfonylurea compounds bearing ortho-alkoxy substitutions

A series of 51 novel sulfonylurea compounds with ortho-alkoxy substituent at phenyl ring were chemically synthesized and spectroscopically characterized. The biological activities of the target compounds were evaluated using the enzyme inhibition against acetohydroxyacid synthase (AHAS; EC 2.2.1.6) from fungal or plant source, as well as cell-based antifungal assay and greenhouse pot herbicidal assay. Among the target compounds, 6e showed desirable antifungal activity against Candida albicans standard isolate sc5314 with minimum inhibition concentration (MIC) of 0.39 mg/L (0.98 μM) after 24 h, and 6a demonstrated promising pre-emergence herbicidal activity against Echinochloacrus-galli at 30 g/ha dosage. Representative compounds 6a, 6e, and 6i showed no cell cytotoxicity even at 40 mg/L concentration. Theoretical DFT calculations indicated HOMO maps should be considered to understand the structure-activity relationships. The present study has hence provided useful information for further discovery of novel antifungal agents or selective herbicides.

Discovery of ortho-Alkoxy Substituted Novel Sulfonylurea Compounds That Display Strong Herbicidal Activity against Monocotyledon Grasses

In the present study, we have designed and synthesized a series of 42 novel sulfonylurea compounds with ortho-alkoxy substitutions at the phenyl ring and evaluated their herbicidal activities. Some target compounds showed excellent herbicidal activity against monocotyledon weed species. When applied at 7.5 g ha^-1, 6-11 exhibited more potent herbicidal activity against barnyard grass (Echinochloa crus-galli) and crab grass (Digitaria sanguinalis) than commercial acetohydroxyacid synthase (AHAS; EC 2.2.1.6) inhibitors triasulfuron, penoxsulam, and nicosulfuron at both pre-emergence and postemergence conditions. 6-11 was safe for peanut for postemergence application at this ultralow dosage, suggesting that it could be considered a potential herbicide candidate for peanut fields. Although 6-11 and triasulfuron share similar chemical structures and have close K_i values for plant AHAS, a significant difference has been observed between their LUMO maps from DFT calculations, which might be a possible factor that leads to their different behaviors toward monocotyledon weed species.

QSAR/QSPR models based on quantum chemistry for risk assessment of pesticides according to current European legislation

In Europe, agencies and official organizations involved in the pesticide control such as the EFSA, ECHA, JRC and ECETOC or even the OECD are pointing out that the software tools based on quantitative structure relationship models, i.e. QSAR and QSPR, have a huge potential to improve the pesticide risk assessment process. In this sense, these non-animal test methods can promote the competitiveness of agriculture in this region: the consumer safety is increased with them due to the possibility of perform an overall better risk assessment of the degradation products and metabolites from pesticides. However, the use of theses computational-based (in silico) tools must be much more systematised and harmonised, improving their validation and including case studies to test them. To open databases, incorporating critical data in an orderly manner for building the models, becomes also necessary. Moreover, quantum chemistry through the Density Functional Theory should be promoted as tool for calculation of quantum descriptors, especially for the study of similar compounds with the same carbon skeleton but differing substitution patterns, e.g. isomers.

Fluorogenic Assay for Acetohydroxyacid Synthase: Design and Applications

Acetohydroxyacid synthase (AHAS) exists in plants and many microorganisms (including gut flora) but not in mammals, making it an attractive drug target. Fluorescent-based methods should be practical for high-throughput screening of inhibitors. Herein, we describe the development of the first AHAS fluorogenic assay based on an intramolecular charge transfer (ICT)-based fluorescent probe. The assay is facile, sensitive, and continuous and can be applied toward various AHASs from different species, AHAS mutants, and crude cell lysates. The fluorogenic assay was successfully applied for (1) high-throughput screening of commerical herbicides toward different AHASs for choosing matching herbicides, (2) identification of a Soybean AHAS gene with broad-spectrum herbicide resistance, and (3) identification of selective inhibitors toward intestinal-bacterial AHASs. Among the AHAS inhibitors, an active agent was found for selective inhibition of obesity-associated Ruminococcus torques growth, implying the possibility of AHAS inhibitors for the ultimate goal toward antiobesity therapeutics. The fluorogenic assay opens the door for high-throughput programs in AHAS-related fields, and the design principle might be applied for development of fluorogenic assays of other synthases.

Chemical preparation, biological evaluation and 3D-QSAR of ethoxysulfuron derivatives as novel antifungal agents targeting acetohydroxyacid synthase

Accetohydroxyacid synthase (AHAS) is the first enzyme involved in the biosynthetic pathway of branched-chain amino acids. Earlier gene mutation of Candida albicans in a mouse model suggested that this enzyme is a promising target of antifungals. Recent studies have demonstrated that some commercial AHAS-inhibiting sulfonylurea herbicides exerted desirable antifungal activity. In this study, we have designed and synthesized 68 novel ethoxysulfulron (ES) derivatives and evaluated their inhibition constants (K_i) against C. albicans AHAS and cell based minimum inhibitory concentration (MIC) values. The target compounds 5-1, 5-10, 5-22, 5-31 and 5-37 displayed stronger AHAS inhibitions than ES did. Compound 5-1 had the best K_i of 6.7 nM against fungal AHAS and MIC values of 2.5 mg/L against Candida albicans and Candica parapsilosis after 72 h. A suitable nematode model was established here and the antifungal activity of 5-1 was further evaluated in vivo. A possible binding mode was simulated via molecular docking and a comparative field analysis (CoMFA) model was constructed to understand the structure-activity relationship. The current study has indicated that some ES derivatives should be considered as promising hits to develop antifungal drugs with novel biological target.

Synthesis, characterization, cytotoxic screening, and density functional theory studies of new derivatives of quinazolin-4(3H)-one Schiff bases

A series of novel derivatives of quinazolinone Schiff bases were synthesized from benzoic acid starting material and evaluated for potential cytotoxic activities against the human breast adenocarcinoma (MCF-7) and the human colon adenocarcinoma (HT-29) cell lines. Compared to the reference drug, these compounds showed good cytotoxic activities against studied cell lines especially compounds 4d and 4e. The ground-state geometries of these compounds (4a-g) were optimized at the B3LYP/6–31G* density functional theory (DFT) level. Then maximum absorptions electron affinity, ionization potential, electronegativity (χ), energy gap (E_gap), hardness (η), softness (S), electrophilicity (ω), and electrophilicity index (ω_i) were calculated and discussed. The quantitative structure-activity relationship (QSAR) properties including the physicochemical parameters were also evaluated and studied. The computed properties of our novel synthesized compounds were compared with erlotinib compound.

Triazolopyrimidines as a New Herbicidal Lead for Combating Weed Resistance Associated with Acetohydroxyacid Synthase Mutation

Acetohydroxyacid synthase (AHAS; also known as acetolactate synthase; EC 2.2.1.6, formerly EC 4.1.3.18) is the first common enzyme in the biosynthetic pathway leading to the branched-chain amino acids in plants and a wide range of microorganisms. Weed resistance to AHAS-inhibiting herbicides, increasing at an exponential rate, is becoming a global problem and leading to an urgent demand of developing novel compounds against both resistant and wild AHAS. In the present work, a series of novel 2-aroxyl-1,2,4-triazolopyrimidine derivatives (a total of 55) were designed and synthesized with the aim to discover an antiresistant lead compound. Fortunately, the screening results indicated that many of the newly synthesized compounds showed a better, even excellent, inhibition effect against both the wild-type Arabidopsis thaliana AHAS and P197L mutants. Among them, compounds 5-3 to 5-17, compounds 5-19 to 5-26, compounds 5-28 to 5-45, and compound 5-48 have the lower values of resistance factor (RF) and display a potential power to overcome resistance associated with the P197L mutation in the enzyme levels. Further greenhouse in vivo assay showed that compounds 5-15 and 5-20 displayed "moderate" to "good" herbicidal activity against both the wild type-and the resistant (P197L mutation) Descurainia sophia, even at a rate as low as 0.9375 (g of ai/ha). The above results indicated that these two compounds could be used as new leads for the future development of antiresistance herbicides.

A quantitative structure-activity relationship approach for assessing toxicity of mixture of organic compounds

Four types of reactivity indices were employed to construct quantitative structure-activity relationships for the assessment of toxicity of organic chemical mixtures. Results of analysis indicated that the maximum positive charge of the hydrogen atom and the inverse of the apolar surface area are the most important descriptors for the toxicity of mixture of benzene and its derivatives to Vibrio fischeri. The toxicity of mixture of aromatic compounds to green alga Scenedesmus obliquus is mainly affected by the electron flow and electrostatic interactions. The electron-acceptance chemical potential and the maximum positive charge of the hydrogen atom are found to be the most important descriptors for the joint toxicity of aromatic compounds.

The enantioselectivity in asymmetric ketone hydrogenation catalyzed by RuH2(diphosphine)(diamine) complexes: insights from a 3D-QSSR and DFT study

The 3D-QSSR method was carried out to investigate the enantioselectivity of the asymmetric ketone hydrogenation (AKH) catalyzed by RuH₂(diphosphine)(diamine) complexes integrating with DFT method, which could provide a way to design homogeneous transition-metal catalysts.

Synthesis, biological evaluation and molecular docking studies of some pyrimidine derivatives

Some novel pyrimidine-5-carbonitrile derivatives bearing various substituent have been synthesized. The structures of target compounds were confirmed by elemental analysis and spectral data. Some selected members of the newly synthesized compounds were investigated for their cytotoxic potency against certain human tumor cell lines. Five representative active anticancer compounds 6a, 6c, 6d, 17a and 18a were subjected to docking using MOE program on the 3D structure of two enzymes, namely; thymidylate synthase and dihydrofolate reductase. The antimicrobial activities of the synthesized compounds were tested against Staphylococcus aureus, Pseudomonas aeruginosa, Shigella flexneri and Candida albicans. Compounds 2c, 7a and 9c showed broad spectrum antimicrobial activity.

The minimum activation peptide from ilvH can activate the catalytic subunit of AHAS from different species

Acetohydroxyacid synthases (AHASs), which catalyze the first step in the biosynthesis of branched-chain amino acids, are composed of a catalytic subunit (CSU) and a regulatory subunit (RSU). The CSU harbors the catalytic site, and the RSU is responsible for the activation and feedback regulation of the CSU. Previous results from Chipman and co-workers and our lab have shown that heterologous activation can be achieved among isozymes of Escherichia coli AHAS. It would be interesting to find the minimum peptide of ilvH (the RSU of E. coli AHAS III) that could activate other E. coli CSUs, or even those of ## species. In this paper, C-terminal, N-terminal, and C- and N-terminal truncation mutants of ilvH were constructed. The minimum peptide to activate ilvI (the CSU of E. coli AHAS III) was found to be ΔN 14-ΔC 89. Moreover, this peptide could not only activate its homologous ilvI and heterologous ilvB (CSU of E. coli AHAS I), but also heterologously activate the CSUs of AHAS from Saccharomyces cerevisiae, Arabidopsis thaliana, and Nicotiana plumbaginifolia. However, this peptide totally lost its ability for feedback regulation by valine, thus suggesting different elements for enzymatic activation and feedback regulation. Additionally, the apparent dissociation constant (Kd ) of ΔN 14-ΔC 89 when binding CSUs of different species was found to be 9.3-66.5 μM by using microscale thermophoresis. The ability of this peptide to activate different CSUs does not correlate well with its binding ability (Kd ) to these CSUs, thus implying that key interactions by specific residues is more important than binding ability in promoting enzymatic reactions. The high sequence similarity of the peptide ΔN 14-ΔC 89 to RSUs across species hints that this peptide represents the minimum activation motif in RSU and that it regulates all AHASs.

A perspective on the role of quantitative structure-activity and structure-property relationships in herbicide discovery

BACKGROUND

For the last 15 years the agrochemical industry has focused on using genetic modification to put genes that confer resistance to existing commercial herbicides into crop plants rather than on discovering new herbicides with novel modes of action. The widespread appearance of weeds resistant to those herbicides is now causing the industry to revive their herbicide discovery programs.

RESULTS

Elucidation of quantitative structure-activity relationships (QSARs) played a major role in the discovery and development of existing commercial herbicides, but the advent of genetically modified crops has caused published work (at least) in the area to drift from the industrial arena into academic studies. The focus has also turned inward, to refining models for established herbicide targets instead of elucidating new ones.

CONCLUSION

This perspective highlights the importance of QSARs and quantitative structure-property relationships (QSPRs) to herbicide discovery in an historical context and provides some guidance as to how they might profitably be applied going forward.

Preliminary X-ray crystallographic studies of the catalytic subunit of Escherichia coli AHAS II with its cofactors

Acetohydroxyacid synthase (AHAS) is the first common enzyme in the branched-chain amino-acid biosynthesis pathway and is the target of several classes of commercial herbicides. In this study, the Escherichia coli ilvG gene that encodes the catalytic subunit of AHAS II was cloned into the pET28a vector and expressed in soluble form at high levels in E. coli strain BL21 (DE3) cells. The protein was purified using Ni(2+)-chelating chromatography followed by size-exclusion chromatography. The catalytic subunit of E. coli AHAS II was cocrystallized with its cofactors Mg(2+), FAD and ThDP using the sitting-drop vapour-diffusion method and the crystals diffracted to 2.80 Å resolution.

4,6-Dimethoxy-2-(methylsulfonyl)pyrimidine

The asymmetric unit of the title compound, C₇H₁₀N₂O₄S, comprises of two independent mol­ecules (A and B) which differ in the orientation of the methyl­sulfonyl group [C—S—C—N = 157.98 (13)° in mol­ecule A and 6.09 (18)° in mol­ecule B]. In the crystal structure, mol­ecules of type A are linked into chains along the a axis by inter­molecular C—H⋯O hydrogen bonds. The B mol­ecules are linked to these chains by C—H⋯O hydrogen bonds.

Synthesis and biological activity of 4-(4,6-disubstituted-pyrimidin-2-yloxy)phenoxy acetates

Ten novel 4-(4,6-dimethoxypyrimidin-2-yloxy)phenoxy acetates and 4-(4,6-dimethylpyrimidin-2-yloxy)phenoxy acetates were synthesized with hydroquinone, 2-methylsulfonyl-4,6-disubstituted-pyrimidine and chloroacetic ester as starting materials. The products were characterized by IR, ¹H-NMR, MS spectra and elemental analyses. Preliminary bioassay indicates that the target compounds possess high herbicidal activity against monocotyledonous plants such as Digitaria sanguinalis L. at concentrations of 100 mg/L and 50 mg/L.

Docking and 3D-QSAR studies of acetohydroxy acid synthase inhibitor sulfonylurea derivatives

Docking and three dimensional quantitative-structure activity relationship (3D-QSAR) studies were performed on acetohydroxy acid synthase (AHAS) inhibitor sulfonylurea analogues with potential herbicidal activity. The 3D-QSAR studies were carried out using shape, spatial and electronic descriptors along with a few structural parameters. Genetic function approximation (GFA) was used as the chemometric tool for this analysis. The whole data set (n = 45) was divided into a training set (75% of the data set) and a test set (remaining 25%) on the basis of the K-means clustering technique on a standardised topological, physicochemical and structural descriptor matrix. Models developed from the training set were used to predict the activity of the test set compounds. All models were validated internally, externally and using the Y-randomisation technique. Docking studies suggested that the molecules bind within a pocket of the enzyme formed by some important amino acid residues (Met351, Asp375, Arg377, Gly509, Met570 and Val571). In QSAR studies, molecular shape analysis showed that bulky substitution at the R(1) position may enhance AHAS inhibitory activity. Charged surface area descriptors suggested that negative charge distributed over a large surface area may enhance this activity. The hydrogen bond acceptor parameter supported the charged surface area descriptors and suggested that, for better activity, the number of electronegative atoms present in the molecule should be high. The spatial descriptors show that, for better activity, the molecules should possess a bulky substituent and a small substitution at the R(2) and R(3) positions, respectively.

Fragment-based quantitative structure-activity relationship (FB-QSAR) for fragment-based drug design

In cooperation with the fragment-based design a new drug design method, the so-called "fragment-based quantitative structure-activity relationship" (FB-QSAR) is proposed. The essence of the new method is that the molecular framework in a family of drug candidates are divided into several fragments according to their substitutes being investigated. The bioactivities of molecules are correlated with the physicochemical properties of the molecular fragments through two sets of coefficients in the linear free energy equations. One coefficient set is for the physicochemical properties and the other for the weight factors of the molecular fragments. Meanwhile, an iterative double least square (IDLS) technique is developed to solve the two sets of coefficients in a training data set alternately and iteratively. The IDLS technique is a feedback procedure with machine learning ability. The standard Two-dimensional quantitative structure-activity relationship (2D-QSAR) is a special case, in the FB-QSAR, when the whole molecule is treated as one entity. The FB-QSAR approach can remarkably enhance the predictive power and provide more structural insights into rational drug design. As an example, the FB-QSAR is applied to build a predictive model of neuraminidase inhibitors for drug development against H5N1 influenza virus.

Application of quantum topological molecular similarity descriptors in QSPR study of the O-methylation of substituted phenols

The usefulness of a novel type of electronic descriptors called quantum topological molecular similarity (QTMS) indices for describing the quantitative effects of molecular electronic environments on the O-methylation kinetic of substituted phenols has been investigated. QTMS theory produces for each molecule a matrix of descriptors, containing bond (or structure) information in one dimension and electronic effects in another dimension, instead of other methods producing a vector of descriptors for each molecule. A collection of chemometrics tools including principal component analysis (PCA), partial least squares (PLS), and genetic algorithms (GA) were used to model the structure-kinetic data. PCA separated the bond and descriptor effects, and PLS modeled the effects of these parameters on the rate constant data, and GA selected the most relevant subset of variables. The model performances were validated by both cross-validation and external validation. The results indicated that the proposed models could explain about 95% of variances in the rate constant data. The significant effects of variables on the reaction kinetic were identified by calculating variable important in projection (VIP). It was found that the rate constant of esterification of phenols is highly influenced by the electronic properties of the C2--C1--O--H fragment of the parent molecule. Indeed, the C2--X and C4--X bonds (corresponding to ortho and para substituents) were found as highly influential parameters. All of the eight calculated QTMS indices were found significant however, lambda1, lambda2, lambda3, epsilon, and K(r) were detected as highly influential parameters.