Multiscale Computational Study on the Catalytic Mechanism of the Nonmetallo Amidase Maleamate Amidohydrolase (NicF).

Maleamate amidohydrolase (NicF) is a key enzyme in vitamin B3 metabolism that catalyzes the hydrolysis of maleamate to produce maleic acid and ammonia. Unlike most members from the amidohydrolase superfamily it does not require a metal ion. Here, we use multiscale computational enzymology to investigate the catalytic mechanism, substrate binding, oxyanion hole, and roles of key active site residues of NicF from Bordetella bronchiseptica. In particular, molecular dynamics (MD) simulations, quantum mechanics/molecular mechanics (QM/MM) and QTAIM methods have been applied. The mechanism of the NicF-catalyzed reaction proceeds by a nucleophilic addition-elimination sequence involving the formation of a thioester enzyme intermediate (IC2 in stage 1) followed by hydrolysis of the thioester bond to form the products (stage 2). Consequently, the formation of IC2 in stage 1 is the rate-limiting step with a barrier of 88.8 kJ·mol-1 relative to the reactant complex, RC. Comparisons with related metal-dependent enzymes, particularly the zinc-dependent nicotinamidase from Streptococcus pneumonia (SpNic), have also been made to further illustrate unique features of the present mechanism. Along with -NH- donor groups of the oxyanion hole (i.e., HN-Thr146, HN-Cys150), the active site ß-hydroxyl of threonine (HO-ßThr146) is concluded to play a role in stabilizing the carbonyl oxygen of maleamate during the mechanism.

Ketamine Inhibition of the Pentameric Ligand-Gated Ion Channel GLIC.

Ketamine inhibits pentameric ligand-gated ion channels (pLGICs), including the bacterial pLGIC from Gloeobacter violaceus (GLIC). The crystal structure of GLIC shows R-ketamine bound to an extracellular intersubunit cavity. Here, we performed molecular dynamics simulations of GLIC in the absence and presence of R- or S-ketamine. No stable binding of S-ketamine in the original cavity was observed in the simulations, largely due to its unfavorable access to residue D154, which provides important electrostatic interactions to stabilize R-ketamine binding. Contrary to the symmetric binding shown in the crystal structure, R-ketamine moved away from some of the binding sites and was bound to GLIC asymmetrically at the end of simulations. The asymmetric binding is consistent with the experimentally measured negative cooperativity of ketamine binding to GLIC. In the presence of R-ketamine, all subunits showed changes in structure and dynamics, irrespective of binding stability; the extracellular intersubunit cavity expanded and intersubunit electrostatic interactions involved in channel activation were altered. R-ketamine binding promoted a conformational shift toward closed GLIC. Conformational changes near the ketamine-binding site were propagated to the interface between the extracellular and transmembrane domains, and further to the pore-lining TM2 through two pathways: pre-TM1 and the ß1-ß2 loop. Both signaling pathways have been predicted previously using the perturbation-based Markovian transmission model. The study provides a structural and dynamics basis for the inhibitory modulation of ketamine on pLGICs.

Computational investigations on the catalytic mechanism of maleate isomerase: the role of the active site cysteine residues.

The maleate isomerase (MI) catalysed isomerization of maleate to fumarate has been investigated using a wide range of computational modelling techniques, including small model DFT calculations, QM-cluster approach, quantum mechanical/molecular mechanical approach (QM/MM in the ONIOM formalism) and molecular dynamics simulations. Several fundamental questions regarding the mechanism were answered in detail, such as the activation and stabilization of the catalytic Cys in a rather hydrophobic active site. The two previously proposed mechanisms were considered, where either enediolate or succinyl-Cys intermediate forms. Small model calculations as well as an ONIOM-based approach suggest that an enediolate intermediate is too unstable. Furthermore, the formation of succinyl-Cys intermediate via the nucleophilic attack of Cys76(-) on the substrate C2 (as proposed experimentally) was found to be energetically unfeasible in both QM-cluster and ONIOM approaches. Instead, our results show that Cys194, upon activation via the substrate, acts as a nucleophile and Cys76 acts as an acid/base catalyst, forming a succinyl-Cys intermediate in a concerted fashion. Indeed, the calculated PA of Cys76 is always higher than that of Cys194 before or upon substrate binding in the active site. Furthermore, the mechanism proceeds via multiple steps by substrate rotation around C2-C3 with the assistance of the now negatively charged Cys76, leading to the formation of fumarate. Finally, our calculated barrier is in good agreement with experiment. These findings represent a novel mechanism in the racemase superfamily.

A multi-scale computational study on the mechanism of Streptococcus pneumoniae Nicotinamidase (SpNic).

Nicotinamidase (Nic) is a key zinc-dependent enzyme in NAD metabolism that catalyzes the hydrolysis of nicotinamide to give nicotinic acid. A multi-scale computational approach has been used to investigate the catalytic mechanism, substrate binding and roles of active site residues of Nic from Streptococcus pneumoniae (SpNic). In particular, density functional theory (DFT), molecular dynamics (MD) and ONIOM quantum mechanics/molecular mechanics (QM/MM) methods have been employed. The overall mechanism occurs in two stages: (i) formation of a thioester enzyme-intermediate (IC2) and (ii) hydrolysis of the thioester bond to give the products. The polar protein environment has a significant effect in stabilizing reaction intermediates and in particular transition states. As a result, both stages effectively occur in one step with Stage 1, formation of IC2, being rate limiting barrier with a cost of 53.5 kJ·mol-1 with respect to the reactant complex, RC. The effects of dispersion interactions on the overall mechanism were also considered but were generally calculated to have less significant effects with the overall mechanism being unchanged. In addition, the active site lysyl (Lys103) is concluded to likely play a role in stabilizing the thiolate of Cys136 during the reaction.

QM/MM investigation of the discriminatory pre-transfer editing mechanism operated by Lysyl-tRNA synthetase.

Aminoacyl-tRNA synthetases (aaRSs) are essential enzymes that remarkable facilitate the aminoacylation process during translation. With a high fidelity, the mischarged tRNA is prevented through implementing pre- and post-transfer proofreading mechanisms. For instance, Lysine-tRNA synthetase charges the native substrate, lysine, to its cognate tRNA. In spite of the great structural similarity between lysine to the noncognate and toxic ornithine, with the side chain of lysine being only one methylene group longer, LysRS is able to achieve this discrimination with a high efficiency. In this work, the hybrid quantum mechanics/molecular mechanics (QM/MM) investigation was applied to probe the pre-transfer editing mechanism catalyzed by lysyl-tRNA synthetase to reject the noncognte aminoacyl, L-ornityl (Orn), compared to the cognate substrate, L-lysyl. Particularly, the self-cyclization pre-transfer editing mechanism was explored for the two substrates. The substrate-assisted self-cyclization editing of Orn-AMP, where its phosphate moiety acts as the catalytic base, is found to be the rate-determining step with an energy barrier of 101.2 kJ mol-1. Meanwhile, the corresponding rate-limiting pathway for the native Lys-AMP lies at 140.2 kJ mol-1. This observation clearly indicated the infeasibility of this catalytic scenario in the presence of the native substrate. Interestingly, a thermodynamically favorable cyclic product of -92.9 kJ mol-1 with respect to the aminoacyl reactant complex demonstrated evidence of a successful pre-transfer editing. This reaction resulted in the discharge of the on-cognate -ornithine derivative from LysU's active site. These valuable mechanistic insights are valuable to enrich our knowledge of this extremely efficient and specific catalytic machinery of LysRS.Communicated by Ramaswamy H. Sarma.

A Molecular Dynamics (MD) and Quantum Mechanics/Molecular Mechanics (QM/MM) study on Ornithine Cyclodeaminase (OCD): a tale of two iminiums.

Ornithine cyclodeaminase (OCD) is an NAD+-dependent deaminase that is found in bacterial species such as Pseudomonas putida. Importantly, it catalyzes the direct conversion of the amino acid L-ornithine to L-proline. Using molecular dynamics (MD) and a hybrid quantum mechanics/molecular mechanics (QM/MM) method in the ONIOM formalism, the catalytic mechanism of OCD has been examined. The rate limiting step is calculated to be the initial step in the overall mechanism: hydride transfer from the L-ornithine's C(α)-H group to the NAD+ cofactor with concomitant formation of a C(α)=NH(2)+ Schiff base with a barrier of 90.6 kJ mol-1. Importantly, no water is observed within the active site during the MD simulations suitably positioned to hydrolyze the C(α)=NH(2)+ intermediate to form the corresponding carbonyl. Instead, the reaction proceeds via a non-hydrolytic mechanism involving direct nucleophilic attack of the δ-amine at the C(α)-position. This is then followed by cleavage and loss of the α-NH(2) group to give the Δ1-pyrroline-2-carboxylate that is subsequently reduced to L-proline.