The widespread use of glyphosate and the high dependence of the agricultural industry on this herbicide cause environmental pollution and pose a threat to living organisms. One of the appropriate solutions in sustainable agriculture to deal with pollution caused by glyphosate and its metabolites is creating a symbiotic relationship between plants and mycorrhizal fungi. Glomalin-related soil protein is a key protein for the bioremediation of glyphosate and its metabolite aminomethyl phosphonic acid in soil. This study uses homology modeling, molecular docking, and molecular dynamic simulation approaches to investigate the binding mechanism of glomalin-related soil protein from arbuscular mycorrhiza (GmHsp60) with glyphosate and its metabolite and the role of soil protein in the removal and sequestering of common agricultural soil pollutants. GmHsp60 protein structure was predicted by homology modeling, and the quality of the generated model was assessed. Then, the interaction between glyphosate and aminomethyl phosphonic acid and the modeled GmHsp60 protein was explored by molecular docking. Based on docking results, GmHsp60 has an efficient role in the bioremediation of glyphosate and aminomethyl phosphonic acid (-6.03 and -5.34 kcal/mol). Glyphosate forms three hydrogen bonds with Lys258, Gly262, and Glu58 of GmHsp60, and aminomethyl phosphonic acid forms three hydrogen bonds with Lys258, Gly261, and Gly262 of GmHsp60. In addition, the glyphosate's and its metabolite's stability was confirmed by molecular docking simulations and binding free energy calculations using MM/PBSA analysis. This study provides a molecular-level understanding of GmHsp60 expression and function for glyphosate bioremediation.
Keywords: Glyphosate bioremediation; arbuscular mycorrhizal fungi; glomalin; molecular docking; molecular dynamic simulations.