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.     

Substrate binding favors enhanced NO binding to P450cam

Ferric cytochrome P450cam from Pseudomonas putida (P450cam) in buffer solution at physiological pH 7.4 reversibly binds NO to yield the nitrosyl complex P450cam(NO). The presence of 1R-camphor affects the dynamics of NO binding to P450cam and enhances the association and dissociation rate constants significantly. In the case of the substrate-free form of P450cam, subconformers are evident and the NO binding kinetics are much slower than in the presence of the substrate. The association and dissociation processes were investigated by both laser flash photolysis and stopped-flow techniques at ambient and high pressure. Large and positive values of S and V observed for NO binding to and release from the substrate-free P450cam complex are consistent with the operation of a limiting dissociative ligand substitution mechanism, where the lability of coordinated water dominates the reactivity of the iron(III)-heme center with NO. In contrast, NO binding to P450cam in the presence of camphor displays negative activation entropy and activation volume values that support a mechanism dominated by a bond formation process. Volume profiles for the binding of NO appear to be a valuable approach to explain the differences observed for P450cam in the absence and presence of the substrate and enable the clarification of the underlying reaction mechanisms at a molecular level. Changes in spin state of the iron center during the binding/release of NO contribute significantly to the observed volume effects. The results are discussed in terms of relevance for the biological function of cytochrome P450 and in context to other investigations of the related reactions between NO and imidazole- and thiolate-ligated iron(III) hemoproteins.     

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.     

Convergence of Free Energy Profile of Coumarin in Lipid Bilayer

Atomistic molecular dynamics (MD) simulations of druglike molecules embedded in lipid bilayers are of considerable interest as models for drug penetration and positioning in biological membranes. Here we analyze partitioning of coumarin in dioleoylphosphatidylcholine (DOPC) bilayer, based on both multiple, unbiased 3 μs MD simulations (total length) and free energy profiles along the bilayer normal calculated by biased MD simulations (～7 μs in total). The convergences in time of free energy profiles calculated by both umbrella sampling and z-constraint techniques are thoroughly analyzed. Two sets of starting structures are also considered, one from unbiased MD simulation and the other from "pulling" coumarin along the bilayer normal. The structures obtained by pulling simulation contain water defects on the lipid bilayer surface, while those acquired from unbiased simulation have no membrane defects. The free energy profiles converge more rapidly when starting frames from unbiased simulations are used. In addition, z-constraint simulation leads to more rapid convergence than umbrella sampling, due to quicker relaxation of membrane defects. Furthermore, we show that the choice of RESP, PRODRG, or Mulliken charges considerably affects the resulting free energy profile of our model drug along the bilayer normal. We recommend using z-constraint biased MD simulations based on starting geometries acquired from unbiased MD simulations for efficient calculation of convergent free energy profiles of druglike molecules along bilayer normals. The calculation of free energy profile should start with an unbiased simulation, though the polar molecules might need a slow pulling afterward. Results obtained with the recommended simulation protocol agree well with available experimental data for two coumarin derivatives.     

Common system setup for the entire catalytic cycle of cytochrome P450(cam) in quantum mechanical/molecular mechanical studies

We describe a system setup that is applicable to all species in the catalytic cycle of cytochrome P450(cam). The chosen procedure starts from the X-ray coordinates of the ferrous dioxygen complex and follows a protocol that includes the careful assignment of protonation states, comparison between different conceivable hydration schemes, and system preparation through a series of classical minimizations and molecular dynamics (MD) simulations. The resulting setup was validated by quantum mechanical/molecular mechanical (QM/MM) calculations on the resting state, the pentacoordinated ferric and ferrous complexes, Compound I, the transition state and hydroxo intermediate of the C--H hydroxylation reaction, and the product complex. The present QM/MM results are generally consistent with those obtained previously with individual setups. Concerning hydration, we find that saturating the protein interior with water is detrimental and leads to higher structural flexibility and catalytically inefficient active-site geometries. The MD simulations favor a low water density around Asp251 that facilitates side chain rotation of protonated Asp251 during the conversion of Compound 0 to Compound I. The QM/MM results for the two preferred hydration schemes (labeled SE-1 and SE-4) are similar, indicating that slight differences in the solvation close to the active site are not critical as long as camphor and the crystallographic water molecules preserve their positions in the experimental X-ray structures.     

Biotransformation of the sesquiterpene (+)-valencene by cytochrome P450cam and P450BM-3

The sesquiterpenoids are a large class of naturally occurring compounds with biological functions and desirable properties. Oxidation of the sesquiterpene (+)-valencene by wild type and mutants of P450cam from Pseudomonas putida, and of P450BM-3 from Bacillus megaterium, have been investigated as a potential route to (+)-nootkatone, a fine fragrance. Wild type P450cam did not oxidise (+)-valencene but the mutants showed activities up to 9.8 nmol (nmol P450)(-1) min(-1), with (+)-trans-nootkatol and (+)-nootkatone constituting >85% of the products. Wild type P450BM-3 and mutants had higher activities (up to 43 min(-1)) than P450cam but were much less selective. Of the many products, cis- and trans-(+)-nootkatol, (+)-nootkatone, cis-(+)-valencene-1,10-epoxide, trans-(+)-nootkaton-9-ol, and (+)-nootkatone-13S,14-epoxide were isolated from whole-cell reactions and characterised. The selectivity patterns suggest that (+)-valencene has one binding orientation in P450cam but multiple orientations in P450BM-3.     

Force field development for organic molecules: Modifying dihedral and 1–n pair interaction parameters

We comprehensively illustrate a general process of fitting all‐atom molecular mechanics force field (FF) parameters based on quantum mechanical calculations and experimental thermodynamic data. For common organic molecules with free dihedral rotations, this FF format is comprised of the usual bond stretching, angle bending, proper and improper dihedral rotation, and 1–4 scaling pair interactions. An extra format of 1–n scaling pair interaction is introduced when a specific intramolecular rotation is strongly hindered. We detail how the preferred order of fitting all intramolecular FF parameters can be determined by systematically generating characteristic configurations. The intermolecular Van der Waals parameters are initially taken from the literature data but adjusted to obtain a better agreement between the molecular dynamics (MD) simulation results and the experimental observations if necessary. By randomly choosing the molecular configurations from MD simulation and comparing their energies computed from FF parameters and quantum mechanics, the FF parameters can be verified self‐consistently. Using an example of a platform chemical 3‐hydroxypropionic acid, we detail the comparison between the new fitting parameters and the existing FF parameters. In particular, the introduced systematic approach has been applied to obtain the dihedral angle potential and 1–n scaling pair interaction parameters for 48 organic molecules with different functionality. We suggest that this procedure might be used to obtain better dihedral and 1–n interaction potentials when they are not available in the current widely used FF. © 2014 Wiley Periodicals, Inc.     

Protein dynamics in cytochrome P450 molecular recognition and substrate specificity using 2D IR vibrational echo spectroscopy

Cytochrome (cyt) P450s hydroxylate a variety of substrates that can differ widely in their chemical structure. The importance of these enzymes in drug metabolism and other biological processes has motivated the study of the factors that enable their activity on diverse classes of molecules. Protein dynamics have been implicated in cyt P450 substrate specificity. Here, 2D IR vibrational echo spectroscopy is employed to measure the dynamics of cyt P450(cam) from Pseudomonas putida on fast time scales using CO bound at the active site as a vibrational probe. The substrate-free enzyme and the enzyme bound to both its natural substrate, camphor, and a series of related substrates are investigated to explicate the role of dynamics in molecular recognition in cyt P450(cam) and to delineate how the motions may contribute to hydroxylation specificity. In substrate-free cyt P450(cam), three conformational states are populated, and the structural fluctuations within a conformational state are relatively slow. Substrate binding selectively stabilizes one conformational state, and the dynamics become faster. Correlations in the observed dynamics with the specificity of hydroxylation of the substrates, the binding affinity, and the substrates' molecular volume suggest that motions on the hundreds of picosecond time scale contribute to the variation in activity of cyt P450(cam) toward different substrates.     

An automated approach for the parameterization of accurate intermolecular force‐fields: Pyridine as a case study

An automated protocol is proposed and validated, which integrates accurate quantum mechanical calculations with classical numerical simulations. Intermolecular force fields, (FF) suitable for molecular dynamics (MD) and Monte Carlo simulations, are parameterized through a novel iterative approach, fully based on quantum mechanical data, which has been automated and coded into the PICKY software, here presented. The whole procedure is tested and validated for pyridine, whose bulk phase, described through MD simulations performed with the specifically parameterized FF, is characterized by computing several of its thermodynamic, structural, and transport properties, comparing them with their experimental counterparts. © 2011 Wiley Periodicals, Inc.     

Kinetic isotope effects implicate the iron-oxene as the sole oxidant in P450-catalyzed N-dealkylation

Multiple oxidants have been implicated as playing a role in cytochrome P450-mediated oxidations. Herein, we report results on N-dealkylation, one of the most facile reactions mediated by P450 enzymes. We have employed the N-oxides of a series of para-substituted 13C2H2-labeled N,N-dimethylanilines to function as both substrates and surrogate oxygen atom donors for P450cam and P4502E1. Kinetic isotope effect profiles obtained using the N-oxide system were found to closely match the profiles produced using the complete NAD(P)H/NAD(P)-P450 reductase/O2 system. The results are consistent with oxidation occurring solely through an iron-oxene species.     

Structure and thermodynamics of alpha-, beta-, and gamma-cyclodextrin dimers. Molecular dynamics studies of the solvent effect and free binding energies

The alpha-, beta-, and gamma-cyclodextrin (CyDs) dimers were studied by molecular dynamics (MD) simulations in water as an explicit solvent. The relative stability of dimers and the involved molecular interactions were determined. Three possible starting orientations were considered for the dimers: head-to-head, head-to-tail, and tail-to-tail. MD simulations were performed over a period of 5 ns to ensure the stability of the system for both the CyD dimers and monomers. The MM-PBSA methodology was used to obtain the free binding energy of the dimers and to determine the most stable arrangement for each solvated CyD. In a vacuum, MD simulations provided the head-to-head orientation as the most stable orientation for the three CyDs, while in aqueous solution the, the head-to-tail orientation was found to be the most stable for the alpha-CyD dimer and the tail-to-tail orientation the most stable for the beta- and gamma-CyD dimers.     

Compound I in heme thiolate enzymes: a comparative QM/MM study

This study directly compares the active species of heme enzymes, so-called Compound I (Cpd I), across the heme-thiolate enzyme family. Thus, sixty-four different Cpd I structures are calculated by hybrid quantum mechanical/molecular mechanical (QM/MM) methods using four different cysteine-ligated heme enzymes (P450(cam), the mutant P450(cam)-L358P, CPO and NOS) with varying QM region sizes in two multiplicities each. The overall result is that these Cpd I species are similar to each other with regard to many characteristic features. Hence, using the more stable CPO Cpd I as a model for P450 Cpd I in experiments should be a reasonable approach. However, systematic differences were also observed, and it is shown that NOS stands out in most comparisons. By analyzing the electrical field generated by the enzyme on the QM region, one can see that (a) the protein exerts a large influence and modifies all the Cpd I species compared with the gas-phase situation and (b) in NOS this field is approximately planar to the heme plane, whereas it is approximately perpendicular in the other enzymes, explaining the deviating results on NOS. The calculations on the P450(cam) mutant L358P show that the effects of removing the hydrogen bond between the heme sulfur and L358 are small at the Cpd I stage. Finally, Mossbauer parameters are calculated for the different Cpd I species, enabling future comparisons with experiments. These results are discussed in the broader context of recent findings of Cpd I species that exhibit large variations in the electronic structure due to the presence of the substrate.     

Complexes of a Zn‐metalloenzyme binding site with hydroxamate‐containing ligands. A case for detailed benchmarkings of polarizable molecular mechanics/dynamics potentials when the experimental binding structure is unknown

Zn‐metalloproteins are a major class of targets for drug design. They constitute a demanding testing ground for polarizable molecular mechanics/dynamics aimed at extending the realm of quantum chemistry (QC) to very long‐duration molecular dynamics (MD). The reliability of such procedures needs to be demonstrated upon comparing the relative stabilities of competing candidate complexes of inhibitors with the recognition site stabilized in the course of MD. This could be necessary when no information is available regarding the experimental structure of the inhibitor–protein complex. Thus, this study bears on the phosphomannose isomerase (PMI) enzyme, considered as a potential therapeutic target for the treatment of several bacterial and parasitic diseases. We consider its complexes with 5‐phospho‐d ‐arabinonohydroxamate and three analog ligands differing by the number and location of their hydroxyl groups. We evaluate the energy accuracy expectable from a polarizable molecular mechanics procedure, SIBFA. This is done by comparisons with ab initio quantum‐chemistry (QC) calculations in the following cases: (a) the complexes of the four ligands in three distinct structures extracted from the entire PMI‐ligand energy‐minimized structures, and totaling up to 264 atoms; (b) the solvation energies of several energy‐minimized complexes of each ligand with a shell of 64 water molecules; (c) the conformational energy differences of each ligand in different conformations characterized in the course of energy‐minimizations; and (d) the continuum solvation energies of the ligands in different conformations. The agreements with the QC results appear convincing. On these bases, we discuss the prospects of applying the procedure to ligand‐macromolecule recognition problems. © 2016 Wiley Periodicals, Inc.     

Fluorescent probes for cytochrome p450 structural characterization and inhibitor screening

We have synthesized two luminescent probes (D-4-Ad and D-8-Ad) that target cytochrome P450cam. D-4-Ad luminescence is quenched by F?rster energy transfer upon binding (Kd = 0.83 muM) but is restored when the probe is displaced from the active site by camphor. In contrast, D-8-Ad (Kd approximately 0.02 muM) is not displaced from the enzyme, even in the presence of a large excess of camphor. The 2.2 A resolution crystal structure of the D-8-Ad:P450cam complex reveals extensive hydrophobic contacts between the probe and the enzyme, which result from the conformational flexibility of the B', F, and G helices. Probes with properties similar to those of D-4-Ad potentially could be useful for screening P450 inhibitors.     

Towards the Accurate Thermodynamic Characterization of Enzyme Reaction Mechanisms

We employed QM/MM molecular dynamics (MD) simulations to characterize the rate-limiting step of the glycosylation reaction of pancreatic α-amylase with combined DFT/molecular dynamics methods (PBE/def2-SVP : AMBER). Upon careful choice of four starting active site conformations based on thorough reactivity criteria, Gibbs energy profiles were calculated with umbrella sampling simulations within a statistical convergence of 1–2 kcal ⋅ mol⁻¹. Nevertheless, Gibbs activation barriers and reaction energies still varied from 11.0 to 16.8 kcal ⋅ mol⁻¹ and −6.3 to +3.8 kcal ⋅ mol⁻¹ depending on the starting conformations, showing that despite significant state-of-the-art QM/MM MD sampling (0.5 ns/profile) the result still depends on the starting structure. The results supported the one step dissociative mechanism of Asp197 glycosylation preceded by an acid-base reaction by the Glu233, which are qualitatively similar to those from multi-PES QM/MM studies, and thus support the use of the latter to determine enzyme reaction mechanisms.     

Thermodynamic and kinetic studies on the binding of nitric oxide to a new enzyme mimic of cytochrome p450

A new model for the P450 enzyme carrying a SO(3)(-) ligand coordinated to iron(III) (complex 2) reversibly binds NO to yield the nitrosyl adduct. The rate constant for NO binding to 2 in toluene is of the same order of magnitude as that found for the nitrosylation of the native, substrate-bound form of P450(cam) (E.S-P450(cam)). Large and negative activation entropy and activation volume values for the binding of NO to complex 2 support a mechanism that is dominated by bond formation with concomitant iron spin change from S = (5)/(2) to S = 0, as proposed for the reaction between NO and E.S-P450(cam). In contrast, the dissociation of NO from 2(NO) was found to be several orders of magnitude faster than the corresponding reaction for the E.S-P450(cam)/NO system. In a coordinating solvent such as methanol, the alcohol coordinates to iron(III) of 2 at the distal position, generating a six-coordinate, high-spin species 5. The reaction of NO with 5 in methanol was found to be much slower in comparison to the nitrosylation reaction of 2 in toluene. This behavior can be explained in terms of a mechanism in which methanol must be displaced during Fe-NO bond formation. The thermodynamic and kinetic data for NO binding to the new model complexes of P450 (2 and 5) are discussed in reference to earlier results obtained for closely related nitrosylation reactions of cytochrome P450(cam) (in the presence and in the absence of the substrate) and a thiolate-ligated iron(III) model complex.     

Understanding the pH-dependent behavior of graphene oxide aqueous solutions: a comparative experimental and molecular dynamics simulation study

Understanding the pH-dependent behavior of graphene oxide (GO) aqueous solutions is important to the production of assembled GO or reduced GO films for electronic, optical, and biological applications. We have carried out a comparative experimental and molecular dynamics (MD) simulation study to uncover the mechanisms behind the aggregation and the surface activity of GO at different pH values. At low pH, the carboxyl groups are protonated such that the GO sheets become less hydrophilic and form aggregates. MD simulations further suggest that the aggregates exhibit a GO-water-GO sandwichlike structure and as a result are stable in water instead of precipitating. However, at high pH, the deprotonated carboxyl groups are very hydrophilic such that individual GO sheets prefer to dissolve in bulk water like a regular salt. The GO aggregates formed at low pH are found to be surface-active and do not exhibit characteristic features of surfactant micelles. Our findings suggest that GO does not behave like conventional surfactants in pH 1 and 14 aqueous solutions. The molecular-level understanding of the solution behavior of GO presented here can facilitate and improve the experimental techniques used to synthesize and sort large, uniform GO dispersions in a solution phase.     

3D-QSAR methods on the basis of ligand–receptor complexes. Application of COMBINE and GRID/GOLPE methodologies to a series of CYP1A2 ligands

Many heterocyclic amines (HCA) present in cooked food exert a genotoxic activity when they are metabolised (N-oxidated) by the human cytochrome P450 1A2 (CYP1A2h). In order to rationalize the observed differences in activity of this enzyme on a series of 12 HCA, 3D-QSAR methods were applied on the basis of models of HCA–CYP1A2h complexes. The CYP1A2h enzyme model has been previously reported and was built by homology modeling based on cytochrome P450 BM3. The complexes were automatically generated applying the AUTODOCK software and refined using AMBER. A COMBINE analysis on the complexes identified the most important enzyme–ligand interactions that account for the differences in activity within the series. A GRID/GOLPE analysis was then performed on just the ligands, in the conformations and orientations found in the modeled complexes. The results from both methods were concordant and confirmed the advantages of incorporating structural information from series of ligand–receptor complexes into 3D-QSAR methodologies.     

Roles of the proximal heme thiolate ligand in cytochrome p450(cam).

To examine the roles of the proximal thiolate iron ligand, the C357H mutant of P450(cam) (CYP101) was characterized by resonance Raman, UV, circular dichroism, and activity measurements. The C357H mutant must be reconstituted with hemin for activity to be observed. The reconstituted enzyme is a mixture of high and low spin species. Low temperature (10 degrees C), low enzyme concentration (1 microM), high camphor concentration (1 mM), and 5--50 mM buffer concentrations increase the high to low spin ratio, but under no conditions examined was the protein more than 60% high spin. The C357H mutant has a poorer K(m) for camphor (23 vs 2 microM) and a poorer K(d) for putidaredoxin (50 vs 20 microM) than wild-type P450(cam). The mutant also exhibits a greatly decreased camphor oxidation rate, elevated uncoupling rate, and much greater peroxidase activity. Electron transfer from putidaredoxin to the mutant is much slower than to the wild-type even though redox potential measurements show that the electron transfer remains thermodynamically favored. These experiments confirm that the thiolate ligand facilitates the O--O bond cleavage by P450 enzymes and also demonstrate that this ligand satisfies important roles in protein folding, substrate binding, and electron transfer.