Chemical accuracy in QM/MM calculations on enzyme-catalysed 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.     

Mechanisms of reaction in cytochrome P450: Hydroxylation of camphor in P450cam

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.     

Calculation of chemical and phase equilibria via simulated annealing

We present a new class of techniques for the solution of the chemical and phase equilibria problem for reacting species in a closed system. The minimisation of the Gibbs free energy for all the species in the system is conducted using the technique of simulated annealing (SA). The SA objective function incorporates non‐ideal equations of state. This new approach is demonstrably able to solve multi‐species and multi‐phase LTCE problems in ideal‐gas solutions, ideal solutions and mixtures of ideal and non‐ideal solutions. This revised version was published online in July 2006 with corrections to the Cover Date.