ABSTRACT
One widely known herbicide target is 4-hydroxyphenylpyruvate dioxygenase (HPPD). Avena sativa HPPD is less sensitive to mesotrione (herbicide) than Arabidopsis thaliana HPPD. HPPD inhibitor-sensitivity is governed by the dynamic behavior of the C-terminal α-helix of HPPD (H11) in closed and open forms. However, the specific relationship between the plant inhibitor sensitivity and H11 dynamic behavior remains unclear. Herein, we determined the conformational changes in H11 to understand the inhibitor-sensitivity mechanism based on free-energy calculations using molecular dynamics simulations. The calculated free-energy landscapes revealed that Arabidopsis thaliana HPPD preferred the open form of H11 in the apo form and the closed-like form in complex with mesotrione, whereas Avena sativa HPPD exhibited the opposite tendency. We also identified some important residues involved in the dynamic behavior of H11. Therefore, inhibitor sensitivity is governed by indirect interactions due to the protein flexibility caused by the conformational changes of H11.
Subject(s)
4-Hydroxyphenylpyruvate Dioxygenase , Arabidopsis , Dioxygenases , Herbicides , 4-Hydroxyphenylpyruvate Dioxygenase/metabolism , Arabidopsis/metabolism , Cyclohexanones/pharmacology , Herbicides/pharmacology , Herbicides/chemistry , Enzyme Inhibitors/chemistryABSTRACT
We proposed a novel QSAR (quantitative structure-activity relationship) procedure called LERE (linear expression by representative energy terms)-QSAR involving molecular calculations such as ab initio fragment molecular orbital and generalized Born/surface area ones. We applied LERE-QSAR to two datasets for the free-energy changes during complex formation between carbonic anhydrase and a series of substituted benzenesulfonamides. The first compound set (Set I) and the second one (Set II) include relatively small substituents and alkyl chains of different lengths in the benzene ring, respectively. Variation of the inhibitory activity in Set I is expressed as the combination of Hammett σ and the hydrophobic substituent constant π in classical QSAR, and variation in Set II only by π. LERE-QSAR analyses clearly revealed that effects of σ and π on the activity variations in Sets I and II are consistently explainable with the energy terms in the LERE formulation, and provide more detailed and direct information as to the binding mechanism. The proposed procedure was demonstrated to provide a quantitative basis for understanding ligand-protein interactions at the electronic and atomic levels.
Subject(s)
Carbonic Anhydrases/chemistry , Sulfonamides/chemistry , Carbonic Anhydrases/metabolism , Protein Binding , Quantitative Structure-Activity Relationship , Thermodynamics , BenzenesulfonamidesABSTRACT
Quantitative structure-activity relationship analyses on the free energy change during complex formation between substituted benzenesulfonamides (BSAs) and bovine carbonic anhydrase II (bCA II) were performed using generilized Born/surface area (GB/SA) and ab initio fragment molecular orbital (FMO) calculations for the whole complex structures. The result shows that the overall free energy change is governed by the contribution from solvation and dissociation free energy changes accompanying by complex formation. The FMO-IFIE (interfragment interaction energy) analysis quantitatively revealed that the intrinsic interaction energy of bCA II with BSAs is mostly from interactions with amino acid residues in the active site of bCA II. The "Zn block" (Zn(2+) and three histidine residues coordinated to Zn(2+)) in the active site shows the lowest interaction energy and the greatest variance of interaction energy with BSAs through their coordination interaction. The proposed procedure was demonstrated to provide a quantitative basis for understanding a ligand-protein interaction at electronic and atomic levels.