Based on the Virtual Screening of Multiple Pharmacophores, Docking and Molecular Dynamics Simulation Approaches toward the Discovery of Novel HPPD Inhibitors

4-Hydroxyphenylpyruvate Dioxygenase Inhibitors: From Molecular Design to Synthesis

Weed resistance is a critical issue in crop production. Among the known herbicides, 4-hydroxyphenylpyruvate dioxygenase (HPPD) inhibitors are crucial for addressing weed resistance. HPPD inhibitors constitute a pivotal aspect of contemporary crop protection strategies. The advantages of these herbicides are their broad weed spectrum, flexible application, and excellent compatibility with other herbicides. They also exhibit satisfactory crop selectivity and low toxicity and are environmentally friendly. An increasing number of new HPPD inhibitors have been designed by combining computer-aided drug design with conventional design approaches. Herein, the molecular design and structural features of innovative HPPD inhibitors are reviewed to guide the development of new HPPD inhibitors possessing an enhanced biological efficacy.

Integrated Virtual Screening and Validation toward Potential HPPD Inhibition Herbicide

4-Hydroxyphenylpyruvate dioxygenase (EC 1.13.11.27, HPPD) is one of the most widely studied herbicide targets and has gained significant attention. To identify potential effective HPPD inhibitors, a rational multistep virtual screening workflow was built, which included CBP models (based on the receptor-ligand interactions in the crystal complex), Hypogen models with activity prediction ability (according to the derivation of structure-activity relationships from a set of molecules with reported activity values), and a consensus docking procedure (consisting of LibDock, Glide, and CDOCKER). About 1 million molecules containing diketone or β-keto-enol substructures were filtered by Lipinski's rules, CBP model, and Hypogen model. A total of 12 compounds with similar docking postures were generated by consensus docking. Eventually, four molecules were screened based on the specific binding pattern and affinity of the HPPD inhibitor. The biological evaluation in vivo displayed that compounds III-1 and III-2 exhibited comparable herbicidal activity to isoxaflutole and possessed superior safety on various crops (wheat, rice, sorghum, and maize). The ADMET prediction (absorption, distribution, metabolism, excretion, and toxicity) showed that compound III possessed relatively good toxicological results. This work provides a theoretical basis and valuable reference for the virtual screening and molecular design of novel HPPD inhibition herbicides.

In-silico bioprospecting of secondary metabolites from endophytic Streptomyces spp. against Magnaporthe oryzae, a cereal killer fungus

Rice blast disease, caused by Magnaporthe oryzae, is the most devastating cereal killer worldwide. Note that melanin pigment is an essential factor of M. oryzae virulence, thus fungicides interfering with melanin biosynthesizing enzymes would reduce the pathogenicity. Scytalone dehydratase (SDH) is the key target for commercial fungicides, like carpropamid, due to its role in the dehydration reaction of the fungal melanin pathway. However, a single-point mutation (V75M) in SDH elicits resistance to carpropamid. A lack of effective fungicides against this resistant strain expedited the quest for novel bioactive inhibitors. Currently, bacterial endophytes like Streptomyces have been heralded for synthesizing bioactive metabolites to protect plants from phytopathogens. The literature search led to the identification of 21 Streptomyces spp. symbionts of paddy that can suppress M. oryzae growth. An antiSMASH server was used to explore Streptomyces spp. gene clusters and found 4463 putative metabolites. Besides, 745 unique metabolites were subjected to a series of virtual screening techniques. Ideally, this process identified five potential SDH inhibitors. The docking result highlights that the metabolite pseudopyronine A interacted hydrophobically with both Val75 of SDHWT and Met75 of SDHV75M targets. Moreover, pseudopyronine A has a higher binding free energy with SDHWT (- 89.94 kcal/mol) and SDHV75M (- 71.95 kcal/mol). Interestingly, the pyranones scaffold of pseudopyronine A was reported for antifungal activity against phytopathogens. Dynamic behavior confirms that pseudopyronine A has excellent conformational states with both SDHWT and SDHV75M. Altogether, we hope that this study creates a new avenue for the discovery of novel phytopathogen inhibitors from endophytes.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03859-7.

Finding potential inhibitors for Main protease (Mpro) of SARS-CoV-2 through virtual screening and MD simulation studies

SARS-CoV-2 is a highly hazardous species that can infect people with Covid-19 disease, dramatically increasing mortality rates worldwide. Plenty of researches have been done to find drugs or inhibitors, with this study aiming to identify an inhibitor within the ChEMBL database using computational approaches. From the ChEMBL library, 19,43,048 compounds which are known type of small compounds and proteins were downloaded and docked with the Main protease (Mpro). After performing compound screening using Lipinski's rule, Qikprop analysis following with virtual Screening, Induced Fit Docking (IFD) and MM-GBSA analysis with the Glide and Prime modules of Schrödinger, the best complex was subjected to MD simulation with Desmond. According to the docking results, small protein 2,371,668 and compound 1,090,395 were docked with Main protease with -12.6, -12.0 kcal/mol dock score and interacted with the functional site residues His 41 and Cys 145, as well as the binding site residues Thr 26, Phe 140, Asn 142, Gly 143, Glu 166, and Gln 189. Complex structures were shown to be steadier by the MD simulation study because both the ligands heavy atoms and the protein Cα atoms' RMSD values fell within acceptable ranges. As a result, this research suggests that the molecule CHEMBL2371668 and the compound CHEMBL1090395 may inhibit the activity of Main protease, and the usefulness of these molecules will be examined further through experimental research.

A Protocol for Activated Bioorthogonal Fluorescence Labeling and Imaging of 4-Hydroxyphenylpyruvate Dioxygenase in Plants

4-Hydroxyphenylpyruvate dioxygenase (HPPD) plays a crucial role in the synthesis of nutrients needed to maintain optimal plant growth. Its level is closely linked to the extent of abiotic stress experienced by plants. Moreover, it is also the target of commercial herbicides. Therefore, labeling of HPPD in plants not only enables visualization of its tissue distribution and cellular uptake, it also facilitates assessment of abiotic stress of plants and provides information needed for the development of effective environmentally friendly herbicides. In this study, we created a method for fluorescence labeling of HPPD that avoids interference with the normal growth of plants. In this strategy, a perylene-linked dibenzyl-cyclooctyne undergoes strain-promoted azide-alkyne cycloaddition with an azide-containing HPPD ligand. The activation-based labeling process results in a significant emission enhancement caused by the change in the fluorescent forms from an excimer to a monomer. Notably, this activated bioorthogonal strategy is applicable to visualizing HPPD in Arabidopsis thaliana, and assessing its response to multiple abiotic stresses. Also, it can be employed to monitor in vivo levels and locations of HPPD in crops. Consequently, the labeling strategy will be a significant tool in investigations of HPPD-related abiotic stress mechanisms, discovering novel herbicides, and uncovering unknown biological functions.

Pharmacophore based virtual screening for identification of effective inhibitors to combat HPV 16 E6 driven cervical cancer

Targeting HPV16 E6 has emerged as an effective drug target for the treatment/management of cervical cancer. We utilized pharmacophore-based virtual screening, molecular docking, absorption, distribution, metabolism and excretion (ADME) prediction, and molecular dynamics simulation approach for identifying potential inhibitors of HPV16 E6. Initially, we generated a ligand-based pharmacophore model based on the features of four known HPV16 E6 inhibitors (CA24, CA25, CA26, and CA27) via the PHASE module implanted in the Schrödinger suite. We constructed four-point pharmacophore features viz., three hydrogen bond acceptors (A) and one aromatic ring (R). The common pharmacophore feature further employed as a query for virtual screening against the ASINEX database via Schrödinger suite. The pharmacophore-based virtual screening filtered out top 2000 hits, based on the fitness score. We then applied the high throughput virtual screening (HTVS), standard precision (SP) and extra precision (XP). 1000 compounds were obtained from HTVS docking. Based on the glide score, they were further filtered to 500 hits by employing docking in standard precision mode. Finally, the best four hits and a negative molecule were identified using docking in XP mode. The four lead compounds and a negative molecule were then further subjected to ADME profile prediction by engaging Qikprop module. The ADME properties of the four lead molecules indicate good pharmacokinetic (PK) properties rather than the negative molecule. The binding stability of the HPV16 E6-hit complexes were investigated at a different time scale (100 ns) by using the desmond package and the results were examined using Root Mean Square Deviation (RMSD) and Root Mean Square Fluctuation (RMSF) and it revealed the stability of the protein-ligand complex throughout the simulation. Key residues, CYS 51 and GLN 107, also play a crucial role in enhancing the stability of the protein-ligand complex during the simulation. Furthermore, the binding free energy of the HPV16 E6-leads complexes was analyzed by prime which revealed that the ΔGbind coulomb and ΔGbind vdW interactions are crucially contributes to the binding affinity. In order to validate the computational findings, the efficacy of benzoimidazole and benzotriazole were ascertained for regulating ME180 cervical cancer cell survival, migration and ability to release MMP-2.

Discovery of novel HPPD inhibitors: Virtual screening, molecular design, structure modification and biological evaluation

4-Hydroxyphenylpyruvate dioxygenase (EC 1.13.11.27, HPPD, a Fe(II)/α-ketoglutarate dependent oxygenases), is a popular herbicide target. In this work, two pharmacophore models based on common molecular characteristics (HipHop) and receptor-ligand complex (CBP) were generated for virtual screening for HPPD inhibitors. About 1,000,000 molecules containing diketone structure from PubChem were filtered by Lipinski's rules to build a 3D database. Then the database was screened through combining HipHop model, CBP model, ADMET (absorption, distribution, metabolism, excretion and toxicity) prediction and molecular docking. Subsequently, based on the specific binding mode and affinity of HPPD inhibitors, 4 molecules with high -CDOCKER energy, good aqueous solubility and human safety predicative properties values were screened. From the screening results and combined with previous work, three novel HPPD inhibitors were designed and synthesized through fragment splicing and bioisosterism strategies. Compound IV-a exhibited similar inhibition of Arabidopsis thaliana HPPD (AtHPPD) and herbicidal activity as mesotrione. Crop selectivity showed that compound IV-a had better crop safety than mesotrione. Comparing the molecular properties, ADMET and molecular docking studies indicated that compounds IV-a exhibited better properties than mesotrione, which could be further modified as novel HPPD inhibitor herbicides.

Computational investigation of potent inhibitors against SARS-CoV-2 2'-O-methyltransferase (nsp16): Structure-based pharmacophore modeling, molecular docking, molecular dynamics simulations and binding free energy calculations

The Coronavirus Disease 2019 (COVID-19), caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, has created unprecedented public health and economic crises around the world. SARS-CoV-2 2'-O-methyltransferase (nsp16) adds a "cap" to viral RNA to maintain the stability of viral RNA, and inhibition of nsp16 activity may reduce viral proliferation, making this protein an attractive drug target. Here, we report the identification of several small molecule inhibitors of nsp16 by virtual screening. First, the nsp16-sinefungin complex (PDB ID: 6WKQ) was selected from the protein data bank. Asp6912, Cys6913, Asp6897 and Asp6928 were determined to be the key amino acids for sinefungin binding in the crystal structure of nsp16-sinefungin complex by molecular dynamics simulation. The complex structures in the stable binding trajectory of nsp16-sinefungin were than clustered through molecular dynamics RMSD analysis. Six clusters were generated, and six representative structures were selected to construct the pharmacophore based on the structure. These six pharmacophores were superimposed on the binding pocket to simplify and pick the common characteristics. The compounds obtained by the pharmacophore screening from Bionet and Chembiv databases were docked into the nsp16 active pocket. The candidate compounds were selected according to the molecular docking score and then screened by MM/GBSA. Finally, four candidate compounds were obtained. Four sets of 150ns molecular dynamics simulations were performed to determine whether candidate compounds could maintain stable interactions with key amino acids. The results of MD and MM/PBSA energy decomposition indicated that C1 and C2 could form a stable complex system with nsp16, and could form strong hydrogen bonds and salt bridges with the key amino acid Asp6897 and Asp6928. This study thus identifies and attempts to validate for the first time the potential inhibitory activities of C1 and C2 against nsp16, allowing the development of potent anti-COVID-19 drugs and unique treatment strategies.

Fludarabine, a Potential DNA-Dependent RNA Polymerase Inhibitor, as a Prospective Drug against Monkeypox Virus: A Computational Approach

Monkeypox is a zoonotic contagious disease that has recently re-emerged in different countries worldwide. Due to the lack of an effective treatment that eliminates the virus, there is an urgent need to find effective drugs to stop the spread of the multi-country outbreak. The current study aimed to use computational methods to quickly identify potentially effective drugs against the Monkeypox virus (MPXV). Three MPXV proteins were targeted in this study due to their essential role in viral replication (a DNA-Dependent RNA Polymerase subunit (A6R)), a protein involved in cell entry (D8L), and a protein catalyzing the envelopment of intracellular mature virus particles (F13L). We virtually screened a library of 1615 FDA-approved compounds, utilizing different in-silico approaches including computational modeling, molecular docking, molecular dynamic (MD) simulation, and MM-GBSA. The compound Fludarabine was found to have the best docking score (−7.53 kcal/mol) in relation to the MPXV A6R protein. Additionally, Fludarabine showed in-silico activity on the D8L and F13L proteins. During the whole period of the 100 ns MD simulation, the complex of A6R and Fludarabine exhibited the best stability. This stability was reflected in a good score of MM-GBSA, with an average value of −44.62 kcal/mole in a range between −53.26 and −35.49 and a low value of standard deviation (3.76). Furthermore, Fludarabine blocked efficiently the Asn175 residue which has an important role in the attachment of the virus to a host cell. The results of this study recommend more in vitro studies on this compound, as a starting point to develop a novel treatment against MPXV.

Computer-Assisted Proposition of Promising Aryloxyacetic Acid Derivatives as HPPD Inhibitors

A series of aryloxyacetic acid derivatives have demonstrated promising herbicidal performance by inhibition of the hydroxyphenylpyruvate deoxygenase (HPPD) enzyme. We hereby applied quantitative structure-activity relationship (QSAR) and docking strategies to model and chemically understand the bioactivities of these compounds and subsequently propose unprecedented analogues aiming at improving the herbicidal and environmental properties. Bulky halogens at the 2-, 3-, 4-, and 6-positions of an aromatic ring, CF₃ in 4-position, and the 2-NO₂ group in a phenyl ring appear to favor the HPPD inhibition. At the same time, Me and OMe substituents contribute to decreasing the pK_i values. Accordingly, a few compounds were proposed and the candidate with 2,4,6-triBr substituents demonstrated an estimated pK_i similar to those of the best library compounds. This finding was corroborated by the docking scores of the ligand-enzyme interactions. In addition, the high calculated lipophilicity of some proposed agrochemicals suggests that they should have low soil mobility and, therefore, are not prone to easily leach out and reach groundwater, despite causing other ecological issues.

Survey on the Recent Advances in 4-Hydroxyphenylpyruvate Dioxygenase (HPPD) Inhibition by Diketone and Triketone Derivatives and Congeneric Compounds: Structural Analysis of HPPD/Inhibitor Complexes and Structure-Activity Relationship Considerations

The serendipitous discovery of the HPPD inhibitors from allelopathic plants opened the way for searching new and effective herbicidal agents by application of classical hit-to-lead optimization approaches. A plethora of active and selective compounds were discovered that belong to three major classes of cyclohexane-based triketones, pyrazole-based diketones, and diketonitriles. In addition, to enhance inhibitory constant and herbicidal activity, many efforts were also made to gain broader weed control, crop safety, and eventual agricultural applicability. Moreover, HPPD inhibitors emerged as therapeutic agents for inherited and metabolic human diseases as well as vector-selective insecticides in the control of hematophagous arthropods. Given the large set of experimental data available, structure-activity relationship analysis could be used to derive suggestions for next generation optimized compounds.

Virtual screening based on pharmacophore model for developing novel HPPD inhibitors

4-Hydroxyphenylpyruvate dioxygenase (HPPD) is an important target for herbicide design. A multilayered virtual screening workflow was constructed by combining two pharmacophore models based on ligand and crystal complexes, molecular docking, molecular dynamics (MD), and biological activity determination to identify novel small-molecule inhibitors of HPPD. About 110, 000 compounds of Bailingwei and traditional Chinese medicine databases were screened. Of these, 333 were analyzed through docking experiments. Five compounds were selected by analyzing the binding pattern of inhibitors with amino acid residues in the active pocket. All five compounds could produce stable coordination with cobalt ion, and form favorable π-π interactions. MD simulation demonstrated that Phe381 and Phe424 made large contributions to the strength of binding. The enzyme activity experiment verified that compound-139 displayed excellent potency against AtHPPD (IC₅₀ = 0.742 μM), however, compound-5222 had inhibitory effect on human HPPD (IC₅₀ = 6 nM). Compound-139 exhibited herbicidal activity to some extent on different gramineous weeds. This work provided a strong insight into the design and development of novel HPPD inhibitor using in silico techniques.

Design, synthesis and biological activity of novel triketone-containing quinoxaline as HPPD inhibitor

BACKGROUND

4-Hydroxyphenyl pyruvate dioxygenase (EC 1.13.11.27, HPPD) is one of the important target enzymes used to address the issue of weed control. HPPD-inhibiting herbicides can reduce the carotenoid content in plants and hinder photosynthesis, eventually causing albinism and death. Exploring novel HPPD-inhibiting herbicides is a significant direction in pesticide research. In the process of exploring new high-efficiency HPPD inhibitors, a series of novel quinoxaline derivatives were designed and synthesized using an active fragment splicing strategy.

RESULTS

The title compounds were unambiguously characterized by infrared, ¹ H NMR, ¹³ C NMR, and high-resolution mass spectroscopy. The results of the in vitro tests indicated that the majority of the title compounds showed potent inhibition of Arabidopsis thaliana HPPD (AtHPPD). Preliminary bioevaluation results revealed that a number of novel compounds displayed better or excellent herbicidal activity against broadleaf and monocotyledonous weeds. Compound III-5 showed herbicidal effects comparable to those of mesotrione at a rate of 150 g of active ingredient (ai)/ha for post-emergence application. The results of molecular dynamics verified that compound III-5 had a more stable protein-binding ability. Molecular docking results showed that compound III-5 and mesotrione shared homologous interplay with the surrounding residues. In addition, the enlarged aromatic ring system adds more force, and the hydrogen bond formed can enhance the synergy with π-π stacking.

CONCLUSIONS

The present work indicates that compound III-5 may be a potential lead structure for the development of new HPPD inhibitors.

DCLK1 and its interaction partners: An effective therapeutic target for colorectal cancer

Doublecortin-like kinase protein 1 (DCLK1) is a microtubule-associated protein with a C-terminal serine/threonine kinase domain. Its expression was first reported in radial glial cells, where it serves an essential role in early neurogenesis, and since then, other functions of the DCLK1 protein have also been identified. Initially considered to be a marker of quiescent gastrointestinal and pancreatic stem cells, DCLK1 has recently been identified in the gastrointestinal tract as a marker of tuft cells. It has also been implicated in different types of cancer, where it regulates several vital pathways, such as Kras signaling. However, its underlying molecular mechanisms remain unclear. The present review discusses the different roles of DCLK1 and its interactions with other proteins that are homologically similar to DCLK1 to develop a novel therapeutic strategy to target cancer cells more accurately.

Phenoxypyridine derivatives containing natural product coumarins with allelopathy as novel and promising proporphyrin IX oxidase-inhibiting herbicides: Design, synthesis and biological activity study

To seek novel and safe protoporphyrinogen oxidase (PPO, EC 1.3.3.4) inhibitors with excellent herbicidal activity. A series of novel phenoxypyridine derivatives containing natural product coumarins with allelopathy were designed and synthesized based on bioisosterism and active subunit combination in this research. Compounds W3.1 and W3.4, with the half-maximal inhibitory concentration (IC₅₀) value of 0.02653 mg/L and 0.01937 mg/L, respectively, displayed excellent herbicidal activity in greenhouse. Their herbicidal activity was similar to commercial herbicide oxyfluorfen (IC₅₀ = 0.04943 mg/L). The best field inhibitory effect of compounds W3.1 and W3.4 recorded was at doses of 450 g ai/ha and 300 g ai/ha, respectively. Compound W3.4 had the best herbicidal activity among all the target compounds in this paper. Molecular docking analysis revealed that compounds W3.1 and W3.4 could form a hydrogen bonds with the amino acid AGR-98 and a π-π superposition with the amino acid PHE-398, respectively, which was similar to the oxyfluorfen. The crop selectivity tests results indicated that maize, cotton and soybean showed high tolerance to compound W3.4. Compound W3.4 reduced the C_a and C_b contents of wheat and rice, but had less effect on maize, cotton and soybean. Selectivity of compound W3.4 in maize, cotton and soybean were appeared to be due to reduced absorption of the herbicide compared to wheat and rice. Compound W3.4 deserves further attention as a candidate structure for new herbicides.

Potential Inhibitors of Fascin From A Database of Marine Natural Products: A Virtual Screening and Molecular Dynamics Study

Marine nature products are unique compounds that are produced by the marine environment including plants, animals, and microorganisms. The wide diversity of marine natural products have great potential and are versatile in terms of drug discovery. In this paper, we use state-of-the-art computational methods to discover inhibitors from marine natural products to block the function of Fascin, an overexpressed protein in various cancers. First, virtual screening (pharmacophore model and molecular docking) was carried out based on a marine natural products database (12015 molecules) and provided eighteen molecules that could potentially inhibit the function of Fascin. Next, molecular mechanics generalized Born surface area (MM/GBSA) calculations were conducted and indicated that four molecules have higher binding affinities than the inhibitor NP-G2-029, which was validated experimentally. ADMET analyses of pharmacokinetics demonstrated that one of the four molecules does not match the criterion. Finally, ligand Gaussian accelerated molecular dynamics (LiGaMD) simulations were carried out to validate the three inhibitors binding to Fascin stably. In addition, dynamic interactions between protein and ligands were analyzed systematically. Our study will accelerate the development of the cancer drugs targeting Fascin.