Modelling Biological Molecules and Simulating Enzyme-Catalysed Reactions

Dr. Adrian Mulholland (Adrian.Mulholland@bris.ac.uk)
 

A postdoctoral research assistant position is available, to start on June 1st 2005 or later. This project will develop accurate QM/MM methods for modelling enzyme catalysis. The position is funded by the EPSRC. The aim of the project is to push forward the state of the art in describing the fundamental mechanisms of biological catalysis. Hybrid methods combining a quantum mechanical (QM) treatment of a small number of atoms within the active site with molecular mechanics (MM) for the surrounding protein have already been shown to yield valuable insight into enzyme reaction mechanisms. The PDRA will develop QM/MM methods using cutting edge quantum chemistry for the QM region and will investigate techniques for exploring reaction paths and for computing activation free energies for enzymes. It is a collaborative project, involving Dr. Fred Manby and Dr. Jeremy Harvey. Please contact one of us for more details.

 

PhD projects for 2005 - these projects will pay an enhanced student stipend of £2000 extra per year above standard rates throughout the project

Two BBSRC Strategic Research Studentships are available to begin in October 2005. The projects are in the areas of:

Examining fundamental principles of enzyme catalysis by computational enzymology

Relating central nervous system enzyme sequence and function through high-performance computing and e-science

These projects will involve the use of modelling methods to understand fundamental aspects of enzymes relevant to biocatalysis, medicine and human health (for example, enzymes involved in the central nervous system). UK citizens are eligible for these studentships. Please contact me for more details. Other projects include studying quantum mechanical effects, such as tunnelling, in enzyme-catalysed proton and hydride transfer reactions.

Other PhD projects to start in 2005: applications from EU citizens are welcome. Collaborative projects with experimental groups are available (see e.g. research into protein structure in the group of Dr. A. Hadfield, Biochemistry) . Start dates can be flexible.

Meet the group...

How a biological reaction happens. This picture shows the transition state in the enzyme para-hydroxybenzoate hydroxylase, as found by QM/MM modelling (refs. 5, 9 and 15 below). An OH group (centre) is transferred from a flavin cofactor (right) to the substrate (left). The transition state is stabilized by a hydrogen bond interaction with a key group in the enzyme (shown as a dotted line). This catalytic effect helps to make the biodegradation of aromatic compounds efficient. Para-hydroxybenzoate hydroxylase (PHBH) is an important enzyme in the microbial biodegradation of a wide variety of aromatic chemicals, including pollutants and lignin, a major component of wood and so among the most abundant of all biopolymers. PHBH acts by oxygenating its substrates (adding a hydroxyl (OH) group to the aromatic ring), and is a key member of an important class of oxygenase enzymes.

Here is a poster (pdf-quite large) describing some of our recent work, presented at the Protein Society meeting in Boston 2003  Work from the group has also recently been highlighted in the newsletter of the CSAR national supercomputing centre, CSAR Focus (Summer 2003, 10, 12-13)  and in the magazine of the BBSRC, BBSRC Business ('Biology Joins the High Performance Computing Club') ; Please contact me for more details.

        Most biochemical reactions depend on enzyme catalysis, and understanding how enzymes 'work' at the molecular level is a fundamental problem.  Enzymes are remarkably efficient and specific catalysts, but despite intensive experimental investigations, the detailed origins of their rate accelerations remain unclear.  This question is of crucial importance in biology, and also for the development of protein catalysts for practical applications.  Better understanding is vital for analysing the activities of mutant or designed proteins, and for the design of inhibitors as pharmaceutical lead compounds.
   

Citrate synthase (refs.1,4&8)
 

        Computer simulations are a good way to study enzyme reactions.  They can provide information which is often inaccessible experimentally, such as details of unstable species (for example transition states and reaction intermediates) and on energetic contributions to catalysis.  Simulating an enzyme reaction is a challenging problem, and requires the use of specialized techniques.  An effective approach is to use combined quantum mechanical/molecular mechanical (QM/MM) methods.  Molecular mechanics methods can be used to study protein conformational changes, dynamics and binding, but generally they can't be applied to processes involving the breaking or making of chemical bonds.  For chemical reactions a quantum mechanical description is needed, which can be achieved by the combined QM/MM approach.  The small QM region contains the groups involved in the reaction (e.g. the catalytic residues and the substrate(s)) and is influenced by the surroundings (represented more simply by standard molecular mechanics), so including the effects of the environment.  In this way the reaction in an enzyme can be studied, and contributions of individual groups can be analysed.

        The central theme of my work is the use of simulation methods to investigate enzyme mechanism and dynamics (in several cases in active collaboration with experimental studies).  Programming expertise is not required for most of the present work. Projects offered include:

Some references (please contact me if you would like reprints):
  1. C.M. Bathelt, L. Ridder, A.J. Mulholland & J.N. Harvey 'Aromatic hydroxylation by cytochrome P450: Model calculations of mechanism and substituent effects' J. Am. Chem. Soc., 125, 15004 -15005 (2003).
  2. K.E. Ranaghan, L. Ridder, B. Szefczyk, W.A. Sokalski, J.C. Hermann, and A.J. Mulholland 'Transition state stabilization and substrate strain in enzyme catalysis: ab initio QM/MM modelling of the chorismate mutase reaction' Organic and Biomolecular Chemistry (2004) 2, 968-980
  3. J. Zurek, A.L. Bowman, W.A. Sokalski & A.J.Mulholland 'MM and QM/MM Modeling of Threonyl-tRNA Synthetase: Model Preparation, Testing and Simulations', Structural Chemistry, 15, 405-414 (2004).
  4. K.E. Ranaghan & A.J. Mulholland 'Conformational effects in enzyme catalysis: QM/MM free energy calculation of the 'NAC' contribution in chorismate mutase' Chem. Commun. (2004) (10), 1238-1239
  5. C.M. Bathelt, L. Ridder, A.J. Mulholland & J.N. Harvey 'Mechanism and Structure-Reactivity Relationships for Aromatic Hydroxylation by Cytochrome P450' Organic and Biomolecular Chemistry, (2004) 2, 2998-3005.
  6. J.C. Hermann, L. Ridder, A.J. Mulholland, and H.-D. Hoeltje, Identification of Glu166 as the General Base in the Acylation Reaction of Class A beta-Lactamses through QM/MM Modeling  J. Am. Chem. Soc. 125, 9590-9591 (2003)
  7. L. Ridder and A.J. Mulholland, Modeling biotransformation reactions by combined quantum mechanical/molecular mechanical approaches: from structure to activity  Curr. Topics in Medicinal Chemistry. 3, 1241-1256 (2003)
  8. K. Ranaghan, L. Ridder, B. Szefczyk, W. A. Sokalski and A.J. Mulholland, Insights into enzyme catalysis from QM/MM modelling: transition state stabilization in chorismate mutase  Mol. Phys. 101, 2695-2714 (2003)
  9. L. Ridder, J.N. Harvey, I.M.C.M. Rietjens, J. Vervoort and A.J. Mulholland, Ab initio QM/MM Modeling of the Hydroxylation Step in p-Hydroxybenzoate Hydroxylase   J. Phys. Chem. B 107, 2118-2126 (2003)
  10. L. Ridder, I.M.C.M. Rietjens, J. Vervoort and A.J. Mulholland, Quantum Mechanical/Molecular Mechanical Free Energy Simulations of the Glutathione S-Transferase (M1-1) Reaction with Phenanthrene 9,10-Oxide  J. Am. Chem. Soc. 124, 9926-9936 (2002)
  11. A.J. Mulholland, P.D. Lyne and M. Karplus, Ab Initio QM/MM Study of the Citrate Synthase Mechanism. A Low-Barrier Hydrogen Bond is not Involved  J. Am. Chem. Soc. 122, 534-535 (2000)
  12. A.T. Hadfield and A.J. Mulholland, Active Site Dynamics of ASADH- A Bacterial Biosynthetic Enzyme  Int. J. Quant. Chem., Biophys. Q. 73, 137-146 (1999)
  13. L. Ridder, A.J. Mulholland, J. Vervoort and I.M.C.M. Rietjens, Correlation of Calculated Activation Energies with Experimental Rate Constants for an Enzyme-Catalyzed Aromatic Hydroxylation  J. Am. Chem. Soc. 120, 7641-7642 (1998)
  14. A.J. Mulholland and W.G. Richards, Modeling Enzyme Reaction Intermediates and Transition States: Citrate Synthase   J.Phys.Chem. B 102, 6635-6646 (1998)
  15. A.J. Mulholland, G.H. Grant and W.G. Richards, Review: Computer modelling of enzyme catalysed reaction mechanisms.  Protein Engineering 6, 133-147 (1993)
  16. P.D. Lyne, A.J. Mulholland and W.G. Richards, Insights into Chorismate Mutase Catalysis from a Combined QM/MM Simulation of the Enzyme ReactionJ. Am. Chem. Soc. 117, 11345-11350 (1995)
  17. A.J. Mulholland and M. Karplus, Simulations of enzymic reactions, Biochem. Soc. Trans. 24, 247-254 (1996)
  18. A.J. Mulholland and W.G. Richards, Acetyl-CoA Enolization in Citrate Synthase: A Quantum Mechanical/Molecular Mechanical (QM/MM) Study, Proteins: Structure, Function, and Genetics 27, 9-25 (1997)
  19. L.Ridder, A.J. Mulholland, I.M.C.M. Rietjens and J. Vervoort, A Quantum Mechanical/Molecular Mechanical Study of the Hydroxylation of Phenol and Halogenated Derivatives by Phenol Hydroxylase, J. Am. Chem. Soc. 122, 8728-8738 (2000)
  20.  


Back to School of Chemistry Home Page