The catalytic role of water in the binding site of l,d-transpeptidase 2 within acylation mechanism: A QM/MM (ONIOM) modelling

Investigating the biological actions of some Schiff bases using density functional theory study

Schiff base ligands have wide varieties of application in several fields. One of which is the biological actions they possess such as anti-fungal, anti-bacterial, anti-malarial, and anti-viral characteristics. In this study, some synthesized phenylimino-based Schiff bases were investigated using density functional theory (DFT) to unravel their biological descriptors. The gas-phase quantum chemical calculation was done on the Schiff base 3-((E)-(phenylimino)methyl)benzene-1,2-diol and other synthesized analogues to evaluate their reactivity and stability properties including the substituent effect on the basic molecule. The Coulomb-attenuating method (CAM-B3LYP) functional was employed for the theoretical calculations. The Nuclear Magnetic Resonance (NMR), Fourier Transform-Infrared (FT-IR), Ultraviolet/visible spectroscopies calculated agrees with the experimental values. The obtained charge transfer and electronic features provide useful information regarding the active sites for biological application in the compounds.

Mechanisms of Proteolytic Enzymes and Their Inhibition in QM/MM Studies

Experimental evidence for enzymatic mechanisms is often scarce, and in many cases inadvertently biased by the employed methods. Thus, apparently contradictory model mechanisms can result in decade long discussions about the correct interpretation of data and the true theory behind it. However, often such opposing views turn out to be special cases of a more comprehensive and superior concept. Molecular dynamics (MD) and the more advanced molecular mechanical and quantum mechanical approach (QM/MM) provide a relatively consistent framework to treat enzymatic mechanisms, in particular, the activity of proteolytic enzymes. In line with this, computational chemistry based on experimental structures came up with studies on all major protease classes in recent years; examples of aspartic, metallo-, cysteine, serine, and threonine protease mechanisms are well founded on corresponding standards. In addition, experimental evidence from enzyme kinetics, structural research, and various other methods supports the described calculated mechanisms. One step beyond is the application of this information to the design of new and powerful inhibitors of disease-related enzymes, such as the HIV protease. In this overview, a few examples demonstrate the high potential of the QM/MM approach for sophisticated pharmaceutical compound design and supporting functions in the analysis of biomolecular structures.

Mechanistic insight on the inhibition of D, D-carboxypeptidase from Mycobacterium tuberculosis by β-lactam antibiotics: an ONIOM acylation study

Mycobacterium tuberculosis cell wall is intricate and impermeable to many agents. A D, D-carboxypeptidase (DacB1) is one of the enzymes involved in the biosynthesis of cell wall peptidoglycan and catalyzes the terminal D-alanine cleavage from pentapeptide precursors. Catalytic activity and mechanism by which DacB1 functions is poorly understood. Herein, we investigated the acylation mechanism of DacB1 by β-lactams using a 6-membered ring transition state model that involves a catalytic water molecule in the reaction pathway. The full transition states (TS) optimization plus frequency were achieved using the ONIOM (B3LYP/6-31 + G(d): AMBER) method. Subsequently, the activation free energies were computed via single-point calculations on fully optimized structures using B3LYP/6-311++(d,p): AMBER and M06-2X/6-311++(d,p): AMBER with an electronic embedding scheme. The 6-membered ring transition state is an effective model to examine the inactivation of DacB1 via acylation by β-lactams antibiotics (imipenem, meropenem, and faropenem) in the presence of the catalytic water. The ΔG^# values obtained suggest that the nucleophilic attack on the carbonyl carbon is the rate-limiting step with 13.62, 19.60 and 30.29 kcal mol^-1 for Imi-DacB1, Mero-DacB1 and Faro-DacB1, respectively. The electrostatic potential (ESP) and natural bond orbital (NBO) analysis provided significant electronic details of the electron-rich region and charge delocalization, respectively, based on the concerted 6-membered ring transition state. The stabilization energies of charge transfer within the catalytic reaction pathway concurred with the obtained activation free energies. The outcomes of this study provide important molecular insight into the inactivation of D, D-carboxypeptidase by β-lactams.

Synthesis and recycling of the mycobacterial cell envelope

•

The unique mycobacterial cell wall is a layered structure of carbohydrates and lipids.

•

The architecture and biosynthesis of the cell wall have largely been elucidated.

•

Increasing evidence indicates that each cell wall layer is remodelled and recycled.

•

There are opportunities to discover new essential enzymes in cell wall metabolism.

•

Cell wall metabolism is a validated source of targets for tuberculosis drug discovery.

Mycobacterium tuberculosis (Mtb), the causative agent of the disease tuberculosis, is a recognised global health concern. The efficacy of the current treatment regime is under threat due to the emergence of antibiotic resistance, directing an urgent requirement for the discovery of new anti-tubercular agents and drug targets. The mycobacterial cell wall is a well-validated drug target for Mtb and is composed of three adaptive macromolecular structures, peptidoglycan, arabinogalactan and mycolic acids, an array of complex lipids and carbohydrates. The majority of the enzymes involved in cell wall synthesis have been established, whilst studies directed towards the mechanisms of remodelling and recycling have been neglected. This review briefly describes mycobacterial cell wall synthesis, and focuses on aspects of remodelling and recycling, thus highlighting opportunities for future research.

Demystifying the catalytic pathway of Mycobacterium tuberculosis isocitrate lyase

Pulmonary tuberculosis, caused by Mycobacterium tuberculosis, is one of the most persistent diseases leading to death in humans. As one of the key targets during the latent/dormant stage of M. tuberculosis, isocitrate lyase (ICL) has been a subject of interest for new tuberculosis therapeutics. In this work, the cleavage of the isocitrate by M. tuberculosis ICL was studied using quantum mechanics/molecular mechanics method at M06-2X/6-31+G(d,p): AMBER level of theory. The electronic embedding approach was applied to provide a better depiction of electrostatic interactions between MM and QM regions. Two possible pathways (pathway I that involves Asp108 and pathway II that involves Glu182) that could lead to the metabolism of isocitrate was studied in this study. The results suggested that the core residues involved in isocitrate catalytic cleavage mechanism are Asp108, Cys191 and Arg228. A water molecule bonded to Mg²⁺ acts as the catalytic base for the deprotonation of isocitrate C(2)–OH group, while Cys191 acts as the catalytic acid. Our observation suggests that the shuttle proton from isocitrate hydroxyl group C(2) atom is favourably transferred to Asp108 instead of Glu182 with a lower activation energy of 6.2 kcal/mol. Natural bond analysis also demonstrated that pathway I involving the transfer of proton to Asp108 has a higher intermolecular interaction and charge transfer that were associated with higher stabilization energy. The QM/MM transition state stepwise catalytic mechanism of ICL agrees with the in vitro enzymatic assay whereby Asp108Ala and Cys191Ser ICL mutants lost their isocitrate cleavage activities.

Concerted hydrolysis mechanism of HIV-1 natural substrate against subtypes B and C-SA PR: insight through molecular dynamics and hybrid QM/MM studies

It is well known that understanding the catalytic mechanism of HIV-1 PR is the rationale on which its inhibitors were developed; therefore, a better understanding of the mechanism of natural substrate hydrolysis is important. Herein, the reaction mechanism of HIV-1 natural substrates with subtypes B and common mutant in South Africa (subtype C-SA) protease were studied through transition state modelling, using a general acid-general base (GA-GB) one-step concerted process. The activation free energies of enzyme-substrate complexes were compared based on their rate of hydrolysis using a two-layered ONIOM (B3LYP/6-31++G(d,p):AMBER) method. We expanded our computational model to obtain a better understanding of the mechanism of hydrolysis as well as how the enzyme recognises or chooses the cleavage site of the scissile bonds. Using this model, a potential substrate-based inhibitor could be developed with better potency. The calculated activation energies of natural substrates in our previous study correlated well with experimental data. A similar trend was observed for the Gag and Gag-Pol natural substrates in the present work for both enzyme complexes except for the PR-RT substrate. Natural bond orbital (NBO) analysis was also applied to determine the extent of charge transfer within the QM part of both enzymes considered and the PR-RT natural substrate. The result of this study shows that the method can be utilized as a dependable computational technique to rationalize lead compounds against specific targets.

Revealing electronic features governing hydrolysis of cephalosporins in the active site of the L1 metallo-β-lactamase

The QM/MM simulations followed by electron density feature analysis are carried out to deepen the understanding of the reaction mechanism of cephalosporin hydrolysis in the active site of the L1 metallo-β-lactamase. The differences in reactivity of ten similar cephalosporin compounds are explained by using an extended set of bonding descriptors. The limiting step of the reaction is characterized by the proton transfer to the nitrogen atom of the cephalosporin thiazine ring accompanied with formation of the C₄ Created by potrace 1.16, written by Peter Selinger 2001-2019 ]]> C₃ double bond in its N–C₄–C₃ fragment. The temporary N⋯H–O_w hydrogen bond, which is formed in the transition state of the limiting step of the reaction was recognized as a key atomic interaction governing the reactivity of various cephalosporins. Non-local real-space bonding descriptors show that different extent of localization of electron lone pair at N atom in the transition state affect the reactivity of compounds: smaller electron localization is typical for the less reactive species. In particular, the Fermi hole analysis shows how exchange electron correlation in the N⋯H–O_w fragment control electron lone pair localization. Delocalization tensor, linear response kernel and source function indicate that features of electron delocalization in the N–C₄–C₃ fragment of cephalosporins in the transition state complexes determine the differences in C₄–C₃ bond for substrates with high and low rate constants. The C₄–C₃ bond of the N–C₄–C₃ fragment at the transition state is similar to that of the preceding intermediate for the less reactive species and resembles the features of the enzyme–product complex for more reactive compounds. The power and limitations of the descriptors applied for solving the problem are discussed and the generality of approach is stressed.

Combination of QM/MM and modern bonding descriptors explains different reactivity of cephalosporins in the active site of the L1 metallo-β-lactamase.

Inhibition mechanism of L,D-transpeptidase 5 in presence of the β-lactams using ONIOM method

Tuberculosis (TB) is one of the world's deadliest diseases resulting from infection by the bacterium, Mycobacterium tuberculosis (M.tb). The L,D-transpeptidase enzymes catalyze the synthesis of 3 → 3 transpeptide linkages which are predominant in the peptidoglycan of the M.tb cell wall. Carbapenems is class of β-lactams that inactivate L,D-transpeptidases by acylation, although differences in antibiotic side chains modulate drug binding and acylation rates. Herein, we used a two-layered our Own N-layer integrated Molecular Mechanics ONIOM method to investigate the catalytic mechanism of L,D-transpeptidase 5 (Ldt_Mt5) by β-lactam derivatives. Ldt_Mt5 complexes with six β-lactams, ZINC03788344 (1), ZINC02462884 (2), ZINC03791246 (3), ZINC03808351 (4), ZINC03784242 (5) and ZINC02475683 (6) were simulated. The QM region (high-level) comprises the β-lactam, one water molecule and the Cys360 catalytic residue, while the rest of the Ldt_Mt5 residues were treated with AMBER force field. The activation energies (ΔG^#) were calculated with B3LYP, M06-2X and ωB97X density functionals with 6-311++G(2d, 2p) basis set. The ΔG^# for the acylation of Ldt_Mt5 by the selected β-lactams were obtained as 13.67, 20.90, 22.88, 24.29, 27.86 and 28.26 kcal mol^-1respectively. Several of the compounds showed an improved ΔG^# when compared to the previously calculated energies for imipenem and meropenem for the acylation step for Ldt_Mt5. This model provides further validation of the catalytic inhibition mechanism of LDTs with atomistic detail.