Key substrate recognition residues in the active site of a plant cytochrome P450, CYP73A1. Homology guided site-directed mutagenesis

Engineering medicinal plant-derived CYPs: a promising strategy for production of high-valued secondary metabolites

Cytochorme P450s (CYPs) play a critical role in the catalysis of secondary metabolite biosynthetic pathways. For their commercial use, various strategies for metabolic pathway engineering using CYP as a potential target have been explored. Plants produce a vast diversity of secondary metabolites which are being used to treat various ailments and diseases. Some of these metabolites are difficult to obtain in large quantities limiting their industrial use. Cytochrome P450 enzymes (CYPs) are important catalysts in the biosynthesis of highly valued secondary metabolites, and are found in all domains of life. With the development of high-throughput sequencing and high-resolution mass spectrometry, new biosynthetic pathways and associated CYPs are being identified. In this review, we present CYPs identified from medicinal plants as a potential game changer in the metabolic engineering of secondary metabolic pathways. We present the achievements made so far in enhancing the production of important bioactivities through pathway engineering, giving some popular examples. At last, current challenges and possible strategies to overcome the limitations associated with CYP engineering to enhance the biosynthesis of target secondary metabolites are also highlighted.

Functional characterization of phenylalanine ammonia-lyase- and cinnamate 4-hydroxylase-encoding genes from Lycoris radiata, a galanthamine-producing plant

Galanthamine (GAL), the well-known Amaryllidaceae alkaloid, is a clinically used drug for the treatment of Alzheimer's disease. L-Phenylalanine (Phe) and trans-cinnamic acid (CA) were enzymatically transformed into the catechol portion of GAL. Herein, a Phe ammonia-lyase-encoding gene LrPAL3 and a cinnamate 4-hydroxylase-encoding gene LrC4H were cloned from Lycoris radiata, a GAL-producing plant. LrPAL3 was overexpressed in Escherichia coli and purified to homogeneity. LrPAL3 catalyzes the forward deamination conversion of L-Phe into trans-CA. The 3-chloro- and 4-fluoro-L-Phe were deaminated to generate the corresponding 3-chloro- and 4-fluoro-trans-CA by LrPAL3. LrPAL3-catalyzed reverse hydroamination was confirmed by the conversion of trans-CA into L-Phe with exceptional regio- and stereo-selectivity. LrC4H was overexpressed in E. coli with tCamCPR, a cytochrome P450 reductase-encoding gene. LrC4H catalyzes the regioselective para-hydroxylation on trans-CA to form p-coumaric acid. The transcriptional levels of both LrPAL3 and LrC4H were positively associated with the GAL contents within the leaves and flowers of L. radiata, which suggested that their expression and function are co-regulated and involved in the biosynthesis of GAL. The present investigations on the biosynthetic genes of GAL will promote the development of synthetic biology platforms for this kind of important drug via metabolic engineering.

The isolation and functional characterization of three liverwort genes encoding cinnamate 4-hydroxylase

The plant phenylpropanoid pathway is responsible for the synthesis of a wide variety of secondary metabolites. The second step in phenylpropanoid synthesis is carried out by the cytochrome P450 monooxygenase enzyme cinnamate 4-hydroxylase (C4H), which catalyzes the p-hydroxylation of trans-cinnamic acid to p-coumarate. Genes encoding C4H have been characterized in many vascular plant species, but as yet not in any bryophyte species. Here, a survey of the transcriptome sequences of four liverwort species was able to identify eight putative C4Hs. The three liverwort C4H genes taken forward for isolation and functional characterization were harbored by Plagiochasma appendiculatum (PaC4H) and Marchantia paleacea (MpC4H1 and MpC4H2). When the genes were heterologously expressed in yeast culture, an assay of enzyme activity indicated that PaC4H and MpC4H1 had a higher level of activity than MpC4H2. The favored substrate (trans-cinnamic acid) of all three liverwort C4Hs was the same as that of higher plant C4Hs. The co-expression of PaC4H in yeast cells harboring PaPAL (a P. appendiculatum ene encoding phenylalanine ammonia lyase) allowed the conversion of L-phenylalanine to p-coumaric acid. Furthermore, the expression level of PaC4H was enhanced after treatment with abiotic stress inducers UV irradiation or salicylic acid in the thallus of P. appendiculatum. The likelihood is that high activity C4Hs evolved in the liverworts and have remained highly conserved across the plant kingdom.

Structural and energetic analysis to provide insight residues of CYP2C9, 2C11 and 2E1 involved in valproic acid dehydrogenation selectivity

Docking and molecular dynamics (MD) simulation have been two computational techniques used to gain insight about the substrate orientation within protein active sites, allowing to identify potential residues involved in the binding and catalytic mechanisms. In this study, both methods were combined to predict the regioselectivity in the binding mode of valproic acid (VPA) on three cytochrome P-450 (CYP) isoforms CYP2C9, CYP2C11, and CYP2E1, which are involved in the biotransformation of VPA yielding reactive hepatotoxic intermediate 2-n-propyl-4-pentenoic acid (4nVPA). There are experimental data about hydrogen atom abstraction of the C4-position of VPA to yield 4nVPA, however, there are not structural evidence about the binding mode of VPA and 4nVPA on CYPs. Therefore, the complexes between these CYP isoforms and VPA or 4nVPA were studied to explore their differences in binding and energetic stabilization. Docking results showed that VPA and 4nVPA are coupled into CYPs binding site in a similar conformation, but it does not explain the VPA hydrogen atom abstraction. On the other hand, MD simulations showed a set of energetic states that reorient VPA at the first ns, then making it susceptible to a dehydrogenation reaction. For 4nVPA, multiple binding modes were observed in which the different states could favor either undergo other reaction mechanism or ligand expulsion from the binding site. Otherwise, the energetic and entropic contribution point out a similar behavior for the three CYP complexes, showing as expected a more energetically favorable binding free energy for the complexes between CYPs and VPA than with 4nVPA.

Analysis of flavonoids and the flavonoid structural genes in brown fiber of upland cotton

Backgroud

As a result of changing consumer preferences, cotton (Gossypium Hirsutum L.) from varieties with naturally colored fibers is becoming increasingly sought after in the textile industry. The molecular mechanisms leading to colored fiber development are still largely unknown, although it is expected that the color is derived from flavanoids.

Experimental Design

Firstly, four key genes of the flavonoid biosynthetic pathway in cotton (GhC4H, GhCHS, GhF3′H, and GhF3′5′H) were cloned and studied their expression profiles during the development of brown- and white cotton fibers by QRT-PCR. And then, the concentrations of four components of the flavonoid biosynthetic pathway, naringenin, quercetin, kaempferol and myricetin in brown- and white fibers were analyzed at different developmental stages by HPLC.

Result

The predicted proteins of the four flavonoid structural genes corresponding to these genes exhibit strong sequence similarity to their counterparts in various plant species. Transcript levels for all four genes were considerably higher in developing brown fibers than in white fibers from a near isogenic line (NIL). The contents of four flavonoids (naringenin, quercetin, kaempferol and myricetin) were significantly higher in brown than in white fibers and corresponding to the biosynthetic gene expression levels.

Conclusions

Flavonoid structural gene expression and flavonoid metabolism are important in the development of pigmentation in brown cotton fibers.

Cinnamate 4-Hydroxylase (C4H) genes from Leucaena leucocephala: a pulp yielding leguminous tree

Leucaena leucocephala is a leguminous tree species accounting for one-fourth of raw material supplied to paper and pulp industry in India. Cinnamate 4-Hydroxylase (C4H, EC 1.14.13.11) is the second gene of phenylpropanoid pathway and a member of cytochrome P450 family. There is currently intense interest to alter or modify lignin content of L. leucocephala. Three highly similar C4H alleles of LlC4H1 gene were isolated and characterized. The alleles shared more than 98 % sequence identity at amino acid level to each other. Binding of partial promoter of another C4H gene LlC4H2, to varying amounts of crude nuclear proteins isolated from leaf and stem tissues of L. leucocephala formed two loose and one strong complex, respectively, suggesting that the abundance of proteins that bind with the partial C4H promoter is higher in stem tissue than in leaf tissue. Quantitative Real Time PCR study suggested that among tissues of same age, root tissues had highest level of C4H transcripts. Maximum transcript level was observed in 30 day old root tissue. Among the tissues investigated, C4H activity was highest in 60 day old root tissues. Tissue specific quantitative comparison of lignin from developing seedling stage to 1 year old tree stage indicated that Klason lignin increased in tissues with age.

Mechanism of the plant cytochrome P450 for herbicide resistance: a modelling study

Plant cytochrome P450 is a key enzyme responsible for the herbicide resistance but the molecular basis of the mechanism is unclear. To understand this, four typical plant P450s and a widely resistant herbicide chlortoluron were analysed by carrying out homology modelling, molecular docking, molecular dynamics simulations and binding free energy analysis. Our results demonstrate that: (i) the putative hydrophobic residues located in the F-helix and polar residues in I-helix are critical in the herbicide resistance; (ii) the binding mode analysis and binding free energy calculation indicate that the distance between catalytic site of chlortoluron and heme of P450, as well as the binding affinity are key elements affecting the resistance for plants. In conclusion, this work provides a new insight into the interactions of plant P450s with herbicide from a molecular level, offering valuable information for the future design of novel effective herbicides which also escape from the P450 metabolism.

Biotransformation of isonitrosoacetophenone (2-keto-2-phenyl-acetaldoxime) in tobacco cell suspensions

Nicotiana tabacum cell suspensions, 2 g wet wt/ml, rapidly took up 1 mM isonitrosoacetophenone (INAP), a plant-derived stress metabolite with anti-oxidative and anti-fungal properties, producing 4'-hexopyranosyloxy-3'-methoxyisonitrosoacetophenone in 54 % yield over 18 h. Unconverted INAP was at 33 μM. UPLC-MS/MS analyses with MassFragment software were used for metabolite identification. INAP had been hydroxylated at its meta- and para-positions as well as undergoing subsequent methoxylation and glycosylation. INAP is thus recognized by the enzymatic machinery of the phenylpropanoid pathway and is converted to a molecule with a substitution pattern similar to ferulic acid.

Homology modeling of the three membrane proteins of the dhurrin metabolon: catalytic sites, membrane surface association and protein-protein interactions

Formation of metabolons (macromolecular enzyme complexes) facilitates the channelling of substrates in biosynthetic pathways. Metabolon formation is a dynamic process in which transient structures mediated by weak protein-protein interactions are formed. In Sorghum, the cyanogenic glucoside dhurrin is derived from l-tyrosine in a pathway involving the two cytochromes P450 (CYPs) CYP79A1 and CYP71E1, a glucosyltransferase (UGT85B1), and the redox partner NADPH-dependent cytochrome P450 reductase (CPR). Experimental evidence suggests that the enzymes of this pathway form a metabolon. Homology modeling of the three membrane bound proteins was carried out using the Sybyl software and available relevant crystal structures. Residues involved in tight positioning of the substrates and intermediates in the active sites of CYP79A1 and CYP71E1 were identified. In both CYPs, hydrophobic surface domains close to the N-terminal trans-membrane anchor and between the F' and G helices were identified as involved in membrane anchoring. The proximal surface of both CYPs showed positively charged patches complementary to a negatively charged bulge on CPR carrying the FMN domain. A patch of surface exposed, positively charged amino acid residues positioned on the opposite face of the membrane anchor was identified in CYP71E1 and might be involved in binding UGT85B1 via a hypervariable negatively charged loop in this protein.

Improved cytochrome P450 3A4 molecular models accurately predict the Phe215 requirement for raloxifene dehydrogenation selectivity

The use of molecular modeling in conjunction with site-directed mutagenesis has been extensively used to study substrate orientation within cytochrome P450 active sites and to identify potential residues involved in the positioning and catalytic mechanisms of these substrates. However, because docking studies utilize static models to simulate dynamic P450 enzymes, the effectiveness of these studies is strongly dependent on accurate enzyme models. This study employed a cytochrome P450 3A4 (CYP3A4) crystal structure (Protein Data Bank entry 1W0E) to predict the sites of metabolism of the known CYP3A4 substrate raloxifene. In addition, partial charges were incorporated into the P450 heme moiety to investigate the effect of the modified CYP3A4 model on metabolite prediction with the ligand docking program Autodock. Dehydrogenation of raloxifene to an electrophilic diquinone methide intermediate has been linked to the potent inactivation of CYP3A4. Active site residues involved in the positioning and/or catalysis of raloxifene supporting dehydrogenation were identified with the two models, and site-directed mutagenesis studies were conducted to validate the models. The addition of partial charges to the heme moiety improved the accuracy of the docking studies, increasing the number of conformations predicting dehydrogenation and facilitating the identification of substrate-active site residue interactions. On the basis of the improved model, the Phe215 residue was hypothesized to play an important role in orienting raloxifene for dehydrogenation through a combination of electrostatic and steric interactions. Substitution of this residue with glycine or glutamine significantly decreased dehydrogenation rates without concurrent changes in the rates of raloxifene oxygenation. Thus, the improved structural model predicted novel enzyme-substrate interactions that control the selective dehydrogenation of raloxifene to its protein-binding intermediate.

Isolation and characterization of a gene encoding cinnamate 4-hydroxylase from Parthenocissus henryana

Cinnamate 4-hydroxylase (C4H, EC 1.14.13.11) plays an important role in the phenylpropanoid pathway, which produces many economically important secondary metabolites. A gene coding for C4H, designated as PhC4H (GenBank accession no. DQ211885) was isolated from Parthenocissus henryana. The full-length PhC4H cDNA is 1,747 bp long with a 1,518-bp open reading frame encoding a protein of 505 amino acids, a 40-bp 5' non-coding region and a 189-bp 3'-untranslated region. Secondary structure of the deduced PhC4H protein consists of 41.78% alpha helix, 15.64% extended strand and 42.57% random coil. The genomic DNA of PhC4H is 2,895 bp long and contains two introns; intron I is 205-bp and intron II is 1,172-bp (GenBank accession no. EU440734). DNA gel blot analysis revealed that there might be a single copy of PhC4H in Parthenocissus henryana genome. By using anchored PCR, a 963-bp promoter sequence was isolated and it contains many responsive elements conserved in the upstream region of PAL, C4H and 4CL including the P-, A-, L- and H-boxes.

Single mutations change CYP2F3 from a dehydrogenase of 3-methylindole to an oxygenase

Pulmonary cytochrome P450 2F3 (CYP2F3) catalyzes the dehydrogenation of the pneumotoxin 3-methylindole (3MI) to an electrophilic intermediate, 3-methyleneindolenine, which is responsible for the toxicity of the parent compound. Members of the CYP2F subfamily are the only enzymes known to exclusively dehydrogenate 3MI, without detectable formation of oxygenation products. Thus, CYP2F3 is an attractive model to study dehydrogenation mechanisms. The purpose of this study was to identify specific residues that could facilitate 3MI dehydrogenation. Both single and double mutations were constructed to study the molecular mechanisms that direct dehydrogenation. Double mutations in substrate recognition sites (SRS) 1 produced an inactive enzyme, while double mutants in SRS 4 did not alter 3MI metabolism. However, double mutations in SRS 5 and SRS 6 successfully introduced oxygenase activity to CYP2F3. Single mutations in SRS 5, SRS 6 and near SRS 2 also introduced 3MI oxygenase activity. Mutants S474H and D361T oxygenated 3MI but also increased dehydrogenation rates, while G214L, E215Q and S475I catalyzed 3MI oxygenation exclusively. A homology model of CYP2F3 was precisely consistent with specific dehydrogenation of 3MI via initial hydrogen atom abstraction from the methyl group. In addition, intramolecular kinetic deuterium isotope studies demonstrated an isotope effect ( K H/ K D) of 6.8. This relatively high intramolecular deuterium isotope effect confirmed the initial hydrogen abstraction step; a mutant (D361T) that retained the dehydrogenation reaction exhibited the same deuterium isotope effect. The results showed that a single alteration, such as a serine to isoleucine change at residue 475, dramatically switched catalytic preference from dehydrogenation to oxygenation.

Molecular cloning of two genes encoding cinnamate 4-hydroxylase (C4H) from oilseed rape (Brassica napus)

Cinnamate 4-hydroxylase (C4H) is a key enzyme of phenylpropanoid pathway, which synthesizes numerous secondary metabolites to participate in development and adaption. Two C4H isoforms, the 2192-bp BnC4H-1 and 2108-bp BnC4H-2, were cloned from oilseed rape (Brassica napus). They both have two introns and a 1518-bp open reading frame encoding a 505-amino-acid polypeptide. BnC4H-1 is 57.73 kDa with an isoelectric point of 9.11, while 57.75 kDa and 9.13 for BnC4H-2. They share only 80.6% identities on nucleotide level but 96.6% identities and 98.4% positives on protein level. Showing highest homologies to Arabidopsis thaliana C4H, they possess a conserved p450 domain and all P450-featured motifs, and are identical to typical C4Hs at substrate-recognition sites and active site residues. They are most probably associated with endoplasmic reticulum by one or both of the N- and C-terminal transmembrane helices. Phosphorylation may be a necessary post-translational modification. Their secondary structures are dominated by alpha helices and random coils. Most helices locate in the central region, while extended strands mainly distribute before and after this region. Southern blot indicated about 9 or more C4H paralogs in B. napus. In hypocotyl, cotyledon, stem, flower, bud, young- and middle-stage seed, they are co-dominantly expressed. In root and old seed, BnC4H-2 is dominant over BnC4H-1, with a reverse trend in leaf and pericarp. Paralogous C4H numbers in Brassicaceae genomes and possible roles of conserved motifs in 5' UTR and the 2nd intron are discussed.

Engineering proteinase K using machine learning and synthetic genes

Background

Altering a protein's function by changing its sequence allows natural proteins to be converted into useful molecular tools. Current protein engineering methods are limited by a lack of high throughput physical or computational tests that can accurately predict protein activity under conditions relevant to its final application. Here we describe a new synthetic biology approach to protein engineering that avoids these limitations by combining high throughput gene synthesis with machine learning-based design algorithms.

Results

We selected 24 amino acid substitutions to make in proteinase K from alignments of homologous sequences. We then designed and synthesized 59 specific proteinase K variants containing different combinations of the selected substitutions. The 59 variants were tested for their ability to hydrolyze a tetrapeptide substrate after the enzyme was first heated to 68°C for 5 minutes. Sequence and activity data was analyzed using machine learning algorithms. This analysis was used to design a new set of variants predicted to have increased activity over the training set, that were then synthesized and tested. By performing two cycles of machine learning analysis and variant design we obtained 20-fold improved proteinase K variants while only testing a total of 95 variant enzymes.

Conclusion

The number of protein variants that must be tested to obtain significant functional improvements determines the type of tests that can be performed. Protein engineers wishing to modify the property of a protein to shrink tumours or catalyze chemical reactions under industrial conditions have until now been forced to accept high throughput surrogate screens to measure protein properties that they hope will correlate with the functionalities that they intend to modify. By reducing the number of variants that must be tested to fewer than 100, machine learning algorithms make it possible to use more complex and expensive tests so that only protein properties that are directly relevant to the desired application need to be measured. Protein design algorithms that only require the testing of a small number of variants represent a significant step towards a generic, resource-optimized protein engineering process.

Non-natural cinnamic acid derivatives as substrates of cinnamate 4-hydroxylase

Cinnamate 4-hydroxylase (C4H), a monooxygenase in the plant phenylpropanoid pathway, was assayed for its ability to hydroxylate 29 substrate analogues. Nine of the tested analogues with various aromatic side chains, including 3-coumaric acid, were metabolized by C4H. Seven products from these reactive analogues were characterized using LC/MS, 1H NMR and 13C NMR spectroscopic analysis. For example, caffeic acid was the product of 3-coumaric acid. The products 4-hydroxy-2-chlorocinnamic acid and 4-hydroxy-2-ethoxycinnamic acid are novel compounds that have not been previously reported. The kinetic parameters of C4H towards these analogues were determined.

Functional interaction of nitrogenous organic bases with cytochrome P450: A critical assessment and update of substrate features and predicted key active-site elements steering the access, binding, and orientation of amines

The widespread use of nitrogenous organic bases as environmental chemicals, food additives, and clinically important drugs necessitates precise knowledge about the molecular principles governing biotransformation of this category of substrates. In this regard, analysis of the topological background of complex formation between amines and P450s, acting as major catalysts in C- and N-oxidative attack, is of paramount importance. Thus, progress in collaborative investigations, combining physico-chemical techniques with chemical-modification as well as genetic engineering experiments, enables substantiation of hypothetical work resulting from the design of pharmacophores or homology modelling of P450s. Based on a general, CYP2D6-related construct, the majority of prospective amine-docking residues was found to cluster near the distal heme face in the six known SRSs, made up by the highly variant helices B', F and G as well as the N-terminal portion of helix C and certain beta-structures. Most of the contact sites examined show a frequency of conservation < 20%, hinting at the requirement of some degree of conformational versatility, while a limited number of amino acids exhibiting a higher level of conservation reside close to the heme core. Some key determinants may have a dual role in amine binding and/or maintenance of protein integrity. Importantly, a series of non-SRS elements are likely to be operative via long-range effects. While hydrophobic mechanisms appear to dominate orientation of the nitrogenous compounds toward the iron-oxene species, polar residues seem to foster binding events through H-bonding or salt-bridge formation. Careful uncovering of structure-function relationships in amine-enzyme association together with recently developed unsupervised machine learning approaches will be helpful in both tailoring of novel amine-type drugs and early elimination of potentially toxic or mutagenic candidates. Also, chimeragenesis might serve in the construction of more efficient P450s for activation of amine drugs and/or bioremediation.

Kinetic and molecular analysis of 5-epiaristolochene 1,3-dihydroxylase, a cytochrome P450 enzyme catalyzing successive hydroxylations of sesquiterpenes

The final step of capsidiol biosynthesis is catalyzed by 5-epiaristolochene dihydroxylase (EAH), a cytochrome P450 enzyme that catalyzes the regio- and stereospecific insertion of two hydroxyl moieties into the bicyclic sesquiterpene 5-epiaristolochene (EA). Detailed kinetic studies using EA and the two possible monohydroxylated intermediates demonstrated the release of 1beta-hydroxy-EA ((OH)EA) at high EA concentrations and a 10-fold catalytic preference for 1beta(OH)EA versus 3alpha(OH)EA, indicative of a preferred reaction order of hydroxylation at C-1, followed by that at C-3. Sequence alignments and homology modeling identified active-site residues tested for their contribution to substrate specificity and overall enzymatic activity. Mutants EAH-S368C and EAH-S368V exhibited wild-type catalytic efficiencies for 1beta(OH)EA biosynthesis, but were devoid of the successive hydroxylation activity for capsidiol biosynthesis. In contrast to EAH-S368C, EAH-S368V catalyzed the relative equal biosynthesis of 1beta(OH)EA, 2beta(OH)EA, and 3beta(OH)EA from EA with wild-type efficiency. Moreover, EAH-S368V converted approximately 1.5% of these monohydroxylated products to their respective ketone forms. Alanine and threonine mutations at position 368 were significantly compromised in their conversion rates of EA to capsidiol and correlated with 3.6- and 5.7-fold increases in their Km values for the 1beta(OH)EA intermediate, respectively. A role for Ile486 in the successive hydroxylations of EA was also suggested by the EAH-I468A mutant, which produced significant amounts 1beta(OH)EA, but negligible amounts of capsidiol from EA. The altered product profile of the EAH-I486A mutant correlated with a 3.6-fold higher Km for EA and a 4.4-fold slower turnover rate (kcat) for 1beta(OH)EA. These kinetic and mutational studies were correlated with substrate docking predictions to suggest how Ser368 and Ile486 might contribute to active-site topology, substrate binding, and substrate presentation to the oxo-Fe-heme reaction center.

Engineering of a water-soluble plant cytochrome P450, CYP73A1, and NMR-based orientation of natural and alternate substrates in the active site

CYP73A1 catalyzes cinnamic acid hydroxylation, a reaction essential for the synthesis of lignin monomers and most phenolic compounds in higher plants. The native CYP73A1, initially isolated from Jerusalem artichoke (Helianthus tuberosus), was engineered to simplify purification from recombinant yeast and improve solublity and stability in the absence of detergent by replacing the hydrophobic N terminus with the peptitergent amphipathic sequence PD1. Optimized expression and purification procedures yielded 4 mg engineered CYP73A1 L(-1) yeast culture. This water-soluble enzyme was suitable for 1H-nuclear magnetic resonance (NMR) investigation of substrate positioning in the active site. The metabolism and interaction with the enzyme of cinnamate and four analogs were compared by UV-visible and 1H-NMR analysis. It was shown that trans-3-thienylacrylic acid, trans-2-thienylacrylic acid, and 4-vinylbenzoic acid are good ligands and substrates, whereas trans-4-fluorocinnamate is a competitive inhibitor. Paramagnetic relaxation effects of CYP73A1-Fe(III) on the 1H-NMR spectra of cinnamate and analogs indicate that their average initial orientation in the active site is parallel to the heme. Initial orientation and distances of ring protons to the iron do not explain the selective hydroxylation of cinnamate in the 4-position or the formation of single products from the thienyl compounds. Position adjustments are thus likely to occur during the later steps of the catalytic cycle.