A survey of active site access channels in cytochromes P450

In cytochrome P450s, the active site is situated deep inside the protein next to the heme cofactor, and is often completely isolated from the surrounding solvent. To identify routes by which substrates may enter into and products exit from the active site, random expulsion molecular dynamics simulations were performed for three cytochrome P450s: CYP101, CYP102A1 and CYP107A1 [J. Mol. Biol. 303 (2000) 797; Proc. Natl. Acad. Sci. USA 99 (2002) 5361]. Amongst the different pathways identified, one pathway was found to be common to all three cytochrome P450s although the mechanism of ligand passage along it was different in each case and apparently adapted to the substrate specificity of the enzyme. Recently, a number of new crystal structures of cytochrome P450s have been solved. Here, we analyse the open channels leading to the active site that these structures reveal. We find that in addition to showing the common pathway, they provide experimental evidence for the existence of three additional channels that were identified by simulation. We also discuss how the location of xenon binding sites in CYP101 suggests a role for one of the pathways identified by molecular dynamics simulations as a route for gaseous species, such as oxygen, to access the active site.     

Crystallographic insights into a cobalt (III) sepulchrate based alternative cofactor system of P450 BM3 monooxygenase

P450 BM3 is a multi-domain heme-containing soluble bacterial monooxygenase. P450 BM3 and variants are known to oxidize structurally diverse substrates. Crystal structures of individual domains of P450 BM3 are available. However, the spatial organization of the full-length protein is unknown. In this study, crystal structures of the P450 BM3 M7 heme domain variant with and without cobalt (III) sepulchrate are reported. Cobalt (III) sepulchrate acts as an electron shuttle in an alternative cofactor system employing zinc dust as the electron source. The crystal structure shows a binding site for the mediator cobalt (III) sepulchrate at the entrance of the substrate access channel. The mediator occupies an unusual position which is far from the active site and distinct from the binding of the natural redox partner (FAD/NADPH binding domain).     

Active-site hydration and water diffusion in cytochrome P450cam: a highly dynamic process

Long-timescale molecular dynamics simulations (300 ns) are performed on both the apo- (i.e., camphor-free) and camphor-bound cytochrome P450cam (CYP101). Water diffusion into and out of the protein active site is observed without biased sampling methods. During the course of the molecular dynamics simulation, an average of 6.4 water molecules is observed in the camphor-binding site of the apo form, compared to zero water molecules in the binding site of the substrate-bound form, in agreement with the number of water molecules observed in crystal structures of the same species. However, as many as 12 water molecules can be present at a given time in the camphor-binding region of the active site in the case of apo-P450cam, revealing a highly dynamic process for hydration of the protein active site, with water molecules exchanging rapidly with the bulk solvent. Water molecules are also found to exchange locations frequently inside the active site, preferentially clustering in regions surrounding the water molecules observed in the crystal structure. Potential-of-mean-force calculations identify thermodynamically favored trans-protein pathways for the diffusion of water molecules between the protein active site and the bulk solvent. Binding of camphor in the active site modifies the free-energy landscape of P450cam channels toward favoring the diffusion of water molecules out of the protein active site.     

The dynamics of camphor in the cytochrome P450 CYP101D2

The recent crystal structures of CYP101D2, a cytochrome P450 protein from the oligotrophic bacterium Novosphingobium aromaticivorans DSM12444 revealed that both the native (substrate‐free) and camphor‐soaked forms have open conformations. Furthermore, two other potential camphor‐binding sites were also identified from electron densities in the camphor‐soaked structure, one being located in the access channel and the other in a cavity on the surface near the F‐helix side of the F‐G loop termed the substrate recognition site. These latter sites may be key intermediate positions on the pathway for substrate access to or product egress from the active site. Here, we show via the use of unbiased atomistic molecular dynamics simulations that despite the open conformation of the native and camphor‐bound crystal structures, the underlying dynamics of CYP101D2 appear to be very similar to other CYP proteins. Simulations of the native structure demonstrated that the protein is capable of sampling many different conformational substates. At the same time, simulations with the camphor positioned at various locations within the access channel or recognition site show that movement towards the active site or towards bulk solvent can readily occur on a short timescale, thus confirming many previously reported in silico studies using steered molecular dynamics. The simulations also demonstrate how the fluctuations of an aromatic gate appear to control access to the active site. Finally, comparison of camphor‐bound simulations with the native simulations suggests that the fluctuations can be of similar level and thus are more representative of the conformational selection model rather than induced fit.     

Conformational states of cytochrome P450cam revealed by trapping of synthetic molecular wires

Members of the ubiquitous cytochrome P450 family catalyze a vast range of biologically significant reactions in mammals, plants, fungi, and bacteria. Some P450s display a remarkable promiscuity in substrate recognition, while others are very specific with respect to substrate binding or regio and stereo-selective catalysis. Recent results have suggested that conformational flexibility in the substrate access channel of many P450s may play an important role in controlling these effects. Here, we report the X-ray crystal structures at 1.8A and 1.5A of cytochrome P450cam complexed with two synthetic molecular wires, D-4-Ad and D-8-Ad, consisting of a dansyl fluorophore linked to an adamantyl substrate analog via an alpha,omega-diaminoalkane chain of varying length. Both wires bind with the adamantyl moiety in similar positions at the camphor-binding site. However, each wire induces a distinct conformational response in the protein that differs from the camphor-bound structure. The changes involve significant movements of the F, G, and I helices, allowing the substrate access channel to adapt to the variable length of the probe. Wire-induced opening of the substrate channel also alters the I helix bulge and Thr252 at the active site with binding of water that has been proposed to assist in peroxy bond cleavage. The structures suggest that the coupling of substrate-induced conformational changes to active-site residues may be different in P450cam and recently described mammalian P450 structures. The wire-induced changes may be representative of the conformational intermediates that must exist transiently during substrate entry and product egress, providing a view of how substrates enter the deeply buried active site. They also support observed examples of conformational plasticity that are believed be responsible for the promiscuity of drug metabolizing P450s. Observation of such large changes in P450cam suggests that substrate channel plasticity is a general property inherent to all P450 structures.     

Structural evidence for a functionally relevant second camphor binding site in P450cam: model for substrate entry into a P450 active site

P450cam has long served as a prototype for the cytochrome P450 (CYP) gene family. But, little is known about how substrate enters its active site pocket, and how access is achieved in a way that minimizes exposure of the reactive heme. We hypothesize that P450cam may first bind substrate transiently near the mobile F-G helix that covers the active site pocket. Such a two-step binding process is kinetically required if P450cam rarely populates an open conformation-as suggested by previous literature and the inability to obtain a crystal structure of P450cam in an open conformation. Such a mechanism would minimize exposure of the heme by allowing P450cam to stay in a closed conformation as long as possible, since only brief flexing into an open conformation would be required to allow substrate entry. To test this model, we have attempted to dock a second camphor molecule into the crystal structure of camphor-bound P450cam. The docking identified only one potential entry site pocket, a well-defined cavity on the F-helix side of the F-G flap, 16 A from the heme iron. Location of this entry site pocket is consistent with our NMR T1 relaxation-based measurements of distances for a camphor that binds in fast exchange (active site camphor is known to bind in slow exchange). Presence of a second camphor binding site is also confirmed with [(1)H-(13)C] HSQC titrations of (13)CH3-threonine labeled P450cam. To confirm that camphor can bind outside of the active site pocket, (13)CH3-S-pyridine was bound to the heme iron to physically block the active site, and to serve as an NMR chemical shift probe. Titration of this P450cam-pyridine complex confirms that camphor can bind to a site outside the active site pocket, with an estimated Kd of 43 microM. The two-site binding model that is proposed based on these data is analogous to that recently proposed for CYP3A4, and is consistent with recent crystal structures of P450cam bound to tethered-substrates, which force a partially opened conformation.     

Substrate Access to Cytochrome P450cam: A Comparison of a Thermal Motion Pathway Analysis with Molecular Dynamics Simulation Data

Cytochrome P450 enzymes are hemoprotein monooxygenases that catalyse the oxidation of a variety of compounds. The mechanism by which camphor, the natural substrate of Cytochrome P450cam (P450cam), accesses the active site is a long-standing puzzle, although putative access channels have been proposed. A thermal motion pathway (TMP) analysis was performed on the crystal structure of P450cam with camphor bound. Hereby, three distinct thermal motion pathway families (TMPFs) were found. Possible substrate access channels obtained by this analysis based on B-factors are compared with exit channels explored by molecular dynamics simulations (MDS) by imposing an artificial expulsion force on the substrate in addition to the standard MD force field. Two out of three TMPFs are supported by results obtained with the random expulsion MDS method. However, the pathway found by the TMP method to have the highest average B-factor could not be observed by MDS. The pathway proposed from crystallographic data, which is a small opening above the active site located near residues 185, 87 and 395 corresponds to the TMPF with the second highest average B-factor.     

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

Cytochrome P450s form a ubiquitous protein family with functions including the synthesis and degradation of many physiologically important compounds and the degradation of xenobiotics. Cytochrome P450cam from Pseudomonas putida has provided a paradigm for the structural understanding of cytochrome P450s. However, the mechanism by which camphor, the natural substrate of cytochrome P450cam, accesses the buried active site is a long-standing puzzle. While there is recent crystallographic and simulation evidence for opening of a substrate-access channel in cytochrome P450BM-3, for cytochrome P450cam, no such conformational changes have been observed either in different crystal structures or by standard molecular dynamics simulations. Here, a novel simulation method, random expulsion molecular dynamics, is presented, in which substrate-exit channels from the buried active site are found by imposing an artificial randomly oriented force on the substrate, in addition to the standard molecular dynamics force field. The random expulsion molecular dynamics method was tested in simulations of the substrate-bound structure of cytochrome P450BM-3, and then applied to complexes of cytochrome P450cam with different substrates and with product. Three pathways were identified, one of which corresponds to a channel proposed earlier on the basis of crystallographic and site-directed mutagenesis data. Exit via the water-filled channel, which was previously suggested to be a product exit channel, was not observed. The pathways obtained by the random expulsion molecular dynamics method match well with thermal motion pathways obtained by an analysis of crystallographic B-factors. In contrast to large backbone motions (up to 4 A) observed in cytochrome P450BM-3 for the exit of palmitoleic acid, passage of camphor through cytochrome P450cam only requires small backbone motions (less than 2.4 A) in conjunction with side-chain rotations. Concomitantly, in almost all the exit trajectories, salt-links that have been proposed to act as ionic tethers between secondary structure elements of the protein, are perturbed.     

Structure of microsomal cytochrome P450 2B4 complexed with the antifungal drug bifonazole: insight into P450 conformational plasticity and membrane interaction

To better understand ligand-induced structural transitions in cytochrome P450 2B4, protein-ligand interactions were investigated using a bulky inhibitor. Bifonazole, a broad spectrum antifungal agent, inhibits monooxygenase activity and induces a type II binding spectrum in 2B4dH(H226Y), a modified enzyme previously crystallized in the presence of 4-(4-chlorophenyl)imidazole (CPI). Isothermal titration calorimetry and tryptophan fluorescence quenching indicate no significant burial of protein apolar surface nor altered accessibility of Trp-121 upon bifonazole binding, in contrast to recent results with CPI. A 2.3 A crystal structure of 2B4-bifonazole reveals a novel open conformation with ligand bound in the active site, which is significantly different from either the U-shaped cleft of ligand-free 2B4 or the small active site pocket of 2B4-CPI. The O-shaped active site cleft of 2B4-bifonazole is widely open in the middle but narrow at the top. A bifonazole molecule occupies the bottom of the active site cleft, where helix I is bent approximately 15 degrees to accommodate the bulky ligand. The structure also defines unanticipated interactions between helix C residues and bifonazole, suggesting an important role of helix C in azole recognition by mammalian P450s. Comparison of the ligand-free 2B4 structure, the 2B4-CPI structure, and the 2B4-bifonazole structure identifies structurally plastic regions that undergo correlated conformational changes in response to ligand binding. The most plastic regions are putative membrane-binding motifs involved in substrate access or substrate binding. The results allow us to model the membrane-associated state of P450 and provide insight into how lipophilic substrates access the buried active site.     

Crystal structure of inhibitor-bound P450BM-3 reveals open conformation of substrate access channel

P450BM-3 is an extensively studied P450 cytochrome that is naturally fused to a cytochrome P450 reductase domain. Crystal structures of the heme domain of this enzyme have previously generated many insights into features of P450 structure, substrate binding specificity, and conformational changes that occur on substrate binding. Although many P450s are inhibited by imidazole, this compound does not effectively inhibit P450BM-3. Omega-imidazolyl fatty acids have previously been found to be weak inhibitors of the enzyme and show some unusual cooperativity with the substrate lauric acid. We set out to improve the properties of these inhibitors by attaching the omega-imidazolyl fatty acid to the nitrogen of an amino acid group, a tactic that we used previously to increase the potency of substrates. The resulting inhibitors were significantly more potent than their parent compounds lacking the amino acid group. A crystal structure of one of the new inhibitors bound to the heme domain of P450BM-3 reveals that the mode of interaction of the amino acid group with the enzyme is different from that previously observed for acyl amino acid substrates. Further, required movements of residues in the active site to accommodate the imidazole group provide an explanation for the low affinity of imidazole itself. Finally, the previously observed cooperativity with lauric acid is explained by a surprisingly open substrate-access channel lined with hydrophobic residues that could potentially accommodate lauric acid in addition to the inhibitor itself.     

Three clusters of conformational states in p450cam reveal a multistep pathway for closing of the substrate access channel

Conformational changes in the substrate access channel have been observed for several forms of cytochrome P450, but the extent of conformational plasticity exhibited by a given isozyme has not been completely characterized. Here we present crystal structures of P450cam bound to a library of 12 active site probes containing a substrate analogue tethered to a variable linker. The structures provide a unique view of the range of protein conformations accessible during substrate binding. Principal component analysis of a total of 30 structures reveals three discrete clusters of conformations: closed (P450cam-C), intermediate (P450cam-I), and fully open (P450cam-O). Relative to P450cam-C, the P450cam-I state results predominantly from a retraction of helix F, while both helices F and G move in concert to reach the fully open P450cam-O state. Both P450cam-C and P450cam-I are well-defined states, while P450cam-O shows evidence of a somewhat broader distribution of conformations and includes the open form recently seen in the absence of substrate. The observed clustering of protein conformations over a wide range of ligand variants suggests a multistep closure of the enzyme around the substrate that begins by conformational selection from an ensemble of open conformations and proceeds through a well-defined intermediate, P450cam-I, before full closure to the P450cam-C state in the presence of small substrates. This multistep pathway may have significant implications for a full understanding of substrate specificity, kinetics, and coupling of substrate binding to P450 function.     

Structural analysis of mammalian cytochrome P450 2B4 covalently bound to the mechanism-based inactivator tert-butylphenylacetylene: insight into partial enzymatic activity

A combined structural and computational analysis of rabbit cytochrome P450 2B4 covalently bound to the mechanism-based inactivator tert-butylphenylacetylene (tBPA) has yielded insight into how the enzyme retains partial activity. Since conjugation to tBPA modifies a highly conserved active site residue, the residual activity of tBPA-labeled 2B4 observed in previous studies was puzzling. Here we describe the first crystal structures of a modified mammalian P450, which show an oxygenated metabolite of tBPA conjugated to Thr 302 of helix I. These results are consistent with previous studies that identified Thr 302 as the site of conjugation. In each structure, the core of 2B4 remains unchanged, but the arrangement of plastic regions differs. This results in one structure that is compact and closed. In this conformation, tBPA points toward helix B', making a 31° angle with the heme plane. This conformation is in agreement with previously performed in silico experiments. However, dimerization of 2B4 in the other structure, which is caused by movement of the B/C loop and helices F through G, alters the position of tBPA. In this case, tBPA lies almost parallel to the heme plane due to the presence of helix F' of the opposite monomer entering the active site to stabilize the dimer. However, docking experiments using this open form show that tBPA is able to rotate upward to give testosterone and 7-ethoxy-4-trifluoromethylcoumarin access to the heme, which could explain the previously observed partial activity.     

What common structural features and variations of mammalian P450s are known to date?

Sufficient structural information on mammalian cytochromes P450 has now been published (including seventeen X-ray structures of these enzymes by June 2006) to allow characteristic features of these enzymes to be identified, including: (i) the presence of a common fold, typical of all P450s, (ii) similarities in the positioning of the heme cofactor, (iii) the spatial arrangement of certain structural elements, and (iv) the access/egress paths for substrates and products, (v) probably common orientation in the membrane, (vi) characteristic properties of the active sites with networks of water molecules, (vii) mode of interaction with redox partners and (viii) a certain degree of flexibility of the structure and active site determining the ease with which the enzyme may bind the substrates. As well as facilitating the identification of common features, comparison of the available structures allows differences among the structures to be identified, including variations in: (i) preferred access/egress paths to/from the active site, (ii) the active site volume and (iii) flexible regions. The availability of crystal structures provides opportunities for molecular dynamic simulations, providing data that are apparently complementary to experimental findings but also allow the dynamic behavior of access/egress paths and other dynamic features of the enzymes to be explored.     

A Multiscale Approach to Modelling Drug Metabolism by Membrane-Bound Cytochrome P450 Enzymes

Cytochrome P450 enzymes are found in all life forms. P450s play an important role in drug metabolism, and have potential uses as biocatalysts. Human P450s are membrane-bound proteins. However, the interactions between P450s and their membrane environment are not well-understood. To date, all P450 crystal structures have been obtained from engineered proteins, from which the transmembrane helix was absent. A significant number of computational studies have been performed on P450s, but the majority of these have been performed on the solubilised forms of P450s. Here we present a multiscale approach for modelling P450s, spanning from coarse-grained and atomistic molecular dynamics simulations to reaction modelling using hybrid quantum mechanics/molecular mechanics (QM/MM) methods. To our knowledge, this is the first application of such an integrated multiscale approach to modelling of a membrane-bound enzyme. We have applied this protocol to a key human P450 involved in drug metabolism: CYP3A4. A biologically realistic model of CYP3A4, complete with its transmembrane helix and a membrane, has been constructed and characterised. The dynamics of this complex have been studied, and the oxidation of the anticoagulant R-warfarin has been modelled in the active site. Calculations have also been performed on the soluble form of the enzyme in aqueous solution. Important differences are observed between the membrane and solution systems, most notably for the gating residues and channels that control access to the active site. The protocol that we describe here is applicable to other membrane-bound enzymes.     

Crystal structure of the Mycobacterium tuberculosis P450 CYP121-fluconazole complex reveals new azole drug-P450 binding mode

Azole and triazole drugs are cytochrome P450 inhibitors widely used as fungal antibiotics and possessing potent antimycobacterial activity. We present here the crystal structure of Mycobacterium tuberculosis cytochrome P450 CYP121 in complex with the triazole drug fluconazole, revealing a new azole heme ligation mode. In contrast to other structurally characterized cytochrome P450 azole complexes, where the azole nitrogen directly coordinates the heme iron, in CYP121 fluconazole does not displace the aqua sixth heme ligand but occupies a position that allows formation of a direct hydrogen bond to the aqua sixth heme ligand. Direct ligation of fluconazole to the heme iron is observed in a minority of CYP121 molecules, albeit with severe deviations from ideal geometry due to close contacts with active site residues. Analysis of both ligand-on and -off structures reveals the relative position of active site residues derived from the I-helix is a key determinant in the relative ratio of on and off states. Regardless, both ligand-bound states lead to P450 inactivation by active site occlusion. This previously unrecognized means of P450 inactivation is consistent with spectroscopic analyses in both solution and in the crystalline form and raises important questions relating to interaction of azoles with both pathogen and human P450s.     

Key Mutations Alter the Cytochrome P450 BM3 Conformational Landscape and Remove Inherent Substrate Bias

Cytochrome P450 monooxygenases (P450s) have enormous potential in the production of oxychemicals, due to their unparalleled regio- and stereoselectivity. The Bacillus megaterium P450 BM3 enzyme is a key model system, with several mutants (many distant from the active site) reported to alter substrate selectivity. It has the highest reported monooxygenase activity of the P450 enzymes, and this catalytic efficiency has inspired protein engineering to enable its exploitation for biotechnologically relevant oxidations with structurally diverse substrates. However, a structural rationale is lacking to explain how these mutations have such effects in the absence of direct change to the active site architecture. Here, we provide the first crystal structures of BM3 mutants in complex with a human drug substrate, the proton pump inhibitor omeprazole. Supported by solution data, these structures reveal how mutation alters the conformational landscape and decreases the free energy barrier for transition to the substrate-bound state. Our data point to the importance of such “gatekeeper” mutations in enabling major changes in substrate recognition. We further demonstrate that these mutants catalyze the same 5-hydroxylation reaction as performed by human CYP2C19, the major human omeprazole-metabolizing P450 enzyme.     

Probing ligand binding modes of human cytochrome P450 2J2 by homology modeling, molecular dynamics simulation, and flexible molecular docking

Cytochrome P450 (P450) 2J2 catalyzes epoxidation of arachidonic acid to eicosatrienoic acids, which are related to a variety of diseases such as coronary artery disease, hypertension, and carcinogenesis. Recent experimental data also suggest that P450 2J2 could be a novel biomarker and a potential target for cancer therapy. However, the active site topology and substrate specificity of this enzyme remain unclear. In this study, a three-dimensional model of human P450 2J2 was first constructed on the basis of the crystal structure of human P450 2C9 in complex with a substrate using homology modeling method, and refined by molecular dynamics simulation. Flexible docking approaches were then employed to dock four ligands into the active site of P450 2J2 in order to probe the ligand-binding modes. By analyzing the results, active site architecture and certain key residues responsible for substrate specificity were identified on the enzyme, which might be very helpful for understanding the enzyme's biological role and providing insights for designing novel inhibitors of P450 2J2.     

The role of cytochrome P450 BM3 phenylalanine-87 and threonine-268 in binding organic hydroperoxides

BackgroundCytochrome P450 (P450) BM3, from Bacillus megaterium, catalyzes a wide range of chemical reactions and is routinely used as a model system to study mammalian P450 reactions and structure.MethodsThe metabolism of 2,6-di-tert-butyl-4-hydroperoxy-4-methyl-2,5-cyclohexadienone (BHTOOH) and 2-tert-butyl-4-hydroperoxy-4-methyl-2,5-cyclohexadien-1-one (BMPOOH) was examined with P450 BM3 and with the conserved T268 and F87 residues mutated to investigate their effects on organic hydroperoxide metabolism. To determine the effects of the mutations on the active site volume and architecture, the X-ray crystal structure of the F87A/T268A P450 BM3 heme domain (BMP) was determined and compared to previous structures. To investigate the interactions of the substrates with the F87 and T268 residues, BHTOOH and BMPOOH were docked into the BMP X-ray crystal structures.ResultsLower metabolism of BHTOOH and BMPOOH was observed in the WT P450 BM3 and the T268A P450 BM3 mutant than in the F87A and F87A/T268A P450 BM3 mutants. Large differences were found in the F–G loop regions and active site cavity volumes for the F87A mutated structures.ConclusionsAnalysis of the metabolism, X-ray crystal structures, and molecular docking simulations suggests that P450 BM3 activity toward BHTOOH and BMPOOH is mediated through substrate recognition by T268 and F87, and the active site cavity volume. Based on this information, a simplified representation is presented with the relative orientation of organic hydroperoxides in the P450 BM3 active site.General significanceThe metabolism results and structural analysis of this model P450 allowed us to rationalize the structural factors that influence organic hydroperoxide metabolism.     

Exploring coumarin egress channels in human cytochrome P450 2A6 by random acceleration and steered molecular dynamics simulations

The kinetic analysis of coumarin oxidation by CYP2A6 suggested that substrate binding and release occurred in the multiple steps and such events proceeded rapidly. However, the crystal structure of the CYP2A6-coumarin complex reveals that no obvious channel is open enough to allow coumarin to pass through. Thus, an intriguing and important question arises: how coumarin enters and exits the active site, which is deeply buried at the center of CYP2A6 fold. In this study, geometric analysis of the potential openings was first performed on all the available crystal structures of CYP2A6. And then, random acceleration molecular dynamics simulations were used to explore the possible substrate egress channels in CYP2A6. Two channels were most frequently observed. Afterwards, steered molecular dynamics simulations were performed and potentials of mean force were constructed to compare the preference of the two channels serving as the substrate egress channel. The results showed that channel 2c, which is located between helices I and G and the helix B'-C region, was the most likely channel for coumarin egress. The opening of channel 2c was characterized by a rotation of Phe111 together with a bending of helix B'. Our findings will not only be helpful for understanding the unbinding mechanism of coumarin and for identifying structural determinants related to the biological function of CYP2A6, but also provide further insight into the channel selectivity of P450s.     

Exploring the binding sites of the haloalkane dehalogenase DhlA from Xanthobacter autotrophicus GJ10

The catalytic site of haloalkane dehalogenase DhlA is buried more than 10 A from the protein surface. While potential access channels to this site have been reported, the precise mechanism of substrate import and product export is still unconfirmed. We used computational methods to examine surface pockets and their putative roles in ligand access to and from the catalytic site. Computational solvent mapping moves small organic molecule as probes over the protein surface in order to identify energetically favorable sites, that is, regions that tend to bind a variety of molecules. The mapping of three DhlA structures identifies seven such regions, some of which have been previously suggested to be involved in the binding and the import/export of substrates or products. These sites are the active site, the putative entrance of the channel leading to the active site, two pockets that bind Br- ions, a pocket in the slot region, and two additional sites between the main domain and the cap of DhlA. We also performed mapping and free energy analysis of the DhlA structures using the substrate, 1,2-dichloroethane, and halide ions as probes. The findings were compared to crystallographic data and to results obtained by CAVER, a program developed for finding routes from protein clefts and cavities to the surface. Solvent mapping precisely reproduced all three Br- binding sites identified by protein crystallography and the openings to four channels found by CAVER. The analyses suggest that (i) the active site has the highest affinity for the substrate molecule, (ii) the substrate initially binds at the entrance of the main tunnel, (iii) the site Br2, close to the entrance, is likely to serve as an intermediate binding site in product export, (iv) the site Br3, induced in the structure at high concentrations of Br-, could be part of an auxiliary route for product release, and (v) three of the identified sites are likely to be entrances of water-access channels leading to the active site. For comparison, we also mapped haloalkane dehalogenases DhaA and LinB, both of which contain significantly larger and more solvent accessible binding sites than DhlA. The mapping of DhaA and LinB places the majority of probes in the active site, but most of the other six regions consistently identified in DhlA were not observed, suggesting that the more open active site eliminates the need for intermediate binding sites for the collision complex seen in DhlA.