Dataset Information

Ab initio quantum mechanical/molecular mechanical molecular dynamics simulation of enzyme catalysis: the case of histone lysine methyltransferase SET7/9.

ABSTRACT: To elucidate enzyme catalysis through computer simulation, a prerequisite is to reliably compute free energy barriers for both enzyme and solution reactions. By employing on-the-fly Born-Oppenheimer molecular dynamics simulations with the ab initio quantum mechanical/molecular mechanical approach and the umbrella sampling method, we have determined free energy profiles for the methyl-transfer reaction catalyzed by the histone lysine methyltransferase SET7/9 and its corresponding uncatalyzed reaction in aqueous solution, respectively. Our calculated activation free energy barrier for the enzyme catalyzed reaction is 22.5 kcal/mol, which agrees very well with the experimental value of 20.9 kcal/mol. The difference in potential of mean force between a corresponding prereaction state and the transition state for the solution reaction is computed to be 30.9 kcal/mol. Thus, our simulations indicate that the enzyme SET7/9 plays an essential catalytic role in significantly lowering the barrier for the methyl-transfer reaction step. For the reaction in solution, it is found that the hydrogen bond network near the reaction center undergoes a significant change, and there is a strong shift in electrostatic field from the prereaction state to the transition state, whereas for the enzyme reaction, such an effect is much smaller and the enzyme SET7/9 is found to provide a preorganized electrostatic environment to facilitate the methyl-transfer reaction. Meanwhile, we find that the transition state in the enzyme reaction is a little more dissociative than that in solution.

SUBMITTER: Wang S

PROVIDER: S-EPMC2527688 | biostudies-literature | 2007 Apr

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Ab initio quantum mechanical/molecular mechanical molecular dynamics simulation of enzyme catalysis: the case of histone lysine methyltransferase SET7/9.

Wang Shenglong S Hu Po P Zhang Yingkai Y

The journal of physical chemistry. B 20070322 14

To elucidate enzyme catalysis through computer simulation, a prerequisite is to reliably compute free energy barriers for both enzyme and solution reactions. By employing on-the-fly Born-Oppenheimer molecular dynamics simulations with the ab initio quantum mechanical/molecular mechanical approach and the umbrella sampling method, we have determined free energy profiles for the methyl-transfer reaction catalyzed by the histone lysine methyltransferase SET7/9 and its corresponding uncatalyzed reac ...[more]

PMID: 17388541

Similar Datasets

Project description:The present study successfully implemented the ab initio quantum mechanical charge field molecular dynamics (QMCF MD) formalism for the investigation of structural and dynamical properties of hydrated cobalt-porphyrin complexes. Considering the significance of cobalt ions in biological systems (for instance, vitamin B12), which reportedly incorporate cobalt ions in a d6, low spin, +3 state chelated in the corrin ring, an analog of porphyrin, the current study is focused on cobalt in the oxidation states +2 and +3 bound to the parent porphyrin lead structures embedded in an aqueous solution. These cobalt-porphyrin complexes were investigated in terms of their structural and dynamical properties at the quantum chemical level. The structural attributes of these hydrated complexes revealed the contrasting features of the water binding to these solutes, including a detailed evaluation of the associated dynamics. The study also yielded notable findings in regard to the respective electronic configurations vs coordination, which suggested that Co(II)-POR possesses a 5-fold square pyramidal coordination geometry in an aqueous solution containing the metal ion coordinating to four nitrogen atoms of the porphyrin ring and one axial water as the fifth ligand. On the other hand, high-spin Co(III)-POR was hypothesized to be more stable due to the smaller size-to-charge ratio of the cobalt ion, but the high-spin complex demonstrated unstable structural and dynamical behavior. However, the corresponding properties of the hydrated Co(III)LS-POR revealed a stable structure in an aqueous solution, thus suggesting the Co(III) ion to be in a low-spin state when bound to the porphyrin ring. Moreover, the structural and dynamical data were augmented by computing the free energy of water binding to the cobalt ions and the solvent-accessible surface area, which provide further information on thermochemical properties of the metal-water interaction and the hydrogen bonding potential of the porphyrin ring in these hydrated systems.

Project description:DNA polymerases play a crucial role in the cell cycle due to their involvement in genome replication and repair. Understanding the reaction mechanism by which these polymerases carry out their function can provide insights into these processes. Recently, the crystal structures of human DNA polymerase lambda (Pollambda) have been reported both for pre- and post-catalytic complexes [García-Díaz et al., DNA Repair 3 (2007), 1333]. Here we employ the pre-catalytic complex as a starting structure for the determination of the catalytic mechanism of Pollambda using ab initio quantum mechanical/molecular mechanical methods. The reaction path has been calculated using Mg(2+) and Mn(2+) as the catalytic metals. In both cases the reaction proceeds through a two-step mechanism where the 3'-OH of the primer sugar ring is deprotonated by one of the conserved Asp residues (D490) in the active site before the incorporation of the nucleotide to the nascent DNA chain. A significant charge transfer is observed between both metals and some residues in the active site as the reaction proceeds. The optimized reactant and product structures agree with the reported crystal structures. In addition, the calculated reaction barriers for both metals are close to experimentally estimated barriers. Energy decomposition analysis to explain individual residue contributions suggests that several amino acids surrounding the active site are important for catalysis. Some of these residues, including R420, R488 and E529, have been implicated in catalysis by previous mutagenesis experiments on the homologous residues on Polbeta. Furthermore, Pollambda residues R420 and E529 found to be important from the energy decomposition analysis, are homologous to residues R183 and E295 in Polbeta, both of which are linked to cancer. In addition, residues R386, E391, K422 and K472 appear to have an important role in catalysis and could be a potential target for mutagenesis experiments. There is partial conservation of these residues across the Pol X family of DNA polymerases.

Dataset Information

Ab initio quantum mechanical/molecular mechanical molecular dynamics simulation of enzyme catalysis: the case of histone lysine methyltransferase SET7/9.

Publications

Ab initio quantum mechanical/molecular mechanical molecular dynamics simulation of enzyme catalysis: the case of histone lysine methyltransferase SET7/9.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets