Identification of novel natural inhibitor for NorM – a multidrug and toxic compound extrusion transporter – an insilico molecular modeling and simulation studies

The emergence of bacterial multidrug resistance is an increasing problem in treatment of infectious diseases. An important cause for the multidrug resistance of bacteria is the expression of multidrug efflux transporters. The multidrug and toxic compound extrusion (MATE) transporters are most recently recognized as unique efflux system for extrusion of antimicrobials and therapeutic drugs due to energy stored in either Na+ or H+ electrochemical gradient. In the present study, high throughput virtual screening of natural compound collections against NorM – a MATE transporter from Neisseria gonorrhea (NorM-NG) has been carried out followed by flexible docking. The molecular simulation in membrane environment has been performed for understanding the stability and binding energetic of top lead compounds. Results identified a compound from the Indian medicinal plant “Terminalia chebula” which has good binding free energy compared to substrates (rhodamine 6 g, ethidium) and more favorable interactions with the central cavity forming active site residues. The compound has restricted movement in TM7, TM8, and TM1, thus blocking the disruption of Na+ – coordination along with equilibrium state bias towards occlude state of NorM transporter. Thus, this compound blocks the effluxing pathway of antimicrobial drugs and provides as a natural bioactive lead inhibitor against NorM transporter in drug-resistant gonorrhea.

[1]  Alexander D. MacKerell,et al.  Update of the CHARMM all-atom additive force field for lipids: validation on six lipid types. , 2010, The journal of physical chemistry. B.

[2]  R. Friesner,et al.  New insights about HERG blockade obtained from protein modeling, potential energy mapping, and docking studies. , 2006, Bioorganic & medicinal chemistry.

[3]  A Srinivas Reddy,et al.  Virtual screening in drug discovery -- a computational perspective. , 2007, Current protein & peptide science.

[4]  K Schulten,et al.  VMD: visual molecular dynamics. , 1996, Journal of molecular graphics.

[5]  Sarita Rajender Potlapally,et al.  Identification of novel leads applying in silico studies for Mycobacterium multidrug resistant (MMR) protein , 2014, Journal of biomolecular structure & dynamics.

[6]  Laxmikant V. Kalé,et al.  Scalable molecular dynamics with NAMD , 2005, J. Comput. Chem..

[7]  Mihaly Mezei,et al.  Simulaid: A simulation facilitator and analysis program , 2010, J. Comput. Chem..

[8]  Ion binding and internal hydration in the multidrug resistance secondary active transporter NorM investigated by molecular dynamics simulations. , 2012, Biochemistry.

[9]  I. Bahar,et al.  Normal mode analysis of biomolecular structures: functional mechanisms of membrane proteins. , 2010, Chemical reviews.

[10]  R. Stein,et al.  Na⁺-substrate coupling in the multidrug antiporter norm probed with a spin-labeled substrate. , 2013, Biochemistry.

[11]  Alexander D. MacKerell,et al.  CHARMM general force field: A force field for drug‐like molecules compatible with the CHARMM all‐atom additive biological force fields , 2009, J. Comput. Chem..

[12]  Mu-Ping Nieh,et al.  Fluid phase lipid areas and bilayer thicknesses of commonly used phosphatidylcholines as a function of temperature. , 2011, Biochimica et biophysica acta.

[13]  Klaus Schulten,et al.  Accelerating Molecular Modeling Applications with GPU Computing , 2009 .

[14]  P. Bag,et al.  Evaluation of the Antimicrobial Activity of Some Medicinal Plants against Enteric Bacteria with Particular Reference to Multi-Drug Resistant Vibrio cholerae , 2009 .

[15]  J. Symerský,et al.  Structures of a Na+-coupled, substrate-bound MATE multidrug transporter , 2013, Proceedings of the National Academy of Sciences.

[16]  T. Piggot,et al.  The NorM MATE transporter from N. gonorrhoeae: insights into drug and ion binding from atomistic molecular dynamics simulations. , 2014, Biophysical journal.

[17]  T. Darden,et al.  A smooth particle mesh Ewald method , 1995 .

[18]  Asad U. Khan,et al.  Structure-based screening of inhibitors against KPC-2: designing potential drug candidates against multidrug-resistant bacteria , 2014, Journal of biomolecular structure & dynamics.

[19]  Sinosh Skariyachan,et al.  Exploring insights for virulent gene inhibition of multidrug resistant Salmonella typhi, Vibrio cholerae, and Staphylococcus areus by potential phytoligands via in silico screening , 2014, Journal of biomolecular structure & dynamics.

[20]  R. Friesner,et al.  Novel procedure for modeling ligand/receptor induced fit effects. , 2006, Journal of medicinal chemistry.

[21]  Calvin Yu-Chian Chen,et al.  TCM Database@Taiwan: The World's Largest Traditional Chinese Medicine Database for Drug Screening In Silico , 2011, PloS one.

[22]  G J Williams,et al.  The Protein Data Bank: a computer-based archival file for macromolecular structures. , 1978, Archives of biochemistry and biophysics.

[23]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[24]  Matthew P. Repasky,et al.  Glide: a new approach for rapid, accurate docking and scoring. 1. Method and assessment of docking accuracy. , 2004, Journal of medicinal chemistry.

[25]  V. Rajendran,et al.  Drug resistance mechanism of PncA in Mycobacterium tuberculosis , 2014, Journal of biomolecular structure & dynamics.

[26]  Jianpeng Ma,et al.  CHARMM: The biomolecular simulation program , 2009, J. Comput. Chem..

[27]  I. Paulsen,et al.  The multidrug efflux protein NorM is a prototype of a new family of transporters , 1999, Molecular microbiology.

[28]  Matthew P. Repasky,et al.  Extra precision glide: docking and scoring incorporating a model of hydrophobic enclosure for protein-ligand complexes. , 2006, Journal of medicinal chemistry.

[29]  Ivet Bahar,et al.  ProDy: Protein Dynamics Inferred from Theory and Experiments , 2011, Bioinform..

[30]  Hege S. Beard,et al.  Glide: a new approach for rapid, accurate docking and scoring. 2. Enrichment factors in database screening. , 2004, Journal of medicinal chemistry.

[31]  G. Chang,et al.  Structure of a cation-bound multidrug and toxic compound extrusion transporter , 2010, Nature.

[32]  M. Cunningham.  Dynamical probing of allosteric control in nuclear receptors , 2012, Journal of Molecular Modeling.

[33]  D. Beglov,et al.  Atomic Radii for Continuum Electrostatics Calculations Based on Molecular Dynamics Free Energy Simulations , 1997 .

[34]  Velin Z. Spassov,et al.  Introducing an Implicit Membrane in Generalized Born/Solvent Accessibility Continuum Solvent Models , 2002 .

[35]  J. Bajorath,et al.  Quo vadis, virtual screening? A comprehensive survey of prospective applications. , 2010, Journal of medicinal chemistry.

[36]  Jianing Song,et al.  Insights on Na+ binding and conformational dynamics in multidrug and toxic compound extrusion transporter NorM , 2014, Proteins.

[37]  Mihaly Mezei,et al.  TRAJELIX: A Computational Tool for the Geometric Characterization of Protein Helices During Molecular Dynamics Simulations , 2006, J. Comput. Aided Mol. Des..

[38]  Yoshiki Tanaka,et al.  Structural basis for the drug extrusion mechanism by a MATE multidrug transporter , 2013, Nature.

[39]  Ioannis G. Tsoulos,et al.  Eucb: A C++ program for molecular dynamics trajectory analysis , 2011, Comput. Phys. Commun..

[40]  Woody Sherman,et al.  Protein and ligand preparation: parameters, protocols, and influence on virtual screening enrichments , 2013, Journal of Computer-Aided Molecular Design.

[41]  Taehoon Kim,et al.  CHARMM‐GUI: A web‐based graphical user interface for CHARMM , 2008, J. Comput. Chem..

[42]  W. Im,et al.  Automated Builder and Database of Protein/Membrane Complexes for Molecular Dynamics Simulations , 2007, PloS one.

[43]  Hyeon Joo,et al.  OPM database and PPM web server: resources for positioning of proteins in membranes , 2011, Nucleic Acids Res..

[44]  Brian K. Shoichet,et al.  Virtual screening of chemical libraries , 2004, Nature.

[45]  William L. Jorgensen,et al.  Temperature dependence of TIP3P, SPC, and TIP4P water from NPT Monte Carlo simulations: Seeking temperatures of maximum density , 1998, J. Comput. Chem..

[46]  Yuji Sugita,et al.  Analysis of lipid surface area in protein–membrane systems combining voronoi tessellation and monte carlo integration methods , 2012, J. Comput. Chem..

[47]  Laxmikant V. Kale,et al.  NAMD2: Greater Scalability for Parallel Molecular Dynamics , 1999 .

[48]  Toni Giorgino,et al.  MEMBPLUGIN: studying membrane complexity in VMD , 2014, Bioinform..

[49]  Conrad C. Huang,et al.  UCSF Chimera—A visualization system for exploratory research and analysis , 2004, J. Comput. Chem..