DNA polymerases are enzymes responsible for the synthesis of DNA from nucleotides. Understanding their molecular fundamentals is a prerequisite for elucidating their aberrant activities in diseases such as cancer. Here we have carried out ab initio quantum mechanical/molecular mechanical (QM/MM) studies on the nucleotidyl-transfer reaction catalyzed by the lesion-bypass DNA polymerase IV (Dpo4) from Sulfolobus solfataricus, with template guanine and Watson-Crick paired dCTP as the nascent base pair. The results suggested a novel water-mediated and substrate-assisted (WMSA) mechanism: the initial proton transfer to the alpha-phosphate of the substrate via a bridging crystal water molecule is the rate-limiting step, the nucleotidyl-transfer step is associative with a metastable pentacovalent phosphorane intermediate, and the pyrophosphate leaving is facilitated by a highly coordinated proton relay mechanism through mediation of water which neutralizes the evolving negative charge. The conserved carboxylates, which retain their liganding to the two Mg2+ ions during the reaction process, are found to be essential in stabilizing transition states. This WMSA mechanism takes specific advantage of the unique structural features of this low-fidelity lesion-bypass Y-family polymerase, which has a more spacious and solvent-exposed active site than replicative and repair polymerases.

Histone lysine methylation is emerging as an important mechanism to regulate chromatin structure and gene activity. To provide theoretical understanding of its reaction mechanism and product specificity, ab initio quantum mechanical/molecular mechanical free energy (QM/MM-FE) calculations and molecular dynamics simulations have been carried out to investigate the histone lysine methyltransferase SET7/9. It is found that the methyl-transfer reaction catalyzed by SET7/9 is a typical in-line S(N)2 nucleophilic substitution reaction with a transition state of 70% dissociative character. The calculated average free energy barrier at the MP2(6-31+G) QM/MM level is 20.4 +/- 1.1 kcal/mol, consistent with the activation barrier of 20.9 kcal/mol estimated from the experimental reaction rate. The barrier fluctuation has a strong correlation with the nucleophilic attack distance and angle in the reactant complex. The calculation results show that the product specificity of SET7/9 as a monomethyltransferase is achieved by disrupting the formation of near-attack conformations for the dimethylation reaction.

To characterize the catalytic mechanism for zinc-dependent histone deacetylases (HDAC), we have carried out density functional theory QM/MM studies on the deacetylation reaction catalyzed by a histone-deacetylase-like protein (HDLP). The calculation results do not support the previous mechanistic hypothesis, but suggest a lower protonation state for the active site as well as a 4-fold zinc coordination during the reaction process. To characterize such mechanistic difference is not only significant for our fundamental understanding of its inner workings but also crucial for the design of HDAC inhibitors.

Sir2 enzymes catalyze the NAD + dependent protein deacetylation and play critical roles in epigenetics, cell death and lifespan regulation. In spite of a current flurry of experimental studies, the catalytic mechanism for this unique and important class of enzymes remains elusive. Employing on-the-fly Born-Oppenheimer molecular dynamics simulations with the B3LYP/6-31G(d) QM/MM potential and the umbrella sampling method, we have characterized the initial step of the Sir2Tm catalyzed reaction, which is also the most controversial portion of its mechanism. Our results indicate that the nicotinamide cleavage reaction employs a highly dissociative and concerted displacement mechanism: the cleavage of the glycosidic bond is facilitated by the nucleophilic participation of the acetyl-lysine, and the dissociative transition state has a significant oxocarbenium ion character. During this step of the reaction, the Sir2Tm enzyme strongly stabilizes the covalent O-alkylamidate intermediate whereas its effect on the transition state is quite minimal. In addition, functional roles of key residues and motifs have been elucidated. This work further demonstrates the feasibility and applicability of the state-of-the-art ab initio QM/MM molecular dynamics approach in simulating enzyme reactions.

AbstracsThe different coordination modes and fast ligand exchange of zinc coordination has been suggested to be one key catalytic feature of the zinc ion which makes it an invaluable metal in biological catalysis. However, partly due to the well known difficulties for zinc to be characterized by spectroscopy methods, evidence for dynamic nature of the catalytic zinc coordination has so far mainly been indirect. In this work, Born-Oppenheimer ab initio QM/MM molecular dynamics simulation has been employed, which allows for a first-principle description of the dynamics of the metal active site while properly including effects of the heterogeneous and fluctuating protein environment. Our simulations have provided direct evidence regarding inherent flexibility of the catalytic zinc coordination shell in Thermolysin (TLN) and Histone Deacetylase 8 (HDAC8). We have observed different coordination modes and fast ligand exchange during the picosecond's timescale. For TLN, the coordination of the carboxylate group of Glu166 to Zinc is found to continuously change between monodentate and bidentate manner dynamically; while for HDAC8, the flexibility mainly comes from the coordination to a non-amino-acid ligand. Such distinct dynamics in the zinc coordination shell between two enzymes suggests that the catalytic role of Zinc in TLN and HDAC8 is likely to be different in spite of the fact that both catalyze the hydrolysis of amide bond. Meanwhile, considering that such Born-Oppenheimer ab initio QM/MM MD simulations are very much desired but are widely considered to be too computationally expensive to be feasible, our current study demonstrates the viability and powerfulness of this state-of-the-art approach in simulating metalloenzymes.Zinc is relatively abundant in biological materials. Approximately 10% of the total human proteome have been identified to bind with zinc in vivo from a bioinformatics investigation1 and they play very crucial roles in all forms of life2 -6 . For mononuclear zinc enzymes, a typical metal coordination environment contains three amino acid side chain ligands (His, Glu, Asp and Cys) and one/two small molecule(s). 3, 7 , 8 The flexibility of zinc coordination, which allows different coordination modes and fast ligand exchange, has been suggested to be one key catalytic feature of the zinc ion which makes it an invaluable metal in biological catalysis.9 However, partly due to the well known difficulties for zinc to be characterized by spectroscopy methods 10,11 , evidence for dynamic nature of the catalytic zinc coordination has so far mainly In order to provide deep insights into the dynamics and flexibility of the zinc catalytic site, which would be essential in characterizing their catalytic mechanisms and rational design of novel inhibitors for zinc enzymes, we have carried out DFT QM/MM Born-Oppenheimer molecular dynamics (BOMD) simulations on TLN and HDAC8. Although semi-empirical QM/MM BOMD simulations of some zinc-dependent enzymes have been carried out 29-32, one main concern is the acc...