Insight into herbicide resistance of W574L mutant Arabidopsis thaliana acetohydroxyacid synthase: molecular dynamics simulations and binding free energy calculations

Molecular docking and mutation sites of CYP57A1 enzyme with Fomesafen

Fomesafen is a diphenyl ether herbicide developed by Zeneca Group PLC (UK), mainly used in soybean and peanut fields to control annual and perennial broad-leaved weeds. Fomesafen has strong persistence in the soil, slow degradation rate, and is prone to harm subsequent sensitive crops. This study utilized Autodock molecular docking technology to investigate the binding and interaction between degradation enzyme CYP57A1 and small molecules of fomesafen herbicides. The CYP57A1 gene cloned from a fomesafen-resistant fungus Fusarium verticilloids, belongs to a fragment of the P450 family, contains 587 bases, encodes 190 amino acids, and has an isoelectric point of 5.16. Visualization of the active surface of the protein receptor reveals that fomesafen is located in the cavity formed by the CYP57A1 protein and the cavity is small and tightly, the proteins are connected to small molecules through hydrogen bonds, halogen atom and π - cation interactions. Molecular modification of CYP57A1 enzyme was carried out using virtual amino acid mutation technology. Four key amino acids, LEU143, MET52, PHE176, and GLU177, were subjected to site-specific mutations. This study successfully constructed mutant engineered bacteria with stable protein expression. Mutations (1) MET52 > TRP showed a a decrease in enzyme activity, and the degradation rate of fomesafen was only 7.8 % of the wild-type. It is believed that MET52 is a key active site for the binding of CYP57A1 enzyme to small molecules of fomesafen, playing a crucial role in the degradation of fomesafen by this enzyme. This provides new insights into the impact on the degradation activity of fomesafen.

Discovery of New 2-[(4,6-Dimethoxy-1,3,5-triazin-2-yl)oxy]-6-(substituted phenoxy)benzoic Acids as Flexible Inhibitors of Arabidopsis thaliana Acetohydroxyacid Synthase and Its P197L Mutant

In the search for new antiresistance acetohydroxyacid synthase (AHAS, EC 2.2.1.6) inhibitors to combat weed resistance associated with AHAS mutations, a series of 2-[(4,6-dimethoxy-1,3,5-triazin-2-yl)oxy]-6-(substituted phenoxy)benzoic acids 11-38 were designed and synthesized via the strategy of conformational flexibility analysis. Compounds 21, 22, 26, 33, 36, and 38 with high potency against both wild-type AtAHAS and its P197L mutant were identified as promising candidates with low resistance factors (RF, defined as the ratio between the k_i values toward P197L mutant and wild-type AHAS) ranging from 0.73 to 6.32. Especially, compound 22 (RF = 0.73) was further identified as the most potent antiresistance AHAS inhibitor because of its significantly reduced resistance level compared with that of tribenuron-methyl (RF = 2650) and bispyribac (RF = 4.57). Furthermore, compounds 26, 33, 36, and 38 also displayed promising herbicidal activities against sensitive and resistant (P197L) Descurainia sophia at the dosage of 75-150 g of active ingredient (ai)/ha. Notably, compounds 33 and 38 still maintained over 60% herbicidal activity toward the resistant weed even at much lower dosages (37.5 g ai/ha). Therefore, the designed scaffold has the great potential to discover new candidate compounds for the control of weed resistance associated with AHAS mutation.

Inhibition of acetylcholinesterase by two genistein derivatives: kinetic analysis, molecular docking and molecular dynamics simulation

In this study two genistein derivatives (G1 and G2) are reported as inhibitors of acetylcholinesterase (AChE) and butyrylcholinesterase (BuChE), and differences in the inhibition of AChE are described. Although they differ in structure by a single methyl group, the inhibitory effect of G1 (IC₅₀=264 nmol/L) on AChE was 80 times stronger than that of G2 (IC₅₀=21,210 nmol/L). Enzyme-kinetic analysis, molecular docking and molecular dynamics (MD) simulations were conducted to better understand the molecular basis for this difference. The results obtained by kinetic analysis demonstrated that G1 can interact with both the catalytic active site and peripheral anionic site of AChE. The predicted binding free energies of two complexes calculated by the molecular mechanics/generalized born surface area (MM/GBSA) method were consistent with the experimental data. The analysis of the individual energy terms suggested that a difference between the net electrostatic contributions (ΔE_ele+ΔG_GB) was responsible for the binding affinities of these two inhibitors. Additionally, analysis of the molecular mechanics and MM/GBSA free energy decomposition revealed that the difference between G1 and G2 originated from interactions with Tyr124, Glu292, Val294 and Phe338 of AChE. In conclusion, the results reveal significant differences at the molecular level in the mechanism of inhibition of AChE by these structurally related compounds.

Molecular docking and molecular dynamics studies on the interactions of hydroxylated polybrominated diphenyl ethers to estrogen receptor alpha

Environmental estrogens have attracted great concerns. Recent studies have indicated that some hydroxylated polybrominated diphenyl ethers (HO-PBDEs) can interact with estrogen receptor (ER), and exhibit estrogenic activity. However, interactions between HO-PBDEs and ER are not well understood. In this work, molecular docking and molecular dynamics (MD) simulations were performed to characterize interactions of two HO-PBDEs (4'-HO-BDE30 and 4'-HO-BDE121) with ERα. Surflex-Dock was employed to reveal the probable binding conformations of the compounds at the active site of ERα; MD simulation was used to determine the detailed binding process. The driving forces of the binding between HO-PBDEs and ERα were van der Waals and electrostatic interactions. The decomposition of the binding free energy indicated that the hydrogen bonds between the residues Glu353, Gly521 and ligands were crucial for anchoring the ligands into the active site of ERα and stabilizing their conformations. The results showed that different interaction modes and different specific interactions with some residues were responsible for the different estrogenic activities of the two HO-PBDEs.

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.