Refined models of New Delhi metallo-beta-lactamase-1 with inhibitors: an QM/MM modeling study

A modified bonded model approach for molecular dynamics simulations of New Delhi Metallo-β-lactamase

Modelling metalloproteins using the classical force fields is challenging. Several methods have been devised to model metalloproteins in force fields. Of these methods, the bonded model, combined with Restrained Electrostatic Potential (RESP) charge fitting, proved its superiority. The latter method was facilitated by the development of the python-based Metal Centre Parameter Builder (MCPB.py) AmberTool. However, the standard bonded model method offered by the MCPB.py tool may not be appropriate for validating and refining the binding modes predicted by docking when crystal structures are lacking. That is because the representation of coordination interactions between any bound ligand and metal ions by covalent bonds can hinder the flexibility of the ligand. Therefore, a new modification to the standard bonded model approach is proposed here. Molecular dynamics (MD) simulations based on the new modified bonded model (MBM) approach avoid the bias caused by coordination bonds and, unlike hybrid QM/MM MD, allow for sufficient sampling of the binding mode given the currently available computational power. The MBM MD approach reproduced the studied crystal structure conformations of New Delhi Metallo-β-lactamase 1 (NDM-1). Furthermore, the MBM approach described the binding interactions of intact β-lactams with NDM-1 reasonably, and predicted a non-productive binding mode for the poor NDM-1 substrate aztreonam whilst predicting productive binding modes for known good substrates. This study presents a useful MD method for metallo-β-lactamases and provides better understanding of β-lactam substrates recognition by NDM-1. The proposed MBM approach might also be useful in the investigation of other metal-containing protein targets.

Crystallography and QM/MM Simulations Identify Preferential Binding of Hydrolyzed Carbapenem and Penem Antibiotics to the L1 Metallo-β-Lactamase in the Imine Form

Widespread bacterial resistance to carbapenem antibiotics is an increasing global health concern. Resistance has emerged due to carbapenem-hydrolyzing enzymes, including metallo-β-lactamases (MβLs), but despite their prevalence and clinical importance, MβL mechanisms are still not fully understood. Carbapenem hydrolysis by MβLs can yield alternative product tautomers with the potential to access different binding modes. Here, we show that a combined approach employing crystallography and quantum mechanics/molecular mechanics (QM/MM) simulations allow tautomer assignment in MβL:hydrolyzed antibiotic complexes. Molecular simulations also examine (meta)stable species of alternative protonation and tautomeric states, providing mechanistic insights into β-lactam hydrolysis. We report the crystal structure of the hydrolyzed carbapenem ertapenem bound to the L1 MβL from Stenotrophomonas maltophilia and model alternative tautomeric and protonation states of both hydrolyzed ertapenem and faropenem (a related penem antibiotic), which display different binding modes with L1. We show how the structures of both complexed β-lactams are best described as the (2S)-imine tautomer with the carboxylate formed after β-lactam ring cleavage deprotonated. Simulations show that enamine tautomer complexes are significantly less stable (e.g., showing partial loss of interactions with the L1 binuclear zinc center) and not consistent with experimental data. Strong interactions of Tyr32 and one zinc ion (Zn1) with ertapenem prevent a C6 group rotation, explaining the different binding modes of the two β-lactams. Our findings establish the relative stability of different hydrolyzed (carba)penem forms in the L1 active site and identify interactions important to stable complex formation, information that should assist inhibitor design for this important antibiotic resistance determinant.

Multiscale Workflow for Modeling Ligand Complexes of Zinc Metalloproteins

Zinc metalloproteins are ubiquitous, with protein zinc centers of structural and functional importance, involved in interactions with ligands and substrates and often of pharmacological interest. Biomolecular simulations are increasingly prominent in investigations of protein structure, dynamics, ligand interactions, and catalysis, but zinc poses a particular challenge, in part because of its versatile, flexible coordination. A computational workflow generating reliable models of ligand complexes of biological zinc centers would find broad application. Here, we evaluate the ability of alternative treatments, using (nonbonded) molecular mechanics (MM) and quantum mechanics/molecular mechanics (QM/MM) at semiempirical (DFTB3) and density functional theory (DFT) levels of theory, to describe the zinc centers of ligand complexes of six metalloenzyme systems differing in coordination geometries, zinc stoichiometries (mono- and dinuclear), and the nature of interacting groups (specifically the presence of zinc-sulfur interactions). MM molecular dynamics (MD) simulations can overfavor octahedral geometries, introducing additional water molecules to the zinc coordination shell, but this can be rectified by subsequent semiempirical (DFTB3) QM/MM MD simulations. B3LYP/MM geometry optimization further improved the accuracy of the description of coordination distances, with the overall effectiveness of the approach depending upon factors, including the presence of zinc-sulfur interactions that are less well described by semiempirical methods. We describe a workflow comprising QM/MM MD using DFTB3 followed by QM/MM geometry optimization using DFT (e.g., B3LYP) that well describes our set of zinc metalloenzyme complexes and is likely to be suitable for creating accurate models of zinc protein complexes when structural information is more limited.

Discovery of the Novel Inhibitor Against New Delhi Metallo-β-Lactamase Based on Virtual Screening and Molecular Modelling

New Delhi metallo-β-lactamase (NDM-1), one of the metallo-β-lactamases (MBLs), leads to antibiotic resistance in clinical treatments due to the strong ability of hydrolysis to almost all kinds of β-lactam antibiotics. Therefore, there is the urgent need for the research and development of the novel drug-resistant inhibitors targeting NDM-1. In this study, ZINC05683641 was screened as potential NDM-1 inhibitor by virtual screening and the inhibitor mechanism of this compound was explored based on molecular dynamics simulation. The nitrocefin assay showed that the IC₅₀ value of ZINC05683641 was 13.59 ± 0.52 μM, indicating that the hydrolytic activity of NDM-1 can be obviously suppressed by ZINC05683641. Further, the binding mode of ZINC05683641 with NDM-1 was obtained by molecular modeling, binding free energy calculation, mutagenesis assays and fluorescence-quenching assays. As results, ILE-35, MET-67, VAL-73, TRP-93, CYS-208, ASN-220 and HIS-250 played the key roles in the binding of NDM-1 with ZINC05683641. Interestingly, these key residues were exactly located in the catalytic activity region of NDM-1, implying that the inhibitor mechanism of ZINC05683641 against NDM-1 was the competitive inhibition. These findings will provide an available approach to research and develop new drug against NDM-1 and treatment for bacterial resistance.

A molecular dynamics study of the complete binding process of meropenem to New Delhi metallo-β-lactamase 1

The mechanism of substrate hydrolysis of New Delhi metallo-β-lactamase 1 (NDM-1) has been reported, but the process in which NDM-1 captures and transports the substrate into its active center remains unknown. In this study, we investigated the process of the substrate entry into the NDM-1 activity center through long unguided molecular dynamics simulations using meropenem as the substrate. A total of 550 individual simulations were performed, each of which for 200 ns, and 110 of them showed enzyme-substrate binding events. The results reveal three categories of relatively persistent and noteworthy enzyme-substrate binding configurations, which we call configurations A, B, and C. We performed binding free energy calculations of the enzyme-substrate complexes of different configurations using the molecular mechanics Poisson-Boltzmann surface area method. The role of each residue of the active site in binding the substrate was investigated using energy decomposition analysis. The simulated trajectories provide a continuous atomic-level view of the entire binding process, revealing potentially valuable regions where the enzyme and the substrate interact persistently and five possible pathways of the substrate entering into the active center, which were validated using well-tempered metadynamics. These findings provide important insights into the binding mechanism of meropenem to NDM-1, which may provide new prospects for the design of novel metallo-β-lactamase inhibitors and enzyme-resistant antibiotics.

Molecular modeling studies and anti-TB activity of trisubstituted indolizine analogues; molecular docking and dynamic inputs

A series of trisubstituted indolizine analogues has been designed as a result of a fragment-based approach to target the inhibition of mycobacterial enoyl-acyl carrier protein reductase. Anti-tuberculosis (TB) screening of the characterized compounds by a resazurin microplate assay method revealed that ethyl group at second position of indolizine nucleus exhibited activity against susceptible and multidrug-resistant strains of Mycobacterium tuberculosis at concentration of 5.5 and 11.3 μg/mL, respectively. A molecular docking study was also conducted to evaluate the stability of the active compounds, and compound with ethyl substitution at second position of indolizine nucleus showed the highest free binding energy of ΔG -24.11 (kcal/mol), a low clash score of 3.04, and high lipo score of -13.33. Indolizine analog with ethyl substitution at second position demonstrated Molecular Mechanics/Generalized Born Surface Area (-23.85 kcal/mol). Two molecular dynamics studies were computed (100 ps and 50 ns) to calculate the relationship between the potential and kinetic energies of the active anti-TB compound with time and temperature. The discovery of this lead may have a positive impact on anti-TB drug discovery.

Designing of inhibitors against CTX-M-15 type β-lactamase: potential drug candidate against β-lactamases-producing multi-drug-resistant bacteria

CTX-M-15 are the most prevalent types of β-lactamases that hydrolyze almost all antibiotics of β-lactam group lead to multiple-antibiotic resistance in bacteria. Three β-lactam inhibitors are available for use in combination with different antibiotics of cephalosporine group against the CTX-M-15-producing strains. Therefore, strategies to identify novel anti β-lactamase agents with specific mechanisms of action are the need of an hour. In this study, we screened three novel non-β-lactam inhibitors against CTX-M-15 by multi-step virtual screening approach. The potential for virtually screened drugs was estimated through in vitro cell assays. Hence, we proposed a study to understand the binding mode of CTX-M-15 with inhibitors using bioinformatics and experimental approach. We calculated the dissociation constants (K_d), association constant (K_a), stoichiometry (n) and binding energies (ΔG) of compounds with the respective targets. Molecular dynamic simulation carried out for 25 ns, revealed that these complexes were found stable throughout the simulation with relative RMSD in acceptable range. Moreover, microbiological and kinetic studies further confirmed high efficacies of these inhibitors by reducing the minimum inhibitory concentration (MIC) and catalysis of antibiotics by β-lactamases in the presence of inhibitors. Therefore, we conclude that these potential inhibitors may be used as a lead molecule for future drug candidates against β-lactamases-producing bacteria.