Proximal Pocket Controls Alkene Oxidation Selectivity of Cytochrome P450 and Chloroperoxidase toward Small, Nonpolar Substrates.

This paper examines the influence of the proximal pockets of cytochrome P450CAM and chloroperoxidase (CPO) on the relative favorability of catalytic epoxidation and allylic hydroxylation of olefins, a type of alkene oxidation selectivity. The study employs quantum mechanical models of the active site to isolate the proximal pocket's influence on the barrier for the selectivity-determining step for each reaction, using cyclohexene and cis-ß-methylstyrene as substrates. The proximal pocket is found to preference epoxidation by 2-5 kcal/mol, the largest value being for CPO, converting the active heme-thiolate moiety from being intrinsically hydroxylation-selective to being intrinsically epoxidation-selective. This theoretical study, the first to correctly predict these enzymes' preference for epoxidation of allylic substrates, strongly suggests that the proximal pocket is the key determinant of alkene oxidation selectivity. The selectivity for epoxidation can be rationalized in terms of the proximal pocket's modulation of the thiolate's electron "push" and consequent influence on the heme redox potential and the basicity of the trans ligand.

How the Proximal Pocket May Influence the Enantiospecificities of Chloroperoxidase-Catalyzed Epoxidations of Olefins.

Chloroperoxidase-catalyzed enantiospecific epoxidations of olefins are of significant biotechnological interest. Typical enantiomeric excesses are in the range of 66%-97% and translate into free energy differences on the order of 1 kcal/mol. These differences are generally attributed to the effect of the distal pocket. In this paper, we show that the influence of the proximal pocket on the electron transfer mechanism in the rate-limiting event may be just as significant for a quantitatively accurate account of the experimentally-measured enantiospecificities.

Chloroperoxidase-Catalyzed Epoxidation of Cis-ß-Methylstyrene: NH-S Hydrogen Bonds and Proximal Helix Dipole Change the Catalytic Mechanism and Significantly Lower the Reaction Barrier.

Proximal hydrogen bonding of the axial sulfur with the backbone amides (NH-S) is a conserved feature of heme-thiolate enzymes such as chloroperoxidase (CPO) and cytochrome P450 (P450). In CPO, the effect of NH-S bonds is amplified by the dipole moment of the proximal helix. Our gas-phase DFT studies show that the proximal pocket effect significantly enhances CPO's reactivity toward the epoxidation of olefinic substrates. Comparison of models with and without proximal pocket residues shows that with them, the barrier for Cß-O bond formation is lowered by about ∼4.6 kcal/mol, while Cα-O-Cß ring closure becomes barrierless. The dipole moment of the proximal helix was estimated to contribute 1/3 of the decrease, while the rest is attributed to the effect of NH-S bonds. The decrease of the reaction barrier correlates with increased electron density transfer to residues of the proximal pocket. The effect is most pronounced on the doublet spin surface and involves a change in the electron-transfer mechanism. A full enzyme QMMM study on the doublet spin surface gives about the same barrier as the gas-phase DFT study. The free-energy barrier was estimated to be in agreement with the experimental results for the CPO-catalyzed epoxidation of styrene.

Proximal Pocket Hydrogen Bonds Significantly Influence the Mechanism of Chloroperoxidase Compound I Formation.

The influence of backbone hydrogen bonds to the sulfur atom of the proximal thiolate (NH···S hydrogen bonds) on the formation of compound I in chloroperoxidase is investigated with DFT calculations. Reaction profiles for the transformation of the ferric resting state into compound I in the presence of a peroxide substrate are calculated for a model system incorporating the heme and key proximal and distal amino acid residues. We find that NH···S hydrogen bonds (1) reduce the barrier for the formation of compound 0 by 7.6 kcal/mol, (2) increase the stability of compound 0 by 5.2 kcal/mol, (3) reduce the stability of compound I relative to compound 0 by 6.2 kcal/mol, and (4) reduce the stability of protonated compound 0, favoring a hybrid homo-heterolytic relative to a classic heterolytic mechanism for O-O bond scission. In general, the influence of the NH···S hydrogen bonds can be traced to a reduction in the pKa of the heme-bound substrate. We find that the hydrogen bond networks on the proximal and distal sides of the heme function together to modulate the mechanism of reaction. These results confirm and extend long-standing theories that the NH···S hydrogen bonds in heme thiolate proteins influence reactivity by tuning the thiolate "push" effect.

A possible mechanism for redox control of human neuroglobin activity.

Neuroglobin (Ngb) promotes neuron survival under hypoxic/ischemic conditions. In vivo and in vitro assays provide evidence for redox-regulated functioning of Ngb. On the basis of X-ray crystal structures and our MD simulations, a mechanism for redox control of human Ngb (hNgb) activity via the influence of the CD loop on the active site is proposed. We provide evidence that the CD loop undergoes a strand-to-helix transition when the external environment becomes sufficiently oxidizing, and that this CD loop conformational transition causes critical restructuring of the active site. We postulate that the strand-to-helix mechanics of the CD loop allows hNgb to utilize the lability of Cys46/Cys55 disulfide bonding and of the Tyr44/His64/heme propionate interaction network for redox-controlled functioning of hNgb.

Theoretical study of HOCl-catalyzed keto-enol tautomerization of ß-cyclopentanedione in an explicit water environment.

The mechanism of acid-catalyzed keto-enol tautomerization of ß-cyclopentanedione (CPD) in solution is studied computationally. Reaction profiles are first calculated for a limited solvation environment using ab initio and density functional methods. Barrier heights for systems including up to one hydration shell of explicit water molecules depend strongly on the number of waters involved in proton transfer and to a lesser but significant extent on the number of waters forming hydrogen bonds with waters in the proton-transfer chain (each such water reduces the barrier by 4.4 kcal/mol on average). Barriers of 8-13 kcal/mol were obtained when a full or nearly full hydration shell was present, consistent with calculations for nonacid-catalyzed keto-enol tautomerization of related molecules. The presence of HOCl reduced the barrier by 4.5 kcal/mol in relation to the gas phase, consistent with the well-known principle that keto-enol tautomerization can be acid- or base-catalyzed. The reaction was also modeled beginning with snapshots of reactant conformations taken from a 300 K molecular dynamics simulation of CPD, HOCl, and 324 explicit waters. Reaction profiles were calculated at a QM/MM level with waters in the first hydration shell either fixed or energy-minimized at each step along the reaction coordinate. A substantial variation in barrier height was observed in both cases, depending primarily on electrostatic interactions (hydrogen bonding) with first-hydration-shell waters and, to a lesser extent, on electrostatic interactions with more distant waters and geometric distortion effects. For the lowest barriers, the extent of barrier reduction by waters involved in proton transfer is consistent with the limited solvation results, but further barrier reduction due to hydrogen bonding to waters involved in proton transfer is not observed. It is postulated that this is because highly flexible structures such as extensive hydrogen bonding networks optimal for reaction are entropically disfavored and so may not contribute significantly to the observed reaction rate.

Paramagnetic nuclear magnetic resonance relaxation and molecular mechanics studies of the chloroperoxidase-indole complex: insights into the mechanism of chloroperoxidase-catalyzed regioselective oxidation of indole.

To unravel the mechanism of chloroperoxidase (CPO)-catalyzed regioselective oxidation of indole, we studied the structure of the CPO-indole complex using nuclear magnetic resonance (NMR) relaxation measurements and computational techniques. The dissociation constant (KD) of the CPO-indole complex was calculated to be approximately 21 mM. The distances (r) between protons of indole and the heme iron calculated via NMR relaxation measurements and molecular docking revealed that the pyrrole ring of indole is oriented toward the heme with its 2-H pointing directly at the heme iron. Both KD and r values are independent of pH in the range of 3.0-6.5. The stability and structure of the CPO-indole complex are also independent of the concentration of chloride or iodide ion. Molecular docking suggests the formation of a hydrogen bond between the NH group of indole and the carboxyl O of Glu 183 in the binding of indole to CPO. Simulated annealing of the CPO-indole complex using r values from NMR experiments as distance restraints reveals that the van der Waals interactions were much stronger than the Coulomb interactions in the binding of indole to CPO, indicating that the association of indole with CPO is primarily governed by hydrophobic rather than electrostatic interactions. This work provides the first experimental and theoretical evidence of the long-sought mechanism that leads to the "unexpected" regioselectivity of the CPO-catalyzed oxidation of indole. The structure of the CPO-indole complex will serve as a lighthouse in guiding the design of CPO mutants with tailor-made activities for biotechnological applications.

Chloroperoxidase-catalyzed epoxidation of cis-ß-methylstyrene: distal pocket flexibility tunes catalytic reactivity.

Chloroperoxidase, the most versatile heme protein, has a hybrid active site pocket that shares structural features with peroxidases and cytochrome P450s. The simulation studies presented here show that the enzyme possesses a remarkable ability to efficiently utilize its hybrid structure, assuming structurally different peroxidase-like and P450-like distal pocket faces and thereby enhancing the inherent catalytic capability of the active center. We find that, during epoxidation of cis-ß-methylstyrene (CBMS), the native peroxidase-like aspect of the distal pocket is diminished as the polar Glu183 side chain is displaced away from the active center and the distal pocket takes on a more hydrophobic, P450-like, aspect. The P450-like distal pocket provides a significant enthalpic stabilization of ∼4 kcal/mol of the 14 kcal/mol reaction barrier for gas-phase epoxidation of CMBS by an oxyferryl heme-thiolate species. This stabilization comes from breathing of the distal pocket. As until recently the active site of chloroperoxidase was postulated to be inflexible, these results suggest a new conceptual understanding of the enzyme's versatility: catalytic reactivity is tuned by flexibility of the distal pocket.

Enantiospecificity of chloroperoxidase-catalyzed epoxidation: biased molecular dynamics study of a cis-ß-methylstyrene/chloroperoxidase-compound I complex.

Molecular dynamics simulations of an explicitly solvated cis-ß-methylstyrene/chloroperoxidase-Compound I complex are performed to determine the cause of the high enantiospecificity of epoxidation. From the simulations, a two-dimensional free energy potential is calculated to distinguish binding potential wells from which reaction to 1S2R and 1R2S epoxide products may occur. Convergence of the free energy potential is accelerated with an adaptive biasing potential. Analysis of binding is followed by analysis of 1S2R and 1R2S reaction precursor structures in which the substrate, having left the binding wells, places its reactive double bond in steric proximity to the oxyferryl heme center. Structural analysis of binding and reaction precursor conformations is presented. We find that 1), a distortion of Glu(183) is important for CPO-catalyzed epoxidation as was postulated previously based on experimental results; 2), the free energy of binding does not provide significant differentiation between structures leading to the respective epoxide enantiomers; and 3), CPO's enantiospecificity toward cis-ß-methylstyrene is likely to be caused by a specific group of residues which form a hydrophobic core surrounding the oxyferryl heme center.

The conformation of end-groups is one determinant of carotenoid topology suitable for high fidelity molecular recognition: a study of beta- and epsilon-end-groups.

Conformation affects a carotenoid's ability to bind selectively to proteins. We calculated adiabatic energy profiles for rotating the ring end-groups around the C6C7 bond and for flexing of the ring with respect to the polyene chain. The choice of computational methods is important. A low, 4.2 kcal/mol barrier to rotation exists for a beta-ring. An 8.3 kcal/mol barrier exists for rotation of an epsilon-ring. Rotation of the epsilon-ring is sensitive to substitution at C3. In the absence of external forces neither beta- nor epsilon-rings are rotationally constrained. The nearly parallel alignment of the beta-ring to the C6C7 bond axis contrasts to the more perpendicular orientation of the epsilon-ring. Flexion of a beta-ring to the minimized epsilon-ring conformation requires approximately 23 kcal/mol; extension of the epsilon-ring to the minimized beta-ring conformation requires approximately 8 kcal/mol. Selectivity associated with beta- versus epsilon-rings is dominated by the inability of the beta-ring to flex to minimize protein/ring steric interactions and maximize van der Waal's attractions with the binding site.

1,3-dipolar cycloadditions of trimethylsilyldiazomethane revisited: steric demand of the dipolarophile and the influence on product distribution.

The 1,3-dipolar cycloaddition of trimethylsilyldiazomethane with alpha,beta-unsaturated esters was examined. The resulting 1-pyrazolines isomerize to regioisomeric 2-pyrazolines (a or b) or undergo desilylation (c). Acrylates yield only b or c. beta-Substituted dipolarophiles may yield all three types of products. This work demonstrates that the distribution of 2-pyrazoline products is highly dependent on the relative configuration of the substituents on the 1-pyrazoline intermediate.

Correlation times and adiabatic barriers for methyl rotation in SNase.

The relation of rotational correlation times to adiabatic rotational barriers for alanine methyl groups in staphylococcal nuclease (SNase) is investigated. The hypothesis that methyl rotational barriers may be useful probes of local packing in proteins is supported by an analysis of ten X-ray crystal structures of SNase mutants. The barrier heights are consistent across a set of ten structures of a native SNase and mutants containing single-point mutations or single or double insertions, most in a ternary SNase complex. The barriers for different methyls have a range of 7.5 kcal/mol, which at 300 K would correspond to a five-order-of-magnitude range in correlation time. It is demonstrated that adiabatic rotational barriers can fluctuate significantly during an MD simulation of hydrated SNase, but that a Boltzmann weighted average is predictive of rotational correlation times determined from correlation functions. Even if a given methyl is on average quite sterically hindered, infrequently sampled low-barrier conformations may dominate the Boltzmann distribution. This result is consistent with the observed uniformity of NMR correlation times for (13)C-labeled methyls. The methyl barriers in simulation fluctuate on multiple time scales, which can make the precise relationship between methyl rotational correlation time and methyl rotation barriers complicated. The implications of these issues for the interpretation of correlation times determined from NMR and simulation are discussed.

Methyl dynamics in crystalline amino acids: MD and NMR.

Correlation times for rotation of deuterated methyls in crystalline leucine, valine, and cyclo-L-alanyl-L-alanine are calculated with molecular dynamics and compared with NMR data. The simulations distinguish between methyls having different steric environments in the crystal, yielding correlation times differing by a factor of up to 30 for methyls within a given crystal. MD and NMR correlation times agree to within a factor of 2. However, averaging over nonequivalent methyls can yield correlation functions that, although actually multiexponential, are well fit by single exponentials. This may have significance for interpreting NMR data; previous NMR data did not distinguish between the methyls in these crystals. Adiabatic rotational barriers calculated with the X-ray structure differ from effective barriers during simulation by up to +/-1 kcal/mol; the difference indicates that dynamical effects have a significant role in determining rotational correlation times.