Structural basis for inhibition of a GH116 β-glucosidase and its missense mutants by GBA2 inhibitors: Crystallographic and quantum chemical study

The crystal structure of the Thermoanaerobacterium xylanolyticum in glycoside hydrolase family 116 (TxGH116) β-glucosidase provides a structural model for human GBA2 glucosylceramidase, an enzyme defective in hereditary spastic paraplegia and a potential therapeutic target for treating Gaucher disease. To assess the therapeutic potential of known inhibitors, the X-ray structure of TxGH116 in complex with isofagomine (IFG) was determined at 2.0 Å resolution and showed the IFG bound in a relaxed chair conformation. The binding of IFG and 7 other iminosugar inhibitors to wild-type and mutant enzymes (Asp508His and Arg786His) mimicking GBA2 pathogenic variants was then evaluated computationally by two-layered ONIOM calculations (at the B3LYP:PM7 level). Calculations showed that six charged residues, Glu441, Asp452, His507, Asp593, Glu777, and Arg786 influence inhibitor binding most. His507, Glu777 and Arg786, form strong hydrogen bonds with the inhibitors (∼1.4-1.6 Å). Thus, the missense mutation of one of these residues in Arg786His has a greater effect on the interaction energies for all inhibitors compared to Asp508His. In line with the experimental data for the inhibitors that have been tested, the favorable interaction energy between the inhibitors and the TxGH116 protein followed the trend: isofagomine > 1-deoxynojirimycin > glucoimidazole > N-butyl-deoxynojirimycin ≈ N-nonyl-deoxynojirimycin > conduritol B epoxide ≈ azepane 1 > azepane 2. The obtained structural and energetic properties and comparison to the GBA2 model can lead to understanding of structural requirement for inhibitor binding in GH116 to aid the design of high potency GBA2 inhibitors.

Systematic Functional and Computational Analysis of Glucose-Binding Residues in Glycoside Hydrolase Family GH116

Glycoside hydrolases (GH) bind tightly to the sugar moiety at the glycosidic bond being hydrolyzed to stabilize its transition state conformation. We endeavored to assess the importance of glucose-binding residues in GH family 116 (GH116) β-glucosidases, which include human β-glucosylceramidase 2 (GBA2), by mutagenesis followed by kinetic characterization, X-ray crystallography, and ONIOM calculations on Thermoanaerobacterium xylanolyticum TxGH116, the structural model for GH116 enzymes. Mutations of residues that bind at the glucose C3OH and C4OH caused 27–196-fold increases in KM for p-nitrophenyl-β-D-glucoside, and significant decreases in the kcat, up to 5000-fold. At the C6OH binding residues, mutations of E777 decreased the kcat/KM by over 60,000-fold, while R786 mutants increased both the KM (40-fold) and kcat (2–4-fold). The crystal structures of R786A and R786K suggested a larger entrance to the active site could facilitate their faster rates. ONIOM binding energy calculations identified D452, H507, E777, and R786, along with the catalytic residues E441 and D593, as strong electrostatic contributors to glucose binding with predicted interaction energies > 15 kcal mol−1, consistent with the effects of the D452, H507, E777 and R786 mutations on enzyme kinetics. The relative importance of GH116 active site residues in substrate binding and catalysis identified in this work improves the prospects for the design of inhibitors for GBA2 and the engineering of GH116 enzymes for hydrolytic and synthetic applications.

Computational Screening of Phenylamino-Phenoxy-Quinoline Derivatives against the Main Protease of SARS-CoV-2 Using Molecular Docking and the ONIOM Method

In the search for new anti-HIV-1 agents, two forms of phenylamino-phenoxy-quinoline derivatives have been synthesized, namely, 2-phenylamino-4-phenoxy-quinoline and 6-phenylamino-4-phenoxy-quinoline. In this study, the binding interactions of phenylamino-phenoxy-quinoline derivatives and six commercially available drugs (hydroxychloroquine, ritonavir, remdesivir, S-217622, N3, and PF-07321332) with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) main protease (Mpro) were investigated using molecular docking and the ONIOM method. The molecular docking showed the hydrogen bonding and hydrophobic interactions of all the compounds in the pocket of SARS-CoV-2 main protease (Mpro), which plays an important role for the division and proliferation of the virus into the cell. The binding free energy values between the ligands and Mpro ranged from -7.06 to -10.61 kcal/mol. The molecular docking and ONIOM results suggested that 4-(2',6'-dimethyl-4'-cyanophenoxy)-2-(4″-cyanophenyl)-aminoquinoline and 4-(4'-cyanophenoxy)-2-(4″-cyanophenyl)-aminoquinoline have low binding energy values and appropriate molecular properties; moreover, both compounds could bind to Mpro via hydrogen bonding and Pi-Pi stacking interactions with amino acid residues, namely, HIS41, GLU166, and GLN192. These amino acids are related to the proteolytic cleavage process of the catalytic triad mechanisms. Therefore, this study provides important information for further studies on synthetic quinoline derivatives as antiviral candidates in the treatment of SARS-CoV-2.

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.

The Driving Force for the Acylation of β-Lactam Antibiotics by L,D-Transpeptidase 2: Quantum Mechanics/Molecular Mechanics (QM/MM) Study

β-lactam antibiotics, which are used to treat infectious diseases, are currently the most widely used class of antibiotics. This study focused on the chemical reactivity of five- and six-membered ring systems attached to the β-lactam ring. The ring strain energy (RSE), force constant (FC) of amide (C-N), acylation transition states and second-order perturbation stabilization energies of 13 basic structural units of β-lactam derivatives were computed using the M06-2X and G3/B3LYP multistep method. In the ring strain calculations, an isodesmic reaction scheme was used to obtain the total energies. RSE is relatively greater in the five-(1a-2c) compared to the six-membered ring systems except for 4b, which gives a RSE that is comparable to five-membered ring lactams. These variations were also observed in the calculated inter-atomic amide bond distances (C-N), which is why the six-membered ring lactams C-N bond are more rigid than those with five-membered ring lactams. The calculated ΔG^# values from the acylation reaction of the lactams (involving the S-H group of the cysteine active residue from L,D transpeptidase 2) revealed a faster rate of C-N cleavage in the five-membered ring lactams especially in the 1-2 derivatives (17.58 kcal mol^-1 ). This observation is also reflected in the calculated amide bond force constant (1.26 mDyn/A) indicating a weaker bond strength, suggesting that electronic factors (electron delocalization) play more of a role on reactivity of the β-lactam ring, than ring strain.

Torsional flexibility of undecorated catechol diether compound as potent NNRTI targeting HIV-1 reverse transcriptase

Conformational adaptation of non-nucleoside reverse transcriptase inhibitor (NNRTI) via torsional flexibility is found to be very significant for targeting human immunodeficiency virus type 1 reverse transcriptase (HIV-1 RT) mutants. Catechol diether derivative including flexible torsions is new potent NNRTI with picomolar activity. Moreover, this derivative also reveals the good solubility, low toxicity and potent inhibition for HIV-1 mutants. In this study, torsional flexibility of an undecorated catechol diether compound in the binding pocket of wild type and mutants (Y181C and K103N/Y181C) HIV-1 RT is investigated by using QM/MM calculations. From the results, the uracil ring is found to exhibit more flexibility in the NNIBP. On the contrary, potential energy surfaces show that high energy is encountered by changing of the corresponding torsion of the cyanovinyl aryl ring indicating the limitation for torsional flexibility. For pointing out the key interaction for the binding, the residual interaction energies are performed by means of QM calculations. Important attractive interactions through hydrogen bonds between the inhibitor and K102, K/N103, V106, and Y188 are observed. The catechol ring is proposed to be modified in order to strengthen interactions with surrounding amino acids. The results may help for the designing of new potent NNRTIs.

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

Mycobacterium tuberculosis is the causative agent of Tuberculosis. Formation of 3 → 3 crosslinks in the peptidoglycan layer of M. tuberculosis is catalyzed by l,d-transpeptidases. These enzymes can confer resistance against classical β-lactams that inhibit enzymes that generate 4 → 3 peptidoglycan crosslinks. The focus of this study is to investigate the catalytic role of water molecules in the acylation mechanism of the β-lactam ring within two models; 4- and 6-membered ring systems using two-layered our Own N-layer integrated Molecular Mechanics ONIOM (B3LYP/6-311++G(2d,2p): AMBER) model. The obtained thermochemical parameters revealed that the 6-membered ring model best describes the inhibition mechanism of acylation which indicates the role of water in the preference of 6-membered ring reaction pathway. This finding is in accordance with experimental data for the rate-limiting step of cysteine protease with the same class of inhibitor and binding affinity for both inhibitors. As expected, the ΔG^# results also reveal that the 6-membered ring reaction pathway is the most favourable. The electrostatic potential (ESP) and the natural bond orbital analysis (NBO) showed stronger interactions in 6-membered ring transition state (TS-6) mechanism involving water in the active site of the enzyme. This study could be helpful in the development of novel antibiotics against l,d-transpeptidase.

A ligand-free strategy for the copper-catalysed direct alkynylation of trifluoromethyl ketones

A superbase-induced copper-catalyzed direct alkynylation method for trifluoromethyl ketones was developed without any additional ligand under mild conditions.