Crystal structure of cytochrome P-450cam complexed with the (1S)-camphor enantiomer

Sniffing out camphor: the fine balance between hydrogen bonding and London dispersion in the chirality recognition with α-fenchol

Binary complexes between the chiral monoterpenoids camphor and α-fenchol were explored with vibrational and rotational jet spectroscopy as well as density functional theory in order to explore how chirality can influence the binding preferences in gas-phase complexes. The global minimum structures of the two diastereomers were assigned. It is found that chirality recognition leads to different compromises in the fine balance between intermolecular interactions. While one isomer features a stronger hydrogen bond, the other one is more tightly arranged and stabilized by larger London dispersion interactions. These new spectroscopic results help understand the influence of chirality in molecular aggregation and unveil the kind of interactions involved between a chiral alcohol and a chiral ketone with large dispersion contributions.

Bioinorganic Chemistry on Electrodes: Methods to Functional Modeling

One of the major goals of bioinorganic chemistry has been to mimic the function of elegant metalloenzymes. Such functional modeling has been difficult to attain in solution, in particular, for reactions that require multiple protons and multiple electrons (nH⁺/ne^-). Using a combination of heterogeneous electrochemistry, electrode and molecule design one may control both electron transfer (ET) and proton transfer (PT) of these nH⁺/ne^- reactions. Such control can allow functional modeling of hydrogenases (H⁺ + e^- → 1/2 H₂), cytochrome c oxidase (O₂ + 4 e^- + 4 H⁺ → 2 H₂O), monooxygenases (RR'CH₂ + O₂ + 2 e^- + 2 H⁺ → RR'CHOH + H₂O) and dioxygenases (S + O₂ → SO₂; S = organic substrate) in aqueous medium and at room temperatures. In addition, these heterogeneous constructs allow probing unnatural bioinspired reactions and estimation of the inner- and outer-sphere reorganization energy of small molecules and proteins.

Rejigging Electron and Proton Transfer to Transition between Dioxygenase, Monooxygenase, Peroxygenase, and Oxygen Reduction Activity: Insights from Bioinspired Constructs of Heme Enzymes

Nature has employed heme proteins to execute a diverse set of vital life processes. Years of research have been devoted to understanding the factors which bias these heme enzymes, with all having a heme cofactor, toward distinct catalytic activity. Among them, axial ligation, distal super structure, and substrate binding pockets are few very vividly recognized ones. Detailed mechanistic investigation of these heme enzymes suggested that several of these enzymes, while functionally divergent, use similar intermediates. Furthermore, the formation and decay of these intermediates depend on proton and electron transfer processes in the enzyme active site. Over the past decade, work in this group, using in situ surface enhanced resonance Raman spectroscopy of synthetic and biosynthetic analogues of heme enzymes, a general idea of how proton and electron transfer rates relate to the lifetime of different O2 derived intermediates has been developed. These findings suggest that the enzymatic activities of all these heme enzymes can be integrated into one general cycle which can be branched out to different catalytic pathways by regulating the lifetime and population of each of these intermediates. This regulation can further be achieved by tuning the electron and proton transfer steps. By strategically populating one of these intermediates during oxygen reduction, one can navigate through different catalytic processes to a desired direction by altering proton and electron transfer steps.

Structure-based mechanisms: On the way to apply alcohol dehydrogenases/reductases to organic-aqueous systems

Alcohol dehydrogenases/reductases catalyze enantioselective syntheses of versatile chiral compounds relying on direct hydride transfer from cofactor to substrates, or to an intermediate and then to substrates. Since most of the substrates catalyzed by alcohol dehydrogenases/reductases are insoluble in aqueous solutions, increasing interest has been turning to organic-aqueous systems. However, alcohol dehydrogenases/reductases are normally instable in organic solvents, leading to the unsatisfied enantioselective synthesis efficiency. The behaviors of these enzymes in organic solvents at an atomic level are unclear, thus it is of great importance to understand its structure-based mechanisms in organic-aqueous systems to improve their relative stability. Here, we summarized the accessible structures of alcohol dehydrogenases/reductases in Protein Data Bank crystallized in organic-aqueous systems, and compared the structures of alcohol dehydrogenases/reductases which have different tolerance towards organic solvents. By understanding the catalytic behaviors and mechanisms of these enzymes in organic-aqueous systems, the efficient enantioselective syntheses mediated by alcohol dehydrogenases/reductases and further challenges are also discussed through solvent engineering and enzyme-immobilization in the last decade.

A Single Iron Porphyrin Shows pH Dependent Switch between "Push" and "Pull" Effects in Electrochemical Oxygen Reduction

The "push-pull" effects associated with heme enzymes manifest themselves through highly evolved distal amino acid environments and axial ligands to the heme. These conserved residues enhance their reactivities by orders of magnitude relative to small molecules that mimic the primary coordination. An instance of a mononuclear iron porphyrin with covalently attached pendent phenanthroline groups is reported which exhibit reactivity indicating a pH dependent "push" to "pull" transition in the same molecule. The pendant phenanthroline residues provide proton transfer pathways into the iron site, ensuring selective 4e^-/4H⁺ reduction of O₂ to water. The protonation of these residues at lower pH mimics the pull effect of peroxidases, and a coordination of an axial hydroxide ligand at high pH emulates the push effect of P450 monooxygenases. Both effects enhance the rate of O₂ reduction by orders of magnitude over its value at neutral pH while maintaining exclusive selectivity for 4e^-/4H⁺ oxygen reduction reaction.

Heme-FeIII Superoxide, Peroxide and Hydroperoxide Thermodynamic Relationships: FeIII-O2•- Complex H-Atom Abstraction Reactivity

Establishing redox and thermodynamic relationships between metal-ion-bound O2 and its reduced (and protonated) derivatives is critically important for a full understanding of (bio)chemical processes involving dioxygen processing. Here, a ferric heme peroxide complex, [(F8)FeIII-(O22-)]- (P) (F8 = tetrakis(2,6-difluorophenyl)porphyrinate), and a superoxide complex, [(F8)FeIII-(O2•-)] (S), are shown to be redox interconvertible. Using Cr(η-C6H6)2, an equilibrium state where S and P are present is established in tetrahydrofuran (THF) at -80 °C, allowing determination of the reduction potential of S as -1.17 V vs Fc+/0. P could be protonated with 2,6-lutidinium triflate, yielding the low-spin ferric hydroperoxide species, [(F8)FeIII-(OOH)] (HP). Partial conversion of HP back to P using a derivatized phosphazene base gave a P/HP equilibrium mixture, leading to the determination of pKa = 28.8 for HP (THF, -80 °C). With the measured reduction potential and pKa, the O-H bond dissociation free energy (BDFE) of hydroperoxide species HP was calculated to be 73.5 kcal/mol, employing the thermodynamic square scheme and Bordwell relationship. This calculated O-H BDFE of HP, in fact, lines up with an experimental demonstration of the oxidizing ability of S via hydrogen atom transfer (HAT) from TEMPO-H (2,2,6,6-tetramethylpiperdine-N-hydroxide, BDFE = 66.5 kcal/mol in THF), forming the hydroperoxide species HP and TEMPO radical. Kinetic studies carried out with TEMPO-H(D) reveal second-order behavior, kH = 0.5, kD = 0.08 M-1 s-1 (THF, -80 °C); thus, the hydrogen/deuterium kinetic isotope effect (KIE) = 6, consistent with H-atom abstraction by S being the rate-determining step. This appears to be the first case where experimentally derived thermodynamics lead to a ferric heme hydroperoxide OO-H BDFE determination, that FeIII-OOH species being formed via HAT reactivity of the partner ferric heme superoxide complex.

Analyses of cobalt-ligand and potassium-ligand bond lengths in metalloproteins: trends and patterns

Cobalt and potassium are biologically important metal elements that are present in a large array of proteins. Cobalt is mostly found in vivo associated with a corrin ring, which represents the core of the vitamin B12 molecule. Potassium is the most abundant metal in the cytosol, and it plays a crucial role in maintaining membrane potential as well as correct protein function. Here, we report a thorough analysis of the geometric properties of cobalt and potassium coordination spheres that was performed with high resolution on a representative set of structures from the Protein Data Bank and complemented by quantum mechanical calculations realized at the DFT level of theory (B3LYP/ SDD) on mononuclear model systems. The results allowed us to draw interesting conclusions on the structural characteristics of both Co and K centers, and to evaluate the importance of effects such as their association energies and intrinsic thermodynamic stabilities. Overall, the results obtained provide useful data for enhancing the atomic models normally applied in theoretical and computational studies of Co or K proteins performed at the quantum mechanical level, and for developing molecular mechanical parameters for treating Co or K coordination spheres in molecular mechanics or molecular dynamics studies.

CO, NO and O2 as Vibrational Probes of Heme Protein Interactions

The gaseous XO molecules (X = C, N or O) bind to the heme prosthetic group of heme proteins, and thereby activate or inhibit key biological processes. These events depend on interactions of the surrounding protein with the FeXO adduct, interactions that can be monitored via the frequencies of the Fe-X and X-O bond stretching modes, νFeX and νXO. The frequencies can be determined by vibrational spectroscopy, especially resonance Raman spectroscopy. Backbonding, the donation of Fe dπ electrons to the XO π* orbitals, is a major bonding feature in all the FeXO adducts. Variations in backbonding produce negative νFeX/νXO correlations, which can be used to gauge electrostatic and H-bonding effects in the protein binding pocket. Backbonding correlations have been established for all the FeXO adducts, using porphyrins with electron donating and withdrawing substituents. However the adducts differ in their response to variations in the nature of the axial ligand, and to specific distal interactions. These variations provide differing vantages for evaluating the nature of protein-heme interactions. We review experimental studies that explore these variations, and DFT computational studies that illuminate the underlying physical mechanisms.

A novel type of allosteric regulation: functional cooperativity in monomeric proteins

Cooperative functional properties and allosteric regulation in cytochromes P450 play an important role in xenobiotic metabolism and define one of the main mechanisms of drug-drug interactions. Recent experimental results suggest that ability to bind simultaneously two or more small organic molecules can be the essential feature of cytochrome P450 fold, and often results in rich and complex pattern of allosteric behavior. Manifestations of non-Michaelis kinetics include homotropic and heterotropic activation and inhibition effects depending on the stoichiometric ratios of substrate and effector, changes in the regio- and stereospecificity of catalytic transformations, and often give rise to the clinically important drug-drug interactions. In addition, functional response of P450 systems is modulated by the presence of specific and non-specific effector molecules, metal ions, membrane incorporation, formation of homo- and hetero-oligomers, and interactions with the protein redox partners. In this article we briefly overview the main factors contributing to the allosteric effects in cytochromes P450 with the main focus on the sources of cooperative behavior in xenobiotic metabolizing monomeric heme enzymes with their conformational flexibility and extremely broad substrate specificity. The novel mechanism of functional cooperativity in P450 enzymes does not require substantial binding cooperativity, rather it implies the presence of one or more binding sites with higher affinity than the single catalytically active site in the vicinity of the heme iron.

A PREDICTIVE CYTOCHROME P450 MONOOXYGENASE FUNCTIONAL MODEL FOR GENERIC HYDROXYLATION BY RHIZOPUS ORYZAE ATCC 11145

An active site model for generic hydroxylation of xenobiotes by the fungus Rhizopus oryzae ATCC 11145 has been developed. The model was constructed using data from known metabolites of R. oryzae, as reported in the chemical literature, including diterpenes and steroids. The model takes the form of an irregular polyhedron of 340 Å3 volume, with two hydrophilic patches and an active site region. It was validated using data for other known substrates, including terpenes, steroids, octalin, and hydrindenone derivatives, drawn from the literature. The interpretive and predictive abilities of the model are discussed.

Changing the substrate specificity of P450cam towards diphenylmethane by semi-rational enzyme engineering

A focused library comprising nine residues of the active site of P450cam monooxygenase resulting in ∼ 300,000 protein variants was screened for activity on diphenylmethane (DPM). The assay was based on the depletion of NADH by an in vitro reconstituted P450cam system in a 96-well scale. The throughput was increased by the parallel cultivation, purification and analysis of 20 variants per well (cluster screening). Thus ∼ 20,000 protein variants could be screened in summary of which five were found to transform DPM with a specific activity of up to 75% of the wild-type activity on d-camphor and a coupling rate of 7-18%. One variant converting DPM to 4-hydroxydiphenylmethane (4HDPM) was subjected to site-directed mutagenesis and saturation mutagenesis, which revealed the particular importance of positions F87, Y96 and L244 for substrate selectivity and the possibility for further improvements of this variant. Moreover, a reduction in size of the amino acid at position 396 decreased specific activity dramatically but increased coupling and switched the main product formation from 4HDPM towards diphenylmethanol.

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.

The rate-limiting step in P450 hydroxylation of hydrocarbons a direct comparison of the "somersault" versus the "consensus" mechanism involving compound I

Model theoretical quantum mechanical (QM) calculations are described for the P-450 hydroxylation of methane, isobutane, and camphor that compare the concerted somersault H-abstraction mechanism with the oxidation step involving Cpd I. Special emphasis has been placed on maintaining a balanced basis set in the oxidation step. QM calculations, employing the 6-311+G(d,p) basis set on the Fe atom and all of the key surrounding atoms involved in the C-H abstraction step, reaffirm a mechanism involving rearrangement of the iron hydroperoxide group (FeO-OH --> FeO...HO(*)) in concert with hydrogen abstraction from the C-H bond of the substrate by the incipient bound hydroxyl radical HO(*). The barrier for the somersault rearrangement of model Cpd 0 (FeO-OH) is calculated to be 21.4 kcal/mol in the absence of substrate. The overall activation energy for the oxidation of camphor involving the somersault motion of the FeO-OH group of P450 model porphyrin iron(III) hydroperoxide [Por(SH)Fe(III)-OOH(-)] --> [Por(SH)Fe(III)-O....HO(-)] in concert with hydrogen abstraction is DeltaE(++) = 12.4 kcal/mol. The corresponding abstraction of the hydrogen atom from the C-H bond of camphor by Cpd I has an activation barrier of 17.6 kcal/mol. Arguments are presented that the somersault rearrangement is induced by steric compression at the active site. Kinetic isotope effect data are discussed that provides compelling evidence for a rate-limiting step involving C-H bond cleavage.

Conformational plasticity and structure/function relationships in cytochromes P450

The cytochrome P450s are a superfamily of enzymes that are found in all kingdoms of living organisms, and typically catalyze the oxidative addition of atomic oxygen to an unactivated C-C or C-H bond. Over 8000 nonredundant sequences of putative and confirmed P450 enzymes have been identified, but three-dimensional structures have been determined for only a small fraction of these. While all P450 enzymes for which structures have been determined share a common global fold, the flexibility and modularity of structure around the active site account for the ability of P450 enzymes to accommodate a vast number of structurally dissimilar substrates and support a wide range of selective oxidations. In this review, known P450 structures are compared, and some structural criteria for prediction of substrate selectivity and reaction type are suggested. The importance of dynamic processes such as redox-dependent and effector-induced conformational changes in determining catalytic competence and regio- and stereoselectivity is discussed, and noncrystallographic methods for characterizing P450 structures and dynamics, in particular, mass spectrometry and nuclear magnetic resonance spectroscopy are reviewed.

S K-edge XAS and DFT calculations on cytochrome P450: covalent and ionic contributions to the cysteine-Fe bond and their contribution to reactivity

Experimental covalencies of the Fe-S bond for the resting low-spin and substrate-bound high-spin active site of cytochrome P450 are reported. DFT calculations on the active site indicate that one H-bonding interaction from the protein backbone is needed to reproduce the experimental values. The H-bonding to the thiolate from the backbone decreases the anisotropic pi covalency of the Fe-S bond lowering the barrier of free rotation of the exchangeable axial ligand, which is important for reactivity. The anionic axial thiolate ligand is calculated to lower the Fe(III/II) reduction potential of the active site by more than 1 V compared to a neutral imidazole ligand. About half of this derives from its covalent bonding and half from its electrostatic interaction with the oxidized Fe. This axial thiolate ligand increases the pK(a) of compound 0 (Fe(III)-hydroperoxo) favoring its protonation which promotes O-O bond heterolysis forming compound I. The reactivity of compound I is calculated to be relatively insensitive to the nature of the axial ligand due to opposing reduction potential and proton affinity contributions to the H-atom abstraction energy.

MolAxis: efficient and accurate identification of channels in macromolecules

Channels and cavities play important roles in macromolecular functions, serving as access/exit routes for substrates/products, cofactor and drug binding, catalytic sites, and ligand/protein. In addition, channels formed by transmembrane (TM) proteins serve as transporters and ion channels. MolAxis is a new sensitive and fast tool for the identification and classification of channels and cavities of various sizes and shapes in macromolecules. MolAxis constructs corridors, which are pathways that represent probable routes taken by small molecules passing through channels. The outer medial axis of the molecule is the collection of points that have more than one closest atom. It is composed of two-dimensional surface patches and can be seen as a skeleton of the complement of the molecule. We have implemented in MolAxis a novel algorithm that uses state-of-the-art computational geometry techniques to approximate and scan a useful subset of the outer medial axis, thereby reducing the dimension of the problem and consequently rendering the algorithm extremely efficient. MolAxis is designed to identify channels that connect buried cavities to the outside of macromolecules and to identify TM channels in proteins. We apply MolAxis to enzyme cavities and TM proteins. We further utilize MolAxis to monitor channel dimensions along Molecular Dynamics trajectories of a human Cytochrome P450. MolAxis constructs high quality corridors for snapshots at picosecond time-scale intervals substantiating the gating mechanism in the 2e substrate access channel. We compare our results with previous tools in terms of accuracy, performance and underlying theoretical guarantees of finding the desired pathways. MolAxis is available on line as a web-server and as a stand alone easy-to-use program (http://bioinfo3d.cs.tau.ac.il/MolAxis/).

Conformational equilibrium of cytochrome P450 BM-3 complexed with N-palmitoylglycine: a replica exchange molecular dynamics study

UV-vis absorbance measurements and associated studies of cytochrome P450 BM-3 in complex with N-palmitoylglycine (NPG) indicate that a conformational change occurs in the active site of the complex where the terminal atoms of the ligand move from a site distant from the heme iron, as seen in the low temperature crystal structure to a site proximal to the heme iron at biological temperatures. We employ replica exchange molecular dynamics simulations to study this conformational change. The population of the proximal state is found to increase with temperature in agreement with UV-vis absorbance and NMR measurements. In addition to the conformations characterized by X-ray crystallography and computer modeling, this study shows that a new conformational state is significantly populated at room temperature. The observed increase in the population of conformations where the terminal atoms of NPG are proximal to the heme iron with increasing temperature indicates that the proximal state is stabilized by conformational entropy. A proposal for the origin of this entropic stabilization is provided on the basis of the structure of the newly identified state. We use the temperature weighted histogram (T-WHAM) method to characterize the transition state regions of the conformational ensemble and propose a mechanism of interconversion between these low free energy conformational states.

Second derivatives in generalized Born theory

Generalized Born solvation models offer a popular method of including electrostatic aspects of solvation free energies within an analytical model that depends only upon atomic coordinates, charges, and dielectric radii. Here, we describe how second derivatives with respect to Cartesian coordinates can be computed in an efficient manner that can be distributed over multiple processors. This approach makes possible a variety of new methods of analysis for these implicit solvation models. We illustrate three of these methods here: the use of Newton-Raphson optimization to obtain precise minima in solution; normal mode analysis to compute solvation effects on the mechanical properties of DNA; and the calculation of configurational entropies in the MM/GBSA model. An implementation of these ideas, using the Amber generalized Born model, is available in the nucleic acid builder (NAB) code, and we present examples for proteins with up to 45,000 atoms. The code has been implemented for parallel computers using both the OpenMP and MPI environments, and good parallel scaling is seen with as many as 144 OpenMP processing threads or MPI processing tasks.

Conformational dynamics of the cytochrome P450 BM3/N-palmitoylglycine complex: the proposed "proximal-distal" transition probed by temperature-jump spectroscopy

The ferric spin state equilibrium of the heme iron was analyzed in wild-type cytochrome P450 BM3 and its F87G mutant by using temperature (T)-jump relaxation spectroscopy in combination with static equilibrium experiments. No relaxation process was measurable in the substrate-free enzyme indicating a relaxation process with a rate constant>10,000 s(-1). In contrast, a slow spin state transition process was observed in the N-palmitoylglycine (NPG)-bound enzyme species. This transition occurred with an observed rate constant (298 K) of approximately 800 s(-1) in the wild-type, and approximately 2500 s(-1) in the F87G mutant, suggesting a significant contribution of the phenylalanine side chain to a reaction step rate limiting the actual spin state transition. These findings are discussed in terms of an equilibrium between different binding modes of the substrate, including a position 7.5 A away from the heme iron ("distal") and the catalytically relevant "proximal" binding site, and are in accordance with results from X-ray crystallography, NMR studies, and molecular dynamics simulations.

The status of high-valent metal oxo complexes in the P450 cytochromes

The oxidative prowess of the P450 cytochromes in physiological reactions is attributed to the production of a high-valent iron-oxo complex, or Compound I intermediate, in the reaction cycle. Despite many years of study, however, the full electronic description of this fleeting intermediate still remains an active area of study. In this manuscript, the current status of the isolation and characterization of the P450 oxo-Fe(IV) is examined and compared to analogous states in related heme enzymes. In addition, the utilization of cofactor exchange to stabilize high-valent oxo-states in the P450 is addressed. Structural and spectroscopic studies on manganese reconstituted P450, and its corresponding oxo-complex, are presented.

Thermal Equilibrium of High- and Low-Spin Forms of Cytochrome P450 BM-3: Repositioning of the Substrate?

We demonstrate that cytochrome P450 BM-3 in complex with N-palmitoylglycine undergoes a spin state change between room temperature, where optimal activity is seen, and low temperatures, where X-ray diffraction characterization has been carried out. On the basis of NMR measurements of the full-length protein, this spin state change is likely to be accompanied by a general structural rearrangement in the enzyme pocket. The substrate remains bound at all temperatures. We propose that the substrate may "slide" from a position directly atop the heme (thus displacing the ligating water) to the more distant position (thus restoring the ligating water) as the temperature is lowered. This proposal is evaluated on the basis of computational modeling of the protein-ligand complex, using a novel induced fit methodology. We thereby generate a structure with the ligand in close contact with the heme, similar in energy to the experimental structure. With this combination of theory and experiment we provide a specific proposal of how ligands may be positioned for chemistry for this enzyme.

Cytochrome P450cam-monoterpene interactions

The use of the molecular mechanics AMBER force field (FF) to predict product profiles for the hydroxylation of the monoterpenes 1R-camphor, 1S-camphor, 1R-norcamphor, 1S-norcamphor and camphane by the enzyme cytochrome P450cam from the soil bacterium Pseudomonas putida was investigated. Predictions were carried out by applying multiple substrate (starting) orientations in the enzyme pocket in two procedures: a procedure based on molecular dynamics (MD) and a procedure based on short MD simulations followed by geometry optimisations. The latter (GO) procedure is faster and enabled the use of more monoterpene starting orientations. Monoterpene orientations were transformed into product profiles by applying both energetic and geometrical criteria appropriate for the (monoterpene) hydrogen abstraction reaction. Good predictions compared to experimental data were obtained for most compounds in both the MD and GO procedures. Prior to the product profile calculations, the FF was calibrated by reproducing the experimental data for the binding energy of 1R-camphor and 1S-camphor to P450cam and the energy of vaporisation of water. Focus of the calibration was on the value for the scaling factor for the electrostatic interactions.

Identification of the amino acid residue of CYP27B1 responsible for binding of 25-hydroxyvitamin D3 whose mutation causes vitamin D-dependent rickets type 1

We previously reported the three-dimensional structure of human CYP27B1 (25-hydroxyvitamin D3 1alpha-hydroxylase) constructed by homology modeling. Using the three-dimensional model we studied the docking of the substrate, 25-hydroxyvitamin D3, into the substrate binding pocket of CYP27B1. In this study, we focused on the amino acid residues whose point mutations cause vitamin D-dependent rickets type 1, especially unconserved residues among mitochondrial CYPs such as Gln65 and Thr409. Recently, we successfully overexpressed mouse CYP27B1 by using a GroEL/ES co-expression system. In a mutation study of mouse CYP27B1 that included spectroscopic analysis, we concluded that in a 1alpha-hydroxylation process, Ser408 of mouse CYP27B1 corresponding to Thr409 of human CYP27B1 forms a hydrogen bond with the 25-hydroxyl group of 25-hydroxyvitamin D3. This is the first report that shows a critical amino acid residue recognizing the 25-hydroxyl group of the vitamin D3.

The influence of substrate on the spectral properties of oxyferrous wild-type and T252A cytochrome P450-CAM

To probe whether the nature of the substrate can directly influence the spectral properties of oxyferrous cytochrome P450-CAM, the complex has been investigated in the absence and in the presence of the natural substrate (1R)-camphor (camphor) and of several camphor analogs. The oxyferrous complex of T252A P450-CAM, a mutant lacking the hydroxyl group that forms a hydrogen bond to the heme iron-coordinated dioxygen, has also been studied to gauge the influence of this hydrogen bond. UV-visible absorption and magnetic circular dichroism (MCD) spectra of these oxyferrous adducts prepared and stabilized at -40 degrees C in 60% (v/v) ethylene glycol are generally similar, exhibiting absorption bands at approximately 355, approximately 420, approximately 554, and approximately 585 nm (shoulder) and a characteristic MCD trough at approximately 585 nm. The MCD spectrum of camphor-bound oxyferrous P450-CAM is similar to that of the substrate-free oxyferrous enzyme, but the spectrum of the oxyferrous enzyme differs detectably in the presence of substrate analogs. The spectra of the oxyferrous T252A mutant and wild-type enzyme are overall similar except for Soret band position blue shifts by 2-6 nm for the mutant. 5-Methylenylcamphor (epoxidation substrate) appears to have an anomalous binding mode for the mutant compared with that for the wild-type enzyme. The present results indicate that the structures of the camphor analogs can sensitively influence the physical (spectroscopic) properties of the P450 dioxygen complex and could also affect its reactivity. The ability of substrate to modulate the reactivity of P450 intermediates could be a relevant factor in explaining the remarkable diversity of reactions catalyzed by the enzyme.

Monooxygenation of an aromatic ring by F43W/H64D/V68I myoglobin mutant and hydrogen peroxide. Myoglobin mutants as a model for P450 hydroxylation chemistry

Myoglobin (Mb) is used as a model system for other heme proteins and the reactions they catalyze. The latest novel function to be proposed for myoglobin is a P450 type hydroxylation activity of aromatic carbons (Watanabe, Y., and Ueno, T. (2003) Bull. Chem. Soc. Jpn. 76, 1309-1322). Because Mb does not contain a specific substrate binding site for aromatic compounds near the heme, an engineered tryptophan in the heme pocket was used to model P450 hydroxylation of aromatic compounds. The monooxygenation product was not previously isolated because of rapid subsequent oxidation steps (Hara, I., Ueno, T., Ozaki, S., Itoh, S., Lee, K., Ueyama, N., and Watanabe, Y. (2001) J. Biol. Chem. 276, 36067-36070). In this work, a Mb variant (F43W/H64D/V68I) is used to characterize the monooxygenated intermediate. A modified (+16 Da) species forms upon the addition of 1 eq of H2O2. This product was digested with chymotrypsin, and the modified peptide fragments were isolated and characterized as 6-hydroxytryptophan using matrix-assisted laser desorption ionization time-of-flight tandem mass spectroscopy and 1H NMR. This engineered Mb variant represents the first enzyme to preferentially hydroxylate the indole side chain of Trp at the C6 position. Finally, heme extraction was used to demonstrate that both the formation of the 6-hydroxytryptophan intermediate (+16 Da) and subsequent oxidation to form the +30 Da final product are catalyzed by the heme cofactor, most probably via the compound I intermediate. These results provide insight into the mechanism of hydroxylation of aromatic carbons by heme proteins, demonstrating that non-thiolate-ligated heme enzymes can perform this function. This establishes Mb compound I as a model for P450 type aromatic hydroxylation chemistry.

Compressibility and uncoupling of cytochrome P450cam: high pressure FTIR and activity studies

The effect of the hydrostatic pressure on the CO ligand stretch vibration in cytochrome P450cam-CO bound with various substrates is studied by FTIR. The vibration frequency is linearily shifted to lower values with increasing pressure. The slope of the shift gives the isothermal compressibility of the heme pocket and is found to be related to the high-spin state content in an opposite direction to that previously observed from the pressure-induced shift of the Soret band. This opposite behaviour is explained by the dual effect of heme pocket water molecules both on the CO ligand and on electrostatic potentials produced by the protein at the distal side. The latter effect disturbs ligand-distal side contacts which are needed for a specific proton transfer in oxygen activation when dioxygen is the ligand. Their loss results in uncoupled H(2)O(2) formation.

Crystallographic Studies on the Complex Behavior of Nicotine Binding to P450cam (CYP101)†

Crystallographic and spectroscopic studies have been undertaken to characterize the binding behavior of the non-native substrate nicotine in the active site of the monooxygenase hemoprotein cytochrome P450cam. Despite the existence of a theoretical model that is consistent with the observed distribution of monooxygenation products, the crystal structure of the complex indicates that the primary binding mode of nicotine is unproductive. The structure is confirmed by spectral data that indicate direct coordination of substrate pyridine nitrogen with the heme iron. This would be the proper structure for evaluating binding affinity and inhibition. Reduction of the heme from Fe(III) to Fe(II) and introduction of carbon monoxide into crystals of the nicotine-P450cam complex, to simulate molecular oxygen binding, produces reorientation of the nicotine. This orientation is the appropriate one for predicting regioselectivity and the kinetic features of substrate oxidation. While it is not clear that such complicated behavior will be exhibited for other enzyme-substrate interactions, it is clear that a single crystal structure for a given substrate-enzyme interaction may not provide a good description of the binding mode responsible for product formation.

Crystal structure of human cytochrome P450 2C9 with bound warfarin

Cytochrome P450 proteins (CYP450s) are membrane-associated haem proteins that metabolize physiologically important compounds in many species of microorganisms, plants and animals. Mammalian CYP450s recognize and metabolize diverse xenobiotics such as drug molecules, environmental compounds and pollutants. Human CYP450 proteins CYP1A2, CYP2C9, CYP2C19, CYP2D6 and CYP3A4 are the major drug-metabolizing isoforms, and contribute to the oxidative metabolism of more than 90% of the drugs in current clinical use. Polymorphic variants have also been reported for some CYP450 isoforms, which has implications for the efficacy of drugs in individuals, and for the co-administration of drugs. The molecular basis of drug recognition by human CYP450s, however, has remained elusive. Here we describe the crystal structure of a human CYP450, CYP2C9, both unliganded and in complex with the anti-coagulant drug warfarin. The structure defines unanticipated interactions between CYP2C9 and warfarin, and reveals a new binding pocket. The binding mode of warfarin suggests that CYP2C9 may undergo an allosteric mechanism during its function. The newly discovered binding pocket also suggests that CYP2C9 may simultaneously accommodate multiple ligands during its biological function, and provides a possible molecular basis for understanding complex drug-drug interactions.

Molecular recognition in (+)-alpha-pinene oxidation by cytochrome P450cam

Oxygenated derivatives of the monoterpene (+)-alpha-pinene are found in plant essential oils and used as fragrances and flavorings. (+)-alpha-Pinene is structurally related to (+)-camphor, the natural substrate of the heme monooxygenase cytochrome P450(cam) from Pseudomonas putida. The aim of the present work was to apply the current understanding of P450 substrate binding and catalysis to engineer P450(cam) for the selective oxidation of (+)-alpha-pinene. Consideration of the structures of (+)-camphor and (+)-alpha-pinene lead to active-site mutants containing combinations of the Y96F, F87A, F87L, F87W, and V247L mutations. All mutants showed greatly enhanced binding and rate of oxidation of (+)-alpha-pinene. Some mutants had tighter (+)-alpha-pinene binding than camphor binding by the wild-type. The most active was the Y96F/V247L mutant, with a (+)-alpha-pinene oxidation rate of 270 nmol (nmol of P450(cam))(-)(1) min(-)(1), which was 70% of the rate of camphor oxidation by wild-type P450(cam). Camphor is oxidized by wild-type P450(cam) exclusively to 5-exo-hydroxycamphor. If the gem dimethyl groups of (+)-alpha-pinene occupied similar positions to those found for camphor in the wild-type structure, (+)-cis-verbenol would be the dominant product. All P450(cam) enzymes studied gave (+)-cis-verbenol as the major product but with much reduced selectivity compared to camphor oxidation by the wild-type. (+)-Verbenone, (+)-myrtenol, and the (+)-alpha-pinene epoxides were among the minor products. The crystal structure of the Y96F/F87W/V247L mutant, the most selective of the P450(cam) mutants initially examined, was determined to provide further insight into P450(cam) substrate binding and catalysis. (+)-alpha-Pinene was bound in two orientations which were related by rotation of the molecule. One orientation was similar to that of camphor in the wild-type enzyme while the other was significantly different. Analysis of the enzyme/substrate contacts suggested rationalizations of the product distribution. In particular competition rather than cooperativity between the F87W and V247L mutations and substrate movement during catalysis were proposed to be major factors. The crystal structure lead to the introduction of the L244A mutation to increase the selectivity of pinene oxidation by further biasing the binding orientation toward that of camphor in the wild-type structure. The F87W/Y96F/L244A mutant gave 86% (+)-cis-verbenol and 5% (+)-verbenone. The Y96F/L244A/V247L mutant gave 55% (+)-cis-verbenol but interestingly also 32% (+)-verbenone, suggesting that it may be possible to engineer a P450(cam) mutant that could oxidize (+)-alpha-pinene directly to (+)-verbenone. Verbenol, verbenone, and myrtenol are naturally occurring plant fragrance and flavorings. The preparation of these compounds by selective enzymatic oxidation of (+)-alpha-pinene, which is readily available in large quantities, could have applications in synthesis. The results also show that the protein engineering of P450(cam) for high selectivity of substrate oxidation is more difficult than achieving high substrate turnover rates because of the subtle and dynamic nature of enzyme-substrate interactions.

Specific and non-specific effects of potassium cations on substrate-protein interactions in cytochromes P450cam and P450lin

Substrate binding to cytochrome P450cam is generally considered to be a two-step process. The first step corresponds to the entrance of the substrate, camphor, into the heme pocket. The second step corresponds to a spin transition (low spin-->high spin) of the iron in the protein-substrate complex. This spin transition is related to the mobility of the substrate inside the active site [Biochim Biophys Acta 1338 (1997) 77]. Potassium cations (K(+)) have a specific effect on the spin equilibrium. This is generally attributed to the K(+) ion-induced conformational change of tyrosine 96, the hydroxyl group of which is hydrogen bonded to the keto group of camphor and results in optimum substrate orientation and reduced mobility of this substrate in the active site. In the present paper, we show that K(+) not only affects the substrate-Tyr 96 couple, but acts more globally since K(+) effects are also observed in the Tyr96Phe mutant as well as in complexes with camphor-analogues. Large compounds, that fit well in the heme pocket and bind with higher affinity than camphor, display high spin contents that are less dependent on the presence of K(+). In contrast, K(+) has a significant effect on the high spin content of substrate-cytochrome P450cam complexes with looser interactions. We conclude that large compounds with higher affinities than camphor have more van der Waals contacts with the active site residues. Their mobilities are then reduced and less dependent on the presence of K(+). In this study, we also explored, for comparison, the K(+) effect on the spin transition state of another member of the P450 superfamily, cytochrome P450lin. This effect is not as strong as those observed for cytochrome P450cam. Even though the spin equilibrium does not change dramatically in the presence of K(+) or Na(+), the value of the dissociation constant (K(d)) for linalool binding is significantly affected by ionic strength. Analysis of the thermodynamic parameters for the linalool binding strongly suggests that, similarly to our previous finding for cytochrome P450cam, electrostatic gates participate in the control of substrate access.

Regulation of the properties of the heme-NO complexes in nitric-oxide synthase by hydrogen bonding to the proximal cysteine

Nitric-oxide synthase (NOS) catalyzes the formation of NO and citrulline from l-arginine and oxygen. However, the NO so formed has been found to auto-inhibit the enzymatic activity significantly. We hypothesized that the NO reactivity is in part controlled by hydrogen bonding between the conserved tryptophan residue (position 409 in the neuronal isoform of NOS (nNOS)) and the cysteine residue that forms the proximal bond to the heme. By using resonance Raman spectroscopy and NO as a probe of the heme environment, we show that in the W409F and W409Y mutants of the oxygenase domain of the neuronal enzyme (nNOSox), the Fe-NO bond in the Fe3+NO complex is weaker than in the wild type enzyme, consistent with the loss of a hydrogen bond on the sulfur atom of the proximal cysteine residue. The weaker Fe-NO bond in the W409F and W409Y mutants might result in a faster rate of NO dissociation from the ferric heme in the Trp-409 mutants as compared with the wild type enzyme, which could contribute to the lower accumulation of the inhibitory NO-bound complexes observed during catalysis with the Trp-409 mutants (Adak, S., Crooks, C., Wang, Q., Crane, B. R., Tainer, J. A., Getzoff, E. D., and Stuehr, D. J. (1999) J. Biol. Chem. 274, 26907-26911). The optical and resonance Raman spectra of the Fe2+NO complexes of the Trp-409 mutants differ from those of the wild type enzyme and indicate that a significant population of a five-coordinate Fe2+NO complex is present. These data show that the hydrogen bond provided by the Trp-409 residue is necessary to maintain the thiolate coordination when NO binds to the ferrous heme. Taken together our results indicate that the heme environment on the proximal side of nNOS is critical for the formation of a stable iron-cysteine bond and for the control of the electronic properties of heme-NO complexes.

Conformational relaxation in hemoproteins: the cytochrome P-450cam case

Photodissociation of (CO)P-450(cam)(substrate) complexes was found to trigger a conformational relaxation process that interferes with ligand rebinding at temperatures as low as 140 K even though the protein conformational substates (CS(1)) remain frozen. To analyze the rebinding and relaxation kinetics, we developed a model that takes the distribution of relaxation rates explicitly into account and in which rebinding and relaxation rates are connected by a linear free energy relation. In all complexes heme relaxation occurs first and is probably faster than 100 ns even at 77 K. This is the only process found in substrate-free P-450(cam). Above 140 K and in the presence of a substrate, this initial, fast rebinding state (P) progressively relaxes to another state (P degrees ) in which rebinding is slower. The relaxation rate is independent of solvent rigidity and is governed by the protein's internal dynamics. Rebinding enthalpies in P and P degrees as well as the enthalpy shift brought about by relaxation correlate with the substrate propensity to block access to the iron site. In P degrees the barrier is higher because the substrate is closer to the heme normal and exerts more steric repulsion for CO binding. The relaxation process implies the return of substrate and heme to their ligand-free positions in which access to the heme is reduced.

Insight into protein structure and protein-ligand recognition by Fourier transform infrared spectroscopy

An overview of the application of Fourier transform infrared spectroscopy for the analysis of the structure of proteins and protein-ligand recognition is given. The principle of the technique and of the spectra analysis is demonstrated. Spectral signal assignments to vibrational modes of the peptide chromophore, amino acid side chains, cofactors and metal ligands are summarized. Several examples for protein-ligand recognition are discussed. A particular focus is heme proteins and, as an example, studies of cytochrome P450 are reviewed. Fourier transform infrared spectroscopy in combination with the various techniques such as time-resolved and low-temperature methods, site-directed mutagenesis and isotope labeling is a helpful approach to studying protein-ligand recognition.

Multicopy molecular dynamics simulations suggest how to reconcile crystallographic and product formation data for camphor enantiomers bound to cytochrome P-450cam

Multiple ligand binding modes are possible in many enzyme active sites; their presence in cytochrome P450cam (P450cam) is evident from crystallographic studies of the binding of thiocamphor and phenylimidazoles. Here, we use multicopy molecular dynamics simulations to compare the binding modes of (1R)- and (1S)-camphor in the active site of P450cam. Simulations with (1R)-camphor, the natural substrate, serve to calibrate our protocol: 19 out of 20 copies of (1R)-camphor converged to coordinates very close to those observed for (1R)-camphor in its crystallographic complex with P450cam during the simulations. Simulations with the (1S)-camphor enantiomer showed greater mobility of the substrate, consistent with spectroscopic data, and resulted in 3 major binding modes. One of these is similar to the major conformation (of the two conformations assigned) in a recently determined crystal structure, but this conformation is not correctly oriented for regiospecific hydroxylation at C-5. The simulations, however, provide evidence for reorientation of (1S)-camphor upon formation of the reactive Fe-O intermediate to an orientation suitable for hydroxylation. The simulations thus permit rationalisation of the apparent inconsistency between the crystal structure and the reaction products.

Autocatalytic nitration of P450CAM by peroxynitrite

Peroxynitrite (PN) gains high selectivity as a physiological oxidizing and nitrating agent through catalysis by metal ions. This was established for the heme-thiolate (P450) enzyme prostacyclin synthase which was tyrosine nitrated and inhibited at low PN levels [FEBS Lett. 382 (1996) 101]. Other P450 proteins reacted in a similar manner and a ferryl species (Compound II) has been identified as an intermediate during reactions with PN [Nitric Oxide 3 (1999) 142]. Here we investigated cytochrome P450CAM and found that it catalyzes the decomposition of PN as well as an increased nitration of phenol. The latter at the expense of phenol hydroxylation is characteristic for the proton-assisted PN action. PN also caused self-nitration of P450CAM at several tyrosine residues. Two of them, Y96 and Y305 were largely protected in the presence of the ligand metyrapone. Unlike other heme-thiolate proteins P450CAM did not form distinct spectral intermediates characteristic for Compound II. We conclude that P450CAM serves as a model for the nitration of prostacyclin synthase with respect to its autocatalytic tyrosine nitration and its prevention by blocking the active site.

Hypersensitive radical probes and the mechanisms of cytochrome P450-catalyzed hydroxylation reactions

The title probes are precursors to kinetically calibrated, aryl-substituted cyclopropylcarbinyl radicals that rearrange with picosecond lifetimes. Applications in studies of cytochrome P450-catalyzed hydroxylation reactions are reviewed. Initially confusing results regarding lifetimes of radicals in the hydroxylation reactions were resolved when second-generation probes that distinguish between radicals and cations were employed. The results indicate that two electrophilic oxidizing species are involved in P450-catalyzed hydroxylations, an iron-oxo species that inserts oxygen and a hydroperoxo-iron species that inserts OH(+). The cationic rearrangement products are ascribed to reactions of the protonated alcohol products formed from the latter.

Nitric oxide binding to ferric cytochrome P450: a computational study

The interaction between nitric oxide (NO) and the active site of ferric cytochrome P450 was studied by means of density functional theory (DFT), at the generalized gradient approximation level, and of the SAM1 semiempirical method. The electrostatic effects of the protein environment were included in our DFT scheme by using a hybrid quantum classical approach. The active-site model consisted of an iron(III) porphyrin, the adjacent cysteine residue, and one coordinated water molecule. For this system, spin populations and relative energies for selected spin states were computed. Interestingly, the unpaired electron density, the HOMO, and the LUMO were found to be highly localized on the iron and in an appreciable extent on the sulfur coordinated to the metal. This provides central information about the reactivity of nitric oxide with the active site. Since the substitution of a molecule of H2O by NO has been proposed as being responsible for the inhibition of the cytochrome in the presence of nitric oxide, we have analyzed the thermodynamic feasibility of the ligand exchange process. The structure of the nitrosylated active site was partially optimized using SAM1. A low-spin ground state was obtained for the nitrosyl complex, with a linear Fe-N-O angle. The trends found in Fe-N-O angles and Fe-N lengths of the higher energy spin states provided a notable insight into the electronic configuration of the complex within the framework of the Enemark and Feltham formalism. In relation to the protein environment, it was assessed that the electrostatic field has significant effects on several computed properties. However, in both vacuum and protein environments, the ligand exchange reaction turned out to be exergonic and the relative orders of spin states of the relevant species were the same.