Identification of New Delhi metallo-β-lactamase gene (NDM-1) from a clinical isolate of Acinetobacter junii in China

New Delhi metallo-β-lactamase-1 (NDM-1) is a novel type of metallo-β-lactamase (MBL) responsible for bacterial resistance to β-lactam antibiotics. Acinetobacter junii was previously shown to possess a MBL phenotype; however, the genes responsible for this phenotype were not identified. In this study, we reported the identification of NDM-1 gene in a clinical isolate of A. junii from a child patient in China, which was resistant to all β-lactams except aztreonam but sensitive to aminoglycosides and quinolones. The cloned NDM-1 gene contained an open reading frame of 813 bp and had a nucleotide sequence 99.9% identical (812/813) to reported NDM-1 genes carried by Acinetobacter baumannii , Enterococcus faecium , Escherichia coli , and Klebsiella pneumoniae . Recombinant NDM-1 protein was successfully expressed in E.?coli BL21, and antibiotic sensitivities of the NDM-1-producing E.?coli were largely similar to the A.?junii 1454 isolate. The findings of this study raise attention to the emergence and spread of NDM-1-carrying bacteria in China.     

Mechanistic and spectroscopic studies of metallo-β-lactamase NDM-1

In an effort to biochemically characterize metallo-β-lactamase NDM-1, we cloned, overexpressed, purified, and characterized several maltose binding protein (MBP)-NDM-1 fusion proteins with different N-termini (full-length, Δ6, Δ21, and Δ36). All MBP-NDM-1 fusion proteins were soluble; however, only one, MBP-NDM-1Δ36, exhibited high activity and bound 2 equiv of Zn(II). Thrombin cleavage of this fusion protein resulted in the truncated NDM-1Δ36 variant, which exhibited a k(cat) of 16 s(-1) and a K(m) of 1.1 μM when using nitrocefin as a substrate, bound 2 equiv of Zn(II), and was monomeric in solution. Extended X-ray absorption fine structure studies of the NDM-1Δ36 variant indicate the average metal binding site for Zn(II) in this variant consists of four N/O donors (two of which are histidines) and 0.5 sulfur donor per zinc, with a Zn-Zn distance of 3.38 ?. This metal binding site is very similar to those of other metallo-β-lactamases that belong to the B1 subclass. Pre-steady-state kinetic studies using nitrocefin and chromacef and the NDM-1Δ36 variant indicate that the enzyme utilizes a kinetic mechanism similar to that used by metallo-β-lactamases L1 and CcrA, in which a reactive nitrogen anion is stabilized and its protonation is rate-limiting. While they are very different in terms of amino acid sequence, these studies demonstrate that NDM-1 is structurally and mechanistically very similar to metallo-β-lactamase CcrA.     

Crystal structure of New Delhi metallo-β-lactamase reveals molecular basis for antibiotic resistance

β‐Lactams are the most commonly prescribed class of antibiotics and have had an enormous impact on human health. Thus, it is disquieting that an enzyme called New Delhi metallo‐β‐lactamase‐1 (NDM‐1) can confer Enterobacteriaceae with nearly complete resistance to all β‐lactam antibiotics including the carbapenams. We have determined the crystal structure of Klebsiella pneumoniae apo‐NDM‐1 to 2.1‐Å resolution. From the structure, we see that NDM‐1 has an expansive active site with a unique electrostatic profile, which we propose leads to a broader substrate specificity. In addition, NDM‐1 undergoes important conformational changes upon substrate binding. These changes have not been previously observed in metallo‐β‐lactamase enzymes and may have a direct influence on substrate recognition and catalysis.     

New Delhi metallo-β-lactamase (NDM-1)-producing Klebsiella pneumoniae of sequence type ST11: first identification in a hospital of central Italy

The emergence of novel resistant markers hampers the efficacy of beta-lactam antibiotics to treat infections caused by micro-organisms carrying such resistances. This study investigated the antimicrobial susceptibility pattern, the carpapenem-associated determinants and the molecular epidemiology of Klebsiella pneumoniae showing a New Delhi (NDM) metallo-β-lactamase phenotype, isolated from a patient admitted to intensive care unit of the main hospital for acute care of Molise region, central Italy. Antimicrobial susceptibility was assessed for nineteen antibiotics by disc diffusion and agar dilution methods. Carbapenem-associated resistance determinants were detected through gene-specific amplifications, targeting bla_NDM-1, bla_SHV and bla_TEM, bla_CTX-M, bla_KPC, bla_VIM, bla_IMP, bla_GES and bla_OXA-48-lixe. Molecular characterization was carried out through multilocus sequence typing. The strain showed a multidrug resistant profile, and PCR and sequencing confirmed the presence of bla_NDM-1 gene. Among the multiple resistance-associated determinants tested, the isolate, which was assigned to the sequence type ST11, only harboured bla_SHV and bla_TEM genes. This is the first report of NDM-1 variant in the regional healthcare setting for acute patients, raising significant concerns about the increase in the antimicrobials resistance spread through a different mechanism from the endemic KPC carbapenemase, and underlining the circulation of a virulent clone never identified before in this area.     

Discovery of a novel covalent non-β-lactam inhibitor of the metallo-β-lactamase NDM-1

The inhibition of metallo-β-lactamases (MBL) can prevent the hydrolysis of β-lactam antibiotics and hence is a promising strategy for the treatment of antibiotic resistant infections. In this study, we present a novel reversible covalent inhibitor of the clinically relevant MBL New Delhi metallo-β-lactamase 1 (NDM-1). Electrospray ionization-mass spectrometry (ESI-MS) and single site directed mutagenesis were used to show that the inhibitor forms a covalent bond with Lys224 in the active site of NDM-1. The inhibitor was further characterized using an enzyme inhibition assay, a surface plasmon resonance (SPR) based biosensor assay and covalent docking. The determined inhibition constant (K_I¹) was 580 nM and the inhibition constant for the initial complex (K_I) was 76 μM. To our knowledge, this inhibitor is the first example for a reversible covalent non-β-lactam inhibitor targeting NDM-1 and a promising starting point for the design of potent covalent inhibitors.     

Molecular dynamics simulation and molecular docking studies of 1,4-Dihydropyridines as P-glycoprotein’s allosteric inhibitors

P-glycoprotein (P-gp) is a main factor contributing to multidrug resistance. The effect of this transporter protein on limiting the effectiveness of chemotherapy has been shown by various studies. In a previous report, we synthesized some 14-dihydropyridine (DHP) derivatives as inhibitors of human P-gp. In the present study, a computational approach has been exploited to reveal the main interactions between DHPs and P-gp. In order to do this, homology modeling was performed to obtain a model of the protein. Then, molecular dynamics simulation was used to refine the constructed model of P-gp in the presence of the lipids bilayer. Model validation was performed with several tools. Finally, molecular docking followed by MD simulation of ligand–protein complex was employed to elucidate the binding mode and the dynamical changes of protein with/without DHPs bound. The results emphasized that interaction of the residues Gln912, Ser909, Arg905, Ser474, Val472 with DHPs play a crucial role in the inhibitory of these ligands and this was in a relatively good accordance with the results reported in the experimental studies.     

3D-QSAR,molecular docking and molecular dynamics studies of a series of RORγt inhibitors

The discovery of clinically relevant inhibitors of retinoic acid receptor-related orphan receptor-gamma-t (RORγt) for autoimmune diseases therapy has proven to be a challenging task. In the present work, to find out the structural features required for the inhibitory activity, we show for the first time a three-dimensional quantitative structure–activity relationship (3D-QSAR), molecular docking and molecular dynamics (MD) simulations for a series of novel thiazole/thiophene ketone amides with inhibitory activity at the RORγt receptor. The optimum CoMFA and CoMSIA models, derived from ligand-based superimposition I, exhibit leave-one-out cross-validated correlation coefficient (R²_cv) of .859 and .805, respectively. Furthermore, the external predictive abilities of the models were evaluated by a test set, producing the predicted correlation coefficient (R²_pred) of .7317 and .7097, respectively. In addition, molecular docking analysis was applied to explore the binding modes between the inhibitors and the receptor. MD simulation and MM/PBSA method were also employed to study the stability and rationality of the derived conformations, and the binding free energies in detail. The QSAR models and the results of molecular docking, MD simulation, binding free energies corroborate well with each other and further provide insights regarding the development of novel RORγt inhibitors with better activity.     

Atom-based 3D-QSAR,molecular docking and molecular dynamics simulation assessment of inhibitors for thyroid hormone receptor α and β

HCl, HNO₃ and H₂SO₄ are implicated in atmospheric processes in areas such as polar stratospheric clouds in the stratosphere. Ternary complexes of HCl, HNO₃ and H₂SO₄ were investigated by ab initio calculations at B3LYP level of theory with aug-cc-pVTZ and aug-cc-pVQZ basis sets, taking into account basis set superposition error (BSSE). The results were assessed in terms of structures (five hexagonal cyclic structures and two quasi-pentagonal cyclic structures), inter-monomeric parameters (all ternary complexes built three hydrogen bonds), energetics (seven minima obtained), infrared harmonic vibrational frequencies (red shifting of complexes from monomers), and relative stability of complexes, which were favorable when the temperature decreases under stratospheric conditions, from 298 K to 188 K, and in concrete, at 210 K, 195 K and 188 K.     

In silico–based combinatorial pharmacophore modelling and docking studies of GSK-3β and GK inhibitors of Hippophae

Type 2 diabetes is an inevitably progressive disease, with irreversible β cell failure. Glycogen synthase kinase and Glukokinase, two important enzymes with diverse biological actions in carbohydrate metabolism, are promising targets for developing novel antidiabetic drugs. A combinatorial structure-based molecular docking and pharmacophore modelling study was performed with the compounds of Hippophae salicifolia and H. rhamnoides as inhibitors. Docking with Discovery Studio 3.5 revealed that two compounds from H. salicifolia, viz Lutein D and an analogue of Zeaxanthin, and two compounds from H. rhamnoides, viz Isorhamnetin-3-rhamnoside and Isorhamnetin-7-glucoside, bind significantly to the GSK-3 β receptor and play a role in its inhibition; whereas in the case of Glucokinase, only one compound from both the plants, i.e. vitamin C, had good binding characteristics capable of activation. The results help to understand the type of interactions that occur between the ligands and the receptors. Toxicity predictions revealed that none of the compounds had hepatotoxic effects and had good absorption as well as solubility characteristics. The compounds did not possess plasma protein-binding, crossing blood–brain barrier ability. Further, in vivo and in vitro studies need to be performed to prove that these compounds can be used effectively as antidiabetic drugs.     

Structure of apo- and monometalated forms of NDM-1--a highly potent carbapenem-hydrolyzing metallo-β-lactamase

The New Delhi Metallo-β-lactamase (NDM-1) gene makes multiple pathogenic microorganisms resistant to all known β-lactam antibiotics. The rapid emergence of NDM-1 has been linked to mobile plasmids that move between different strains resulting in world-wide dissemination. Biochemical studies revealed that NDM-1 is capable of efficiently hydrolyzing a wide range of β-lactams, including many carbapenems considered as "last resort" antibiotics. The crystal structures of metal-free apo- and monozinc forms of NDM-1 presented here revealed an enlarged and flexible active site of class B1 metallo-β-lactamase. This site is capable of accommodating many β-lactam substrates by having many of the catalytic residues on flexible loops, which explains the observed extended spectrum activity of this zinc dependent β-lactamase. Indeed, five loops contribute "keg" residues in the active site including side chains involved in metal binding. Loop 1 in particular, shows conformational flexibility, apparently related to the acceptance and positioning of substrates for cleavage by a zinc-activated water molecule.     

Role of invariant water molecules in retaining the active site geometry of β-lactamase: a molecular dynamics simulation study

Water molecules play a critical role in stabilising the three-dimensional architecture, dynamics and function of biological macromolecules. Comparative analysis of structurally similar proteins has shown that there are water molecules conserved in the same relative positions and make similar hydrogen bonds with proteins in all crystal structures. These invariant water molecules are essential for the maintenance of the native structure of proteins. The present study explores the role of invariant water molecules to maintain the active site geometry of β-lactamase enzyme. Thirteen crystal structures of class-A β-lactamase from Staphylococcus aureus have been used in this study. Molecular dynamics simulations of the protein structures were performed in hydrated as well as in dehydrated conditions. The analysis showed that significant changes occur in the active site geometry due to dehydration. These changes can be attributed to the removal of water molecules at the active site.     

Synthesis of oxadiazole-coupled-thiadiazole derivatives as a potent β-glucuronidase inhibitors and their molecular docking study

A new series of oxadiazole with thiadiazole moiety (6–27) were synthesized, characterized by different spectroscopic techniques and evaluated for β-glucuronidase inhibitory potential. Sixteen analogs such as 6, 7, 8, 9, 10, 12, 13, 14, 17, 18, 20, 23, 24, 25, 26 and 27 showed IC₅₀ values in the range of 0.96 ± 0.01 to 46.46 ± 1.10 μM, and hence were found to have excellent inhibitory potential in comparison to standard d-saccharic acid 1,4-lactone (IC₅₀ = 48.4 ± 1.25 μM). Two analogs such as 16 and 19 showed moderate inhibitory potential while analogs 11, 15, 21 and 22 were found inactive. Our study identifies new series of potent β-glucuronidase inhibitors for further investigation. Structure activity relationships were established for all compounds which showed that the activity is varied due to different substituents on benzene ring. The interaction of the compounds with enzyme active site were confirmed with the help of docking studies, which reveals that the electron withdrawing group and hydroxy group make the molecules more favorable for enzyme inhibition.     

Stepwise development of structure–activity relationship of diverse PARP-1 inhibitors through comparative and validated in silico modeling techniques and molecular dynamics simulation

Inhibitors of poly (ADP-ribose) polymerase-1 (PARP-1) enzyme are useful for the treatment of various diseases including cancer. Comparative in silico studies were performed on different ligand-based (2D-QSAR, Kernel-based partial least square (KPLS) analysis, Pharmacophore Search Engine (PHASE) pharmacophore mapping), and structure-based (molecular docking, MM-GBSA analyses, Gaussian-based 3D-QSAR analyses on docked poses) modeling techniques to explore the structure–activity relationship of a diverse set of PARP-1 inhibitors. Two-dimensional (2D)-QSAR highlighted the importance of charge topological index (JGI7), fractional polar surface area (Jurs_FPSA3), and connectivity index (CIC2) along with different molecular fragments. Favorable and unfavorable fingerprints were demonstrated in KPLS analysis, whereas important pharmacophore features (one acceptor, one donor, and two ring aromatic) along with favorable and unfavorable field effects were demonstrated in PHASE-based pharmacophore model. MM-GBSA analyses revealed significance of different polar, non-polar, and solvation energies. Docking-based alignment of ligands was used to perform Gaussian-based 3D-QSAR study that further demonstrated importance of different field effects. Overall, it was found that polar interactions (hydrogen bonding, bridged hydrogen bonding, and pi–cation) play major roles for higher activity. Steric groups increase the total contact surface area but it should have higher fractional polar surface area to adjust solvation energy. Structure-based pharmacophore mapping spotted the positive ionizable feature of ligands as the most important feature for discriminating highly active compounds from inactives. Molecular dynamics simulation, conducted on highly active ligands, described the dynamic behaviors of the protein complexes and supported the interpretations obtained from other modeling analyses. The current study may be useful for designing PARP-1 inhibitors.     

Potent and novel 11β-HSD1 inhibitors identified from shape and docking based virtual screening

Several potent and novel 11β-hydroxysteroid dehydrogenase 1 (11β-HSD1) inhibitors were discovered from in silico screening the commercially available Maybridge database. Among them, seven hit compounds showed good affinity, with IC(50) values lower than 100 nM and the best one 3.7 nM. To select the lead for further optimization, computational ADME/T prediction, the CYP3A4 inhibition and 11β-HSD1 over 11β-HSD2 selectivity test were also performed. Taking all of the above factors into consideration, two promising compounds were selected as lead structures for further development. The employed hierarchical virtual screening protocol not only demonstrates its efficiency, but also provides novel and selective compounds for developing 11β-HSD1 inhibitors to protect against metabolic syndrome.     

3D-QSAR pharmacophore-based virtual screening,molecular docking and molecular dynamics simulation toward identifying lead compounds for NS2B–NS3 protease inhibitors

NS2B–NS3 protease has been identified to serve as lead drug design target due to its significant role in West Nile viral (WNV) and dengue virus (DENV) reproduction and replication. There are currently no approved chemotherapeutic drugs and effective vaccines to inhibit DENV and WNV infections. In this work, 3D-QSAR pharmacophore model has been developed to discover potential inhibitory candidates. Validation through Fischer’s model and decoy test indicate that the developed 3D pharmacophore model is highly predictive for DENV inhibitors, which was then employed to screen ZINC chemical library to obtain reasonable hits. Following ADMET filtering, 15 hits were subjected to further filter through molecular docking and CoMFA modeling. Finally, top three hits were identified as lead compounds or potential inhibitory candidates with IC₅₀ values of ～0.4637?µM and fitness of ～57.73. It is implied from CoMFA modeling that substituents at the side site of benzotriazole such as a p-nitro group (e.g. biphenyl head) and a carbonyl (e.g. carboxylate function) at the side site of furan or amino group may improve bioactivity of ZINC85645245, respectively. Molecular dynamics simulations (MDS) were performed to discover new interactions and reinforce the binding modes from docking for the hits also. The QSAR and MDS results obtained from this work should be useful in determining structural requirements for inhibitor development as well as in designing more potential inhibitors for NS2B–NS3 protease.     

Why tazobactam and sulbactam have different intermediates population with SHV-1 β-lactamase: a molecular dynamics study

The imine intermediates of tazobactam and sulbactam bound to SHV-1 β-lactamase were investigated by molecular dynamics (MD) simulation respectively. Hydrogen bond networks around active site were found different between tazobactam and sulbactam acyl-enzymes. In tazobactam imine intermediate, it was observed that the triazolyl ring formed stable hydrogen bonds with Asn170 and Thr167. The results suggest that conformation of imine determined the population of intermediates. In imine intermediate of tazobactam, the triazolyl ring is trapped in Thr_Asn pocket, and it restricts the rotation of C5-C6 bond so that tazobactam can only generate trans enamine intermediate. Further, conformational cluster analyses are performed to substantiate the results. These findings provide an explanation for the corresponding experimental results, and will be potentially useful in the development of new inhibitors.