Molecular dynamics simulation studies of novel β-lactamase inhibitor

Abstract

New Delhi Metallo-β-Lactamase-1 (NDM-1) has drawn great attention due to its diverse antibiotic resistant activity. It can hydrolyze almost all clinically available β-lactam antibiotics. To inhibit the activity of NDM-1 a new strategy is proposed using computational methods. Molecular dynamics (MD) simulations are used to analyze the molecular interactions between selected inhibitor candidates and NDM-1 structure. The enzyme-ligand complex is subject to binding free energy calculations using MM(PB/GB)SA methods. The role of each residue of the active site contributing in ligand binding affinity is explored using energy decomposition analysis. Furthermore, a hydrogen bonding network between ligand and enzyme active site is observed and key residues are identified ensuring that the ligand stays inside the active site and maintains its movement towards the active site pocket. A production run of 150ns is carried out and results are analyzed using root mean square deviation (RMSD), root mean square fluctuation (RMSF), and radius of gyration (Rg) to explain the stability of enzyme ligand complex. Important active site residue e.g. PHE70, VAL73, TRP93, HIS122, GLN123, ASP124, HIS189, LYS216, CYS208, LYS211, ALA215, HIS250, and SER251 were observed to be involved in ligand attachemet inside the active site pocket, hence depicting its inhibitor potential. Hydrogen bonds involved in structural stability are analyzed through radial distribution function (RDF) and contribution of important residues involved in ligand movement is explained using a novel analytical tool, axial frequency distribution (AFD) to observe the role of important hydrogen bonding partners between ligand atoms and active site residues.

Keywords: Antibiotic resistance; Molecular docking; Molecular dynamics simulations; Multi Drug Resistant Organisms (MDROs); New Delhi metallo-β-lactamase-1 (NDM-1).