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

Acetohydroxyacid synthase (AHAS) is the target enzyme of several classes of herbicides, such as sulfonylureas and imidazolinones. Now many mutant AHASs with herbicide resistance have emerged along with extensive use of herbicides, therefore it is imperative to understand the detailed interaction mechanism and resistance mechanism so as to develop new potent inhibitors for wild-type or resistant AHAS. With the aid of available crystal structures of the Arabidopsis thaliana (At) AHAS-inhibitor complex, molecular dynamics (MD) simulations were used to investigate the interaction and resistance mechanism directly and dynamically at the atomic level. Nanosecond-level MD simulations were performed on six systems consisting of wild-type or W574L mutant AtAHAS in the complex with three sulfonylurea inhibitors, separately, and binding free energy was calculated for each system using the MM-GBSA method. Comprehensive analyses from structural and energetic aspects confirmed the importance of residue W574, and also indicated that W574L mutation might alert the structural charactersistic of the substrate access channel and decrease the binding affinity of inhibitors, which cooperatively weaken the effective channel-blocked effect and finally result in weaker inhibitory effect of inhibitor and corresponding herbicide resistance of W574L mutant. To our knowledge, it is the first report about MD simulations study on the AHAS-related system, which will pave the way to study the interactions between herbicides and wild-type or mutant AHAS dynamically, and decipher the resistance mechanism at the atomic level for better designing new potent anti-resistance herbicides.