Elucidating the mechanism for the reduction of nitrite by copper nitrite reductase--a contribution from quantum chemical studies

Density functional methods have been applied to investigate the properties of the active site of copper-containing nitrite reductases and possible reaction mechanisms for the enzyme catalysis. The results for a model of the active site indicate that a hydroxyl intermediate is not formed during the catalytic cycle, but rather a state with a protonated nitrite bound to the reduced copper. Electron affinity calculations indicate that reduction of the T2 copper site does not occur immediately after nitrite binding. Proton affinity calculations are indicative of substantial pK(a) differences between different states of the T2 site. The calculations further suggest that the reaction does not proceed until uptake of a second proton from the bulk solution. They also indicate that Asp-92 may play both a key role as a proton donor to the substrate, and a structural role in promoting catalysis. In the D92N mutant another base, presumably a nearby histidine (His-249) may take the role as the proton donor. On the basis of these model calculations and available experimental evidence, an ordered reaction mechanism for the reduction of nitrite is suggested. An investigation of the binding modes of the nitric oxide product and the nitrite substrate to the model site has also been made, indicating that nitric oxide prefers to bind in an end-on fashion to the reduced T2 site.     

Electrostatic study of the proton pumping mechanism in bovine heart cytochrome C oxidase

Cytochrome c oxidase (CcO) is the terminal enzyme of the cell respiratory chain in mitochondria and aerobic bacteria. It catalyzes the reduction of oxygen to water and utilizes the free energy of the reduction reaction for proton pumping across the inner-mitochondrial membrane, a process that results in a membrane electrochemical proton gradient. Although the structure of the enzyme has been solved for several organisms, the molecular mechanism of proton pumping remains unknown. In the present paper, continuum electrostatic calculations were employed to evaluate the electrostatic potential, energies, and protonation state of bovine heart cytochrome c oxidase for different redox states of the enzyme along its catalytic cycle. Three different computational models of the enzyme were employed to test the stability of the results. The energetics and pH dependence of the P-->F, F-->O, and O-->E steps of the cycle have been investigated. On the basis of electrostatic calculations, two possible schemes of redox-linked proton pumping are discussed. The first scheme involves His291 as a pump element, whereas the second scheme involves a group linked to propionate D of heme a(3). In both schemes, loading of the pump site is coupled to ET between the two hemes of the enzyme, while transfer of a chemical proton is accompanied by ejection of the pumped H(+). The two models, as well as the energetics results are compared with recent experimental kinetic data. The proton pumping across the membrane is an endergonic process, which requires a sufficient amount of energy to be provided by the chemical reaction in the active site. In our calculations, the conversion of OH(-) to H(2)O provides 520 meV of energy to displace pump protons from a loading site and overall about 635 meV for each electron passing through the system. Assuming that the two charges are translocated per electron against the membrane potential of 200 meV, the model predicts an overall efficiency of 63%.     

Quantum chemistry applied to the mechanisms of transition metal containing enzymes -- cytochrome c oxidase, a particularly challenging case

The Density functional theory (B3LYP) has been used to study the mechanisms of O--O bond cleavage and proton pumping in cytochrome c oxidase. To understand how the energy from the exergonic reduction of molecular oxygen is used to pump protons across the mitochondrial membrane, the energetics of all steps in the catalytic cycle have to be evaluated. For this purpose, models have to be designed that can accurately reproduce relative redox potentials and pKa values within the active site. The present study shows that it is possible to construct such models and to calculate energy profiles which, to a large extent, agree with experimental information. However, the energy profiles point out a problem with an unbalanced partitioning of the energy between the reductive and oxidative half cycles, which is in disagreement with the experimental observation that the proton pumping is evenly distributed between the two half cycles. A conclusion from the present study is, therefore, that something is probably still missing in the modeling of the active site.     

Water-assisted reaction mechanism of monozinc beta-lactamases

Hybrid Car-Parrinello QM/MM calculations are used to investigate the reaction mechanism of hydrolysis of a common beta-lactam substrate (cefotaxime) by the monozinc beta-lactamase from Bacillus cereus (BcII). The calculations suggest a fundamental role for an active site water in the catalytic mechanism. This water molecule binds the zinc ion in the first step of the reaction, expanding the zinc coordination number and providing a proton donor adequately oriented for the second step. The free energy barriers of the two reaction steps are similar and consistent with the available experimental data. The conserved hydrogen bond network in the active site, defined by Asp120, Cys221, and His263, not only contributes to orient the nucleophile (as already proposed), but it also guides the second catalytic water molecule to the zinc ion after the substrate is bound. The hydrolysis reaction in water has a relatively high free energy barrier, which is consistent with the stability of cefotaxime in water solution. The modeled Michaelis complexes for other substrates are also characterized by the presence of an ordered water molecule in the same position, suggesting that this mechanism might be general for the hydrolysis of different beta-lactam substrates.     

Theoretical perspectives on the reaction mechanism of serine proteases: the reaction free energy profiles of the acylation process

The reaction mechanism of serine proteases (trypsin), which catalyze peptide hydrolysis, is studied theoretically by ab initio QM/MM electronic structure calculations combined with Molecular Dynamics-Free Energy Perturbation calculations. We have calculated the entire reaction free energy profiles of the first reaction step of this enzyme (acylation process). The present calculations show that the rate-determining step of the acylation is the formation of the tetrahedral intermediate, and the breakdown of this intermediate has a small energy barrier. The calculated activation free energy for the acylation is approximately 17.8 kcal/mol at QM/MM MP2/(aug)-cc-pVDZ//HF/6-31(+)G/AMBER level, and this reaction is an exothermic process. MD simulations of the enzyme-substrate (ES) complex and the free enzyme in aqueous phase show that the substrate binding induces slight conformational changes around the active site, which favor the alignment of the reactive fragments (His57, Asp102, and Ser195) together in a reactive orientation. It is also shown that the proton transfer from Ser195 to His57 and the nucleophilic attack of Ser195 to the carbonyl carbon of the scissile bond of the substrate occur in a concerted manner. In this reaction, protein environment plays a crucial role to lowering the activation free energy by stabilizing the tetrahedral intermediate compared to the ES complex. The polarization energy calculations show that the enzyme active site is in a very polar environment because of the polar main chain contributions of protein. Also, the ground-state destabilization effect (steric strain) is not a major catalytic factor. The most important catalytic factor of stabilizing the tetrahedral intermediate is the electrostatic interaction between the active site and particular regions of protein: the main chain NH groups in Gly193 and Ser195 (so-called oxyanion hole region) stabilize negative charge generated on the carbonyl oxygen of the scissile bond, and the main chain carbonyl groups in Ile212 approximately Ser214 stabilize a positive charge generated on the imidazole ring of His57.     

Is the bound substrate in nitric oxide synthase protonated or neutral and what is the active oxidant that performs substrate hydroxylation?

We present here results of a series of density functional theory (DFT) studies on enzyme active site models of nitric oxide synthase (NOS) and address the key steps in the catalytic cycle whereby the substrate (L-arginine) is hydroxylated to N(omega)-hydroxo-arginine. It has been proposed that the mechanism follows a cytochrome P450-type catalytic cycle; however, our calculations find an alternative low energy pathway whereby the bound L-arginine substrate has two important functions in the catalytic cycle, namely first as a proton donor and later as the substrate in the reaction mechanism. Thus, the DFT studies show that the oxo-iron active species (compound I) cannot abstract a proton and neither a hydrogen atom from protonated L-arginine due to the strength of the N-H bonds of the substrate. However, the hydroxylation of neutral arginine by compound I and its one electron reduced form (compound II) requires much lower barriers and is highly exothermic. Detailed analysis of proton transfer mechanisms shows that the basicity of the dioxo dianion and the hydroperoxo-iron (compound 0) intermediates in the catalytic cycle are larger than that of arginine, which makes it likely that protonated arginine donates one of the two protons needed during the first catalytic cycle of NOS. Therefore, DFT predicts that in NOS enzymes arginine binds to the active site in its protonated form, but is deprotonated during the oxygen activation process in the catalytic cycle by either the dioxo dianion species or compound 0. As a result of the low ionization potential of neutral arginine, the actual hydroxylation reaction starts with an initial electron transfer from the substrate to compound I to create compound II followed by a concerted hydrogen abstraction/radical rebound from the substrate. These studies indicate that compound II is the actual oxidant in NOS enzymes that performs the hydroxylation reaction of arginine, which is in sharp contrast with the cytochromes P450 where compound II was shown to be a sluggish oxidant. This is the first example of an enzyme where compound II is able to participate in the reaction mechanism. Moreover, arginine hydroxylation by NOS enzymes is catalyzed in a significantly different way from the cytochromes P450 although the active sites of the two enzyme classes are very similar in structure. Detailed studies of environmental effects on the reaction mechanism show that environmental perturbations as appear in the protein have little effect and do not change the energies of the reaction. Finally, a valence bond curve crossing model has been set up to explain the obtained reaction mechanisms for the hydrogen abstraction processes in P450 and NOS enzymes.     

Catalytic mechanism of 6-phosphogluconate dehydrogenase: a theoretical investigation

Density functional calculations are employed to theoretically explore the mechanism of all elementary reaction steps involved in the catalytic reaction of 6-phosphogluconate dehydrogenase (6PGDH). The model systems we choose for the enzyme contain the essential parts of the cofactor (NADP+), the substrate 6-phosphogluconate (6PG), and some key residues (Lys183 and Glu190) in the active site of sheep liver 6PGDH. The effect of the apoenzyme electrostatic environment on the studied reaction is treated by the self-consistent reaction-field method. Our calculations demonstrate that the first step of the catalytic reaction is the formation of a 3-keto 6PG intermediate, which proceeds through a concerted transition state involving a hydride transfer from 6PG to NADP+, and a proton transfer from 6PG to Lys183. The second step is the elimination of a CO2 molecule from 6-PG, concomitant with a proton transfer from Lys183 to 6-PG. In the final step, a concerted double proton transfer (one from Glu190 to the substrate, another from the substrate to Lys183) results in the final product, the keto form of ribulose 5-phosphate (Ru5P). The rate-limiting step is the formation of a 3-keto 6PG intermediate, with a free energy barrier of 22.7 kcal/mol at room temperature in the protein environment, and all three steps are calculated to be thermodynamically favorable. These results are in good agreement with the general acid/general base mechanism suggested from previous experiments for the 6PGDH reaction.     

New insights into the mechanism of the Schiff base hydrolysis catalyzed by type I dehydroquinate dehydratase from S. enterica: a theoretical study

The reaction pathway of Schiff base hydrolysis catalyzed by type I dehydroquinate dehydratase (DHQD) from S. enterica has been studied by performing molecular dynamics (MD) simulations and density functional theory (DFT) calculations and the corresponding potential energy profile has also been identified. On the basis of the results, the catalytic hydrolysis process for the wild-type enzyme consists of three major reaction steps, including nucleophilic attack on the carbon atom involved in the carbon-nitrogen double bond of the Schiff base intermediate by a water molecule, deprotonation of the His143 residue, and dissociation between the product and the Lys170 residue of the enzyme. The remarkable difference between this and the previously proposed reaction mechanism is that the second step here, absent in the previously proposed reaction mechanism, plays an important role in facilitating the reaction through a key proton transfer by the His143 residue, resulting in a lower energy barrier. Comparison with our recently reported results on the Schiff base formation and dehydration processes clearly shows that the Schiff base hydrolysis is rate-determining in the overall reaction catalyzed by type I DHQD, consistent with the experimental prediction, and the calculated energy barrier of ～16.0 kcal mol(-1) is in good agreement with the experimentally derived activation free energy of ～14.3 kcal mol(-1). When the imidazole group of His143 residue is missing, the Schiff base hydrolysis is initiated by a hydroxide ion in the solution, rather than a water molecule, and both the reaction mechanism and the kinetics of Schiff base hydrolysis have been remarkably changed, clearly elucidating the catalytic role of the His143 residue in the reaction. The new mechanistic insights obtained here will be valuable for the rational design of high-activity inhibitors of type I DHQD as non-toxic antimicrobials, anti-fungals, and herbicides.     

Computer simulation of the chemical catalysis of DNA polymerases: discriminating between alternative nucleotide insertion mechanisms for T7 DNA polymerase

Understanding the chemical step in the catalytic reaction of DNA polymerases is essential for elucidating the molecular basis of the fidelity of DNA replication. The present work evaluates the free energy surface for the nucleotide transfer reaction of T7 polymerase by free energy perturbation/empirical valence bond (FEP/EVB) calculations. A key aspect of the enzyme simulation is a comparison of enzymatic free energy profiles with the corresponding reference reactions in water using the same computational methodology, thereby enabling a quantitative estimate for the free energy of the nucleotide insertion reaction. The reaction is driven by the FEP/EVB methodology between valence bond structures representing the reactant, pentacovalent intermediate, and the product states. This pathway corresponds to three microscopic chemical steps, deprotonation of the attacking group, a nucleophilic attack on the P(alpha) atom of the dNTP substrate, and departure of the leaving group. Three different mechanisms for the first microscopic step, the generation of the RO(-) nucleophile from the 3'-OH hydroxyl of the primer, are examined: (i) proton transfer to the bulk solvent, (ii) proton transfer to one of the ionic oxygens of the P(alpha) phosphate group, and (iii) proton transfer to the ionized Asp654 residue. The most favorable reaction mechanism in T7 pol is predicted to involve the proton transfer to Asp654. This finding sheds light on the long standing issue of the actual role of conserved aspartates. The structural preorganization that helps to catalyze the reaction is also considered and analyzed. The overall calculated mechanism consists of three subsequent steps with a similar activation free energy of about 12 kcal/mol. The similarity of the activation barriers of the three microscopic chemical steps indicates that the T7 polymerase may select against the incorrect dNTP substrate by raising any of these barriers. The relative height of these barriers comparing right and wrong dNTP substrates should therefore be a primary focus of future computational studies of the fidelity of DNA polymerases.     

Simulating the proton transfer in gramicidin A by a sequential dynamical Monte Carlo method

The large interest in long-range proton transfer in biomolecules is triggered by its importance for many biochemical processes such as biological energy transduction and drug detoxification. Since long-range proton transfer occurs on a microsecond time scale, simulating this process on a molecular level is still a challenging task and not possible with standard simulation methods. In general, the dynamics of a reactive system can be described by a master equation. A natural way to describe long-range charge transfer in biomolecules is to decompose the process into elementary steps which are transitions between microstates. Each microstate has a defined protonation pattern. Although such a master equation can in principle be solved analytically, it is often too demanding to solve this equation because of the large number of microstates. In this paper, we describe a new method which solves the master equation by a sequential dynamical Monte Carlo algorithm. Starting from one microstate, the evolution of the system is simulated as a stochastic process. The energetic parameters required for these simulations are determined by continuum electrostatic calculations. We apply this method to simulate the proton transfer through gramicidin A, a transmembrane proton channel, in dependence on the applied membrane potential and the pH value of the solution. As elementary steps in our reaction, we consider proton uptake and release, proton transfer along a hydrogen bond, and rotations of water molecules that constitute a proton wire through the channel. A simulation of 8 mus length took about 5 min on an Intel Pentium 4 CPU with 3.2 GHz. We obtained good agreement with experimental data for the proton flux through gramicidin A over a wide range of pH values and membrane potentials. We find that proton desolvation as well as water rotations are equally important for the proton transfer through gramicidin A at physiological membrane potentials. Our method allows to simulate long-range charge transfer in biological systems at time scales, which are not accessible by other methods.     

A theoretical evaluation of possible transition metal electro-catalysts for N2 reduction

Theoretical studies of the possibility of forming ammonia electrochemically at ambient temperature and pressure are presented. Density functional theory calculations were used in combination with the computational standard hydrogen electrode to calculate the free energy profile for the reduction of N(2) admolecules and N adatoms on several close-packed and stepped transition metal surfaces in contact with an acidic electrolyte. Trends in the catalytic activity were calculated for a range of transition metal surfaces and applied potentials under the assumption that the activation energy barrier scales with the free energy difference in each elementary step. The most active surfaces, on top of the volcano diagrams, are Mo, Fe, Rh, and Ru, but hydrogen gas formation will be a competing reaction reducing the faradaic efficiency for ammonia production. Since the early transition metal surfaces such as Sc, Y, Ti, and Zr bind N-adatoms more strongly than H-adatoms, a significant production of ammonia compared with hydrogen gas can be expected on those metal electrodes when a bias of -1 V to -1.5 V vs. SHE is applied. Defect-free surfaces of the early transition metals are catalytically more active than their stepped counterparts.     

Theoretical study toward understanding the catalytic mechanism of pyruvate decarboxylase

Density functional calculations are employed to explore the mechanisms of all elementary reaction steps involved in the catalytic cycle of pyruvate decarboxylase (PDC). Different models are constructed for mimicking the involvement of some key residues in a certain step. The effect of the protein framework on the potential energy profiles of active site models is approximately modeled by fixing some freedoms, based on the crystal structure of the PDC enzyme from Saccharomyces cerevisiae (ScPDC). Our calculations confirm that Glu51 is the most important residue in the formation of the ylide and the release of acetaldehyde via the proton relay between Glu51, N1', and the 4'-amino group of thiamine diphosphate. The presence of Glu477 and Asp28 residues makes the decarboxylation of lactylthiamin diphosphate (LThDP) an endothermic process with a significant free energy barrier. The protonation of the alpha-carbanion to form 2-(1-hydroxyethyl)-thiamin diphosphate is found to go through a concerted double proton transfer transition state involving both Asp28 and His115 residues. The final step, acetaldehyde release, is likely to proceed through a concerted transition state involving carbon-carbon bond-breaking and the deprotonation of the alpha-hydroxyl group. The decarboxylation of LThDP and the protonation of the alpha-carbanion are two rate-limiting steps, relative to the facile occurrence of the ylide formation and acetaldehyde release. The catalytic roles of residues Glu51, Glu477, Asp28, and Gly417 in the active site of ScPDC in individual steps elucidated from the present study are in good agreement with those derived from site-directed mutagenesis.     

Role of Asp102 in the catalytic relay system of serine proteases: a theoretical study

The role of Asp102 in the catalytic relay system of serine proteases is studied theoretically by calculating the free energy profiles of the single proton-transfer reaction by the Asn102 mutant trypsin and the concerted double proton-transfer reaction (so-called the charge-relay mechanism) of the wild-type trypsin. For each reaction, the reaction free energy profile of the rate-determining step (the tetrahedral intermediate formation step) is calculated by using ab initio QM/MM electronic structure calculations combined with molecular dynamics-free energy perturbation method. In the mutant reaction, the free energy monotonically increases along the reaction path. The rate-determining step of the mutant reaction is the formation of tetrahedral intermediate complex, not the base (His57) abstraction of the proton from Ser195. In contrast to the single proton-transfer reaction of the wild-type, MD simulations of the enzyme-substrate complex show that the catalytically favorable alignment of the relay system (the hydrogen bonding network between the mutant triad, His57, Asn102, and Ser195) is rarely observed even in the presence of a substrate at the active site. In the double proton-transfer reaction, the energy barrier is observed at the proton abstraction step, which corresponds to the rate-determining step of the single proton-transfer reaction of the wild-type. Although both reaction profiles show an increase of the activation barrier by several kcals/mol, these increases have different energetic origins: a large energetic loss of the electrostatic stabilization between His57 and Asn102 in the mutant reaction, while the lack of stabilization by the protein environment in the double proton-transfer reaction. Comparing the present results with the single proton transfer of the wild-type, Asp102 is proven to play two important roles in the catalytic process. One is to stabilize the protonated His57, or ionic intermediate, formed during the acylation, and the other is to fix the configuration around the active site, which is favorable to promote the catalytic process. These two factors are closely related to each other and are indispensable for the efficient catalysis. Also the present calculations suggest the importance of the remote site interaction between His57 and Val213-Ser214 at the catalytic transition state.     

Electrocatalytic reductions of nitrite, nitric oxide, and nitrous oxide by thermophilic cytochrome P450 CYP119 in film-modified electrodes and an analytical comparison of its catalytic activities with myoglobin

Previous investigations of nitrite and nitric oxide reduction by myoglobin in surfactant film modified electrodes characterized several distinct steps in the denitrification pathway, including isolation of a nitroxyl adduct similar to that proposed in the P450nor catalytic cycle. To investigate the effect of the axial ligand on these biomimetic reductions, we report here a comparison of the electrocatalytic activity of myoglobin (Mb) with a thermophilic cytochrome P450 CYP119. Electrocatalytic nitrite reduction by CYP119 is very similar to that by Mb: two catalytic waves at analogous potentials are observed, the first corresponding to the reduction of nitric oxide, the second to the production of ammonia. CYP119 is a much more selective catalyst, giving almost exclusively ammonia during the initial half-hour of reductive electrolysis of nitrite. More careful investigations of specific steps in the catalytic cycle show comparable rates of nitrite dehydration and almost identical potentials and lifetimes for ferrous nitroxyl intermediate (Fe(II)-NO(-)) in CYP119 and Mb. The catalytic efficiency of nitric oxide reduction is reduced for CYP119 as compared to Mb, attributable to both a lower affinity of the protein for NO and a decreased rate of N-N coupling. Isotopic labeling studies show ammonia incorporation into nitrous oxide produced during nitrite reduction, as has been termed co-denitrification for certain bacterial and fungal nitrite reductases. Mb has a much higher co-denitrification activity than CYP119. Conversely, CYP119 is shown to be slightly more efficient at the two-electron reduction of N(2)O to N(2). These results suggest that thiolate ligation does not significantly alter the catalytic reactivity, but the dramatic difference in product distribution may suggest an important role for protein stability in the selectivity of biocatalysts.     

Proton exit channels in bovine cytochrome c oxidase

Cytochrome c oxidase (CcO) is the terminal transmembrane enzyme of the respiratory electron transport chain in aerobic cells. It catalyzes the reduction of oxygen to water and utilizes the free energy of the reduction reaction for proton pumping, a process which results in a membrane electrochemical proton gradient. Although the structure of the enzyme has been solved for several organisms, the molecular mechanism of proton pumping and proton exit pathways remain unknown. In our previous work, the continuum electrostatic calculations were employed to evaluate the electrostatic potential, energies, and protonation state of bovine cytochrome c oxidase for different redox states of the enzyme. A possible mechanism of oxygen reduction and proton pumping via His291 was proposed. In this paper, using electrostatic calculations, we examine the proton exit pathways in the enzyme. By monitoring the changes of the protonation states, proton affinities, and energies of electrostatic interactions between the titratable groups in different redox states of CcO, we identified the clusters of strongly interacting residues. Using these data, we detected four possible proton exit points on the periplasmic side of the membrane (Lys171B/Asp173B, His24B/Asp25B, Asp51, and Asp300). We then were able to trace the proton exit pathways and to evaluate the energy profiles along the paths. On the basis of energetic considerations and the conservation of the residues in a protein sequence, the most likely exit pathway is one via the Lys171B/Asp173B site. The obtained results are fully consistent with our His291 model of proton pumping, and provide a rationale for the absence of proton leaking in CcO between the pumping strokes.     

The reaction mechanism of xanthine oxidase: evidence for two-electron chemistry rather than sequential one-electron steps

Current research on xanthine oxidase has favored a mechanism involving base-catalyzed proton abstraction from a Mo-OH group, allowing nucleophilic attack on the substrate and hydride transfer from the substrate to Mo=S group in the active site. During the course of this reaction mechanism, the molybdenum redox cycles from MoVI to MoIV, with reoxidation of the MoIV speices to form the EPR active MoV intermediate. However, it has also been suggested that the reaction occurs in two subsequent one-electron steps. We have determined kinetic parameters kred and kred/Kd for a variety of plausible substrates as well as the one-electron reduction potentials for these substrates. Our data indicate no correlation between these kinetic parameters and their one-electron reduction potentials, as would be expected if the enzyme were using two subsequent one-electron reduction steps. Our results provide additional support to current evidence for the favored two-electron reduction mechanism.     

Co9S8 as a catalyst for electroreduction of O2: quantum chemistry predictions

The slab band quantum computational approach in the Vienna ab initio simulation package (VASP) is used to calculate the adsorption energies of reactants, reaction intermediates, and products in O2 reduction and in water oxidation in acid on three crystallographic surfaces of pentlandite structure Co9S8. Reversible potentials for the reaction steps involving electron and proton transfer are determined by using the energies in a linear Gibbs free energy relationship. On the basis of these results, we find that the partially OH-covered (202) surface is active toward O2 reduction and should have overpotential behavior similar to that observed for platinum electrodes. One structure in the predicted four-electron reduction mechanism is novel: S2- provides an adsorption site for O following O-O bond scission, which, unlike the case of platinum electrodes, takes place prior to the first reduction step.     

铜催化CO2和1-苯基丙炔的氢羧基化反应机理及区域选择性

采用密度泛函理论(DFT)方法研究了在还原剂(EtO)₃SiH存在下, 铜(I) (Cl₂IPrCuF)催化CO₂插入1-苯基丙炔生成α,β不饱和羧酸的反应机理. 计算结果表明, Cl₂IPrCuF 首先与(EtO)₃SiH 生成活性催化剂Cl₂IPrCuH,然后经历三个步骤完成催化反应: (1) Cl₂IPrCuH 与1-苯基丙炔加成生成烯基铜中间体. 由于炔烃的不对称性,烯基铜中间体有两种同分异构体, 最后可导致生成两种对应的α,β不饱和羧酸衍生物; (2) CO₂插入烯基铜中间体得到羧基铜中间体; (3) (EtO)₃SiH 与羧基铜中间体发生σ转位反应形成最终产物, 同时重新生成催化剂Cl₂IPrCuH. 理论研究还表明, 生成两种α,β不饱和羧酸衍生物的反应路径所对应的决速步骤不同, 在Path a 中炔烃插入反应和CO₂插入反应都可能是整个催化反应的决速步骤, 自由能垒分别为68.6 和67.8 kJ·mol^-1, 而在Path b中, 仅炔烃插入反应是整个催化反应的决速步骤, 自由能垒为78.7 kJ·mol^-1. 此结果很好地给出了实验上两种α,β不饱和羧酸衍生物收率不同的原因. 炔烃与Cl₂IPrCuH的加成决定了反应的区域选择性, 其中电子效应是影响反应区域选择性的主要原因.     

Peptide hydrolysis catalyzed by matrix metalloproteinase 2: a computational study

The MMP-2 reaction mechanism is investigated by using different computational methodologies. First, quantum mechanical (QM) calculations are carried out on a cluster model of the active site bound to an Ace-Gly approximately Ile-Nme peptide. Along the QM reaction path, a Zn-bound water molecule attacks the Gly carbonyl group to give a tetrahedral intermediate. The breaking of the C-N bond is completed thanks to the Glu 404 residue that shuttles a proton from the water molecule to Ile-N atom. The gas-phase QM energy barrier is quite low ( approximately 14 kcal/mol), thus suggesting that the essential catalytic machinery is included in the cluster model. A similar reaction path occurs in the MMP-2 catalytic domain bound to an octapeptide substrate according to hybrid QM and molecular mechanical (QM/MM) geometry optimizations. However, the rupture of the Gly( P 1) approximately Ile( P 1') amide bond is destabilized in the static QM/MM calculations, owing to the positioning of the Ile( P 1') side chain inside the MMP-2 S 1' pocket and to the inability of simple energy miminization methodologies to properly relax complex systems. Molecular dynamics simulations show that these steric limitations are overcome easily through structural fluctuations. The energetic effect of structural fluctuations is taken into account by combining QM energies with average MM Poisson-Boltzmann free energies, resulting in a total free energy barrier of 14.8 kcal/mol in good agreement with experimental data. The rate-determining event in the MMP-2 mechanism corresponds to a H-bond rearrangement involving the Glu 404 residue and/or the Glu 404-COOH --> N-Ile( P 1') proton transfer. Overall, the present computational results and previous experimental data complement each other well in order to provide a detailed view of the MMPs catalytic mechanism.