Computational characterization of homologous ligands binding to a deep hydrophobic pocket in Shigella flexneri pilot protein MxiM

The type III secretion system (T3SS) is an important molecular machinery in gram-negative bacteria Shigella flexneri as it provides ways for translocating virulence factors from the bacteria into host cells, eventually leading to severe disease symptoms such as bacillary dysentery. Due to the rising concerns of antibiotics resistance in bactericidal strategy, the anti-virulence strategy that primarily targets the T3SS components becomes an attractive alternative. MxiM, the secretin pilot protein of Shigella flexneri, binds the secretin MxiD and facilitates the formation of the secretin ring in outer membrane in T3SS assembly. MxiM harbors a large hydrophobic pocket that has been shown to be important in MxiM-MxiD interaction. In this work, I examined the ligand binding property of MxiM by performing molecular dynamics (MD) simulations of the association between MxiM and a series of hydrophobic ligands, with simulation time amounted to 30 mu s. MD simulations successfully captured spontaneous ligand binding events in 153 of the 300 trajectories. The ligand binding can be categorized into two types: a fast type, in which the ligand binds quickly into the hydrophobic pocket and a slow type, in which the ligand forms an encounter complex with the protein before binding into the hydrophobic pocket. Using the MxiM-ligand binding poses captured in MD simulations, I additionally performed umbrella-sampling MD simulations with total simulation time amounted to 63 mu s to obtain protein-ligand binding free energies. The relationship between the ligand binding free energy and ligand size appears to be nonlinear and exhibits an exponential decay pattern. In summary, I performed computational characterization of MxiM-hydrophobic ligand binding capabilities and properties, which may provide valuable insights into designing anti-bacterial medicine against antibiotics resistance in Shigella flexneri.