Dataset Information

The QM/MM molecular dynamics and free energy simulations of the acylation reaction catalyzed by the serine-carboxyl peptidase kumamolisin-As.

ABSTRACT: Quantum mechanical/molecular mechanical molecular dynamics and free energy simulations are performed to study the acylation reaction catalyzed by kumamolisin-As, a serine-carboxyl peptidase, and to elucidate the catalytic mechanism and the origin of substrate specificity. It is demonstrated that the nucleophilic attack by the serine residue on the substrate may not be the rate-limiting step for the acylation of the GPH*FF substrate. The present study also confirms the earlier suggestions that Asp164 acts as a general acid during the catalysis and that the electrostatic oxyanion hole interactions may not be sufficient to lead a stable tetrahedral intermediate along the reaction pathway. Moreover, Asp164 is found to act as a general base during the formation of the acyl-enzyme from the tetrahedral intermediate. The role of dynamic substrate assisted catalysis (DSAC) involving His at the P1 site of the substrate is examined for the acylation reaction. It is demonstrated that the bond-breaking and -making events at each stage of the reaction trigger a change of the position for the His side chain and lead to the formation of the alternative hydrogen bonds. The back and forth movements of the His side chain between the C=O group of Pro at P2 and Odelta2 of Asp164 in a ping-pong-like mechanism and the formation of the alternative hydrogen bonds effectively lower the free energy barriers for both the nucleophilic attack and the acyl-enzyme formation and may therefore contribute to the relatively high activity of kumamolisin-As toward the substrates with His at the P1 site.

SUBMITTER: Xu Q

PROVIDER: S-EPMC2533263 | biostudies-literature | 2007 Mar

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

The QM/MM molecular dynamics and free energy simulations of the acylation reaction catalyzed by the serine-carboxyl peptidase kumamolisin-As.

Xu Qin Q Guo Hao-Bo HB Wlodawer Alexander A Nakayama Toru T Guo Hong H

Biochemistry 20070228 12

Quantum mechanical/molecular mechanical molecular dynamics and free energy simulations are performed to study the acylation reaction catalyzed by kumamolisin-As, a serine-carboxyl peptidase, and to elucidate the catalytic mechanism and the origin of substrate specificity. It is demonstrated that the nucleophilic attack by the serine residue on the substrate may not be the rate-limiting step for the acylation of the GPH*FF substrate. The present study also confirms the earlier suggestions that As ...[more]

PMID: 17326662

Similar Datasets

Project description:Quantum mechanical/molecular mechanical (QM/MM) free energy simulations are applied for understanding the mechanism of the acylation reaction catalyzed by sedolisin, a representative serine-carboxyl peptidase, leading to the acyl-enzyme (AE) and first product from the enzyme-catalyzed reaction. One of the interesting questions to be addressed in this work is the origin of the substrate specificity of sedolisin that shows a relatively high activity on the substrates with Glu at P(1) site. It is shown that the bond making and breaking events of the acylation reaction involving a peptide substrate (LLE*FL) seem to be accompanied by local conformational changes, proton transfers as well as the formation of alternative hydrogen bonds. The results of the simulations indicate that the conformational change of Glu at P(1) site and its formation of a low barrier hydrogen bond with Asp-170 (along with the transient proton transfer) during the acylation reaction might play a role in the relatively high specificity for the substrate with Glu at P(1) site. The role of some key residues in the catalysis is confirmed through free energy simulations. Glu-80 is found to act as a general base to accept a proton from Ser-287 during the nucleophilic attack and then as a general acid to protonate the leaving group (N-H of P(1')-Phe) during the cleavage of the scissile peptide bond. Another acidic residue, Asp-170, acts as a general acid catalyst to protonate the carbonyl of P(1)-Glu during the formation of the tetrahedral intermediate and as a general base for the formation of the acyl-enzyme. The energetic results from the free energy simulations support the importance of proton transfer from Asp-170 to the carbonyl of P(1)-Glu in the stabilization of the tetrahedral intermediate and the formation of a low-barrier hydrogen bond between the carboxyl group of P(1)-Glu and Asp-170 in the lowering of the free energy barrier for the cleavage of the peptide bond. Detailed analyses of the proton transfers during acylation are also given.

Project description:Human arginase I (HARGI) is a metalloprotein highly expressed in the liver cytosol and catalyzes the hydrolysis of l-arginine to form l-ornithine and urea. Understanding the reaction mechanism would be highly helpful to design new inhibitor molecules for HARGI as it is a target for heart- and blood-related diseases. In this study, we explored the hydrolysis reaction mechanism of HARGI with antiferromagnetic and ferromagnetic coupling between two Mn(II) ions at the catalytic site by employing molecular dynamics simulations coupled with quantum mechanics and molecular mechanics (QM/MM). The spin states, high-spin ferromagnetic couple (S Mn1 = 5/2, S Mn2 = 5/2), low-spin ferromagnetic couple (S Mn1 = 1/2, S Mn2 = 1/2), high-spin antiferromagnetic couple (S Mn1 = 5/2, S Mn2 = -5/2), and low-spin antiferromagnetic couple (S Mn1 = 1/2, S Mn2 = -1/2) are considered, and the calculated energetics for the complex of the substrate and HARGI are compared. The results show that the high-spin antiferromagnetic couple (S Mn1 = 5/2, S Mn2 = -5/2) is more stable than other spin states. The low-spin ferromagnetic and antiferromagnetic coupled states are highly unstable compared with the corresponding high-spin states. The high-spin antiferromagnetic couple (S Mn1 = 5/2, S Mn2 = -5/2) is stabilized by 0.39 kcal/mol compared with the ferromagnetic couple (S Mn1 = 5/2, S Mn2 = 5/2). The reaction mechanism is independent of spin states; however, the energetics of transition states and intermediates are more stable in the case of the high-spin antiferromagnetic couple (S Mn1 = 5/2, S Mn2 = -5/2) than the corresponding ferromagnetic state. It is evident that the calculated coupling constants are higher for antiferromagnetic states and, interestingly, superexchange coupling is found to occur between Mn(II) ions via hydroxide ions in a reactant. The hydroxide ion enhances the coupling interaction and initiates the catalytic reaction. It is also noted that the first intermediate structure where there is no superexchange coupling is similar to the known inhibitor 2(S)-amino-6-boronohexanoic acid.

Project description:THE RELIABILITY OF FREE ENERGY SIMULATIONS (FES) IS LIMITED BY TWO FACTORS: (a) the need for correct sampling and (b) the accuracy of the computational method employed. Classical methods (e.g., force fields) are typically used for FES and present a myriad of challenges, with parametrization being a principle one. On the other hand, parameter-free quantum mechanical (QM) methods tend to be too computationally expensive for adequate sampling. One widely used approach is a combination of methods, where the free energy difference between the two end states is computed by, e.g., molecular mechanics (MM), and the end states are corrected by more accurate methods, such as QM or hybrid QM/MM techniques. Here we report two new approaches that significantly improve the aforementioned scheme; with a focus on how to compute corrections between, e.g., the MM and the more accurate QM calculations. First, a molecular dynamics trajectory that properly samples relevant conformational degrees of freedom is generated. Next, potential energies of each trajectory frame are generated with a QM or QM/MM Hamiltonian. Free energy differences are then calculated based on the QM or QM/MM energies using either a non-Boltzmann Bennett approach (QM-NBB) or non-Boltzmann free energy perturbation (NB-FEP). Both approaches are applied to calculate relative and absolute solvation free energies in explicit and implicit solvent environments. Solvation free energy differences (relative and absolute) between ethane and methanol in explicit solvent are used as the initial test case for QM-NBB. Next, implicit solvent methods are employed in conjunction with both QM-NBB and NB-FEP to compute absolute solvation free energies for 21 compounds. These compounds range from small molecules such as ethane and methanol to fairly large, flexible solutes, such as triacetyl glycerol. Several technical aspects were investigated. Ultimately some best practices are suggested for improving methods that seek to connect MM to QM (or QM/MM) levels of theory in FES.

Dataset Information

The QM/MM molecular dynamics and free energy simulations of the acylation reaction catalyzed by the serine-carboxyl peptidase kumamolisin-As.

Publications

The QM/MM molecular dynamics and free energy simulations of the acylation reaction catalyzed by the serine-carboxyl peptidase kumamolisin-As.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets