An in silico approach for understanding the molecular evolution of clinically important metallo-beta-lactamases

Abstract

Resistance to carbapenem, considered to be the drug of last resort for treating serious enterobacterial infections, is a growing problem. Metallo-beta-lactamase (MBL) mediated carbapenem resistance is considered to be clinically most important, since no traditional inhibitors work against them. Hence, we have investigated the changes in drug profile, affinity and binding stability of meropenem with clinically important MBLs viz., IMP, VIM and NDM during the course of molecular evolution. Phylogenetic trees were constructed and amino acids positions, presumed to be exposed to positive selection pressure were analyzed. Based on sequence diversity and selection pressure, IMP-1, IMP-8, IMP-9, IMP-21, IMP-27, IMP-20 and IMP-26 among IMP genes; VIM-1, VIM-2, VIM-13, VIM-29, VIM-18 and VIM-7 among VIM genes and NDM-1 had been selected for in silico analysis. Mode of interaction of selected MBL variants with meropenem were analyzed by Autodock4.0 and LIGPLOT analysis. In all the trajectories, we had found an increase in mouth opening/solvent accessible volume/area of the catalytic pocket and decrease in Gibbs' free energy (ΔG°) for binding with meropenem and Michealis-Menten constant (Km) - indicating an increase in choice of drugs, binding stability and meropenem affinity from primitive to advanced MBL variants, with exceptions of IMP-20, IMP-26, VIM-13 and VIM-18 which might be due to sign epistasis. Intergenic comparison revealed NDM-1 might have greater drug profile and catalytic efficiency than IMP-1 and VIM-2 due to largest pocket opening and least distance between the Zn-I ion and β-lactam oxygen of meropenem.

Keywords: Binding stability; Evolutionary trajectory; IMP/VIM/NDM; Meropenem affinity; Metallo-beta-lactamase.