Dataset Information

Reaction pathways and free energy profiles for spontaneous hydrolysis of urea and tetramethylurea: unexpected substituent effects.

ABSTRACT: It has been difficult to directly measure the spontaneous hydrolysis rate of urea and, thus, 1,1,3,3-tetramethylurea (Me4U) was used as a model to determine the "experimental" rate constant for urea hydrolysis. The use of Me4U was based on an assumption that the rate of urea hydrolysis should be 2.8 times that of Me4U hydrolysis because the rate of acetamide hydrolysis is 2.8 times that of N,N-dimethyl-acetamide hydrolysis. The present first-principles electronic-structure calculations on the competing non-enzymatic hydrolysis pathways have demonstrated that the dominant pathway is the neutral hydrolysis via the CN addition for both urea (when pH < ~11.6) and Me4U (regardless of pH), unlike the non-enzymatic hydrolysis of amides where alkaline hydrolysis is dominant. Based on the computational data, the substituent shift of the free energy barrier calculated for the neutral hydrolysis is remarkably different from that for the alkaline hydrolysis, and the rate constant for the urea hydrolysis should be ~1.3 × 10(9)-fold lower than that (4.2 × 10(-12) s(-1)) measured for the Me4U hydrolysis. As a result, the rate enhancement and catalytic proficiency of urease should be 1.2 × 10(25) and 3 × 10(27) M(-1), respectively, suggesting that urease surpasses proteases and all other enzymes in its power to enhance the rate of reaction. All of the computational results are consistent with available experimental data for Me4U, suggesting that the computational prediction for urea is reliable.

SUBMITTER: Yao M

PROVIDER: S-EPMC3870011 | biostudies-literature | 2013 Nov

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Reaction pathways and free energy profiles for spontaneous hydrolysis of urea and tetramethylurea: unexpected substituent effects.

Yao Min M Tu Wenlong W Chen Xi X Zhan Chang-Guo CG

Organic & biomolecular chemistry 20131101 43

It has been difficult to directly measure the spontaneous hydrolysis rate of urea and, thus, 1,1,3,3-tetramethylurea (Me4U) was used as a model to determine the "experimental" rate constant for urea hydrolysis. The use of Me4U was based on an assumption that the rate of urea hydrolysis should be 2.8 times that of Me4U hydrolysis because the rate of acetamide hydrolysis is 2.8 times that of N,N-dimethyl-acetamide hydrolysis. The present first-principles electronic-structure calculations on the co ...[more]

PMID: 24097048

Similar Datasets

Project description:The nature and mechanism of phosphate hydrolysis reactions are of great interest in view of the crucial role of these reactions in key biological processes. Although it is becoming clearer that the ultimate way of resolving mechanistic controversies must involve reliable theoretical studies, it is not widely realized that such studies cannot be performed at present by using most existing automated ways and that only careful systematic studies can lead to meaningful conclusions. The present work clarifies the above point by considering the hydrolysis of phosphate monoesters. The clarification starts by defining the actual issues that should be addressed in careful studies and by highlighting the problems with studies that ignore the need for unique mechanistic definitions (e.g., works that confuse associative and dissociative pathways). We then focus on the analysis of the proton transfer (PT) pathways in phosphate hydrolysis and on recent suggestions that PT involves more than one water molecule. Here we point out that most of the studies that found a proton transfer through several water molecules have not involved a sufficient systematic search of the relevant reaction coordinates. This includes both energy minimization approaches as well as a recent metadynamics (MTD) simulation study. To illustrate the crucial need of exploring the potential surfaces reliably, rather than relying on automated approaches, we present here a very careful study of the free energy landscape along a 3D reaction coordinate (RC) exploring both the standard 2D RC, comprised of the attacking and leaving group reaction coordinates, as well as of the proton transfer (PT) coordinate. Our study points out that QM/MM minimization or MTD studies that concluded that the hydrolysis of phosphate monoesters involves a PT through several water molecules, have not explored carefully the single water (1W) path (that involves a direct PT form the attacking water molecule to the phosphate oxygen). Furthermore, we identified the most likely reason for the difficulty in finding the 1W path by QM/MM minimization methods, as well as by the current MTD simulations. We also discuss the problems with current studies that challenge the phosphate as a base mechanism and emphasize that all recent studies found associative/concerted paths (although many have not realized the meaning of their results). Finally, although we clearly do not have the last word about the 1W versus 2W paths we believe that we illustrated that the crucial mechanistic problems with alternative pathways should not be resolved by just running black box search approaches.

Project description:We report the results of calculations of the Gibbs energy profiles of the guanosine triphosphate (GTP) hydrolysis by the Arl3-RP2 protein complex using molecular dynamics (MD) simulations with ab initio type QM/MM potentials. The chemical reaction of GTP hydrolysis to guanosine diphosphate (GDP) and inorganic phosphate (Pi) is catalyzed by GTPases, the enzymes, which are responsible for signal transduction in live cells. A small GTPase Arl3, catalyzing the GTP → GDP reaction in complex with the activating protein RP2, constitute an essential part of the human vision cycle. To simulate the reaction mechanism, a model system is constructed by motifs of the crystal structure of the Arl3-RP2 complexed with a substrate analog. After selection of reaction coordinates, energy profiles for elementary steps along the reaction pathway GTP + H2O → GDP + Pi are computed using the umbrella sampling and umbrella integration procedures. QM/MM MD calculations are carried out, interfacing the molecular dynamics program NAMD and the quantum chemistry program TeraChem. Ab initio type QM(DFT)/MM potentials are computed with atom-centered basis sets 6-31G** and two hybrid functionals (PBE0-D3 and ωB97x-D3) of the density functional theory, describing a large QM subsystem. Results of these simulations of the reaction mechanism are compared to those obtained with QM/MM calculations on the potential energy surface using a similar description of the QM part. We find that both approaches, QM/MM and QM/MM MD, support the mechanism of GTP hydrolysis by GTPases, according to which the catalytic glutamine side chain (Gln71, in this system) actively participates in the reaction. Both approaches distinguish two parts of the reaction: the cleavage of the phosphorus-oxygen bond in GTP coupled with the formation of Pi, and the enzyme regeneration. Newly performed QM/MM MD simulations confirmed the profile predicted in the QM/MM minimum energy calculations, called here the pathway-I, and corrected its relief at the first elementary step from the enzyme-substrate complex. The QM/MM MD simulations also revealed another mechanism at the part of enzyme regeneration leading to pathway-II. Pathway-II is more consistent with the experimental kinetic data of the wild-type complex Arl3-RP2, whereas pathway-I explains the role of the mutation Glu138Gly in RP2 slowing down the hydrolysis rate.

Dataset Information

Reaction pathways and free energy profiles for spontaneous hydrolysis of urea and tetramethylurea: unexpected substituent effects.

Publications

Reaction pathways and free energy profiles for spontaneous hydrolysis of urea and tetramethylurea: unexpected substituent effects.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets