Molecular dynamics simulations on interactions of five antibiotics with luciferase of Vibrio Qinghaiensis sp.-Q67

Abstract

A large number of antibiotics have been used in the medical industry, agriculture, and animal husbandry industry in recent years. It may cause pollution to the aquatic environment and ultimately threaten to human health due to their prolonged exposure to the environment. We aim to study the toxicity mechanism of enrofloxacin (ENR), chlortetracycline hydrochloride (CTC), trimethoprim (TMP), chloramphenicol (CMP), and erythromycin (ETM) to luciferase of Vibrio Qinghaiensis sp.-Q67 (Q67) by using toxicity testing combined with molecular docking, molecular dynamics, and binding free energy analysis. The curve categories for ENR were different from the other four antibiotics, with ENR being J-type and the rest being S-type, and the toxicity of these five antibiotics (pEC50) followed the order of ENR (7.281) > ETM (6.814) > CMP (6.672) > CTC (6.400) > TMP (6.123), the order of toxicity value is consistent with the the magnitude of the binding free energy (ENR (-47.759 kcal/mol), ETM (-46.821 kcal/mol), CMP (-42.905 kcal/mol), CTC (-40.946 kcal/mol), TMP (-28.251 kcal/mol)). The van der Waals force provided the most important contribution to the binding free energy of the five antibiotics in the binding system with Q67 luciferase. Therefore, the dominant factor for the binding of antibiotics to luciferase was shape compensation. The face-to-face π-π stacking interaction between the diazohexane structure outside the active pocket region and the indoles structure of Phe194 and Phe250 in the molecular structure was the main reason for the highest toxicity value of antibiotic ENR. The hormesis effect of ENR has a competitive binding relationship with the α and β subunits of luciferase. Homology modeling, molecular docking, molecular dynamics simulations and binding free energy calculations were used to derive the toxicity magnitude of different antibiotics against Q67, and insights at the molecular level. The conclusion of toxicological experiments verified the correctness of the simulation results. This study contributes to the understanding of toxicity mechanisms of five antibiotics and facilitates risk assessment of antibiotic contaminants in the aquatic environment.

Keywords: Antibiotic; Homology modeling; Molecular docking; Molecular dynamics simulation.