Pathways and mechanisms for product release in the engineered haloalkane dehalogenases explored using classical and random acceleration molecular dynamics simulations.

Eight mutants of the DhaA haloalkane dehalogenase carrying mutations at the residues lining two tunnels, previously observed by protein X-ray crystallography, were constructed and biochemically characterized. The mutants showed distinct catalytic efficiencies with the halogenated substrate 1,2,3-trichloropropane. Release pathways for the two dehalogenation products, 2,3-dichloropropane-1-ol and the chloride ion, and exchange pathways for water molecules, were studied using classical and random acceleration molecular dynamics simulations. Five different pathways, denoted p1, p2a, p2b, p2c, and p3, were identified. The individual pathways showed differing selectivity for the products: the chloride ion releases solely through p1, whereas the alcohol releases through all five pathways. Water molecules play a crucial role for release of both products by breakage of their hydrogen-bonding interactions with the active-site residues and shielding the charged chloride ion during its passage through a hydrophobic tunnel. Exchange of the chloride ions, the alcohol product, and the waters between the buried active site and the bulk solvent can be realized by three different mechanisms: (i) passage through a permanent tunnel, (ii) passage through a transient tunnel, and (iii) migration through a protein matrix. We demonstrate that the accessibility of the pathways and the mechanisms of ligand exchange were modified by mutations. Insertion of bulky aromatic residues in the tunnel corresponding to pathway p1 leads to reduced accessibility to the ligands and a change in mechanism of opening from permanent to transient. We propose that engineering the accessibility of tunnels and the mechanisms of ligand exchange is a powerful strategy for modification of the functional properties of enzymes with buried active sites.

[1]  J. Koča,et al.  Repositioning the catalytic triad aspartic acid of haloalkane dehalogenase: effects on stability, kinetics, and structure. , 1997, Biochemistry.

[2]  Robert Zwanzig,et al.  Dynamical disorder: Passage through a fluctuating bottleneck , 1992 .

[3]  F. Raushel,et al.  An engineered blockage within the ammonia tunnel of carbamoyl phosphate synthetase prevents the use of glutamine as a substrate but not ammonia. , 2000, Biochemistry.

[4]  M. Fukuda,et al.  Isolation and characterization of Tn5-induced mutants of Pseudomonas paucimobilis UT26 defective in gamma-hexachlorocyclohexane dehydrochlorinase (LinA). , 1993, Bioscience, biotechnology, and biochemistry.

[5]  M. Maté,et al.  Role of the lateral channel in catalase HPII of Escherichia coli , 2008, Protein science : a publication of the Protein Society.

[6]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[7]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[8]  J. Hamilton,et al.  Haloalkane degradation and assimilation by Rhodococcus rhodochrous NCIMB 13064. , 1994, Microbiology.

[9]  W. Kabsch,et al.  Dictionary of protein secondary structure: Pattern recognition of hydrogen‐bonded and geometrical features , 1983, Biopolymers.

[10]  Christian Obinger,et al.  Role of the Main Access Channel of Catalase-Peroxidase in Catalysis* , 2005, Journal of Biological Chemistry.

[11]  J. Sambrook,et al.  Molecular Cloning: A Laboratory Manual , 2001 .

[12]  Jaroslav Koca,et al.  CAVER: a new tool to explore routes from protein clefts, pockets and cavities , 2006, BMC Bioinformatics.

[13]  Ploeg,et al.  Cloning of 1,2-dichloroethane degradation genes of Xanthobacter autotrophicus GJ10 and expression and sequencing of the dhlA gene , 1989, Journal of bacteriology.

[14]  Yuji Nagata,et al.  Haloalkane dehalogenase LinB from Sphingomonas paucimobilis UT26: X-ray crystallographic studies of dehalogenation of brominated substrates. , 2003, Biochemistry.

[15]  Dick B. Janssen,et al.  Specificity and Kinetics of Haloalkane Dehalogenase* , 1996, The Journal of Biological Chemistry.

[16]  Richard H. Henchman,et al.  The dynamics of ligand barrier crossing inside the acetylcholinesterase gorge. , 2003, Biophysical journal.

[17]  J. Damborský,et al.  Crystals of DhaA mutants from Rhodococcus rhodochrous NCIMB 13064 diffracted to ultrahigh resolution: crystallization and preliminary diffraction analysis. , 2008, Acta crystallographica. Section F, Structural biology and crystallization communications.

[18]  J. Newman,et al.  Crystal structure of the haloalkane dehalogenase from Sphingomonas paucimobilis UT26. , 2000, Biochemistry.

[19]  P. Loewen,et al.  An Electrical Potential in the Access Channel of Catalases Enhances Catalysis* , 2003, Journal of Biological Chemistry.

[20]  N. Sreerama,et al.  A self-consistent method for the analysis of protein secondary structure from circular dichroism. , 1993, Analytical biochemistry.

[21]  B. Brooks,et al.  Alternate Pathways for Acetic Acid and Acetate ion Release from Acetylcholinesterase: a Molecular Dynamics Study , 1998 .

[22]  T C Terwilliger,et al.  Dialysis against a Solution Containing 50 Mm Na 2 so 4 and 1 Mm Edta, with the Ph Adjusted to 8.2, Using 1 N H 2 so 4 , 1999 .

[23]  Yuji Nagata,et al.  Halide-stabilizing residues of haloalkane dehalogenases studied by quantum mechanic calculations and site-directed mutagenesis. , 2002, Biochemistry.

[24]  S. Sligar,et al.  Electrostatic control of the substrate access channel in cytochrome P-450cam. , 1994, Biochemistry.

[25]  Jürgen Pleiss,et al.  Multiple molecular dynamics simulations of human p450 monooxygenase CYP2C9: The molecular basis of substrate binding and regioselectivity toward warfarin , 2006, Proteins.

[26]  W. Delano The PyMOL Molecular Graphics System , 2002 .

[27]  J. Damborský,et al.  Weak Activity of Haloalkane Dehalogenase LinB with 1,2,3-Trichloropropane Revealed by X-Ray Crystallography and Microcalorimetry , 2007, Applied and Environmental Microbiology.

[28]  M. A. Andrade,et al.  Evaluation of secondary structure of proteins from UV circular dichroism spectra using an unsupervised learning neural network. , 1993, Protein engineering.

[29]  R. Wade,et al.  Exceptionally stable salt bridges in cytochrome P450cam have functional roles. , 1997, Biochemistry.

[30]  A. Costache,et al.  Structural evidence for a functionally relevant second camphor binding site in P450cam: Model for substrate entry into a P450 active site , 2007, Proteins.

[31]  J. Halpert,et al.  Substrate routes to the buried active site may vary among cytochromes P450: mutagenesis of the F-G region in P450 2B1. , 2002, Chemical research in toxicology.

[32]  Hans-Joachim Wittmann,et al.  LigPath: a module for predictive calculation of a ligand’s pathway into a receptor-application to the gpH1 - receptor , 2006, Journal of molecular modeling.

[33]  J. F. Schindler,et al.  Haloalkane dehalogenases: steady-state kinetics and halide inhibition. , 1999, Biochemistry.

[34]  C. Kennes,et al.  REPLACEMENT OF TRYPTOPHAN RESIDUES IN HALOALKANE DEHALOGENASE REDUCES HALIDE BINDING AND CATALYTIC ACTIVITY , 1995 .

[35]  J. Bujnicki,et al.  Phylogenetic analysis of haloalkane dehalogenases , 2007, Proteins.

[36]  Yuji Nagata,et al.  Catalytic Mechanism of the Haloalkane Dehalogenase LinB from Sphingomonas paucimobilis UT26* , 2003, Journal of Biological Chemistry.

[37]  J. Koča,et al.  Analysis of the reaction mechanism and substrate specificity of haloalkane dehalogenases by sequential and structural comparisons. , 1999, Protein engineering.

[38]  R. Wade,et al.  Comparison of the dynamics of substrate access channels in three cytochrome P450s reveals different opening mechanisms and a novel functional role for a buried arginine , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[39]  C. Di Primo,et al.  Origin of the photoacoustic signal in cytochrome P-450cam: role of the Arg186-Asp251-Lys178 bifurcated salt bridge. , 1997, Biochemistry.

[40]  Sudarko,et al.  A survey of active site access channels in cytochromes P450. , 2004, Journal of inorganic biochemistry.

[41]  Yuji Nagata,et al.  Modification of Activity and Specificity of Haloalkane Dehalogenase from Sphingomonas paucimobilis UT26 by Engineering of Its Entrance Tunnel* , 2003, Journal of Biological Chemistry.

[42]  E. Sedlák,et al.  Cofactor assisted gating mechanism in the active site of NADH oxidase from Thermus thermophilus , 2006, Proteins.

[43]  H. Zhou,et al.  On the human CYP2C9*13 variant activity reduction: a molecular dynamics simulation and docking study. , 2006, Biochimie.

[44]  S. Prasad,et al.  Binding of camphor to Pseudomonas putida cytochrome P450cam: steady‐state and picosecond time‐resolved fluorescence studies , 2000, FEBS letters.

[45]  K. H. Kalk,et al.  Crystallographic analysis of the catalytic mechanism of haloalkane dehalogenase , 1994, Nature.

[46]  A. Dunn,et al.  Probing the open state of cytochrome P450cam with ruthenium-linker substrates , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[47]  Jaroslav Koca,et al.  A Molecular Modeling Study of the Catalytic Mechanism of Haloalkane Dehalogenase: 1. Quantum Chemical Study of the First Reaction Step , 1997, J. Chem. Inf. Comput. Sci..

[48]  M. Zámocký,et al.  Site‐directed mutagenesis of the lower parts of the major substrate channel of yeast catalase A leads to highly increased peroxidatic activity , 1995, FEBS letters.

[49]  P. Kollman,et al.  A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules , 1995 .

[50]  P. Kollman,et al.  A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules J. Am. Chem. Soc. 1995, 117, 5179−5197 , 1996 .

[51]  Yuji Nagata,et al.  Redesigning dehalogenase access tunnels as a strategy for degrading an anthropogenic substrate. , 2009, Nature chemical biology.

[52]  P Argos,et al.  Accessibility to internal cavities and ligand binding sites monitored by protein crystallographic thermal factors , 1998, Proteins.

[53]  Yong Duan,et al.  Chromophore channeling in the G-protein coupled receptor rhodopsin. , 2007, Journal of the American Chemical Society.

[54]  Jerome Baudry,et al.  Ile115Leu mutation in the SRS1 region of an insect cytochrome P450 (CYP6B1) compromises substrate turnover via changes in a predicted product release channel. , 2005, Protein engineering, design & selection : PEDS.

[55]  G Vriend,et al.  WHAT IF: a molecular modeling and drug design program. , 1990, Journal of molecular graphics.

[56]  Lennart Nilsson,et al.  Unbinding of retinoic acid from the retinoic acid receptor by random expulsion molecular dynamics. , 2006, Biophysical journal.

[57]  P. Champion,et al.  Dynamics of nitric oxide rebinding and escape in horseradish peroxidase. , 2006, Journal of the American Chemical Society.

[58]  K. Schulten,et al.  Imaging the migration pathways for O2, CO, NO, and Xe inside myoglobin. , 2006, Biophysical journal.

[59]  Kaihsu Tai,et al.  Properties of water molecules in the active site gorge of acetylcholinesterase from computer simulation. , 2002, Biophysical journal.

[60]  T. Poulos,et al.  Crystal structure of cytochrome P450 14α-sterol demethylase (CYP51) from Mycobacterium tuberculosis in complex with azole inhibitors , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[61]  J. Kingma,et al.  Steady-state and pre-steady-state kinetic analysis of halopropane conversion by a Rhodococcus haloalkane dehalogenase , 2022 .

[62]  F. Gago,et al.  Stepwise dissection and visualization of the catalytic mechanism of haloalkane dehalogenase LinB using molecular dynamics simulations and computer graphics. , 2007, Journal of molecular graphics & modelling.

[63]  Yuji Nagata,et al.  Crystal structure of haloalkane dehalogenase LinB from Sphingomonas paucimobilis UT26 at 0.95 A resolution: dynamics of catalytic residues. , 2004, Biochemistry.

[64]  James J. P. Stewart,et al.  MOPAC: A semiempirical molecular orbital program , 1990, J. Comput. Aided Mol. Des..

[65]  Joel L. Sussman,et al.  The α/β hydrolase fold , 1992 .

[66]  R C Wade,et al.  How do substrates enter and products exit the buried active site of cytochrome P450cam? 2. Steered molecular dynamics and adiabatic mapping of substrate pathways. , 2000, Journal of molecular biology.

[67]  Iwaji Iwasaki,et al.  New Colorimetric Determination of Chloride using Mercuric Thiocyanate and Ferric Ion , 1952 .

[68]  K. Sharp,et al.  Accessibility of oxygen with respect to the heme pocket in horseradish peroxidase , 2003, Proteins.

[69]  L. Mouawad,et al.  Ligand migration and escape pathways in haem proteins. , 2006, Biochemical Society transactions.

[70]  Weihua Li,et al.  POSSIBLE PATHWAY(S) OF TESTOSTERONE EGRESS FROM THE ACTIVE SITE OF CYTOCHROME P450 2B1: A STEERED MOLECULAR DYNAMICS SIMULATION , 2005, Drug Metabolism and Disposition.

[71]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[72]  Michal Otyepka,et al.  Mechanism of enhanced conversion of 1,2,3-trichloropropane by mutant haloalkane dehalogenase revealed by molecular modeling , 2006, J. Comput. Aided Mol. Des..

[73]  I. S. Ridder,et al.  Crystallographic and kinetic evidence of a collision complex formed during halide import in haloalkane dehalogenase. , 1999, Biochemistry.

[74]  M. Silberstein,et al.  Exploring the binding sites of the haloalkane dehalogenase DhlA from Xanthobacter autotrophicus GJ10. , 2007, Biochemistry.

[75]  Dick B. Janssen,et al.  Biodegradation of 1,2,3-Trichloropropane through Directed Evolution and Heterologous Expression of a Haloalkane Dehalogenase Gene , 2002, Applied and Environmental Microbiology.

[76]  Dick B Janssen,et al.  Evolving haloalkane dehalogenases. , 2004, Current opinion in chemical biology.

[77]  K. Sharp,et al.  Water channel of horseradish peroxidase studied by the charge-transfer absorption band of ferric heme , 2004 .

[78]  Pavel Kyslík,et al.  Cloning of an epoxide hydrolase-encoding gene from Aspergillus niger M200, overexpression in E. coli, and modification of activity and enantioselectivity of the enzyme by protein engineering. , 2007, Journal of biotechnology.

[79]  K Schulten,et al.  VMD: visual molecular dynamics. , 1996, Journal of molecular graphics.

[80]  J A McCammon,et al.  Molecular dynamics of mouse acetylcholinesterase complexed with huperzine A. , 1999, Biopolymers.

[81]  R. Wade,et al.  How do substrates enter and products exit the buried active site of cytochrome P450cam? 1. Random expulsion molecular dynamics investigation of ligand access channels and mechanisms. , 2000, Journal of molecular biology.

[82]  Christophe Geourjon,et al.  An interactive graphic program for calculating the secondary structure content of proteins from circular dichroism spectrum , 1993, Comput. Appl. Biosci..

[83]  D. Case,et al.  Exploring protein native states and large‐scale conformational changes with a modified generalized born model , 2004, Proteins.

[84]  J. Sussman,et al.  An electrostatic mechanism for substrate guidance down the aromatic gorge of acetylcholinesterase. , 1993, Proceedings of the National Academy of Sciences of the United States of America.

[85]  Jürgen Pleiss,et al.  Blocking the tunnel: engineering of Candida rugosa lipase mutants with short chain length specificity. , 2002, Protein engineering.

[86]  M. Teeter,et al.  Myoglobin cavities provide interior ligand pathway , 2004, Protein science : a publication of the Protein Society.

[87]  I M Kovach,et al.  Efficient product clearance through exit channels in substrate hydrolysis by acetylcholinesterase , 1994, FEBS letters.

[88]  Joel L Sussman,et al.  How does huperzine A enter and leave the binding gorge of acetylcholinesterase? Steered molecular dynamics simulations. , 2003, Journal of the American Chemical Society.

[89]  Tudor I. Oprea,et al.  Identification of a functional water channel in cytochrome P450 enzymes. , 1997, Proceedings of the National Academy of Sciences of the United States of America.

[90]  Rebecca C Wade,et al.  The ins and outs of cytochrome P450s. , 2007, Biochimica et biophysica acta.

[91]  David S. Goodsell,et al.  Automated docking using a Lamarckian genetic algorithm and an empirical binding free energy function , 1998, J. Comput. Chem..

[92]  Michal Otyepka,et al.  Functionally relevant motions of haloalkane dehalogenases occur in the specificity‐modulating cap domains , 2002, Protein science : a publication of the Protein Society.

[93]  S J Wodak,et al.  Pathways of ligand clearance in acetylcholinesterase by multiple copy sampling. , 2000, Journal of molecular biology.

[94]  William E Bentley,et al.  Protein engineering of toluene 4‐monooxygenase of Pseudomonas mendocina KR1 for synthesizing 4‐nitrocatechol from nitrobenzene , 2004, Biotechnology and bioengineering.