Molecular Docking and Dynamics Simulation of Protein β-Tubulin and Antifungal Cyclic Lipopeptides

To elucidate interactions between the antifungal cyclic lipopeptides iturin A, fengycin, and surfactin produced by Bacillus bacteria and the microtubular protein β-tubulin in plant pathogenic fungi (Fusarium oxysporum, Colletrotrichum gloeosporioides, Alternaria alternata, and Fusarium solani) in molecular docking and molecular dynamics simulations, we retrieved the structure of tubulin co-crystallized with taxol from the Protein Data Bank (PDB) (ID: 1JFF) and the structure of the cyclic lipopeptides from PubChem (Compound CID: 102287549, 100977820, 10129764). Similarity and homology analyses of the retrieved β-tubulin structure with those of the fungi showed that the conserved domains shared 84% similarity, and the root mean square deviation (RMSD) was less than 2 Å. In the molecular docking studies, within the binding pocket, residues Pro274, Thr276, and Glu27 of β-tubulin were responsible for the interaction with the cyclic lipopeptides. In the molecular dynamics analysis, two groups of ligands were formed based on the number of poses analyzed with respect to the RMSD. Group 1 was made up of 10, 100, and 500 poses with distances 0.080 to 0.092 nm and RMSDs of 0.10 to 0.15 nm. For group 2, consisting of 1000 poses, the initial and final distance was 0.1 nm and the RMSDs were in the range of 0.10 to 0.30 nm. These results suggest that iturin A and fengycin bind with higher affinity than surfactin to β-tubulin. These two lipopeptides may be used as lead compounds to develop new antifungal agents or employed directly as biorational products to control plant pathogenic fungi.