Trinuclear Metal Clusters in Catalysis by Terpenoid Synthases

Terpenoid synthases are ubiquitous enzymes that catalyze the formation of structurally and stereochemically diverse isoprenoid natural products. Many isoprenoid coupling enzymes and terpenoid cyclases from bacteria, fungi, protists, plants, and animals share the class I terpenoid synthase fold. Despite generally low amino acid sequence identity among these examples, class I terpenoid synthases contain conserved metal binding motifs that coordinate to a trinuclear metal cluster. This cluster not only serves to bind and orient the flexible isoprenoid substrate in the precatalytic Michaelis complex, but it also triggers the departure of the diphosphate leaving group to generate a carbocation that initiates catalysis. Additional conserved hydrogen bond donors assist the metal cluster in this function. Crystal structure analysis reveals that the constellation of three metal ions required for terpenoid synthase catalysis is generally identical among all class I terpenoid synthases of known structure.     

Quantum chemical calculations of the NHA bound nitric oxide synthase active site: O2 binding and implications for the catalytic mechanism

Density functional theory has been employed to model the binding of the intermediate substrate NHA, by nitric oxide synthases. In particular, the orientation and interactions of possibly catalytically important substrate hydrogens, with and without molecular oxygen bound to the active site heme group, are considered. Without O(2), three possible conformers have been found, with the energetically most favored structure being that in which both protons of the -NHOH moiety of NHA are directed toward the heme group. With oxygen bound, four different structures were found. The energetically lowest structure is again found to have both hydrogens of the -NHOH group pointing toward the heme group, thus forming hydrogen bonds between -NH- and the terminal oxygen, and between -OH and the inner oxygen of the heme-O(2) group. In addition, unprotonated structures of the substrate bound to the active site are considered and the proton affinity calculated.     

Competitive inhibition of aristolochene synthase by phenyl-substituted farnesyl diphosphates: evidence of active site plasticity

Analogues of farnesyl diphosphate (FPP, ) containing phenyl substituents in place of methyl groups have been prepared in syntheses that feature use of a Suzuki-Miyaura reaction as a key step. These analogues were found not to act as substrates of the sesquiterpene cyclase aristolochene synthase from Penicillium roqueforti (AS). However, they were potent competitive inhibitors of AS with K(I)-values ranging from 0.8 to 1.2 microM. These results indicate that the diphosphate group contributes the largest part to the binding of the substrate to AS and that the active sites of terpene synthases are sufficiently flexible to accommodate even substrate analogues with large substituents suggesting a potential way for the generation of non-natural terpenoids. Molecular mechanics simulations of the enzyme bound inhibitors suggested that small changes in orientations of active site residues and subtle alterations of the conformation of the backbones of the inhibitors are sufficient to accommodate the phenyl-farnesyl-diphosphates.     

Prostacyclin and Thromboxane Synthase: New Aspects of Hemethiolate Catalysis

The arachidonic acid metabolites thromboxane A₂ and prostacyclin are highly potent regulators of cell physiology. They are both formed by enzymatic rearrangement of the 9,11-epidioxy prostaglandin H₂ catalyzed, however, by thromboxane and prostacyclin synthase, respectively. The two enzymes have been isolated, sequenced, and characterized as hemethiolate (“P450”) enzymes. The different isomerization products can be explained on the same catalytic principle by a different ligation of the heme centers with the two epidioxy oxygens atoms. This requires different conformations for substrate binding at the active site, which is substantiated by the different inhibitors and amino acid sequences of the enzymes. In a hypothesis which has mechanistic principles in common with the P450-monooxygenases and the allene oxide synthases, oxy radicals are formed first and rearrange to carbon radicals. These could then rapidly be converted into carbocations by the ferrylthiolate or iron(III )thiyl structures formed as intermediates.     

Competition of ester, amide, ether, carbonate, alcohol and epoxide ligands in the dirhodium experiment (chiral discrimination by NMR spectroscopy)

Thirty-eight derivatives of 3-hydroxy-2-methylpropanoic acid, each with two different oxygen functionalities, were synthesized and subjected to the standard dirhodium experiment (1H NMR in the presence of an equimolar amount of the chiral dirhodium tetracarboxylate complex Rh*). Their structures represent ester, amide, carbonate, ether, alcohol and/or epoxy groups. Significant selectivity in the binding of those oxygen groups to the complex were determined. From these results, a priority list in binding to a rhodium atom of Rh* was established: epoxides > primary alcohols > ethers > or = esters > or = amides > carbonates > tertiary alcohols. This sequence allows the prediction of the preferred binding site of oxygen-containing groups in polyfunctional compounds, which frequently occur among natural products, and, particularly, in asymmetric synthesis of such compounds. Differentiation of the enantiomers by the dirhodium experiment is easily accomplished due to numerous signal dispersions in nearly all cases.     

Sesquiterpene synthases: passive catalysts or active players?

Sesquiterpene synthases catalyse the metal dependent turnover of farnesyl diphosphate to generate a class of natural products characterised by an enormous diversity in structure, stereochemistry, biological function and application. It has been proposed that these enzymes take a passive role in the reactions they catalyse and that they serve mostly as stereochemical templates, within which the reactions take place. Here, recent research into the structure and function of sesquiterpene synthases and the mechanisms of the reactions that they catalyse will be reviewed to suggest that these fascinating enzymes play multifaceted active roles in what are arguably the most complex biosynthetic reactions.     

Targeting the Substrate Binding Site of E. coli Nitrile Reductase QueF by Modeling,Substrate and Enzyme Engineering

Nitrile reductase QueF catalyzes the reduction of 2‐amino‐5‐cyanopyrrolo[2,3‐d]pyrimidin‐4‐one (preQ₀) to 2‐amino‐5‐aminomethylpyrrolo[2,3‐d]pyrimidin‐4‐one (preQ₁) in the biosynthetic pathway of the hypermodified nucleoside queuosine. It is the only enzyme known to catalyze a reduction of a nitrile to its corresponding primary amine and could therefore expand the toolbox of biocatalytic reactions of nitriles. To evaluate this new oxidoreductase for application in biocatalytic reactions, investigation of its substrate scope is prerequisite. We report here an investigation of the active site binding properties and the substrate scope of nitrile reductase QueF from Escherichia coli. Screenings with simple nitrile structures revealed high substrate specificity. Consequently, binding interactions of the substrate to the active site were identified based on a new homology model of E. coli QueF and modeled complex structures of the natural and non‐natural substrates. Various structural analogues of the natural substrate preQ₀ were synthesized and screened with wild‐type QueF from E. coli and several active site mutants. Two amino acid residues Cys190 and Asp197 were shown to play an essential role in the catalytic mechanism. Three non‐natural substrates were identified and compared to the natural substrate regarding their specific activities by using wild‐type and mutant nitrile reductase.     

Increasing complexity of a diterpene synthase reaction with a single residue switch

Terpene synthases often catalyze complex reactions involving intricate series of carbocation intermediates. The resulting, generally cyclical, structures provide initial hydrocarbon frameworks that underlie the astonishing structural diversity of the enormous class of terpenoid natural products (>50,000 known), and these enzymes often mediate the committed step in their particular biosynthetic pathway. Accordingly, how terpene synthases specify product outcome has drawn a great deal of attention. In previous work, we have shown that mutational introduction of a hydroxyl group at specific positions within diterpene synthase active sites can "short circuit" complex cyclization and/or rearrangement reactions, resulting in the production of "simpler"' diterpenes. Here we demonstrate that the converse change, substitution of an Ile for Thr at the relevant position in a native pimaradiene synthase, leads to a dramatic increase in reaction complexity. Product outcome is shifted from the tricyclic pimaradiene to a rearranged tetracycle, aphidicol-15-ene. Thus, the nature of the residue at this position acts as a true switch for product outcome. In addition, the ability of aliphatic residue substitution to enable a more complex reaction emphasizes the importance of substrate conformation imposed by a largely inert active site. Furthermore, the profound plasticity of diterpene synthases exemplified by this single residue switch for product outcome is consistent with the screening/diversity-oriented hypothesis of natural products metabolism.     

Study of IspH, a key enzyme in the methylerythritol phosphate pathway using fluoro-substituted substrate analogues

IspH, a [4Fe-4S]-cluster-containing enzyme, catalyzes the reductive dehydroxylation of 4-hydroxy-3-methyl-butenyl diphosphate (HMBPP) to isopentenyl diphosphate (IPP) and dimethylallyl diphosphate (DMAPP) in the methylerythritol phosphate pathway. Studies of IspH using fluoro-substituted substrate analogues to dissect the contributions of several factors to IspH catalysis, including the coordination of the HMBPP C(4)-OH group to the iron-sulfur cluster, the H-bonding network in the active site, and the electronic properties of the substrates, are reported.     

A negative cooperativity mechanism of human CYP2E1 inferred from molecular dynamics simulations and free energy calculations

Human cytochrome P450 2E1 (CYP2E1) participates in the metabolism of over 2% of all the oral drugs. A hallmark peculiar feature of this enzyme is that it exhibits a pronounced negative cooperativity in substrate binding. However the mechanism by which the negative cooperativity occurs is unclear. Here, we performed molecular dynamics simulations and free energy calculations on human CYP2E1 to examine the structural differences between the substrate-free and the enzymes with one and two aniline molecules bound. Our results indicate that although the effector substrate does not bind in the active site cavity, it still can directly interact with the active site residues of human CYP2E1. The interaction of the effector substrate with the active site leads to a reorientation of active site residues, which thereby weakens the interactions of the active substrate with this site. We also identify a conserved residue T303 that plays a crucial role in the negative cooperative binding on the short-range effects. This residue is a key factor in the positioning of substrates and in proton delivery to the active site. Additionally, a long-range effect of the effector substrate is identified in which F478 is proposed to play a key role. As located in the interface between the active and effector sites, this residue structurally links the active and effector sites and is found to play a significant role in affecting substrate access and ligand positioning within the active site. In the negative cooperative binding, this residue can decrease the interactions of the active substrate with the active site by π-π stacking which then lowers the hydroxylation activity for the active substrate. These findings are in agreement with previous experimental observations and thus provide detailed atomistic insight into the poorly understood mechanism of the negative cooperativity in human CYP2E1.     

Density functional study of the interaction between small Au clusters, Au(n) (n=1-7) and the rutile TiO2 surface. II. Adsorption on a partially reduced surface

We use density functional theory to examine the electronic structure of small Au(n) (n=1-7) clusters, supported on a rutile TiO(2)(110) surface having oxygen vacancies on the surface (a partially reduced surface). Except for the monomer, the binding energy of all Au clusters to the partially reduced surface is larger by approximately 0.25 eV than the binding energy to a stoichiometric surface. The bonding site and the orientation of the cluster are controlled by the shape of the highest occupied molecular orbitals (HOMOs) of the free cluster (free cluster means a gas-phase cluster with the same geometry as the supported one). The bond is strong when the lobes of the HOMOs overlap with those of the high-energy states of the clean oxide surface (i.e., with no gold) that have lobes on the bridging and the in-plane oxygen atoms. In other words, the cluster takes a shape and a location that optimizes the contact of its HOMOs with the oxygen atoms. Fivefold coordinated Ti atoms located at a defect site (5c-Ti(*)) participate in the binding only when a protruding lobe of the singly occupied molecular orbital (for odd n) or the lowest unoccupied molecular orbital (for even n) of the free Au(n) cluster points toward a 5c-Ti(*) atom. The oxygen vacancy influences the binding energy of the clusters (except for Au(1)) only when they are in direct contact with the defect. The desorption energy and the total charge on clusters that are close to, but do not overlap with, the vacancy differ little from the values they have when the cluster is adsorbed on a stoichiometric surface. The behavior of Au(1) is rather remarkable. The atom prefers to bind directly to the vacancy site with a binding energy of 1.81 eV. However, it also makes a strong bond (1.21 eV) with any 5c-Ti atom even if that atom is far from the vacancy site. In contrast, the binding of a Au monomer to the 5c-Ti atom of a surface without vacancies is weak (0.45 eV). The presence of the vacancy activates the 5c-Ti atoms by populating states at the bottom of the conduction band. These states are delocalized and have lobes protruding out of the surface at the location of the 5c-Ti atoms. It is the overlap of these lobes with the highest orbital of the Au atom that is the major reason for the bonding to the 5c-Ti atom, no matter how far the latter is from the vacancy. The energy for breaking an adsorbed cluster into two adsorbed fragments is smaller than the kinetic energy of the mass-selected clusters deposited on the surface in experiments. However, this is not sufficient for breaking the cluster upon impact with the surface, since only a fraction of the available energy will go into the reaction coordinate for breakup.     

Effects of slow substrate binding and release in redox enzymes: theory and application to periplasmic nitrate reductase

For redox enzymes, the technique called protein film voltammetry makes it possible to determine the entire profile of activity against driving force by having the enzyme exchanging directly electrons with the rotating-disc electrode onto which it is adsorbed. Both the potential location of the catalytic response and its detailed shape report on the sequence of catalytic events, electron transfers and chemical steps, but the models that have been used so far to decipher this signal lack generality. For example, it was often proposed that substrate binding to multiple redox states of the active site may explain that turnover is greater in a certain window of electrode potential, but no fully analytical treatment has been given. Here, we derive (i) the general current equation for the case of reversible substrate binding to any redox states of a two-electron active site (as exemplified by flavins and Mo cofactors), (ii) the quantitative conditions for an extremum in activity to occur, and (iii) the expressions from which the substrate-concentration dependence of the catalytic potential can be interpreted to learn about the kinetics of substrate binding and how this affects the reduction potential of the active site. Not only does slow substrate binding and release make the catalytic wave shape highly complex, but we also show that it can have important consequences which will escape detection in traditional experiments: the position of the wave (this is the driving force that is required to elicit catalysis) departs from the reduction potential of the active site even at the lowest substrate concentration, and this deviation may be large if substrate binding is irreversible. This occurs in the reductive half-cycle of periplasmic nitrate reductase where irreversibility lowers the driving force required to reduce the active site under turnover conditions and favors intramolecular electron transfer from the proximal [4Fe4S]+ cluster to the active site Mo(V).     

Genome-wide identification,functional and evolutionary analysis of terpene synthases in pineapple

Terpene synthases (TPSs) are vital for the biosynthesis of active terpenoids, which have important physiological, ecological and medicinal value. Although terpenoids have been reported in pineapple (Ananas comosus), genome-wide investigations of the TPS genes responsible for pineapple terpenoid synthesis are still lacking. By integrating pineapple genome and proteome data, twenty-one putative terpene synthase genes were found in pineapple and divided into five subfamilies. Tandem duplication is the cause of TPS gene family duplication. Furthermore, functional differentiation between each TPS subfamily may have occurred for several reasons. Sixty-two key amino acid sites were identified as being type-II functionally divergence between TPS-a and TPS-c subfamily. Finally, coevolution analysis indicated that multiple amino acid residues are involved in coevolutionary processes. In addition, the enzyme activity of two TPSs were tested. This genome-wide identification, functional and evolutionary analysis of pineapple TPS genes provide a new insight into understanding the roles of TPS family and lay the basis for further characterizing the function and evolution of TPS gene family.     

Replacement of Water Molecules in a Phosphate Binding Site by Furanoside‐Appended lin‐Benzoguanine Ligands of tRNA‐Guanine Transglycosylase (TGT)

The enzyme tRNA‐guanine transglycosylase has been identified as a drug target for the foodborne illness shigellosis. A key challenge in structure‐based design for this enzyme is the filling of the polar ribose‐34 pocket. Herein, we describe a novel series of ligands consisting of furanoside‐appended lin‐benzoguanines. They were designed to replace a conserved water cluster and differ by the functional groups at C(2) and C(3) of the furanosyl moiety being either OH or OMe. The unfavorable desolvation of Asp102 and Asp280, which are located close to the ribose‐34 pocket, had a significant impact on binding affinity. While the enzyme has tRNA as its natural substrate, X‐ray co‐crystal structures revealed that the furanosyl moieties of the ligands are not accommodated in the tRNA ribose‐34 site, but at the location of the adjacent phosphate group. A remarkable similarity of the position of the oxygen atoms in these two structures suggests furanosides as a potential phosphate isoster.     

Novel Product Chemistry from Mechanistic Analysis of ent‐Copalyl Diphosphate Synthases from Plant Hormone Biosynthesis

An active‐site water molecule coordinated by conserved histidine and asparagine residues seems to serve as the catalytic base in all ent‐copalyl diphosphate synthases (CPSs). When these residues are substituted by alanine, the mutant CPSs produce stereochemically novel ent‐8‐hydroxy‐CPP. Given the requisite presence of CPSs in all land plants for gibberellin phytohormone biosynthesis, such plasticity presumably underlies the observed extensive diversification of the resulting labdane‐related diterpenoids.     

Ruptenes: A Family of Terpene Analogs Give Insight into Cyclisation Mechanisms by Cascade Disruption

The terpenoid substrate analogs (7R)-6,7-dihydrogeranylgeranyl diphosphate (6,7-dihydro-GGPP) and (7R)-6,7-dihydrogeranylfarnesyl diphosphate (6,7-dihydro-GFPP) were synthesised from (S)-citronellol and enzymatically converted with nine diterpene and two sesterterpene synthases, respectively. In two cases the substrate analogs were converted into diterpenes in cyclisation reactions corresponding to those observed for the native substrate GGPP, while the cyclisation cascade was disrupted or redirected in the other nine cases, leading to products that were named ruptenes. Several of the isolated ruptenes represent deprotonation products of cationic intermediates that are analogs of the intermediates proposed along the cyclisation cascades for the native substrates GGPP or GFPP, thus giving insights into the complex reaction mechanisms of terpene synthase mediated biosynthesis.     

What Determines the Selectivity of Arginine Dihydroxylation by the Nonheme Iron Enzyme OrfP?

The nonheme iron enzyme OrfP reacts with l -Arg selectively to form the 3R,4R-dihydroxyarginine product, which in mammals can inhibit the nitric oxide synthase enzymes involved in blood pressure control. To understand the mechanisms of dioxygen activation of l -Arg by OrfP and how it enables two sequential oxidation cycles on the same substrate, we performed a density functional theory study on a large active site cluster model. We show that substrate binding and positioning in the active site guides a highly selective reaction through C³−H hydrogen atom abstraction. This happens despite the fact that the C³−H and C⁴−H bond strengths of l -Arg are very similar. Electronic differences in the two hydrogen atom abstraction pathways drive the reaction with an initial C³−H activation to a low-energy ⁵σ-pathway, while substrate positioning destabilizes the C⁴−H abstraction and sends it over the higher-lying ⁵π-pathway. We show that substrate and monohydroxylated products are strongly bound in the substrate binding pocket and hence product release is difficult and consequently its lifetime will be long enough to trigger a second oxygenation cycle.     

Molecular interactions of nitrogen-containing bisphosphonates within farnesyl diphosphate synthase

Bisphosphonates, known for their effectiveness in the treatment of osteoporosis, inhibit bone resorption via mechanisms that involve binding to bone mineral and cellular effects on osteoclasts. The major molecular target of nitrogen-containing bisphosphonates (N-BPs) in osteoclasts is farnesyl diphosphate synthase (FPPS). N-BPs likely inhibit this enzyme by mimicking one or more of the natural isoprenoid lipid substrates (GPP/DMAPP and IPP) but the mode of inhibition is not established. The active site of FPPS comprises a subsite for each substrate. Kinetic studies with recombinant human FPPS indicate that both potent (risedronate) and weak (NE-58051) enzyme inhibitors compete with GPP for binding to FPPS, however, binding to this site does not completely explain the difference in potency of the two inhibitors, suggesting that a second binding site may also be a target of bisphosphonate inhibition. Using the docking software suite Autodock, we explored a dual inhibitor binding mode for recombinant human FPPS. Experimental support for dual binding is suggested by Dixon plots for the inhibitors. N-BPs may inhibit by binding to both the GPP and a second site with differences in potency at least partly arising from inhibition at the second site.     

Modulation of substrate binding to naphthalene 1,2-dioxygenase by rieske cluster reduction/oxidation

The active site of the oxygenase component of naphthalene 1,2-dioxygenase (NDO) contains a Rieske Fe-S cluster and a mononuclear non-heme iron, which are contributed by different alpha-subunits in the (alphabeta)(3) structure. The enzyme catalyzes cis-dihydroxylation of aromatic substrates in addition to numerous other adventitious oxidation reactions. High-resolution Mims (2)H-ENDOR spectra have been recorded for the NO-ferrous center of NDO bound with d(8)-naphthalene and d(2)-naphthalene; spectra were collected for the enzyme with the Rieske diiron center both in its oxidized and in its reduced states. A sharp quartet ENDOR pattern from a nearby deuteron of substrate was detected for each substrate. Examination of the sample prepared with 1,4-dideutero-naphthalene shows that the signal arises from D1. The ENDOR data place D1 at a distance of ca. 4.4 A from the mononuclear Fe and with the Fe-D vector being roughly along the Fe-N(O) direction. Because reduction of the Rieske cluster is required for O(2) binding and subsequent catalysis, the effect of its oxidation state on substrate binding was examined. The spectra from the NDO-naphthalene complex reveal two different binding conformations, which change in relative population when the oxidation state of the Rieske cluster is changed. This shift, and the conformational coupling it implies, may hold the key to both oxygen gating and oxygen reactivity for Rieske aromatic dioxygenases.