Transition States in a Protein Environment − ONIOM QM:MM Modeling of Isopenicillin N Synthesis

Comparative reactivity of ferric-superoxo and ferryl-oxo species in heme and non-heme complexes

Ferryl-oxo species have been recognized as a key oxidant in many heme and non-heme enzymes. Recently, less-characterized ferric-superoxo species have been found or suggested to be another electrophilic oxidant. Reactivity of several vital ferryl-oxo and ferric-superoxo model complexes was examined by DFT calculations. Reactivity is found to correlate well with thermodynamic driving force and can increase with higher electrophilicity of the oxidant. Reactivity of the ferric-superoxo oxidants generally is not "superior" to the ferryl-oxo ones. Compared to the high-spin non-heme ferric-superoxo, the lower reactivity of low-spin heme ferric-superoxo, seldom utilized in nature, can be attributed to lower electrophilicity and more pronounced quenching of anti-ferromagnetic coupling between the ferric and superoxo parts. The present comparison should shed some light on mechanistic strategies in heme and non-heme enzymes and provide clues to rational design of ferric-superoxo oxidants.

[Recent progress in the theoretical studies of structure, function, and reaction of biological molecules]

Essential biomolecular functions often involve electron-related events such as chemical reactions and photoluminescence phenomena. Theoretical description of such electronic processes requires the use of quantum mechanics (QM), but the number of atoms that can be handled with QM is usually smaller than the number of atoms present in a single protein. A reasonable strategy is therefore to give priority to a few tens or hundreds of atoms in the system and deal with them quantum mechanically. Lower-priority atoms influence the event occurring in the higher-priority area; therefore, their effect should also be taken into account. Under these circumstances, a reasonable approach is to apply two or more different theoretical methods to differently prioritized subsystems. QM can be combined, for example, with less accurate yet much less demanding molecular mechanics (MM). Our own N-layered integrated molecular orbital and molecular mechanics (ONIOM) method allows for such hybrid calculations, and our group has been applying it to a wide range of biology-related problems. In this paper, we briefly explain the theoretical background and the procedure for the theoretical investigation of biological systems. Subsequently, we provide an overview of some of our recent studies of metalloenzymes and photobiology-related problems.

A QM/QM investigation of the hUNG2 reaction surface: the untold tale of a catalytic residue

Human uracil-DNA glycosylase (hUNG2) is a base excision repair enzyme that removes the damaged base uracil from DNA through hydrolytic deglycosylation of the nucleotide. In the present study, the mechanism of hUNG2 is thoroughly investigated using ONIOM(MPWB1K/6-31G(d):PM3) active-site models to generate reaction potential energy surfaces. Active-site models that differ in the hydrogen-bonding arrangement of the nucleophilic water molecule and/or protonation state of His148 are considered. The large barrier calculated using the model with a cationic His148 verifies that this residue is neutral in the early stages of the reaction. The reaction pathways predicted by two models with a neutral His148 are consistent with a wealth of experimental data on the enzyme, including mutational studies, which supports our approach. On the basis of our calculations, we propose a complete mechanism for the chemical step of hUNG2. In the first part of the reaction, His268, Asn204, and a water molecule work together to stabilize the negative charge forming on the uracil moiety. Subsequently, either Asp145 or His148 can act as the general base that activates the water nucleophile depending on the binding orientation of the water molecule in the active site. However, we propose that His148 preferentially acts as the general base. Therefore, in agreement with previous proposals, we assign the primary function of Asp145 to electrostatic stabilization of the positive charge developing on the sugar moiety during the reaction, which is also consistent with a growing theory that the primary function of active-site carboxylate groups present in many glycosylases is transition state stabilization. Most importantly, our work explains, for the first time, the role of His148 in the chemical step and provides additional support for the inclusion of this amino acid in the list of residues (Asp145 and His268) essential to the chemical step of the hUNG2 mechanism.

Double-lock ratchet mechanism revealing the role of alphaSER-344 in FoF1 ATP synthase

In a majority of living organisms, FoF1 ATP synthase performs the fundamental process of ATP synthesis. Despite the simple net reaction formula, ADP+Pi→ATP+H2O, the detailed step-by-step mechanism of the reaction yet remains to be resolved owing to the complexity of this multisubunit enzyme. Based on quantum mechanical computations using recent high resolution X-ray structures, we propose that during ATP synthesis the enzyme first prepares the inorganic phosphate for the γP-OADP bond-forming step via a double-proton transfer. At this step, the highly conserved αS344 side chain plays a catalytic role. The reaction thereafter progresses through another transition state (TS) having a planar ion configuration to finally form ATP. These two TSs are concluded crucial for ATP synthesis. Using stepwise scans and several models of the nucleotide-bound active site, some of the most important conformational changes were traced toward direction of synthesis. Interestingly, as the active site geometry progresses toward the ATP-favoring tight binding site, at both of these TSs, a dramatic increase in barrier heights is observed for the reverse direction, i.e., hydrolysis of ATP. This change could indicate a "ratchet" mechanism for the enzyme to ensure efficacy of ATP synthesis by shifting residue conformation and thus locking access to the crucial TSs.

QM/MM Studies of the Matrix Metalloproteinase 2 (MMP2) Inhibition Mechanism of (S)-SB-3CT and its Oxirane Analogue

SB-3CT, (4-phenoxyphenylsulfonyl)methylthiirane, is a potent, mechanism-based inhibitor of the gelatinase sub-class of the matrix metalloproteinase (MMP) family of zinc proteases. The gelatinase MMPs are unusual in that there are several examples where both enantiomers of a racemic inhibitor have comparable inhibitory abilities. SB-3CT is one such example. Here, the inhibition mechanism of the MMP2 gelatinase by the (S)-SB-3CT enantiomer and its oxirane analogue is examined computationally, and compared to the mechanism of (R)-SB-3CT. Inhibition of MMP2 by (R)-SB-3CT was shown previously to involve enzyme-catalyzed C-H deprotonation adjacent to the sulfone, with concomitant opening by β-elimination of the sulfur of the three-membered thiirane ring. Similarly to the R enantiomer, (S)-SB-3CT was docked into the active site of MMP2, followed by molecular dynamics simulation to prepare the complex for combined quantum mechanics and molecular mechanics (QM/MM) calculations. QM/MM calculations with B3LYP/6-311+G(d,p) for the QM part (46 atoms) and the AMBER force field for the MM part were used to compare the reaction of (S)-SB-3CT and its oxirane analogue in the active site of MMP2 (9208 atoms). These calculations show that the barrier for the proton abstraction coupled ring opening reaction of (S)-SB-3CT in the MMP2 active site is 4.4 kcal/mol lower than its oxirane analogue, and the ring opening reaction energy of (S)-SB-3CT is only 1.6 kcal/mol less exothermic than its oxirane analogue. Calculations also show that the protonation of the ring-opened products by water is thermodynamically much more favorable for the alkoxide obtained from the oxirane, than for the thiolate obtained from the thiirane. In contrast to (R)-SB-3CT and the R-oxirane analogue, the double bonds of the ring-opened products of (S)-SB-3CT and its S-oxirane analogue have the cis-configuration. Vibrational frequency and intrinsic reaction path calculations on a reduced size QM/MM model (2747 atoms) provide additional insight into the mechanism. These calculations yield 5.9 and 6.7 for the deuterium kinetic isotope effect for C-H bond cleavage in the transition state for the R and S enantiomers of SB-3CT, in good agreement with the experimental results.

The assessment of density functionals for DNA–protein stacked and T-shaped complexes

The present work uses 129 nucleobase – amino acid CCSD(T)/CBS stacking and T-shaped interaction energies as reference data to test the ability of various density functionals with double-zeta quality basis sets, as well as some semi-empirical and molecular mechanics methods, to accurately describe noncovalent DNA–protein π–π and π+–π interactions. The goal of this work is to identify methods that can be used in hybrid approaches (QM/MM, ONIOM) for large-scale modeling of enzymatic systems involving active-site (substrate) π–π contacts. Our results indicate that AMBER is a more appropriate choice for the lower-level method in hybrid techniques than popular semi-empirical methods (AM1, PM3), and suggest that AMBER accurately describes the π–π interactions found throughout DNA–protein complexes. The M06–2X and PBE-D density functionals were found to provide very promising descriptions of the 129 nucleobase – amino acid interaction energies, which suggests that these may be the most suitable methods for describing high-level regions. Therefore, M06–2X and PBE-D with both the 6–31G(d) and 6–31+G(d,p) basis sets were further examined through potential-energy surface scans to better understand how these techniques describe DNA–protein π–π interactions in both minimum and nonminimum regions of the potential-energy surfaces, which is critical information when modeling enzymatic reaction pathways. Our results suggest that studies of stacked nucleobase – amino acid systems should implement the PBE-D/6–31+G(d,p) method. However, if T-shaped contacts are involved and (or) smaller basis sets must be considered due to limitations in computational resources, then M06–2X/6–31G(d) provides an overall excellent description of both nucleobase – amino acid stacking and T-shaped interactions for a range of DNA–protein π–π and π+–π interactions.

Reaction coordinate of isopenicillin N synthase: oxidase versus oxygenase activity

Isopenicillin N synthase (IPNS) can have both oxidase and oxygenase activity depending on the substrate. For the native substrate, ACV, oxidase activity exists; however, for the substrate analogue ACOV, which lacks an amide nitrogen, IPNS exhibits oxygenase activity. The potential energy surfaces for the O-O bond elongation and cleavage were calculated for three different reactions: homolytic cleavage via traditional Fenton chemistry, heterolytic cleavage, and nucleophilic attack. These surfaces show that the hydroperoxide-ferrous intermediate, formed by O(2)-activated H atom abstraction from the substrate, can exploit different reaction pathways and that interactions with the substrate govern the pathway. The hydrogen bonds from hydroperoxide to the amide nitrogen of ACV polarize the sigma* orbital of the peroxide toward the proximal oxygen, facilitating heterolytic cleavage. For the substrate analogue ACOV, this hydrogen bond is no longer present, leading to nucleophilic attack on the substrate intermediate C-S bond. After cleavage of the hydroperoxide, the two reaction pathways proceed with minimal barriers, resulting in the closure of the beta-lactam ring for the oxidase activity (ACV) or formation of the thiocarboxylate for oxygenase activity (ACOV).

Insights into the (superoxo)Fe(III)Fe(III) intermediate and reaction mechanism of myo-inositol oxygenase: DFT and ONIOM(DFT:MM) study

The (superoxo)Fe(III)Fe(III) reactive species and the catalytic reaction mechanism of a diiron enzyme, myo-inositol oxygenase (MIOX), were theoretically investigated by means of density functional theory (DFT) and ONIOM quantum mechanical/molecular mechanical (QM/MM) approaches. The ground state of the (superoxo)Fe(III)Fe(III) intermediate was shown to have a side-on coordination geometry and an S = 1/2 spin state, wherein the two iron sites are antiferromagnetically coupled while the superoxide site and the nearest iron are ferromagnetically coupled. A full reaction pathway leading to a D-glucuronate product from myo-inositol was proposed based on ONIOM computational results. Two major roles of the enzyme surrounding during the catalytic reaction were identified. One is to facilitate the initial H-abstraction step, and the other is to restrict the movement of the substrate via H-bonding interactions in order to avoid unwanted side reactions. In our proposed mechanism, O-O bond cleavage has the highest barrier, thus constituting the rate-limiting step of the reaction. The unique role of the bridging hydroxide ligand as a catalytic base was also identified.

Matrix metalloproteinase 2 inhibition: combined quantum mechanics and molecular mechanics studies of the inhibition mechanism of (4-phenoxyphenylsulfonyl)methylthiirane and its oxirane analogue

The inhibition mechanism of matrix metalloproteinase 2 (MMP2) by the selective inhibitor (4-phenoxyphenylsulfonyl)methylthiirane (SB-3CT) and its oxirane analogue is investigated computationally. The inhibition mechanism involves C-H deprotonation with concomitant opening of the three-membered heterocycle. SB-3CT was docked into the active site of MMP2, followed by molecular dynamics simulation to prepare the complex for combined quantum mechanics and molecular mechanics (QM/MM) calculations. QM/MM calculations with B3LYP/6-311+G(d,p) for the QM part and the AMBER force field for the MM part were used to examine the reaction of these two inhibitors in the active site of MMP2. The calculations show that the reaction barrier for transformation of SB-3CT is 1.6 kcal/mol lower than its oxirane analogue, and the ring-opening reaction energy of SB-3CT is 8.0 kcal/mol more exothermic than that of its oxirane analogue. Calculations also show that protonation of the ring-opened product by water is thermodynamically much more favorable for the alkoxide obtained from the oxirane than for the thiolate obtained from the thiirane. A six-step partial charge fitting procedure is introduced for the QM/MM calculations to update atomic partial charges of the quantum mechanics region and to ensure consistent electrostatic energies for reactants, transition states, and products.

Quantum mechanics/molecular mechanics investigation of the mechanism of phosphate transfer in human uridine-cytidine kinase 2

The mechanisms of enzyme-catalyzed phosphate transfer and hydrolysis reactions are of great interest due to their importance and abundance in biochemistry. The reaction may proceed in a stepwise fashion, with either a pentavalent phosphorane or a metaphosphate anion intermediate, or by a concerted SN2 mechanism. Despite much theoretical work focused on a few key enzymes, a consensus for the mechanism has not been reached, and examples of all three possibilities have been demonstrated. We have investigated the mechanism of human uridine-cytidine kinase 2 (UCK2, EC 2.7.1.48), which catalyzes the transfer of a phosphate group from ATP to the ribose 5'-hydroxyl of cytidine and uridine. UCK2 is normally expressed in human placenta, but is overexpressed in certain cancer cells, where it is responsible for activating a class of antitumor prodrugs. The UCK2 mechanism was investigated by generating a 2D potential energy surface as a function of the P-O bonds forming and breaking, with energies calculated using a quantum mechanics/molecular mechanics potential (B3LYP/6-31G(d):AMBER). The mechanism of phosphate transfer is shown to be concerted, and is accompanied by concerted proton transfer from the 5'-hydroxyl to a conserved active site aspartic acid that serves as a catalytic base. The calculated barrier for this reaction is 15.1 kcal/mol, in relatively good agreement with the experimental barrier of 17.5 kcal/mol. The interactions of the enzyme active site with the reactant, transition state, and product are examined for their implications on the design of anticancer prodrugs or positron emission tomography (PET) reporter probes for this enzyme.