Density functional theory study of the reaction mechanism of the DNA repairing enzyme alkylguanine alkyltransferase

Eliminating Imaginary Vibrational Frequencies in Quantum-Chemical Cluster Models of Enzymatic Active Sites

In constructing finite models of enzyme active sites for quantum-chemical calculations, atoms at the periphery of the model must be constrained to prevent unphysical rearrangements during geometry relaxation. A simple fixed-atom or "coordinate-lock" approach is commonly employed but leads to undesirable artifacts in the form of small imaginary frequencies. These preclude evaluation of finite-temperature free-energy corrections, limiting thermochemical calculations to enthalpies only. Full-dimensional vibrational frequency calculations are possible by replacing the fixed-atom constraints with harmonic confining potentials. Here, we compare that approach to an alternative strategy in which fixed-atom contributions to the Hessian are simply omitted. While the latter strategy does eliminate imaginary frequencies, it tends to underestimate both the zero-point energy and the vibrational entropy while introducing artificial rigidity. Harmonic confining potentials eliminate imaginary frequencies and provide a flexible means to construct active-site models that can be used in unconstrained geometry relaxations, affording better convergence of reaction energies and barrier heights with respect to the model size, as compared to models with fixed-atom constraints.

Mechanism of AGT-Mediated Repair of dG-dC Cross-Links in the Drug Resistance to Chloroethylnitrosoureas: Molecular Docking, MD Simulation, and ONIOM (QM/MM) Investigation

Chloroethylnitrosoureas (CENUs) are important chemotherapies applied in the treatment of cancer. They exert anticancer activity by inducing DNA interstrand cross-links (ICLs) via the formation of two O6-alkylguanine intermediates, O6-chloroethylguanine (O6-ClEtG) and N1,O6-ethanoguanine (N1,O6-EtG). However, O6-alkylguanine-DNA alkyltransferase (AGT), a DNA-repair enzyme, can restore the O6-alkylguanine damages and thereby obstruct the formation of ICLs (dG-dC cross-link). In this study, the inhibitory mechanism of ICL formation was investigated to elucidate the drug resistance of CENUs mediated by AGT in detail. Based on the structures of the substrate-enzyme complexes obtained from docking and MD simulations, two ONIOM (QM/MM) models with different sizes of the QM region were constructed. The model with a larger QM region, which included the substrate (O6-ClEtG or N1,O6-EtG), a water molecule, and five residues (Tyr114, Cys145, His146, Lys165, and Glu172) in the active pocket of AGT, accurately described the repairing reaction and generated the results coinciding with the experimental outcomes. The repair process consists of two sequential steps: hydrogen transfer to form a thiolate anion on Cys145 and alkyl transfer from the O6 site of guanine (the rate-limiting step). The repair of N1,O6-EtG was more favorable than that of O6-ClEtG from both kinetics and thermodynamics aspects. Moreover, the comparison of the repairing process with the formation of dG-dC cross-link and the inhibition of AGT by O6-benzylguanine (O6-BG) showed that the presence of AGT could effectively interrupt the formation of ICLs leading to drug resistance, and the inhibition of AGT by O6-BG that was energetically more favorable than the repair of O6-ClEtG could not prevent the repair of N1,O6-EtG. Therefore, it is necessary to completely eliminate AGT activity before CENUs medication to enhance the chemotherapeutic effectiveness. This work provides reasonable explanations for the supposed mechanism of AGT-mediated drug resistance of CENUs and will assist in the development of novel CENU chemotherapies and their medication strategies.

Catalytic mechanism and pH dependence of a methyltransferase ribozyme (MTR1) from computational enzymology

A methyltransferase ribozyme (MTR1) was selected in vitro to catalyze alkyl transfer from exogenous O6-methylguanine (O6mG) to a target adenine N1, and recently, high-resolution crystal structures have become available. We use a combination of classical molecular dynamics, ab initio quantum mechanical/molecular mechanical (QM/MM) and alchemical free energy (AFE) simulations to elucidate the atomic-level solution mechanism of MTR1. Simulations identify an active reactant state involving protonation of C10 that hydrogen bonds with O6mG:N1. The deduced mechanism involves a stepwise mechanism with two transition states corresponding to proton transfer from C10:N3 to O6mG:N1 and rate-controlling methyl transfer (19.4  kcal·mol-1 barrier). AFE simulations predict the pKa for C10 to be 6.3, close to the experimental apparent pKa of 6.2, further implicating it as a critical general acid. The intrinsic rate derived from QM/MM simulations, together with pKa calculations, enables us to predict an activity-pH profile that agrees well with experiment. The insights gained provide further support for a putative RNA world and establish new design principles for RNA-based biochemical tools.

Using Atomic Confining Potentials for Geometry Optimization and Vibrational Frequency Calculations in Quantum-Chemical Models of Enzyme Active Sites

Quantum-chemical studies of enzymatic reaction mechanisms sometimes use truncated active-site models as simplified alternatives to mixed quantum mechanics molecular mechanics (QM/MM) procedures. Eliminating the MM degrees of freedom reduces the complexity of the sampling problem, but the trade-off is the need to introduce geometric constraints in order to prevent structural collapse of the model system during geometry optimizations that do not contain a full protein backbone. These constraints may impair the efficiency of the optimization, and care must be taken to avoid artifacts such as imaginary vibrational frequencies. We introduce a simple alternative in which terminal atoms of the model system are placed in soft harmonic confining potentials rather than being rigidly constrained. This modification is simple to implement and straightforward to use in vibrational frequency calculations, unlike iterative constraint-satisfaction algorithms, and allows the optimization to proceed without constraint even though the practical result is to fix the anchor atoms in space. The new approach is more efficient for optimizing minima and transition states, as compared to the use of fixed-atom constraints, and also more robust against unwanted imaginary frequencies. We illustrate the method by application to several enzymatic reaction pathways where entropy makes a significant contribution to the relevant reaction barriers. The use of confining potentials correctly describes reaction paths and facilitates calculation of both vibrational zero-point and finite-temperature entropic corrections to barrier heights.

Proton-induced accelerated decay of the fungicide, vinclozolin, on TiO2 surface under solar irradiation: Experimental and DFT study

The photochemical degradation of vinclozolin by addition of titanium dioxide on silica support has been examined both experimental and quantum-chemically. Solar irradiation of vinclozolin on silica with and without addition of titanium dioxide for 6 h resulted in 21% and 97.8% vinclozolin residues, respectively. In both these cases, phototransformation leads to the formation of (3,5-dichlorophenyl isocyanate) and (3,5-dichloroaniline). The presence of the intermediary product resulted from opening of the 2,4-oxazolidine-dione ring is also confirmed by GS-MS and LC-MS chromatography. The proton-induced mechanism of vinclozolin decay at the above experimental conditions is clarified on the base of DFT calculations.

Influence of the Expression Level of O6-Alkylguanine-DNA Alkyltransferase on the Formation of DNA Interstrand Crosslinks Induced by Chloroethylnitrosoureas in Cells: A Quantitation Using High-Performance Liquid Chromatography-Mass Spectrometry

Chloroethylnitrosoureas (CENUs), which are bifunctional alkylating agents widely used in the clinical treatment of cancer, exert anticancer activity by inducing crosslink within a guanine-cytosine DNA base pair. However, the formation of dG-dC crosslinks can be prevented by O⁶-alkylguanine-DNA alkyltransferase (AGT), ultimately leading to drug resistance. Therefore, the level of AGT expression is related to the formation of dG-dC crosslinks and the sensitivity of cells to CENUs. In this work, we determined the CENU-induced dG-dC crosslink in mouse L1210 leukemia cells and in human glioblastoma cells (SF-763, SF-767 and SF-126) containing different levels of AGT using high-performance liquid chromatography coupled with electrospray ionization tandem mass spectrometry. The results indicate that nimustine (ACNU) induced more dG-dC crosslinks in L1210 leukemia cells than those induced by carmustine (BCNU), lomustine (CCNU) and fotemustine (FTMS). This result was consistent with a previously reported cohort study, which demonstrated that ACNU had a better survival gain than BCNU, CCNU and FTMS for patients with high-grade glioma. Moreover, we compared the crosslinking levels and the cytotoxicity in SF-763, SF-767 and SF-126 cells with different AGT expression levels after exposure to ACNU. The levels of dG-dC crosslink in SF-126 cells (low AGT expression) were significantly higher than those in SF-767 (medium AGT expression) and SF-763 (high AGT expression) cells at each time point. Correspondingly, the cytotoxicity of SF-126 was the highest followed by SF-767 and SF-763. The results obtained in this work provided unequivocal evidence for drug resistance to CENUs induced by AGT-mediated repair of DNA ICLs. We postulate that the level of dG-dC crosslink has the potential to be employed as a biomarker for estimating drug resistance and anticancer efficiencies of novel CENU chemotherapies.

Multifaceted roles of alkyltransferase and related proteins in DNA repair, DNA damage, resistance to chemotherapy, and research tools

O(6)-Alkylguanine-DNA alkyltransferase (AGT) is a widely distributed, unique DNA repair protein that acts as a single agent to directly remove alkyl groups located on the O(6)-position of guanine from DNA restoring the DNA in one step. The protein acts only once, and its alkylated form is degraded rapidly. It is a major factor in counteracting the mutagenic, carcinogenic, and cytotoxic effects of agents that form such adducts including N-nitroso-compounds and a number of cancer chemotherapeutics. This review describes the structure, function, and mechanism of action of AGTs and of a family of related alkyltransferase-like proteins, which do not act alone to repair O(6)-alkylguanines in DNA but link repair to other pathways. The paradoxical ability of AGTs to stimulate the DNA-damaging ability of dihaloalkanes and other bis-electrophiles via the formation of AGT-DNA cross-links is also described. Other important properties of AGTs include the ability to provide resistance to cancer therapeutic alkylating agents, and the availability of AGT inhibitors such as O(6)-benzylguanine that might overcome this resistance is discussed. Finally, the properties of fusion proteins in which AGT sequences are linked to other proteins are outlined. Such proteins occur naturally, and synthetic variants engineered to react specifically with derivatives of O(6)-benzylguanine are the basis of a valuable research technique for tagging proteins with specific reagents.

Quantum theoretical study of cleavage of the glycosidic bond of 2'-deoxyadenosine: base excision-repair mechanism of DNA by MutY

The enzyme adenine DNA glycosylase, also called MutY, is known to catalyze base excision repair by removal of adenine from the abnormal 2'-deoxyadenosine:8-oxo-2'-deoxyguanosine pair in DNA. The active site of the enzyme was considered to consist of a glutamic acid residue along with two water molecules. The relevant reaction mechanism involving different barrier energies was studied theoretically. Molecular geometries of the various molecules and complexes involved in the reaction, e.g., the reactant, intermediate, and product complexes as well as transition states, were optimized employing density functional theory at the B3LYP/6-31G(d,p) level in the gas phase. It was followed by single-point energy calculations at the B3LYP/AUG-cc-pVDZ, BHandHLYP/AUG-cc-pVDZ, and MP2/AUG-cc-pVDZ levels in the gas phase. Single-point energy calculations were also carried out at the B3LYP/AUG-cc-pVDZ and BHandHLYP/AUG-cc-pVDZ levels in aqueous media as well as in the solvents chlorobenzene and dichloroethane. For the solvation calculations, the integral equation formalism of the polarizable continuum model (IEF-PCM) was employed. It is found that glutamic acid along with two water molecules would effectively cleave the glycosidic bond of adenosine by a new two-step reaction mechanism proposed here which is different from the three-step mechanism proposed by other authors earlier regarding the working mechanism of MutY.

Mechanism of mycolic acid cyclopropane synthase: a theoretical study

The reaction mechanism of mycolic acid cyclopropane synthase is investigated using hybrid density functional theory. The direct methylation mechanism is examined with a large model of the active site constructed on the basis of the crystal structure of the native enzyme. The important active site residue Glu140 is modeled in both ionized and neutral forms. We demonstrate that the reaction starts via the transfer of a methyl to the substrate double bond, followed by the transfer of a proton from the methyl cation to the bicarbonate present in the active site. The first step is calculated to be rate-limiting, in agreement with experimental kinetic results. The protonation state of Glu140 has a rather weak influence on the reaction energetics. In addition to the natural reaction, a possible side reaction, namely a carbocation rearrangement, is also considered and is shown to have a low barrier. Finally, the energetics for the sulfur ylide proposal, which has already been ruled out, is also estimated, showing a large energetic penalty for ylide formation.

QM/MM study on the reaction mechanism of O6-alkylguanine-DNA alkyltransferase

Combined quantum-mechanical/molecular-mechanical (QM/MM) approaches have been applied to investigate the detailed reaction mechanism of human O(6)-alkylguanine-DNA alkyltransferase (AGT). AGT is a direct DNA repair protein that is capable of repairing alkylated DNA by transferring the methyl group to the thiol group of a cysteine residue (Cys145) in the active site in an irreversible and stoichiometric reaction. Our QM/MM calculations reveal that the methyl group transferring step is expected to occur through two steps, in which the methyl carbocation generating step is the rate-determining step with an energy barrier of 14.4 kcal/mol at the QM/MM B3LYP/6-31G(d,p)//CHARMM22 level of theory. It is different from the previous theoretical studies based on QM calculations by using a cluster model in which the methyl group transferring step is a one-step process with a higher energy barrier.

Quantum chemical modeling of enzymatic reactions: the case of histone lysine methyltransferase

Quantum chemical cluster models of enzyme active sites are today an important and powerful tool in the study of various aspects of enzymatic reactivity. This methodology has been applied to a wide spectrum of reactions and many important mechanistic problems have been solved. Herein, we report a systematic study of the reaction mechanism of the histone lysine methyltransferase (HKMT) SET7/9 enzyme, which catalyzes the methylation of the N-terminal histone tail of the chromatin structure. In this study, HKMT SET7/9 serves as a representative case to examine the modeling approach for the important class of methyl transfer enzymes. Active site models of different sizes are used to evaluate the methodology. In particular, the dependence of the calculated energies on the model size, the influence of the dielectric medium, and the particular choice of the dielectric constant are discussed. In addition, we examine the validity of some technical aspects, such as geometry optimization in solvent or with a large basis set, and the use of different density functional methods.

Repair of O6-methylguanine to guanine by cysteine in the absence and presence of histidine and by cysteine thiolate anion: a quantum chemical study

O6-methylguanine (O6mG) is known to be a potential mutagenic modification of guanine as it mispairs with thymine in DNA and causes GC to AT transversion mutation. It is experimentally known that O6mG can be repaired to guanine by the protein O6-alkylguanine-DNA alkyltransferase (AGT), a cysteine residue being the main active site. In the present work, the mechanisms of repair of cis-O6-methylguanine (O6mG) to guanine due to its reaction with cysteine in the absence and presence of histidine and with cysteine thiolate anion were investigated theoretically using the B3LYP hybrid functional of density functional theory and the second order Møller-Plesset perturbation (MP2) theory. Reactant, intermediate and product complexes as well as transition states involved in these reactions were fully optimized at the B3LYP/6-31 + G* level of theory in the gas phase. The solvent effect of water was treated using the polarizable continuum model (PCM). Single point energy calculations were performed at the B3LYP/AUG-cc-pVDZ and MP2/6-31 + G* levels of theory in the gas phase and aqueous media. It is found that cysteine alone can repair the cis-O6mG to guanine, but the involvement of histidine along with cysteine lowers down the barrier energy significantly. However, when cysteine thiolate anion is used in place of cysteine, the barrier energy is strongly reduced. These results broadly support the suggestions based on experimental studies.

A quantum chemical study of repair of O6-methylguanine to guanine by tyrosine: evaluation of the winged helix-turn-helix model

The winged helix-turn-helix model for the repair of O6-MeG to guanine involving the reaction of O6-MeG with a tyrosine residue of the protein O6-alkylguanine-DNA alkyltransferase (AGT) was examined by studying the reaction mechanism and barrier energies. Molecular geometries of the species and complexes involved in the reaction, i.e. the reactant, intermediate and product complexes as well as transition states, were optimized employing density functional theory in gas phase. It was followed by single point energy calculations using density functional theory along with a higher basis set and second order M(phi)ller-Plesset perturbation theory (MP2) along with two different basis sets in gas phase and aqueous media. For the solvation calculations in aqueous media, the integral equation formalism of the polarizable continuum model (IEF-PCM) was employed. Vibrational frequency analysis was performed for each optimized structure and genuineness of transition states was ensured by visualizing the vibrational modes. It is found that tyrosine can repair O6-MeG to guanine by a two-step reaction. The present results have been compared with those obtained considering the helix-turn-helix model where the repair reaction primarily involves cysteine and occurs in a single-step. It is concluded that the repair through tyrosine envisaged in the winged helix-turn-helix model would be less efficient than that through cysteine envisaged in the helix-turn-helix model.

Progress in ab initio QM/MM free-energy simulations of electrostatic energies in proteins: accelerated QM/MM studies of pKa, redox reactions and solvation free energies

Hybrid quantum mechanical/molecular mechanical (QM/MM) approaches have been used to provide a general scheme for chemical reactions in proteins. However, such approaches still present a major challenge to computational chemists, not only because of the need for very large computer time in order to evaluate the QM energy but also because of the need for proper computational sampling. This review focuses on the sampling issue in QM/MM evaluations of electrostatic energies in proteins. We chose this example since electrostatic energies play a major role in controlling the function of proteins and are key to the structure-function correlation of biological molecules. Thus, the correct treatment of electrostatics is essential for the accurate simulation of biological systems. Although we will be presenting different types of QM/MM calculations of electrostatic energies (and related properties) here, our focus will be on pKa calculations. This reflects the fact that pKa's of ionizable groups in proteins provide one of the most direct benchmarks for the accuracy of electrostatic models of macromolecules. While pKa calculations by semimacroscopic models have given reasonable results in many cases, existing attempts to perform pKa calculations using QM/MM-FEP have led to discrepancies between calculated and experimental values. In this work, we accelerate our QM/MM calculations using an updated mean charge distribution and a classical reference potential. We examine both a surface residue (Asp3) of the bovine pancreatic trypsin inhibitor and a residue buried in a hydrophobic pocket (Lys102) of the T4-lysozyme mutant. We demonstrate that, by using this approach, we are able to reproduce the relevant side chain pKa's with an accuracy of 3 kcal/mol. This is well within the 7 kcal/mol energy difference observed in studies of enzymatic catalysis, and is thus sufficient accuracy to determine the main contributions to the catalytic energies of enzymes. We also provide an overall perspective of the potential of QM/MM calculations in general evaluations of electrostatic free energies, pointing out that our approach should provide a very powerful and accurate tool to predict the electrostatics of not only solution but also enzymatic reactions, as well as the solvation free energies of even larger systems, such as nucleic acid bases incorporated into DNA.