Search

Decapping is a central step in eukaryotic mRNA turnover and in gene expression regulation. The human scavenger decapping enzyme, DcpS, catalyses cap hydrolysis following mRNA degradation. DcpS is a dimeric enzyme, with two active sites. Crystal structures suggest that DcpS must undergo significant conformational changes upon ligand binding, but the mechanism of this transition is unknown. Here, we report two long timescale (20 ns) molecular dynamics simulations of the apo-form of DcpS. The dimer is observed to undergo a strikingly cooperative motion, with one active site closing while the other opens. The amplitude of the conformational change is 6-21 A and the apparent timescale is 4-13 ns. These findings indicate that the crystallographically observed symmetric conformation of apo-form of DcpS is only a minor conformation in solution. The simulations also show that active sites are structurally connected via the domain-swapped dimer structure of the N-terminal domain, even in the absence of a bound ligand. These findings suggest a functional reason for the enzyme existing as a dimer, and may be widely relevant, also for other dimeric proteins., ((c) 2007 Wiley-Liss, Inc.)

Simulations and modelling [e.g. with combined QM/MM (quantum mechanics/molecular mechanics) methods] are increasingly important in investigations of enzyme-catalysed reaction mechanisms. Calculations offer the potential of uniquely detailed, atomic-level insight into the fundamental processes of biological catalysis. Highly accurate methods promise quantitative comparison with experiments, and reliable predictions of mechanisms, revolutionizing enzymology.

Modelling of the mechanism of covalent adduct formation by the inhibitor O-arylcarbamate URB524 in FAAH shows that only one of the two possible inhibitor binding orientations is consistent with the experimentally observed irreversible carbamoylation of the nucleophile serine: this is a potentially crucial insight for designing new covalent inhibitors of this promising drug target.

The combination of quantum mechanics (QM) with molecular mechanics (MM) offers a route to improved accuracy in the study of biological systems, and there is now significant research effort being spent to develop QM/MM methods that can be applied to the calculation of relative free energies. Currently, the computational expense of the QM part of the calculation means that there is no single method that achieves both efficiency and rigor; either the QM/MM free energy method is rigorous and computationally expensive, or the method introduces efficiency-led assumptions that can lead to errors in the result, or a lack of generality of application. In this paper we demonstrate a combined approach to form a single, efficient, and, in principle, exact QM/MM free energy method. We demonstrate the application of this method by using it to explore the difference in hydration of water and methane. We demonstrate that it is possible to calculate highly converged QM/MM relative free energies at the MP2/aug-cc-pVDZ/OPLS level within just two days of computation, using commodity processors, and show how the method allows consistent, high-quality sampling of complex solvent configurational change, both when perturbing hydrophilic water into hydrophobic methane, and also when moving from a MM Hamiltonian to a QM/MM Hamiltonian. The results demonstrate the validity and power of this methodology, and raise important questions regarding the compatibility of MM and QM/MM forcefields, and offer a potential route to improved compatibility.

This paper investigates the use of magnetically active materials in the matching layer of a piezoelectric transducer. This then allows the performance of the device to be dynamically altered by applying an external field. The effect that this new matching layer has on the performance of a typical device is theoretically investigated here. It transpires that the additional flexibility of an active matching layer can be used to maintain the efficiency of the device as the external load is varied.

The large number of degrees of freedom in the design of piezoelectric transducers requires a theoretical model that is computationally efficient so that a large number of iterations can be performed in the design optimisation. The materials used are often lossy, and indeed loss can be used to enhance the operational characteristics of these designs. Motivated by these needs, this paper extends the one-dimensional linear systems model to incorporate frequency dependent elastic loss. The reception sensitivity, electrical impedance and electromechanical coupling coefficient of a 1-3 composite transducer, with frequency dependent loss in the polymer filler, are investigated. By plotting these operating characteristics as a function of the volume fraction of piezoelectric ceramic an optimum design is obtained. A device with a non-standard, high shear attenuation polymer is also simulated and this leads to an increase in the electromechanical coupling coefficient. A comparison with finite element simulations is then performed. This shows that the two methods are in reasonable agreement in their electrical impedance profiles in all the cases considered. The plots are almost identical away from the main resonant peak where the frequency location of the peaks are comparable but there is in some cases a 20% discrepancy in the magnitude of the peak value and in its bandwidth. The finite element model also shows that the use of a high shear attenuation polymer filler damps out the unwanted, low frequency modes whilst maintaining a reasonable impedance magnitude.

Citrate synthase is an archetypal carbon-carbon bond forming enzyme. It promotes the conversion of oxaloacetate (OAA) to citrate by catalyzing the deprotonation (enolization) of acetyl-CoA, followed by nucleophilic attack of the enolate form of this substrate on OAA to form a citryl-CoA intermediate and subsequent hydrolysis. OAA is strongly bound to the active site and its alpha-carbonyl group is polarized. This polarization has been demonstrated spectroscopically, [(Kurz et al., Biochemistry 1985;24:452-457; Kurz and Drysdale, Biochemistry 1987;26:2623-2627)] and has been suggested to be an important catalytic strategy. Substrate polarization is believed to be important in many enzymes. The first step, formation of the acetyl-CoA enolate intermediate, is thought to be rate-limiting in the mesophilic (pig/chicken) enzyme. We have examined the effects of substrate polarization on this key step using quantum mechanical/molecular mechanical (QM/MM) methods. Free energy profiles have been calculated by AM1/CHARMM27 umbrella sampling molecular dynamics (MD) simulations, together with potential energy profiles. To study the influence of OAA polarization, profiles were calculated with different polarization of the OAA alpha-carbonyl group. The results indicate that OAA polarization influences catalysis only marginally but has a larger effect on intermediate stabilization. Different levels of treatment of OAA are compared (MM or QM), and its polarization in the protein and in water analyzed at the B3LYP/6-31+G(d)/CHARMM27 level. Analysis of stabilization by individual residues shows that the enzyme mainly stabilizes the enolate intermediate (not the transition state) through electrostatic (including hydrogen bond) interactions: these contribute much more than polarization of OAA., ((c) 2007 Wiley-Liss, Inc.)

The link atom (LA) and the generalized hybrid orbital (GHO) quantum mechanical/molecular mechanics (QM/MM) boundary treatment methods are compared, in the context of the acylation process (the rate-limiting step) involving the NS3/NS4A HCV serine protease and its NS5A/5B natural substrate. The potential energy surface was calculated, and the free energy along the selected reaction coordinate was obtained from umbrella sampling molecular dynamics simulations, at the AM1/CHARMM27 level. The LA and GHO methods, when properly applied, lead to similar chemical behavior, although the agreement is not quantitative. The choice of QM/MM partitioning is dictated to some extent by the nature of the two different methods, and this influences the results. The free energy profiles obtained by umbrella sampling suggest that the GHO approach is better suited for this system, because it provides a consistent description of the reaction in both the forward and backward directions. This is probably a consequence of the different QM/MM partitioning required by the two different methods (i.e., different numbers of atoms have to be included in the QM region). This finding is therefore likely to be system dependent, so careful testing should be considered for each enzyme application.

The first step of the reaction catalysed by the enzyme citrate synthase is studied here with high level combined quantum mechanical/molecular mechanical (QM/MM) methods (up to the MP2/6-31+G(d)//6-31G(d)/CHARMM level). In the first step of the reaction, acetyl-CoA is deprotonated by Asp375, producing an intermediate, which is the nucleophile for attack on the second substrate, oxaloacetate, prior to hydrolysis of the thioester bond of acetyl-CoA and release of the products. A central question has been whether the nucleophilic intermediate is the enolate of acetyl-CoA, the enol, or an 'enolic' intermediate stabilized by a 'low-barrier' hydrogen bond with His274 at the active site. The imidazole sidechain of His274 is neutral, and donates a hydrogen bond to the carbonyl oxygen of acetyl-CoA in substrate complexes. We have investigated the identity of the nucleophilic intermediate by QM/MM calculations on the substrate (keto), enolate, enol and enolic forms of acetyl-CoA at the active site of citrate synthase. The transition states for proton abstraction from acetyl-CoA by Asp375, and for transfer of the hydrogen bonded proton between His274 and acetyl-CoA have been modelled approximately. The effects of electron correlation are included by MP2/6-31G(d) and MP2/6-31+G(d) calculations on active site geometries produced by QM/MM energy minimization. The results do not support the hypothesis that a low-barrier hydrogen bond is involved in catalysis in citrate synthase, in agreement with earlier calculations. The acetyl-CoA enolate is identified as the only intermediate consistent with the experimental barrier for condensation, stabilized by conventional hydrogen bonds from His274 and a water molecule.

The hydrogen-transfer reaction catalysed by methylamine dehydrogenase (MADH) with methylamine (MA) as substrate is a good model system for studies of proton tunnelling in enzyme reactions--an area of great current interest--for which atomistic simulations will be vital. Here, we present a detailed analysis of the key deprotonation step of the MADH/MA reaction and compare the results with experimental observations. Moreover, we compare this reaction with the related aromatic amine dehydrogenase (AADH) reaction with tryptamine, recently studied by us, and identify possible causes for the differences observed in the measured kinetic isotope effects (KIEs) of the two systems. We have used combined quantum mechanics/molecular mechanics (QM/MM) techniques in molecular dynamics simulations and variational transition state theory with multidimensional tunnelling calculations averaged over an ensemble of paths. The results reveal important mechanistic complexity. We calculate activation barriers and KIEs for the two possible proton transfers identified-to either of the carboxylate oxygen atoms of the catalytic base (Asp428beta)-and analyse the contributions of quantum effects. The activation barriers and tunnelling contributions for the two possible proton transfers are similar and lead to a phenomenological activation free energy of 16.5+/-0.9 kcal mol(-1) for transfer to either oxygen (PM3-CHARMM calculations applying PM3-SRP specific reaction parameters), in good agreement with the experimental value of 14.4 kcal mol(-1). In contrast, for the AADH system, transfer to the equivalent OD1 was found to be preferred. The structures of the enzyme complexes during reaction are analysed in detail. The hydrogen bond of Thr474beta(MADH)/Thr172beta(AADH) to the catalytic carboxylate group and the nonconserved active site residue Tyr471beta(MADH)/Phe169beta(AADH) are identified as important factors in determining the preferred oxygen acceptor. The protein environment has a significant effect on the reaction energetics and hence on tunnelling contributions and KIEs. These environmental effects, and the related clearly different preferences for the two carboxylate oxygen atoms (with different KIEs) in MADH/MA and AADH/tryptamine, are possible causes of the differences observed in the KIEs between these two important enzyme reactions.

Combined quantum mechanics/molecular mechanics (QM/MM) modelling has the potential to answer fundamental questions about enzyme mechanisms and catalysis. Calculations using QM/MM methods can now predict barriers for enzyme-catalysed reactions with unprecedented, near chemical accuracy, i.e. to within 1 kcal/mol in the best cases. Quantitative predictions from first-principles calculations were only previously possible for very small molecules. At this level, quantitative, reliable predictions can be made about the mechanisms of enzyme-catalysed reactions. This development signals a new era of computational biochemistry.

We have investigated the structure and dynamics of the HIV-1 reverse transcriptase (HIV-RT) active site, by modelling the active conformation of the HIV-1 RT/DNA/deoxythymidine triphosphate (dTTP) ternary complex. This has included molecular dynamics simulations with the CHARMM27 force field, and combined quantum mechanics/molecular mechanics (QM/MM) calculations. Three different ternary systems were studied to investigate the effects of different protonation states of the dTTP substrate (a deprotonated and two different mono-protonated triphosphate forms of dTTP at the active site), and the effects of different possible protonation state of potentially catalytic aspartate residues (Asp185 and Asp186) were tested. Several potentially important hydrogen-bonding interactions with amino acids and bound water molecules in the deoxyribonucleoside triphosphate (dNTP) binding pocket were examined. The model of the deprotonated form of the dTTP substrate with the two aspartates in their charged (basic) form seemed to be the most stable and its orientation was in good agreement with X-ray crystallographic structure. In addition, two different semiempirical (AM1/CHARMM and PM3/CHARMM) QM/MM methods were tested for the HIV-RT system, in structural optimizations. Both methods provided conformations of the triphosphate moiety in either fully deprotonated or mono-protonated forms, which agreed well with the experimental structure of dTTP. The only significant difference between the AM1/CHARMM and PM3/CHARMM minimized structures is that the PM3/CHARMM Palpha-O3' optimized distance (important for nucleotide addition) is longer by 0.66 A in the deprotonated system but shorter by 0.37 A in the mono-protonated triphosphate system as compared with those obtained from AM1/CHARMM minimized structure. The obtained results suggest that both of these QM/MM methods, and the stochastic boundary molecular dynamics approach applied in this work, can give reasonable results for modelling the catalytically active complex of this important enzyme. The results provide insight into the structure and interactions of the active site of this important enzyme, with implications for its mechanism, which may be useful in inhibitor design.

Modeling methods allow the identification and analysis of determinants of reactivity and specificity in enzymes. The reaction between glutathione and 1-chloro-2,4-dinitrobenzene (CDNB) is widely used as a standard activity assay for glutathione S-transferases (GSTs). It is important to understand the causes of differences between catalytic GST isoenzymes and the effects of mutations and genetic polymorphisms. Quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulations have been performed here to investigate the addition of the glutathione anion to CDNB in the wild-type M1-1 GST isoenzyme from rat and in three single point mutant (Tyr6Phe, Tyr115Phe, and Met108Ala) M1-1 GST enzymes. We have developed a specifically parameterized QM/MM method (AM1-SRP/CHARMM22) to model this reaction by fitting to experimental heats of formation and ionization potentials. Free energy profiles were obtained from molecular dynamics simulations of the reaction using umbrella sampling and weighted histogram analysis techniques. The reaction in solution has also been simulated and is compared to the enzymatic reaction. The free energies are in excellent agreement with experimental results. Overall the results of the present study show that QM/MM reaction pathway analysis provides detailed insight into the chemistry of GST and can be used to obtain mechanistic insight into the effects of specific mutations on this catalytic process.

Proton tunneling dominates the oxidative deamination of tryptamine catalyzed by the enzyme aromatic amine dehydrogenase. For reaction with the fast substrate tryptamine, a H/D kinetic isotope effect (KIE) of 55 +/- 6 has been reported-one of the largest observed in an enzyme reaction. We present here a computational analysis of this proton-transfer reaction, applying combined quantum mechanics/molecular mechanics (QM/MM) methods (PM3-SRP//PM3/CHARMM22). In particular, we extend our previous computational study (Masgrau et al. Science 2006, 312, 237) by using improved energy corrections, high-level QM/MM methods, and an ensemble of paths to estimate the tunneling contributions. We have carried out QM/MM molecular dynamics simulations and variational transition state theory calculations with small-curvature tunneling corrections. The results provide detailed insight into the processes involved in the reaction. Transfer to the O2 oxygen of the catalytic base, Asp128beta, is found to be the favored reaction both thermodynamically and kinetically, even though O1 is closer in the reactant complex. Comparison of quantum and classical models of proton transfer allows estimation of the contribution of hydrogen tunneling in lowering the barrier to reaction in the enzyme. A reduction of the activation free energy due to tunneling of 3.1 kcal mol-1 is found, which represents a rate enhancement due to tunneling by 2 orders of magnitude. The calculated KIE of 30 is significantly elevated over the semiclassical limit, in agreement with the experimental observations; a semiclassical value of 6 is obtained when tunneling is omitted. A polarization of the C-H bond to be broken is observed due to the close proximity of the catalytic aspartate and the (formally) positively charged imine nitrogen. A comparison is also made with the related quinoprotein methylamine dehydrogenase (MADH)-the much lower KIE of 11 that we obtain for the MADH/methylamine system is found to arise from a more endothermic potential energy surface for the MADH reaction.

A recent paper in this journal uses molecular dynamics methods to study hydrolysis of guanosine triphosphate (GTP). The author reports that cleavage of the molecule occurs in less than 5 ps, and leads to a number of fragments including a free oxygen atom and a reduced magnesium ion. This conclusion is not in agreement with the known biochemistry and chemical reactivity of GTP or with previous computational studies of its hydrolysis reaction.

The fundamental nature of reactivity in cytochrome P450 enzymes is currently controversial. Modelling of bacterial P450cam has suggested an important role for the haem propionates in the catalysis, though this finding has been questioned. Understanding the mechanisms of this enzyme family is important both in terms of basic biochemistry and potentially in the prediction of drug metabolism. We have modelled the hydroxylation of camphor by P450cam, using combined quantum mechanics/molecular mechanics (QM/MM) methods. A set of reaction pathways in the enzyme was determined. We were able to pinpoint the source of the discrepancies in the previous results. We show that when a correct ionization state is assigned to Asp297, no spin density appears on the haem propionates and the protein structure in this region remains preserved. These results indicate that the haem propionates are not involved in catalysis.

QM/MM calculations provide a means for predicting the electronic structure of the metal center in metalloproteins. Two heme peroxidases, Cytochrome c Peroxidase (CcP) and Ascorbate Peroxidase (APX), have a structurally very similar active site, yet have active intermediates with very different electronic structures. We review our recent QM/MM calculations on these systems, and present new computational data. Our results are in good agreement with experiment, and suggest that the difference in electronic structure is due to a large number of small differences in structure from one protein to another. We also discuss recent QM/MM calculations on the active species of cytochrome P450, in which a similar sensitivity of the electronic structure to the environment is found. However, this does not appear to explain different catalytic profiles of the different drug-metabolizing isoforms of this class of enzyme.

It is now widely accepted that enzyme-catalysed C-H bond breakage occurs by quantum mechanical tunnelling. This paradigm shift in the conceptual framework for these reactions away from semi-classical transition state theory (TST, i.e. including zero-point energy, but with no tunnelling correction) has been driven over the recent years by experimental studies of the temperature dependence of kinetic isotope effects (KIEs) for these reactions in a range of enzymes, including the tryptophan tryptophylquinone-dependent enzymes such as methylamine dehydrogenase and aromatic amine dehydrogenase, and the flavoenzymes such as morphinone reductase and pentaerythritol tetranitrate reductase, which produced observations that are also inconsistent with the simple Bell-correction model of tunnelling. However, these data-especially, the strong temperature dependence of reaction rates and the variable temperature dependence of KIEs-are consistent with other tunnelling models (termed full tunnelling models), in which protein and/or substrate fluctuations generate a configuration compatible with tunnelling. These models accommodate substrate/protein (environment) fluctuations required to attain a configuration with degenerate nuclear quantum states and, when necessary, motion required to increase the probability of tunnelling in these states. Furthermore, tunnelling mechanisms in enzymes are supported by atomistic computational studies performed within the framework of modern TST, which incorporates quantum nuclear effects.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Citrate synthase forms citrate by deprotonation of acetyl-CoA followed by nucleophilic attack of this substrate on oxaloacetate, and subsequent hydrolysis. The rapid reaction rate is puzzling because of the instability of the postulated nucleophilic intermediate, the enolate of acetyl-CoA. As alternatives, the enol of acetyl-CoA, or an enolic intermediate sharing a proton with His-274 in a "low-barrier" hydrogen bond have been suggested. Similar problems of intermediate instability have been noted in other enzymic carbon acid deprotonation reactions. Quantum mechanical/molecular mechanical calculations of the pathway of acetyl-CoA enolization within citrate synthase support the identification of Asp-375 as the catalytic base. His-274, the proposed general acid, is found to be neutral. The acetyl-CoA enolate is more stable at the active site than the enol, and is stabilized by hydrogen bonds from His-274 and a water molecule. The conditions for formation of a low-barrier hydrogen bond do not appear to be met, and the calculated hydrogen bond stabilization in the reaction is less than the gas-phase energy, due to interactions with Asp-375 at the active site. The enolate character of the intermediate is apparently necessary for the condensation reaction to proceed efficiently.