P. Ross Wilderman

Identification of syn-pimara-7,15-diene synthase reveals functional clustering of terpene synthases involved in rice phytoalexin/allelochemical biosynthesis.

Rice (Oryza sativa) produces momilactone diterpenoids as both phytoalexins and allelochemicals. Accordingly, the committed step in biosynthesis of these natural products is catalyzed by the class I terpene synthase that converts syn-copalyl diphosphate to the corresponding polycyclic hydrocarbon intermediate syn-pimara-7,15-diene. Here, a functional genomics approach was utilized to identify a syn-copalyl diphosphate specific 9β-pimara-7,15-diene synthase (OsDTS2). To our knowledge, this is the first identified terpene synthase with this particular substrate stereoselectivity and, by comparison with the previously described and closely related ent-copalyl diphosphate specific cassa-12,15-diene synthase (OsDTC1), provides a model system for investigating the enzymatic determinants underlying the observed difference in substrate specificity. Further, OsDTS2 mRNA in leaves is up-regulated by conditions that stimulate phytoalexin biosynthesis but is constitutively expressed in roots, where momilactones are constantly synthesized as allelochemicals. Therefore, transcription of OsDTS2 seems to be an important regulatory point for controlling production of these defensive compounds. Finally, the gene identified here as OsDTS2 has previously been mapped at 14.3 cM on chromosome 4. The class II terpene synthase producing syn-copalyl diphosphate from the universal diterpenoid precursor geranylgeranyl diphosphate was also mapped to this same region. These genes catalyze sequential cyclization steps in momilactone biosynthesis and seem to have been evolutionarily coupled by physical linkage and resulting cosegregation. Further, the observed correlation between physical proximity and common metabolic function indicates that other such class I and class II terpene synthase gene clusters may similarly catalyze consecutive reactions in shared biosynthetic pathways.

Rice Contains Two Disparate ent -Copalyl Diphosphate Synthases with Distinct Metabolic Functions

Rice (Oryza sativa) produces ent-copalyl diphosphate for both gibberellin (GA) phytohormone and defensive phytoalexin biosynthesis, raising the question of how this initial biosynthetic step is carried out for these distinct metabolic processes. Here, a functional genomics approach has been utilized to identify two disparate ent-copalyl diphosphate synthases from rice (OsCPS1ent and OsCPS2ent). Notably, it was very recently demonstrated that only one of these (OsCPS1ent) normally operates in GA biosynthesis as mutations in this gene result in severely impaired growth. Evidence is presented here strongly indicating that the other (OsCPS2ent) is involved in related secondary metabolism producing defensive phytochemicals. In particular, under appropriate conditions, OsCPS2ent mRNA is specifically induced in leaves prior to production of the corresponding phytoalexins. Thus, transcriptional control of OsCPS2ent seems to be an important means of regulating defensive phytochemical biosynthesis. Finally, OsCPS1ent is significantly more similar to the likewise GA-specific gene An1/ZmCPS1ent in maize (Zea mays) than its class II terpene synthase paralogs involved in rice secondary metabolism. Hence, we speculate that this cross-species conservation by biosynthetic process reflects derivation of related secondary metabolism from the GA primary biosynthetic pathway prior to the early divergence between the separate lineages within the cereal/grass family (Poaceae) resulting in modern rice and maize.

Following evolution's lead to a single residue switch for diterpene synthase product outcome

There have been few insights into the biochemical origins of natural product biosynthesis from primary metabolism. Of particular interest are terpene synthases, which often mediate the committed step in particular biosynthetic pathways so that alteration of their product outcome is a key step in the derivation of novel natural products. These enzymes also catalyze complex reactions of significant mechanistic interest. Following an evolutionary lead from two recently diverged, functionally distinct diterpene synthase orthologs from different subspecies of rice, we have identified a single residue that can switch product outcome. Specifically, the mutation of a conserved isoleucine to threonine that acts to convert not only the originally targeted isokaurene synthase into a specific pimaradiene synthase but also has a much broader effect, which includes conversion of the ent-kaurene synthases found in all higher plants for gibberellin phytohormone biosynthesis to the production of pimaradiene. This surprisingly facile switch for diterpene synthase catalytic specificity indicates the ease with which primary (gibberellin) metabolism can be subverted to secondary biosynthesis and may underlie the widespread occurrence of pimaradiene-derived natural products. In addition, because this isoleucine is required for the mechanistically more complex cyclization to tetracyclic kaurene, whereas substitution with threonine “short-circuits” this mechanism to produce the “simpler” tricyclic pimaradiene, our results have some implications regarding the means by which terpene synthases specify product outcome.

Crystal structure of a cytochrome P450 2B6 genetic variant in complex with the inhibitor 4-(4-chlorophenyl)imidazole at 2.0-A resolution.

The structure of the K262R genetic variant of human cytochrome P450 2B6 in complex with the inhibitor 4-(4-chlorophenyl)imidazole (4-CPI) has been determined using X-ray crystallography to 2.0-Å resolution. Production of diffraction quality crystals was enabled through a combination of protein engineering, chaperone coexpression, modifications to the purification protocol, and the use of unique facial amphiphiles during crystallization. The 2B6-4-CPI complex is virtually identical to the rabbit 2B4 structure bound to the same inhibitor with respect to the arrangement of secondary structural elements and the placement of active site residues. The structure supports prior P450 2B6 homology models based on other mammalian cytochromes P450 and is consistent with the limited site-directed mutagenesis studies on 2B6 and extensive studies on P450 2B4 and 2B1. Although the K262R genetic variant shows unaltered binding of 4-CPI, altered binding affinity, kinetics, and/or product profiles have been previously shown with several other ligands. On the basis of new P450 2B6 crystal structure and previous 2B4 structures, substitutions at residue 262 affect a hydrogen-bonding network connecting the G and H helices, where subtle differences could be transduced to the active site. Docking experiments indicate that the closed protein conformation allows smaller ligands such as ticlopidine to bind to the 2B6 active site in the expected orientation. However, it is unknown whether 2B6 undergoes structural reorganization to accommodate bulkier molecules, as previously inferred from multiple P450 2B4 crystal structures.

Plasticity of Cytochrome P450 2B4 as Investigated by Hydrogen-Deuterium Exchange Mass Spectrometry and X-ray Crystallography

Crystal structures of the xenobiotic metabolizing cytochrome P450 2B4 have demonstrated markedly different conformations in the presence of imidazole inhibitors or in the absence of ligand. However, knowledge of the plasticity of the enzyme in solution has remained scant. Thus, hydrogen-deuterium exchange mass spectrometry (DXMS) was utilized to probe the conformations of ligand-free P450 2B4 and the complex with 4-(4-chlorophenyl)imidazole (4-CPI) or 1-biphenyl-4-methyl-1H-imidazole (1-PBI). The results of DXMS indicate that the binding of 4-CPI slowed the hydrogen-deuterium exchange rate over the B′- and C-helices and portions of the F-G-helix cassette compared with P450 2B4 in the absence of ligands. In contrast, there was little difference between the ligand-free and 1-PBI-bound exchange sets. In addition, DXMS suggests that the ligand-free P450 2B4 is predominantly open in solution. Interestingly, a new high resolution structure of ligand-free P450 2B4 was obtained in a closed conformation very similar to the 4-CPI complex. Molecular dynamics simulations performed with the closed ligand-free structure as the starting point were used to probe the energetically accessible conformations of P450 2B4. The simulations were found to equilibrate to a conformation resembling the 1-PBI-bound P450 2B4 crystal structure. The results indicate that conformational changes observed in available crystal structures of the promiscuous xenobiotic metabolizing cytochrome P450 2B4 are consistent with its solution structural behavior.

Plasticity of CYP2B Enzymes: Structural and Solution Biophysical Methods

In the past three years, major advances in understanding cytochrome P450 2B (CYP2B) structure-function relationships have been made through determination of multiple ligand-bound and one ligand-free X-ray crystal structure of CYP2B4 and one ligand-bound X-ray crystal structure of CYP2B6. These structures have provided insight into the features that provide the high degree of plasticity of the enzymes. A combination of a phenylalanine cluster that allows for concerted movement of helices F through G and a conserved set of electrostatic interactions involving Arg(262) facilitates movement of this region to accommodate binding of ligands of various sizes without perturbing most of the P450 fold. Integrating solution based techniques such as NMR or deuterium exchange mass spectrometry (DXMS) with computational methods including molecular docking has provided further insight into enzyme behavior upon ligand binding. In addition, extended molecular dynamics simulations have provided a link between an open and a closed conformation of ligand-free CYP2B4 found in crystal structures. Other studies revealed the utility of rational engineering in improving stability of P450s to facilitate structural studies. The solution and computational results combined with the X-ray crystal structures yield a comprehensive picture of how these enzymes adopt different conformations to bind various ligands.

Investigation by site-directed mutagenesis of the role of cytochrome P450 2B4 non-active-site residues in protein–ligand interactions based on crystal structures of the ligand-bound enzyme

Residues located outside the active site of cytochromes P450 2B have exhibited importance in ligand binding, structural stability and drug metabolism. However, contributions of non‐active‐site residues to the plasticity of these enzymes are not known. Thus, a systematic investigation was undertaken of unique residue–residue interactions found in crystal structures of P450 2B4 in complex with 4‐(4‐chlorophenyl)imidazole (4‐CPI), a closed conformation, or in complex with bifonazole, an expanded conformation. Nineteen mutants distributed over 11 sites were constructed, expressed in Escherichia coli and purified. Most mutants showed significantly decreased expression, especially in the case of interactions found in the 4‐CPI structure. Six mutants (H172A, H172F, H172Q, L437A, E474D and E474Q) were chosen for detailed functional analysis. Among these, the Ks of H172F for bifonazole was ∼ 20 times higher than for wild‐type 2B4, and the Ks of L437A for 4‐CPI was ∼ 50 times higher than for wild‐type, leading to significantly altered inhibitor selectivity. Enzyme function was tested with the substrates 7‐ethoxy‐4‐(trifluoromethyl)coumarin, 7‐methoxy‐4‐(trifluoromethyl)coumarin and 7‐benzyloxyresorufin (7‐BR). H172F was inactive with all three substrates, and L437A did not turn over 7‐BR. Furthermore, H172A, H172Q, E474D and E474Q showed large changes in kcat/KM for each of the three substrates, in some cases up to 50‐fold. Concurrent molecular dynamics simulations yielded distances between some of the residues in these putative interaction pairs that are not consistent with contact. The results indicate that small changes in the protein scaffold lead to large differences in solution behavior and enzyme function.

Structural and thermodynamic basis of (+)-α-pinene binding to human cytochrome P450 2B6.

Despite recent advances in atomic-level understanding of drug and inhibitor interactions with human cytochromes P450, the decades-old questions of chemical and structural determinants of hydrocarbon binding are still unanswered. (+)-α-Pinene is a monoterpene hydrocarbon that is widely distributed in the environment and a potent P450 2B inhibitor. Therefore, a combined biophysical and structural analysis of human P450 2B6 interactions with (+)-α-pinene was undertaken to elucidate the basis of the very high affinity binding. Binding of (+)-α-pinene to the P450 active site was demonstrated by a Type I spectral shift. Thermodynamics of ligand binding were explored using isothermal titration calorimetry and compared to those of P450 2A6, which is much less flexible than 2B6 based on comparison of multiple X-ray crystal structures. Consistent with expectation, entropy is the major driving force for hydrocarbon binding to P450 2A6, as evidenced by the calorimetric results. However, formation of the 2B6-(+)-α-pinene complex has a significant enthalpic component. A 2.0 Å resolution crystal structure of this enzyme-ligand complex reveals that the highly plastic 2B6 utilizes previously unrecognized rearrangements of protein motifs. The results indicate that the specific components of enthalpic contribution to ligand binding are closely tied to the degree of enzyme flexibility.

Rational Engineering of Cytochromes P450 2B6 and 2B11 for Enhanced Stability: Insights Into Structural Importance of Residue 334

Rational mutagenesis was used to improve the thermal stability of human cytochrome P450 2B6 and canine P450 2B11. Comparison of the amino acid sequences revealed seven sites that are conserved between the stable 2B1 and 2B4 but different from those found in the less stable 2B6 and 2B11. P334S was the only mutant that showed increased heterologous expression levels and thermal stability in both 2B6 and 2B11. The mechanism of this effect was explored with pressure-perturbation spectroscopy. Compressibility of the heme pocket in variants of all four CYP2B enzymes containing proline at position 334 are characterized by lower compressibility than their more stable serine 334 counterpart. Therefore, the stabilizing effect of P334S is associated with increased conformational flexibility in the region of the heme pocket. Improved stability of P334S 2B6 and 2B11 may facilitate the studies of these enzymes by X-ray crystallography and biophysical techniques.

Structure-Function Analysis of Mammalian CYP2B Enzymes Using 7-Substituted Coumarin Derivatives as Probes: Utility of Crystal Structures and Molecular Modeling in Understanding Xenobiotic Metabolism

Crystal structures of CYP2B35 and CYP2B37 from the desert woodrat were solved in complex with 4-(4-chlorophenyl)imidazole (4-CPI). The closed conformation of CYP2B35 contained two molecules of 4-CPI within the active site, whereas the CYP2B37 structure demonstrated an open conformation with three 4-CPI molecules, one within the active site and the other two in the substrate access channel. To probe structure-function relationships of CYP2B35, CYP2B37, and the related CYP2B36, we tested the O-dealkylation of three series of related substrates—namely, 7-alkoxycoumarins, 7-alkoxy-4-(trifluoromethyl)coumarins, and 7-alkoxy-4-methylcoumarins—with a C1-C7 side chain. CYP2B35 showed the highest catalytic efficiency (kcat/KM) with 7-heptoxycoumarin as a substrate, followed by 7-hexoxycoumarin. In contrast, CYP2B37 showed the highest catalytic efficiency with 7-ethoxy-4-(trifluoromethyl)coumarin (7-EFC), followed by 7-methoxy-4-(trifluoromethyl)coumarin (7-MFC). CYP2B35 had no dealkylation activity with 7-MFC or 7-EFC. Furthermore, the new CYP2B–4-CPI–bound structures were used as templates for docking the 7-substituted coumarin derivatives, which revealed orientations consistent with the functional studies. In addition, the observation of multiple –Cl and –NH-π interactions of 4-CPI with the aromatic side chains in the CYP2B35 and CYP2B37 structures provides insight into the influence of such functional groups on CYP2B ligand binding affinity and specificity. To conclude, structural, computational, and functional analysis revealed striking differences between the active sites of CYP2B35 and CYP2B37 that will aid in the elucidation of new structure-activity relationships.