Atomic Description of an Enzyme Reaction Dominated by Proton Tunneling

We present an atomic-level description of the reaction chemistry of an enzyme-catalyzed reaction dominated by proton tunneling. By solving structures of reaction intermediates at near-atomic resolution, we have identified the reaction pathway for tryptamine oxidation by aromatic amine dehydrogenase. Combining experiment and computer simulation, we show proton transfer occurs predominantly to oxygen O2 of Asp(128)beta in a reaction dominated by tunneling over approximately 0.6 angstroms. The role of long-range coupled motions in promoting tunneling is controversial. We show that, in this enzyme system, tunneling is promoted by a short-range motion modulating proton-acceptor distance and no long-range coupled motion is required.

The Fe-CO bond energy in myoglobin: a QM/MM study of the effect of tertiary structure

The Fe-CO bond dissociation energy (BDE) in myoglobin (Mb) has been calculated with B3LYP quantum mechanics/molecular mechanics methods for 22 different Mb conformations, generated from molecular dynamics simulations. Our average BDE of 8.1 kcal/mol agrees well with experiment and shows that Mb weakens the Fe-CO bond by 5.8 kcal/mol; the calculations provide detailed atomistic insight into the origin of this effect. BDEs for Mb conformations with the R carbonmonoxy tertiary structure are on average 2.6 kcal/mol larger than those with the T deoxy tertiary structure, suggesting two functionally distinct allosteric states. This allostery is partly explained by the reduction in distal cavity steric crowding as Mb moves from its T to R tertiary structure.

Molecular mechanisms of antibiotic resistance: QM/MM modelling of deacylation in a class A beta-lactamase

Modelling of the first step of the deacylation reaction of benzylpenicillin in the E. coli TEM1 beta-lactamase (with B3LYP/6-31G + (d)//AM1-CHARMM22 quantum mechanics/molecular mechanics methods) shows that a mechanism in which Glu166 acts as the base to deprotonate a conserved water molecule is both energetically and structurally consistent with experimental data; the results may assist the design of new antibiotics and beta-lactamase inhibitors.

Electronic structure of compound I in human isoforms of cytochrome P450 from QM/MM modeling

Human cytochromes P450 play a vital role in drug metabolism. The key step in substrate oxidation involves hydrogen atom abstraction or C=C bond addition by the oxygen atom of the Compound I intermediate. The latter has three unpaired electrons, two on the Fe-O center and one shared between the porphyrin ring and the proximal cysteinyl sulfur atom. Changes in its electronic structure have been suggested to affect reactivity. The electronic and geometric structure of Compound I in three important human subfamilies of cytochrome P450 (P450, 2C, 2B, and 3A) that are major contributors to drug metabolism is characterized here using combined quantum mechanical/molecular mechanical (QM/MM) calculations at the B3LYP:CHARMM27 level. Compound I is remarkably similar in all isoforms, with the third unpaired electron located mainly on the porphyrin ring, and this prediction is not very sensitive to details of the QM/MM methodology, such as the DFT functional, the basis set, or the size of the QM region. The presence of substrate also has no effect. The main source of variability in spin density on the cysteinyl sulfur (from 26 to 50%) is the details of the system setup, such as the starting protein geometry used for QM/MM minimization. This conformational effect is larger than the differences between human isoforms, which are therefore not distinguishable on electronic grounds, so it is unlikely that observed large differences in substrate selectivity can be explained to a large extent in these terms.

Modelling enzyme reaction mechanisms, specificity and catalysis

Modern modelling methods can now give uniquely detailed understanding of enzyme-catalyzed reactions, including the analysis of mechanisms and the identification of determinants of specificity and catalytic efficiency. A new field of computational enzymology has emerged that has the potential to contribute significantly to structure-based design and to develop predictive models of drug metabolism and, for example, of the effects of genetic polymorphisms. This review outlines important techniques in this area, including quantum-chemical model studies and combined quantum-mechanics and molecular-mechanics (QM/MM) methods. Some recent applications to enzymes of pharmacological interest are also covered, showing the types of problems that can be tackled and the insight they can give.

Multiple high-level QM/MM reaction paths demonstrate transition-state stabilization in chorismate mutase: correlation of barrier height with transition-state stabilization

Multiple profiles for the reaction from chorismate to prephenate in the enzyme chorismate mutase calculated with hybrid density functional combined quantum mechanics/molecular mechanics methods (B3LYP/6-31G(d)-CHARMM27) agree well with experiment, and provide direct evidence of transition-state stabilization by this important enzyme, which is at the centre of current debates about the nature of enzyme catalysis.

Mechanisms of Antibiotic Resistance: QM/MM Modeling of the Acylation Reaction of a Class A β-Lactamase with Benzylpenicillin

Understanding the mechanisms by which beta-lactamases destroy beta-lactam antibiotics is potentially vital in developing effective therapies to overcome bacterial antibiotic resistance. Class A beta-lactamases are the most important and common type of these enzymes. A key process in the reaction mechanism of class A beta-lactamases is the acylation of the active site serine by the antibiotic. We have modeled the complete mechanism of acylation with benzylpenicillin, using a combined quantum mechanical and molecular mechanical (QM/MM) method (B3LYP/6-31G+(d)//AM1-CHARMM22). All active site residues directly involved in the reaction, and the substrate, were treated at the QM level, with reaction energies calculated at the hybrid density functional (B3LYP/6-31+Gd) level. Structures and interactions with the protein were modeled by the AM1-CHARMM22 QM/MM approach. Alternative reaction coordinates and mechanisms have been tested by calculating a number of potential energy surfaces for each step of the acylation mechanism. The results support a mechanism in which Glu166 acts as the general base. Glu166 deprotonates an intervening conserved water molecule, which in turn activates Ser70 for nucleophilic attack on the antibiotic. This formation of the tetrahedral intermediate is calculated to have the highest barrier of the chemical steps in acylation. Subsequently, the acylenzyme is formed with Ser130 as the proton donor to the antibiotic thiazolidine ring, and Lys73 as a proton shuttle residue. The presented mechanism is both structurally and energetically consistent with experimental data. The QM/MM energy barrier (B3LYP/ 6-31G+(d)//AM1-CHARMM22) for the enzymatic reaction of 9 kcal mol(-1) is consistent with the experimental activation energy of about 12 kcal mol(-1). The effects of essential catalytic residues have been investigated by decomposition analysis. The results demonstrate the importance of the "oxyanion hole" in stabilizing the transition state and the tetrahedral intermediate. In addition, Asn132 and a number of charged residues in the active site have been identified as being central to the stabilizing effect of the enzyme. These results will be potentially useful in the development of stable beta-lactam antibiotics and for the design of new inhibitors.

QM/MM studies of the electronic structure of the compound I intermediate in cytochrome c peroxidase and ascorbate peroxidase

Cytochrome c peroxidase (CcP) and ascorbate peroxidase (APX) both involve reactive haem oxoferryl intermediates known as 'compound I' species. These two enzymes also have a very similar structure, especially in the vicinity of the haem group. Despite this similarity, the electronic structure of compound I in the two enzymes is known to be very different. Compound I intermediates have three unpaired electrons, two of which are always situated on the Fe-O core, whilst the third is located in a porphyrin orbital in APX and many other compound I species. In CcP, however, this third unpaired electron is positioned on a tryptophan residue lying close to the haem ring. The same residue is present in the same position in APX, yet it is not oxidized in that case. We report QM/MM calculations, using accurate B3LYP density functional theory for the QM region, on the active intermediate for both enzymes. We reproduce the observed difference in electronic structure, and show that it arises as a result of subtle electrostatic effects which affect the ionization potential of both the tryptophan and porphyrin groups. The computed structures of both enzymes do not involve deprotonation of the tryptophan group, or protonation of the oxoferryl oxygen.

QM/MM modelling of oleamide hydrolysis in fatty acid amide hydrolase (FAAH) reveals a new mechanism of nucleophile activation

Fatty acid amide hydrolase (FAAH), a promising target for the treatment of several central and peripheral nervous system disorders, such as anxiety, pain and hypertension, has an unusual catalytic site, and its mechanism has been uncertain; hybrid quantum mechanics/molecular mechanics (QM/MM) calculations reveal a new mechanism of nucleophile activation (involving a Lys-Ser-Ser catalytic triad), with potentially crucial insights for the design of potent and selective inhibitors.

Differential Transition-State Stabilization in Enzyme Catalysis: Quantum Chemical Analysis of Interactions in the Chorismate Mutase Reaction and Prediction of the Optimal Catalytic Field

Chorismate mutase is a key model system in the development of theories of enzyme catalysis. To analyze the physical nature of catalytic interactions within the enzyme active site and to estimate the stabilization of the transition state (TS) relative to the substrate (differential transition state stabilization, DTSS), we have carried out nonempirical variation-perturbation analysis of the electrostatic, exchange, delocalization, and correlation interactions of the enzyme-bound substrate and transition-state structures derived from ab initio QM/MM modeling of Bacillus subtilis chorismate mutase. Significant TS stabilization by approximately -23 kcal/mol [MP2/6-31G(d)] relative to the bound substrate is in agreement with that of previous QM/MM modeling and contrasts with suggestions that catalysis by this enzyme arises purely from conformational selection effects. The most important contributions to DTSS come from the residues, Arg90, Arg7, Glu78, a crystallographic water molecule, Arg116, and Arg63, and are dominated by electrostatic effects. Analysis of the differential electrostatic potential of the TS and substrate allows calculation of the catalytic field, predicting the optimal location of charged groups to achieve maximal DTSS. Comparison with the active site of the enzyme from those of several species shows that the positions of charged active site residues correspond closely to the optimal catalytic field, showing that the enzyme has evolved specifically to stabilize the TS relative to the substrate.

Mechanism and structure–reactivity relationships for aromatic hydroxylation by cytochrome P450

Cytochrome P450 enzymes play a central role in drug metabolism, and models of their mechanism could contribute significantly to pharmaceutical research and development of new drugs. The mechanism of cytochrome P450 mediated hydroxylation of aromatics and the effects of substituents on reactivity have been investigated using B3LYP density functional theory computations in a realistic porphyrin model system. Two different orientations of substrate approach for addition of Compound I to benzene, and also possible subsequent rearrangement pathways have been explored. The rate-limiting Compound I addition to an aromatic carbon atom proceeds on the doublet potential energy surface via a transition state with mixed radical and cationic character. Subsequent formation of epoxide, ketone and phenol products is shown to occur with low barriers, especially starting from a cation-like rather than a radical-like tetrahedral adduct of Compound I with benzene. Effects of ring substituents were explored by calculating the activation barriers for Compound I addition in the meta and para-position for a range of monosubstituted benzenes and for more complex polysubstituted benzenes. Two structure-reactivity relationships including 8 and 10 different substituted benzenes have been determined using (i) experimentally derived Hammett sigma-constants and (ii) a theoretical scale based on bond dissociation energies of hydroxyl adducts of the substrates, respectively. In both cases a dual-parameter approach that employs a combination of radical and cationic electronic descriptors gave good relationships with correlation coefficients R2 of 0.96 and 0.82, respectively. These relationships can be extended to predict the reactivity of other substituted aromatics, and thus can potentially be used in predictive drug metabolism models.

Transition state stabilization and substrate strain in enzyme catalysis: ab initio QM/MM modelling of the chorismate mutase reaction

To investigate fundamental features of enzyme catalysis, there is a need for high-level calculations capable of modelling crucial, unstable species such as transition states as they are formed within enzymes. We have modelled an important model enzyme reaction, the Claisen rearrangement of chorismate to prephenate in chorismate mutase, by combined ab initio quantum mechanics/molecular mechanics (QM/MM) methods. The best estimates of the potential energy barrier in the enzyme are 7.4-11.0 kcal mol(-1)(MP2/6-31+G(d)//6-31G(d)/CHARMM22) and 12.7-16.1 kcal mol(-1)(B3LYP/6-311+G(2d,p)//6-31G(d)/CHARMM22), comparable to the experimental estimate of Delta H(++)= 12.7 +/- 0.4 kcal mol(-1). The results provide unequivocal evidence of transition state (TS) stabilization by the enzyme, with contributions from residues Arg90, Arg7, and Arg63. Glu78 stabilizes the prephenate product (relative to substrate), and can also stabilize the TS. Examination of the same pathway in solution (with a variety of continuum models), at the same ab initio levels, allows comparison of the catalyzed and uncatalyzed reactions. Calculated barriers in solution are 28.0 kcal mol(-1)(MP2/6-31+G(d)/PCM) and 24.6 kcal mol(-1)(B3LYP/6-311+G(2d,p)/PCM), comparable to the experimental finding of Delta G(++)= 25.4 kcal mol(-1) and consistent with the experimentally-deduced 10(6)-fold rate acceleration by the enzyme. The substrate is found to be significantly distorted in the enzyme, adopting a structure closer to the transition state, although the degree of compression is less than predicted by lower-level calculations. This apparent substrate strain, or compression, is potentially also catalytically relevant. Solution calculations, however, suggest that the catalytic contribution of this compression may be relatively small. Consideration of the same reaction pathway in solution and in the enzyme, involving reaction from a 'near-attack conformer' of the substrate, indicates that adoption of this conformation is not in itself a major contribution to catalysis. Transition state stabilization (by electrostatic interactions, including hydrogen bonds) is found to be central to catalysis by the enzyme. Several hydrogen bonds are observed to shorten at the TS. The active site is clearly complementary to the transition state for the reaction, stabilizing it more than the substrate, so reducing the barrier to reaction.

Conformational effects in enzyme catalysis: QM/MM free energy calculation of the ‘NAC’ contribution in chorismate mutase

The controversial 'near attack conformation'(NAC) effect in the important model enzyme chorismate mutase is calculated to be 3.8-4.6 kcal mol(-1) by QM/MM free energy perturbation molecular dynamics methods, showing that the NAC effect by itself does not account for catalysis in this enzyme.

Aromatic Hydroxylation by Cytochrome P450: Model Calculations of Mechanism and Substituent Effects

The mechanism and selectivity of aromatic hydroxylation by cytochrome P450 enzymes is explored using new B3LYP density functional theory computations. The calculations, using a realistic porphyrin model system, show that rate-determining addition of compound I to an aromatic carbon atom proceeds via a transition state with partial radical and cationic character. Reactivity is shown to depend strongly on ring substituents, with both electron-withdrawing and -donating groups strongly decreasing the addition barrier in the para position, and it is shown that the calculated barrier heights can be reproduced by a new dual-parameter equation based on radical and cationic Hammett sigma parameters.

Modeling biotransformation reactions by combined quantum mechanical/molecular mechanical approaches: from structure to activity

An overview of the combined quantum mechanical/molecular mechanical (QM/MM) approach and its application to studies of biotransformation enzymes and drug metabolism is given. Theoretical methods to simulate enzymatic reactions have rapidly developed during the last decade. In particular, QM/MM methods provide detailed insights into enzyme catalyzed reactions, which can be extremely valuable in complementing experimental research. QM/MM methods allow the reacting groups in the active site of an enzyme to be studied at a quantum mechanical level, while the surrounding protein and solvent is included at a classical (and computationally less expensive) molecular mechanical level. Existing QM/MM implementations vary in the level of interaction between the QM and MM regions and in the way the partitioning into QM and MM regions is setup. Some general considerations concerning reaction modeling are discussed and a number of QM/MM studies related to drug metabolism are described. These studies illustrate that theoretical modeling of important metabolic reactions provides detailed insights into mechanisms of reaction and specific catalytic effects of enzyme residues as well as explaining variation in rates of conversion of different metabolites. Such information is essential in the development of methods to predict metabolism of drugs and to understand metabolic effects of genetic polymorphism in biotransformation enzymes.

Identification of Glu166 as the general base in the acylation reaction of class A beta-lactamases through QM/MM modeling

Bacterial class A beta-lactamases are responsible for the most known resistance against beta-lactam antibiotics. With the continuing rise in antibiotic resistance, improved knowledge of the mechanisms of action of these enzymes is needed in the development of effective therapeutic agents and strategies. The mechanism of the deacylation step in class A beta-lactamases is well accepted. In contrast, the mechanism of the acylation step has been uncertain, with several conflicting proposals put forward. We have modeled the acylation step in a class A beta-lactamase, using a combined quantum mechanics/molecular mechanics approach. The results provide an atomic level description of the reaction and show that Glu166 acts as the general base in the reaction, deprotonating Ser70 via an intervening water molecule. Ser70 acts as the nucleophile for attack on the lactam ring in a concerted reaction. The results do not rule out the importance of Lys73 in catalysis, in agreement with experimental data.

Quantum mechanical/molecular mechanical free energy simulations of the glutathione S-transferase (M1-1) reaction with phenanthrene 9,10-oxide

Glutathione S-transferases (GSTs) play an important role in the detoxification of xenobiotics in mammals. They catalyze the conjugation of glutathione to a wide range of electrophilic compounds. Phenanthrene 9,10-oxide is a model substrate for GSTs, representing an important group of epoxide substrates. In the present study, combined quantum mechanical/molecular mechanical (QM/MM) simulations of the conjugation of glutathione to phenanthrene 9,10-oxide, catalyzed by the M1-1 isoenzyme from rat, have been carried out to obtain insight into details of the reaction mechanism and the role of solvent present in the highly solvent accessible active site. Reaction-specific AM1 parameters for sulfur have been developed to obtain an accurate modeling of the reaction, and QM/MM solvent interactions in the model have been calibrated. Free energy profiles for the formation of two diastereomeric products were obtained from molecular dynamics simulations of the enzyme, using umbrella sampling and weighted histogram analysis techniques. The barriers (20 kcal/mol) are in good agreement with the overall experimental rate constant and with the formation of equal amounts of the two diastereomeric products, as experimentally observed. Along the reaction pathway, desolvation of the thiolate sulfur of glutathione is observed, in agreement with solvent isotope experiments, as well as increased solvation of the epoxide oxygen of phenanthrene 9,10-oxide, illustrating an important stabilizing role for active site solvent molecules. Important active site interactions have been identified and analyzed. The catalytic effect of Tyr115 through a direct hydrogen bond with the epoxide oxygen of the substrate, which was proposed on the basis of the crystal structure of the (9S,10S) product complex, is supported by the simulations. The indirect interaction through a mediating water molecule, observed in the crystal structure of the (9R,10R) product complex, cannot be confirmed to play a role in the conjugation step. A selection of mutations is modeled. The Asn8Asp mutation, representing one of the differences between the M1-1 and M2-2 isoenzymes, is identified as a possible factor contributing to the difference in the ratio of product formation by these two isoenzymes. The QM/MM reaction pathway simulations provide new and detailed insight into the reaction mechanism of this important class of detoxifying enzymes and illustrate the potential of QM/MM modeling to complement experimental data on enzyme reaction mechanisms.

Wave propagation in 0-3/3-3 connectivity composites with complex microstructure

This work presents a study of the properties of particulate composites. The whole range of particle volume fraction (0-1) and ideal 0-3, 3-3 and intermediate 0-3/3-3 connectivities are analysed. Two different approaches to produce a realistic model of the complex microstructure of the composites are considered. The first one is based on a random location of mono-dispersed particles in the matrix; while the second incorporates a size distribution of the particles based on experimental measurements. Different particle shapes are also considered. A commercial finite element package was used to study the propagation of acoustic plane waves through the composite materials. Due to the complexity of the problem, and as a first step, a two-dimensional model was adopted. The results obtained for the velocity of sound propagation from the finite element technique are compared with those from other theoretical approaches and with experimental data. The study validates the use of this technique to model acoustic wave propagation in 0-3/3-3 connectivity composites. In addition, the finite element calculations, along with the detailed description of the microstructure of the composite, provide valuable information about the micromechanics of the sample and the influence of the microstructure on macroscopic properties.

Combined quantum mechanical and molecular mechanical reaction pathway calculation for aromatic hydroxylation by p-hydroxybenzoate-3-hydroxylase

The reaction pathway for the aromatic 3-hydroxylation of p-hydroxybenzoate by the reactive C4a-hydroperoxyflavin cofactor intermediate in p-hydroxybenzoate hydroxylase (PHBH) has been investigated by a combined quantum mechanical and molecular mechanical (QM/MM) method. A structural model for the C4a-hydroperoxyflavin intermediate in the PHBH reaction cycle was built on the basis of the crystal structure coordinates of the enzyme-substrate complex. A reaction pathway for the subsequent hydroxylation step was calculated by imposing a reaction coordinate that involves cleavage of the peroxide oxygen-oxygen bond and formation of the carbon-oxygen bond between the C3 atom of the substrate and the distal oxygen of the peroxide moiety of the cofactor. The geometric changes and the Mulliken charge distributions along the calculated reaction pathway are in line with an electrophilic aromatic substitution type of mechanism. The energy barrier of the calculated reaction is considerably lower when the substrate hydroxyl moiety is deprotonated, in comparison with the barrier found with a protonated hydroxyl moiety. This effect of the protonation state of the substrate on the calculated energy barrier supports experimental observations that deprotonation is required for hydroxylation of the substrate. A notable event in the calculated reaction pathway is a lengthening of the peroxide oxygen-oxygen bond at an intermediate stage. Further analysis of the reaction pathway indicates that this oxygen-oxygen bond elongation is accompanied by an increase in electrophilic reactivity on the distal oxygen of the peroxide moiety, which may assist the C-O bond formation in the reaction of the C4a-hydroperoxyflavin intermediate with the substrate. Analysis of the effect of individual active site residues on the reaction reveals a specific transition state stabilization by the backbone carbonyl moiety of Pro293. The crystal water 717 appears to drive the hydroxylation step through a stabilizing hydrogen bond interaction to the proximal oxygen of the C4a-hydroperoxyflavin intermediate, which increases in strength as the hydroperoxyflavin cofactor converts to the anionic (deprotonated) hydroxyflavin.