Structure-activity correlations in pentachlorobenzene oxidation by engineered cytochrome P450cam

We had reported engineering of the heme monooxygenase cytochrome P450cam from Pseudomonas putida with the F87W/Y96F/L244A/V247L mutations for the oxidation of pentachlorobenzene (PeCB), a recalcitrant environmental contaminant, to pentachlorophenol. In order to provide further insights into P450 structure, function and substrate recognition, we have determined the crystal structure of this 4-mutant without a substrate and its complex with PeCB. PeCB is bound face-on to the heme, with a weak Fe--Cl interaction. One PeCB chlorine is located in the cavity generated by the L244A mutation, in striking illustration of the role of this mutation in promoting PeCB binding. The structures also show that the P450(cam) oxygen-binding groove between G248 and T252 is flexible and can tolerate significant deviations from their conformations in the wild type without loss of enzyme activity. Analysis of the PeCB binding interactions led to introduction of the T101A mutation to enable the substrate to reorient during the catalytic cycle for more efficient oxidation. The resultant 5-mutant F87W/Y96F/T101A/L244A/V247L is 3-fold more active for PeCB oxidation than the 4-mutant. Polychlorinated benzene binding by the mutants and the partitioning between substrate oxidation and non-productive (uncoupling) side reactions are correlated with the structural data.     

Comparison of the construction of a 3-D model for human thromboxane synthase using P450cam and BM-3 as templates: implications for the substrate binding pocket

A 3-D model of human thromboxane A₂ synthase (TXAS) was constructedusing a homology modeling approach based on information fromthe 2.0 crystal structure of the hemoprotein domains of cytochromeP450BM-3 and P450cam. P450BM-3 is a bacterial fatty acid monooxygenaseresembling eukaryotic microsomal cytochrome P450s in primarystructure and function. TXAS shares 26.4% residue identity and48.4% residue similarity with the P450BM-3 hemoprotein domain.The homology score between TXAS and P450BM-3 is much higherthan that between TXAS and P450cam. Alignment between TXAS andthe P450BM-3 hemoprotein domain or P450cam was determined throughsequence searches. The P450BM-3 or P450cam main-chain coordinateswere spplied to the TXAS main chain in those sements where thetwo sequences were well aligned. These segments were linkedto one another using a fragment search method, and the sidechains were added to produce a 3-D model for TXAS. A TXAS substrate,prostaglandin H₂ (PGH₂) was docked into the TXAS cavity correspondingto the arachidonic acid binding pocket in P450BM-3 or camphorbinding site in P450cam. Regions of the heme and putative PGH2binding cavities in the TXAS model were identified and analyzed.The segments and residues involved in the active-site pocketof the TXAS model provide reasonable candidates for TXAS proteinengineering and inhibitor design. Comparison of the TXAS modelbased on P450BM-3 with another TXAS model based on the P450BM-3with another TXAS model based on the P450cam structure indicatedthat P450BM-3 is a more suitable template for homology modelingof TXAS.     

Enzyme-catalyzed dehalogenation of pentachloroethane: why F87W-cytochrome P450cam is faster than wild type

Under anaerobic conditions, cytochromes P450 can reductivelydehalogenate heavily halogenated hydrocarbons, such as one-and two-carbon organic solvents. This catalytic capacity hasdrawn attention to the potential use of engineered forms ofP450s in the remediation of contaminated deep subsurface ecosystems.Loida (1994, PhD Thesis, University of Illinois at Urbana-Champaign,IL) and S.G.Sligar (personal communication) have observedrecentlythat an active-site variant of cytochrome P450cam (F87W) dechlorinatespentachloroethane approximately three times faster than thewild-type enzyme. Molecular dynamics simulations have revealedthat the mutant enzyme binding pocket remains smaller, and thatpentachloroethane assumes configurations closer to the heme-Fein the F87W mutant twice as often as in the wild-type enzyme.This result is consistent with a collisional model of dehalogenation,which agrees with experimental observations [Li and Wackett(1993) Biochemistry, 32, 9355–9361] that solutions containingwild-type P450cam dehalogenate pentachloroethane 100 times fasterthan those containing free heme. The simulations suggest thatit is unlikely that Trp87 significantly stabilizes the developingnegative charge on the substrate during carbon-halogen bondreduction. The design of improved microbiai enzymes that incorporateboth steric and electronic effects continues for use in remediatinghalogenated contaminants in situ     

A predicted three-dimensional structure of human cytochrome P450: implications for substrate specificity

A three-dimensional structure for human cytochrome P450IA1 waspredicted based on the crystal coordinates of cytochrome P450camfrom Pseudomonas putida. As there was only 15% residue identitybetween the two enzymes, additional information was used toestablish an accurate sequence alignment that is a prerequisitefor model building. Twelve representative eukaryotic sequenceswere aligned and a net prediction of secondary structure wasmatched against the known -helices and ß-sheets ofP450cam. The cam secondary structure provided a fixed main-chainframework onto which loops of appropriate length from the humanP450IA1 structure were added. The model-built structure of thehuman cytochrome conformed to the requirements for the segregationof polar and nonpolar residues between the core and the surface.The first 44 residues of human cytochrome P450 could not bebuilt into the model and sequence analysis suggested that residues1–26 formed a single membrane-spanning segment. Examinationof the sequences of cytochrome P450s from distinct gene familiessuggested specific residues that could account for the differencesin substrate specificity. A major substrate for P450IA1, 3-methyl-cholanthrene,was fitted into the proposed active site and this planar aromaticmolecule could be accommodated into the available cavity. Residuesthat are likely to interact with the haem were identified. Thesequence similarity between 59 eukaryotic enzymes was representedas a dendrogram that in general clustered according to genefamily. Until a crystallographic structure is available, thismodel-building study identifies potential residues in cytochromeP450s important in the function of these enzymes and these residuesare candidates for site-directed mutagenesis.     

Relative binding free energy calculations of inhibitors to two mutants (Glu46-- Ala/Gln) of ribonuclease T1 using molecular dynamics/free energy perturbation approaches

We present free energy perturbation calculations on the complexesof Glu46— Ala46 (E46A) and Glu46— Gln46 (E46Q) mutantsof ribonuclease T1 (RNaseT1) with inhibitors 2‘-guanosinemonophosphate (GMP) and 2’adenosine monophosphate (AMP)by a thermodynamic perturbation method implemented with moleculardynamics (MD). Using the available crystal structure of theRNaseT1–GMP complex, the structures of E46A-GMP and E46Q-GMPwere model built and equilibrated with MD simulations. The structuresof E46A-AMP and E46Q-AMP were obtained as a final structureof the GMP—AMP perturbation calculation respectively.The calculated difference in the free energy of binding (Gbind)was 0.31 kcal/mol for the E46A system and —1.04 kcal/molfor the E46Q system. The resultant free energies are much smallerthan the experimental and calculated value of 3 kcal/mol forthe native RNase T1, which suggests that both mutants have greaterrelative adenine affinities than native RNaseT1. EspeciallyE46Q is calculated to have a larger affinity for adenine thanguanine, as we suggested previously from the calculation onthe native RNaseT1. Thus, the molecular dynamics/free energyperturbation method may be helpful in protein engineering, directedtoward increasing or changing the substrate specificity of enzymes.     

Free energy calculations of the mutation of Ile96 -> Ala in barnase: contributions to the difference in stability

Free energy calculations were carried out to determine the relativeunfolding free energy of the Ile96 wild type and Ala96 mutantbarnases. The total calculated free energies suggest that substitutionof Ile96 with Ala destabilizes barnase by 3.9 kcal/mol, whichis in good agreement with the independently determined experimentalvalues of 4.0 and 3.3 kcal/mol and a previous simulation. However,a decomposition of the free energy finds the dominant contributionsto this free energy arising from the noncovalent Interactionsbetween the perturbed group and distant residues of barnasein the sequence and water molecules and only a very small contributionfrom covalent interactions. This is in contrast to the previoussimulation, using the dual topology methodology, which produceda decomposition with an {small tilde}60% free energy contributionfrom changes in covalent interactions. The use of the singletopology employed in the present calculations and the dual topologyemployed in the previous study are analyzed in order to understandthe contrast between the present results and the results ofthe previous study.     

Cytochrome P450 active site plasticity: attenuation of imidazole binding in cytochrome P450(cam) by an L244A mutation

We have identified a P450(cam) mutation, L244A, that mitigates the affinity for imidazole and substituted imidazoles while maintaining a high affinity for the natural substrate camphor. The P450(cam) L244A crystal structure solved in the absence of any ligand reveals that the I-helix is displaced inwards by over 1 A in response to the cavity created by the change from leucine to alanine. Furthermore, the crystal structures of imidazole-bound P450(cam) and the 1-methylimidazole-bound P450(cam) L244A mutant reveal that the ligands have distinct binding modes in the two proteins. Whereas in wild-type P450(cam) the imidazole coordinates to the iron in an orientation roughly perpendicular to the plane of the heme, in the L244A mutant the rearranged I helix, and specifically residue Val247, forces the imidazole into an orientation almost parallel to the heme that impairs its ability to coordinate to the heme iron. As a result, the imidazole is much more weakly bound to the mutant than it is to the wild-type enzyme. Despite the constriction of the active site by the mutation, previous work with the L244A mutant has shown that it oxidizes larger substrates than the wild-type enzyme. This paradoxical situation, in which a mutation that nominally increases the active site cavity appears to decrease it, suggests that the mutation actually increases the active site maleability, allowing it to better expand to oxidize larger substrates.     

Application of free energy simulations to the binding of a transition-state-analogue inhibitor to HTV protease

Free energy simulations (slow-change method) have been usedto estimate quantitatively the ratio of the binding constantsof (S) and (R) isomers of a novel HIV protease inhibitor, JG365.As a starting geometry, we used the X-ray crystallographic structureof a complex of HTV protease and JG365 provided by A.Wlodawer.According to our results the (S) configuration, i.e. the formpreviously identified experimentally, binds considerably moretightly to the protease (G° = 2.9 kcal/mol). When the (S)inhibitor is bound, there is a very strong preference for protonationof the Aspl25 (rather than the Asp25) residue of the protease.This study is the first to apply a new method for quantitativelyassessing the precision of free energies calculated by the slow-changemethod     

Intramolecular electron transfer in a cytochrome P450cam system with a site-specific branched structure

Cytochrome P450 (P450) is an attractive oxygenase due to the diverse catalytic reactions and the broad substrate specificity. Class I P450s require an excess concentration (more than 10 times) of iron-sulfur proteins, which transfer electrons to P450s, to attain the maximum catalytic activity and this requirement is a critical bottleneck for practical applications. Here, we show a site-specific branched fusion protein of P450 with its electron transfer proteins using enzymatic cross-linking with transglutaminase. A branched fusion protein of P450 from Pseudomonas putida (P450cam), which was composed of one molecule each of P450cam, putidaredoxin (Pdx) and Pdx reductase, showed higher catalytic activity (306 min(-1)) and coupling efficiency (99%) than the equimolar reconstitution system due to the intramolecular electron transfer. The unique site-specific branched structure simply increased local concentration of proteins without denaturation of each protein. Therefore, enzymatic post-translational protein manipulation can be a powerful alternative to conventional strategies for the creation of multicomponent enzyme systems with novel proteinaceous architecture.     

Molecular mechanics analysis of inhibitor binding to HIV-1 protease

Crystallographic structures of HIV protease with three differentpeptide-mimetic inhibitors were subjected to energy minimizationusing molecular mechanics, the minimized structures analyzedand the inhibitor binding energies calculated. Partial chargeassignment for the hydrogen bonded catalytic aspartk acids,Asp25 and -25', was in good agreement with charge calculationsusing semi-empirical molecular orbital methods. Root mean squaredeviations on minimization were small and similar for both subunitsin the protease dimer. The surface loops, which had the largestB factors, changed most on minimization; the hydrophobic coreand the inhibitor binding site showed little change. The distance-dependentdielectric of D(r) = 4r was found to be preferable to D(r) =r. Distance restraints were applied for the intermolecular hydrogenbonds to maintain the conformation of the inhibitor bindingsite. Using the dielectric of D(r) = 4r, the calculated interactionenergy of the three inhibitors with the protease ranged from–53 to –56 kcal/mol. The groups of the inhibitorswere changed to add or remove a ‘transition state analogue’hydroxyl group, and the loss in energy on the removal of thisgroup was calculated to be 0.9–1.7 kcal/mol. This wouldrepresent 19–36% of the total measured difference in bindingenergy between the inhibitors JG365 and MVT-101.     

Contribution of a disulfide bridge to the stability of 1,3-1,4-P-D-glucan 4-glucanohydrolase from Bacillus licheniformis

Bacillus 1,3-1,4-ß-glucanases possess a highly conserveddisulfide bridge connecting a ß-strand with a solventexposedloop lying on top of the extended binding site cleft The contributionof the disulfide bond and of both individual cysteines (Cys61and Cys90) in the Bacillus licheniformis enzyme to stabilityand activity has been evaluated by protein engineering methods.Reduction of the disulfide bond has no effect on kinetic parameters,has only a minor effect on the activity-temperature profileat high temperatures, and destabilizes the protein by less than0.7 kcal/mol as measured by equilibrium urea denatu ration at37°C. Replacing either of the Cys residues with Ala destabilizesthe protein and lowers the specific activity. C90A retains 70%of wild-type (wt) activity (in terms of Vmax), whereas C61Aand the double mutant C61A–C90A have 10% of wt V_max. Alarger change in free energy of unfolding is seen by equilibriumurea denaturation for the C61A mutation (loop residue, 3.2 kcal/molrelative to reduced wt) as compared with the C90A mutation (ß-strandresidue, 1.8 kcal/mol relative to reduced wt), while the doublemutant C61A–C90A is 0.8 kcal/mol less stable than thesingle C61A mutant. The effects on stability are interpretedas a result of the change in hydrophobic packing that occursupon removal of the sulfur atoms in the Cys to Ala mutations     

The nature of the ion binding interactions in EF-hand peptide analogs: free energy simulation of Asp to Asn mutations

The binding of the La³⁺ ion to a tridecapeptide, which is amodel for the EF-hand in calcium-binding proteins, is studiedhi solution by free energy simulations. The calculations analyzethe effect on the La³⁺ ion binding of the mutation of Asp toAsn for side chains that interact directly with the ion. Theresults are compared with the measurements of Marsden.B.J.,Hodges, R.S. and Sykes, B.D. (1989) Biochemistry, 28,8839, onthe same system. They found that the Asp to Asn mutation hasonly a small effect on the binding; the observed differencesin the free energies on changing one Asp to an Asn are between-0.3 and 1.8 kcal/ mol. This result is analyzed by alchemicalsimulations for the tridecapeptide in the bound Qoop) structureand free (extended) form. The free energy changes due to themutation of an Asp to an Asn are large and positive for boththe bound and free forms. However, since the values of the freeenergy changes are calculated to be similar hi the two forms,the difference in the binding free energy of Asp and Asn peptidesis found to be small, in agreement with experiment. By use ofthermodynamic integration, the various contributions to thefree energy changes are estimated. In the com-plexed form, theAsp to Asn mutation is favored by the reduction in the repulsiveinteraction with other charged residues of the peptide; it isdisfavored by the reduction of the stabilization of the ionand the surrounding water has a small effect. When the peptideadopts an extended conformation in the absence of the ion, themutation Asp to Asn is strongly disfavored by the interactionswith the water and is favored by the interactions within thepeptide. The results demonstrate the essential role of contributionsto the binding of EF-hands from interactions other than thosebetween the ion and the charged amino acid side chains. Theresults obtained from the simulations suggest, in accord withcrystal structures of La³⁺ bound to various ligands, that thecalcium-binding loop complexed with La³⁺ in solution has a significantlydifferent structure from that observed hi proteins.     

Free energy simulations of the HyHEL-10/HEL antibody-antigen complex

Free energy simulations are reported for the N31^L-D mutation,both in the HyHEL-10-HEL antibody-lysozyme complex and in theunliganded antibody, using the thermo-dynamic-cycle perturbationmethod. The present study suggests that the mutation would changethe free energy of binding of the complex by –5.6 kcal/mol(unrestrained free energy simulations), by –0.5 kcal/mol(free energy simulations with a restrained backbone) and by1.8 kcal/ mol (Poisson-Boltzmann calculations, which also usea restrained geometry model). A detailed structural analysishelps in estimating the contributions from various residuesand regions of the system. Enhanced recognition of HEL by themutant HyHEL-10 would arise from the combination of thermodynamicallymore favorable conformational changes of the CDR loops uponassociation and subsequent charge pairing with Lys96 in theantigen.     

Hydrophobic effects on protein/nucleic acid interaction: enhancement of substrate binding by mutating tyrosine 45 to tryptophan in ribonuclease Tl

Hydrophobic effects on binding of ribonuclease Tl to guaninebases of several ribonucleotides have been proved by mutatinga hydrophobic residue at the recognition site and by measuringthe effect on binding. Mutation of a hydrophobic surface residueto a more hydrophobic residue (Tyr45 – Trp) enhances thebinding to ribonucleotides, including mononucleotide inhibitorand product, and a synthetic substrate-analog trinudeotide aswell as the binding to dinucleotide substrates and RNA. Enhancementson binding to non-substrate ribonucleotides by the mutationhave been observed with free energy changes ranging from –2.2 to – 3 .9 kJ/mol. These changes are in good agreementwith that of substrate binding, –2.3 kJ/mol, which iscalculated from Michaelis constants obtained from kinetic studies.It is shown, by comparing the observed and calculated changesin binding free energy with differences in the observed transferfree energy changes of the amino acid side chains from organicsolvents to water, that the enhancement observed on guaninebinding comes from the difference in the hydrophobic effectsof the side chains of tyrosine and tryptophan. Furthermore,a linear relationship between nucleolytic activities and hydrophobicityof the residues (Ala, Phe, Tyr, Trp) at position 45 is observed.The mutation could not change substantially the base specificityof RNase Tl, which exhibits a prime requirement for guaninebases of substrates.     

Free energy perturbation studies on binding of A-74704 and its diester analog to HIV-1 protease

Free energy simulations have been employed to rationalize thebinding differences between A-74704, a pseudo C₂- symmetricinhibitor of HIV-1 protease and its diester analog. The diesteranalog inhibitor, which misses two hydrogen bonds with the enzymeactive site, is surprisingly only 10-fold weaker. The calculatedfree energy difference of 1.7 ± 0.6 kcal/mol is in agreementwith the experimental result. Further, the simulations showthat such a small difference in binding free energies is dueto (1) weaker hydrogen bond interactions between the two (P₁and P₁) NH groups of A-74704 with Gly27/Gly27' carbonyls ofthe enzyme and (2) the higher desolvation free energy of A-74704compared with its ester analog. The results of these calculationsand their implications for design of HIV-1 protease inhibitorsare discussed.     

Substrate mobility in thiocamphor-bound cytochrome P450cam: an explanation of the conflict between the observed product profile and the X-ray structure

Thiocamphor is an unusual substrate for P450_cam in that in theX-ray structure it binds in the active site pocket in two distinctorientations and neither of these orientations are consistentwith the 5-alcohol being the primary product. Other camphoranalogs such as norcamphor or camphane bind in a single orientationconsistent with the 5-alcohol being a major product. We presentan analysis of four 175 ps molecular dynamics trajectories ofthiocamphor-bound cytochrome P450_cam. The first two trajectorieswere calculated for cytochrome P450_cam with thiocamphor boundin both its major and minor crystallographic orientations. Inthe second set of simulations, a single oxygen atom was addedas a distal ligand to the heme group in order to model the putativeferryl oxygen reaction intermediate. Trajectories were againcalculated starting with thiocamphor in its major and minororientations. While the protein dynamics were quite similarin all four trajectories, the substrate showed distinctly differentmotions in each of the trajectories. In particular, the preferredsubstrate orientations were very different in the presence ofthe ferryl oxygen than in the absence of that oxygen. The preferredorientations in the absence of the distal oxygen were consistentwith the 3-akohol being the major product, while the preferredorientations in the presence of the distal oxygen were consistentwith the 5-alcohol being a major product. These simulationsoffer an explanation for the inconsistency between the X-raydata and the product profile.     

Construction and evaluation of a three-dimensional structure of cytochrome P450choP enzyme (CYP105C1)

The purpose of this work was to develop and carefully evaluateimproved strategies for constructing reliable 3-D models ofP450 isozymes. To this end, a unique combination of steps forbuilding and evaluating a model structure was used to builda homology model of the P450choP isozyme, based on knowledgeof the X-ray structures of P450cam, P450terp, P450BM-3 and P450eryF.Specifically, the reliability of this model was examined bysystematic comparisons of its conformational, energetic, environmentaland packing properties and those of the four reference proteinswith corresponding properties from the database of proteinswith known structures. The results showed that the examinedproperties of this model structure are well within the criteriaestablished for reliable structures and are of nearly as goodquality as those of the reference proteins. In addition, theresult from a 120 ps unconstrained MD simulation of the modelwith structural waters provided evidence that the model is stableat room temperature. This 3-D model can now be reliably usedfor explicit characterization of substrate and inhibitor complexes.Most importantly, although it is envisioned that building modelsfor mammalian P450s will be even more challenging, the stepsdescribed here should be very useful in future constructionof 3-D models of mammalian P450 isozymes.     

Ile115Leu mutation in the SRS1 region of an insect cytochrome P450 (CYP6B1) compromises substrate turnover via changes in a predicted product release channel

CYP6B1 represents the principal cytochrome P450 monooxygenase responsible for metabolizing furanocoumarins in Papilio polyxenes, an insect that specializes on host plants containing these toxins. Investigations of the amino acids responsible for the efficient metabolism of these plant toxins has identified Ile115 as one that modulates the rate of furanocoumarin metabolism even though it is predicted to be positioned at the edge of the heme plane and outside substrate contact regions. In contrast to previous expression studies conducted under conditions of limiting P450 reductase showing that the Ile115-to-Leu replacement enhances turnover of xanthotoxin and other furanocoumarins, studies conducted at high P450 reductase indicate that the Ile115-to-Leu replacement reduces turnover of these substrates. Further analysis of substrate binding affinities, heme spin state and NADPH consumption rates indicate that, whereas the I115L replacement mutant displays higher substrate affinity and heme spin state than the wild-type CYP6B1 protein, it utilizes NADPH more slowly than the wild-type CYP6B1 protein at high P450 reductase levels. Molecular models developed for the wild-type CYP6B1 and mutant protein suggest that more constricted channels extending from the catalytic site in the I115L mutant to the P450 surface limit the rate of product release from this mutant catalytic site under conditions not limited by the rate of electron transfer from NADPH.     

Contribution of a proline residue and a salt bridge to the stability of a type I reverse turn in chymotrypsin inhibitor-2

The contributions of the components of a type I reverse turnto the stability of chymotrypsin inhibitor-2 (Lys43-Pro44-Gly45)have been determined by protein engineering methods. A double-mutantcycle was used to determine the interaction between Lys43 andGlu45 by replacing them with alanine. We also mutated Pro44,which gives the geometry of the turn, to alanine and analysedthe stability of the resulting mutants compared with wild-typechymotrypsin inhibitor-2, using equilibrium denaturation inducedby guanldinium chloride. There are decreases in stability (inkcal/mol) of 0.64 = 0.06 for Lys43 - Ala, 0.57 ± 0.15for Glu45 - Ala, 0.95 ± 0.06 for Lys43 - Ala/Glu45 -Ala and 1.93 ± 0.09 for Pro44 - Ala. The free energyof interaction between Lys43 and Glu45 is calculated to be only0.25 ± 0.09 kcal/mol. From the changes in denaturationmidpoint, T_m measured by circular dkhroism, we estimate theenergy of interaction between Lys43 and Glu45 to be 0.36 ±0.07 kcal/mol whereas the contribution of Pro44 is -2.0 kcal/mol.The contribution of the salt bridge to the stability of theprotein is very small and the residue Pro44 plays the key rolein stabilizing the turn     

A Modified Arrhenius Approach to Thermodynamically Study Regioselectivity in Cytochrome P450-Catalyzed Substrate Conversion

The regio- (and stereo-)selectivity and specific activity of cytochrome P450s are determined by the accessibility of potential sites of metabolism (SOMs) of the bound substrate relative to the heme, and the activation barrier of the regioselective oxidation reaction(s). The accessibility of potential SOMs depends on the relative binding free energy (ΔΔG_bind) of the catalytically active substrate-binding poses, and the probability of the substrate to adopt a transition-state geometry. An established experimental method to measure activation energies of enzymatic reactions is the analysis of reaction rate constants at different temperatures and the construction of Arrhenius plots. This is a challenge for multistep P450-catalyzed processes that involve redox partners. We introduce a modified Arrhenius approach to overcome the limitations in studying P450 selectivity, which can be applied in multiproduct enzyme catalysis. Our approach gives combined information on relative activation energies, ΔΔG_bind values, and collision entropies, yielding direct insight into the basis of selectivity in substrate conversion.