Mechanistic Studies of the Protonation-Deprotonation Reactions for Type 1 and Type 2 Isopentenyl Diphosphate:Dimethylallyl Diphosphate Isomerase.

Type 1 and type 2 isopentenyl diphosphate:dimethylallyl diphosphate isomerase (IDI-1 and IDI-2) catalyze the interconversion of isopentenyl diphosphate (IPP) and dimethylallyl diphosphate (DMAPP), the fundamental building blocks for biosynthesis of isoprenoid compounds. Previous studies indicate that both isoforms of IDI catalyze isomerization by a protonation-deprotonation mechanism. IDI-1 and IDI-2 are "sluggish" enzymes with turnover times of ∼10 s-1 and ∼1 s-1, respectively. We measured incorporation of deuterium into IPP and DMAPP in D2O buffer for IDI-1 and IDI-2 under conditions where newly synthesized DMAPP is immediately and irreversibly removed by coupling its release to condensation with l-tryptophan catalyzed by dimethylallyltrytophan synthase. During the course of the reactions, we detected formation of d1, d2, and d3 isotopologues of IPP and DMAPP, which were formed during up to five isomerizations between IPP and DMAPP during each turnover. The patterns for deuterium incorporation into IPP show that d2-IPP is formed in preference to d1-IPP for both enzymes. Analysis of the patterns of deuterium incorporation are consistent with a mechanism involving addition and removal of protons by a concerted asynchronous process, where addition substantially precedes removal, or a stepwise process through a short-lived (<3 ps) tertiary carbocationic intermediate. Previous work with mechanism-based inhibitors and related model studies supports a concerted asynchronous mechanism for the enzyme-catalyzed isomerizations.

Structure-Function Studies of Artemisia tridentata Farnesyl Diphosphate Synthase and Chrysanthemyl Diphosphate Synthase by Site-Directed Mutagenesis and Morphogenesis.

The amino acid sequences of farnesyl diphosphate synthase (FPPase) and chrysanthemyl diphosphate synthase (CPPase) from Artemisia tridentata ssp. Spiciformis, minus their chloroplast targeting regions, are 71% identical and 90% similar. FPPase efficiently and selectively synthesizes the "regular" sesquiterpenoid farnesyl diphosphate (FPP) by coupling isopentenyl diphosphate (IPP) to dimethylallyl diphosphate (DMAPP) and then to geranyl diphosphate (GPP). In contrast, CPPase is an inefficient promiscuous enzyme, which synthesizes the "irregular" monoterpenes chrysanthemyl diphosphate (CPP), lavandulyl diphosphate (LPP), and trace quantities of maconelliyl diphosphate (MPP) from two molecules of DMAPP, and couples IPP to DMAPP to give GPP. A. tridentata FPPase and CPPase belong to the chain elongation protein family (PF00348), a subgroup of the terpenoid synthase superfamily (CL0613) whose members have a characteristic α terpene synthase α-helical fold. The active sites of A. tridentata FPPase and CPPase are located within a six-helix bundle containing amino acids 53 to 241. The two enzymes were metamorphosed into one another by sequentially replacing the loops and helices of the six-helix bundle from enzyme with those from the other. Chain elongation was the dominant activity during the N-terminal to C-terminal metamorphosis of FPPase to CPPase, with product selectivity gradually switching from FPP to GPP, until replacement of the final α-helix, whereupon cyclopropanation and branching activity competed with chain elongation. During the corresponding metamorphosis of CPPase to FPPase, cyclopropanation and branching activities were lost upon replacement of the first helix in the six-helix bundle. Mutations of active site residues in CPPase to the corresponding amino acids in FPPase enhanced chain-elongation activity, while similar mutations in the active site of FPPase failed to significantly promote formation of significant amounts of irregular monoterpenes. Our results indicate that CPPase, a promiscuous enzyme, is more plastic toward acquiring new activities, whereas FPPase is more resistant. Mutations of residues outside of the α terpene synthase fold are important for acquisition of FPPase activity for synthesis of CPP, LPP, and MPP.

Absolute Configuration of Hydroxysqualene. An Intermediate in Bacterial Hopanoid Biosynthesis.

Squalene (SQ) is a key intermediate in hopanoid biosynthesis. Many bacteria synthesize SQ from farnesyl diphosphate (FPP) in three steps: FPP to (1R,2R,3R)-presqualene diphosphate (PSPP), (1R,2R,3R)-PSPP to hydroxysqualene (HSQ), and HSQ to SQ. Chemical, biochemical, and spectroscopic methods were used to establish that HSQ synthase synthesizes (S)-HSQ. In contrast, eukaryotic squalene synthase catalyzes solvolysis of (1R,2R,3R)-PSPP to give (R)-HSQ. The bacterial enzyme that reduces HSQ to SQ does not accept (R)-HSQ as a substrate.

Biosynthesis of Squalene from Farnesyl Diphosphate in Bacteria: Three Steps Catalyzed by Three Enzymes.

Squalene (SQ) is an intermediate in the biosynthesis of sterols in eukaryotes and a few bacteria and of hopanoids in bacteria where they promote membrane stability and the formation of lipid rafts in their hosts. The genes for hopanoid biosynthesis are typically located on clusters that consist of four highly conserved genes-hpnC, hpnD, hpnE, and hpnF-for conversion of farnesyl diphosphate (FPP) to hopene or related pentacyclic metabolites. While hpnF is known to encode a squalene cyclase, the functions for hpnC, hpnD, and hpnE are not rigorously established. The hpnC, hpnD, and hpnE genes from Zymomonas mobilis and Rhodopseudomonas palustris were cloned into Escherichia coli, a bacterium that does not contain genes homologous to hpnC, hpnD, and hpnE, and their functions were established in vitro and in vivo. HpnD catalyzes formation of presqualene diphosphate (PSPP) from two molecules of FPP; HpnC converts PSPP to hydroxysqualene (HSQ); and HpnE, a member of the amine oxidoreductase family, reduces HSQ to SQ. Collectively the reactions catalyzed by these three enzymes constitute a new pathway for biosynthesis of SQ in bacteria.

Dependence of the product chain-length on detergents for long-chain E-polyprenyl diphosphate synthases.

Long-chain E-polyprenyl diphosphate synthases (E-PDS) catalyze repetitive addition of isopentenyl diphosphate (IPP) to the growing prenyl chain of an allylic diphosphate. The polyprenyl diphosphate products are required for the biosynthesis of ubiquinones and menaquinones required for electron transport during oxidative phosphorylation to generate ATP. In vitro, the long-chain PDSs require addition of phospholipids or detergents to the assay buffer to enhance product release and maintain efficient turnover. During preliminary assays of product chain-length with anionic, zwitterionic, and nonionic detergents, we discovered considerable variability. Examination of a series of nonionic PEG detergents with several long-chain E-PDSs from different organisms revealed that in vitro incubations with nonaethylene glycol monododecyl ether or Triton X-100 typically gave chain-lengths that corresponded to those of the isoprenoid moieties in respiratory quinones synthesized in vivo. In contrast, incubations in buffer with n-butanol, CHAPS, DMSO, n-octyl-ß-glucopyranoside, or ß-cyclodextrin or in buffer without detergent typically proceeded more slowly and gave a broad range of chain-lengths.

Prediction of function for the polyprenyl transferase subgroup in the isoprenoid synthase superfamily.

The number of available protein sequences has increased exponentially with the advent of high-throughput genomic sequencing, creating a significant challenge for functional annotation. Here, we describe a large-scale study on assigning function to unknown members of the trans-polyprenyl transferase (E-PTS) subgroup in the isoprenoid synthase superfamily, which provides substrates for the biosynthesis of the more than 55,000 isoprenoid metabolites. Although the mechanism for determining the product chain length for these enzymes is known, there is no simple relationship between function and primary sequence, so that assigning function is challenging. We addressed this challenge through large-scale bioinformatics analysis of >5,000 putative polyprenyl transferases; experimental characterization of the chain-length specificity of 79 diverse members of this group; determination of 27 structures of 19 of these enzymes, including seven cocrystallized with substrate analogs or products; and the development and successful application of a computational approach to predict function that leverages available structural data through homology modeling and docking of possible products into the active site. The crystallographic structures and computational structural models of the enzyme-ligand complexes elucidate the structural basis of specificity. As a result of this study, the percentage of E-PTS sequences similar to functionally annotated ones (BLAST e-value ≤ 1e(-70)) increased from 40.6 to 68.8%, and the percentage of sequences similar to available crystal structures increased from 28.9 to 47.4%. The high accuracy of our blind prediction of newly characterized enzymes indicates the potential to predict function to the complete polyprenyl transferase subgroup of the isoprenoid synthase superfamily computationally.

Mutagenesis of isopentenyl phosphate kinase to enhance geranyl phosphate kinase activity.

Isopentenyl phosphate kinase (IPK) catalyzes the ATP-dependent phosphorylation of isopentenyl phosphate (IP) to form isopentenyl diphosphate (IPP) during biosynthesis of isoprenoid metabolites in Archaea. The structure of IPK from the archeaon Thermoplasma acidophilum (THA) was recently reported and guided the reconstruction of the IP binding site to accommodate the longer chain isoprenoid monophosphates geranyl phosphate (GP) and farnesyl phosphate (FP). We created four mutants of THA IPK with different combinations of alanine substitutions for Tyr70, Val73, Val130, and Ile140, amino acids with bulky side chains that limited the size of the side chain of the isoprenoid phosphate substrate that could be accommodated in the active site. The mutants had substantially increased GP kinase activity, with 20-200-fold increases in k(cat)(GP) and 30-130-fold increases in k(cat)(GP)/K(M)(GP) relative to those of wild-type THA IPK. The mutations also resulted in a 10(6)-fold decrease in k(cat)(IP)/K(M)(IP) compared to that of wild-type IPK. No significant change in the kinetic parameters for the cosubstrate ATP was observed, signifying that binding between the nucleotide binding site and the IP binding site was not cooperative. The shift in substrate selectivity from IP to GP, and to a lesser extent, FP, in the mutants could act as a starting point for the creation of more efficient GP or FP kinases whose products could be exploited for the chemoenzymatic synthesis of radiolabeled isoprenoid diphosphates.

Type II isopentenyl diphosphate isomerase: probing the mechanism with alkyne/allene diphosphate substrate analogues.

Isopentenyl diphosphate isomerase (IDI) catalyzes the interconversion of isopentenyl diphosphate (IPP) and dimethylallyl diphosphate (DMAPP), the basic five-carbon building blocks of isoprenoid molecules. Two structurally unrelated classes of IDIs are known. Type I IPP isomerase (IDI-1) utilizes a divalent metal in a protonation-deprotonation reaction. In contrast, the type II enzyme (IDI-2) requires reduced flavin, raising the possibility that the reaction catalyzed by IDI-2 involves the net addition or abstraction of a hydrogen atom. As part of our studies of the mechanism of isomerization for IDI-2, we synthesized allene and alkyne substrate analogues for the enzyme. These molecules are predicted to be substantially less reactive toward proton addition than IPP and DMAPP but have similar reactivities toward hydrogen atom addition. This prediction was verified by calculations of gas-phase heats of reaction for addition of a proton and of a hydrogen atom to 1-butyne (3) and 1,2-butadiene (4) to form the 1-buten-2-yl carbocation and radical, respectively, and related affinities for 2-methyl-1-butene (5) and 2-methyl-2-butene (6) using G3MP2B3 and CBS-QB3 protocols. Alkyne 1-OPP and allene 2-OPP were not substrates for Thermus thermophilus IDI-2 or Escherichia coli IDI-1 but instead were competitive inhibitors. The experimental and computational results are consistent with a protonation-deprotonation mechanism for the enzyme-catalyzed isomerization of IPP and DMAPP.

Electronic circular dichroism of monomethyl [16O,17O,18O]-phosphate and [16O,17O,18O]-thiophosphate revisited.

Phosphoryl-transfer reactions have long been of interest due to their importance in maintaining numerous cellular functions. A phosphoryl-transfer reaction results in two possible stereochemical outcomes: either retention or inversion of configuration at the transferred phosphorus atom. When the product is phosphate, isotopically-labeled [16O,17O,18O]-phosphate derivatives can be used to distinguish these outcomes; one oxygen must be replaced by sulfur or esterified to achieve isotopic chirality. Conventionally, stereochemical analysis of isotopically chiral phosphate has been based on 31P NMR spectroscopy and involves complex chemical or enzymatic transformations. An attractive alternative would be direct determination of the enantiomeric excess using chiroptical spectroscopy. (S)-Methyl-[16O,17O,18O]-phosphate (MePi*), 7 and enantiomeric [16O,17O,18O]-thiophosphate (TPi*), 10, were previously reported to exhibit weak electronic circular dichroism (ECD), although with 10 the result was considered to be uncertain. We have now re-examined the possibility that excesses of 7 and 10 enantiomers can be detected by ECD spectrometry, using both experimental and theoretical approaches. 7 and both the (R) and (S) enantiomers of 10 (10a,10b) were synthesized by the 'Oxford route' and characterized by 1H, 31P and 17O NMR, and by MS analysis. Weak ECD could be found for 7, with suboptimal S/N. No significant ECD could be detected for the 10 enantiomers. Time-dependent DFT (TDDFT) calculations of the electronic excitation energies and rotational strengths of the same three enantiomers were carried out using the functional B3LYP and the basis set 6-311G**. The isotopically-perturbed geometries were predicted using the anharmonic vibrational frequency calculational code in GAUSSIAN 03. In the case of 10, calculations were also carried out for the hexahydrated complex to investigate the influence of the aqueous solvent. The predicted excitation wavelengths are greater than the observed wavelengths, a not unusual result of TDDFT calculations. The predicted anisotropy ratios are 2.9 x 10(-5) for 7, -5.3 x 10(-6) for 10a/b, and 1.7 x 10(-6) for 10a/b.(H2O)6. For 7 the predicted anisotropy ratio approximates that observed in this work, 4.5 x 10(-5) at 208 nm. For 10a/b, the upper limits of the experimental anisotropy ratios (<5 x 10(-6) at 225 nm, pH 9; <5 x 10(-6) at 236 nm, pH 12) are comparable to the predicted magnitude of the value for 10a/b. The lower predicted value for 10a/b.(H2O)6 suggests that the aqueous environment affects the ECD significantly. Altogether, the TDDFT calculations together with a stereochemical analysis based on NMR and the MS data support the conclusion that the experimental ECD results for MePi* and TPi* may be reliable in order of magnitude.

Recombinant squalene synthase. synthesis of cyclopentyl non-head-to-tail triterpenes.

Incubation of farnesyl diphosphate with recombinant yeast squalene synthase in the absence of NADPH gives a mixture of triterpene hydrocarbons and alcohols, including botryococcene-like compounds with 1'-3 linkages between the farnesyl units. One of these molecules has an unusual cyclopentane structure similar to those recently reported in plant extracts and lakebed sediments.

The determination of the absolute configurations of chiral molecules using vibrational circular dichroism (VCD) spectroscopy.

The vibrational circular dichroism (VCD) spectra of the two enantiomers of a chiral molecule are of equal magnitude and opposite sign: i.e. mirror-image enantiomers give mirror-image VCD spectra. In principle, the absolute configuration (AC) of a chiral molecule can therefore be determined from its VCD spectrum. In practice, the determination of the AC of a chiral molecule from its experimental VCD spectrum requires a methodology which reliably predicts the VCD spectra of its enantiomers. The only reliable methodology developed to date uses the Stephens quantum-mechanical theory of the rotational strengths of fundamental vibrational transitions, developed in the early 1980s, implemented using ab initio density functional theory in the GAUSSIAN program in the mid 1990s. This methodology has by now been widely used in determining ACs from experimental VCD spectra. In this article we discuss the protocol for determining the ACs of chiral molecules with optimum reliability and its implementation for a variety of molecules, including the D3 symmetry perhydrotriphenylene, a thiazino-oxadiazolone recently shown to be a highly active calcium entry channel blocker, the alkaloid natural products schizozygine, iso-schizogaline, and iso-schizogamine, and the iridoid natural products plumericin, iso-plumericin, and prismatomerin. The power of VCD spectroscopy in determining ACs, even for large organic molecules and for substantially conformationally-flexible organic molecules is clearly documented.

Determination of the absolute configurations of natural products via density functional theory calculations of vibrational circular dichroism, electronic circular dichroism, and optical rotation: the iso-schizozygane alkaloids isoschizogaline and isoschizogamine.

The development of density functional theory (DFT) methods for the calculation of vibrational circular dichroism (VCD), electronic circular dichroism (ECD), and transparent spectral region optical rotation (OR) has revolutionized the determination of the absolute configurations (ACs) of chiral molecules using these chiroptical properties. We report the concerted application of DFT calculations of VCD, ECD, and OR to the determination of the ACs of the isoschizozygane alkaloid natural products, isoschizogaline, and isochizogamine, whose ACs have not previously been determined. The ACs of naturally occurring (-)-isoschizogaline and (-)-isoschizogamine, are both determined definitively to be 2R, 7R, 20S, 21S.

Determination of the absolute configurations of pharmacological natural products via density functional theory calculations of vibrational circular dichroism: the new cytotoxic iridoid prismatomerin.

A new highly cytotoxic iridoid has very recently been isolated from Prismatomeris tetrandra and shown to have the structure 3, similar to that of the iridoid oruwacin, 2. We report the determination of the absolute configuration (AC) of the new iridoid, prismatomerin, using vibrational circular dichroism (VCD) spectroscopy. The VCD spectrum of the acetate derivative of 3, 4, is analyzed using the Stephens theory of VCD and density functional theory (DFT). The AC of the naturally occurring 3 is shown to be 1R,5S,8S,9S,10S, identical to that of the naturally occurring iridoid plumericin, 1, also determined using VCD spectroscopy. The [alpha]D values of the natural products 3 and 1 are negative and positive, respectively. Since the ACs of 3 and 1 are identical, it follows that the AC of 3 cannot be correctly determined by empirical comparison of the signs of the [alpha]D values of 3 and 1.

Prismatomerin, a new iridoid from Prismatomeris tetrandra. Structure elucidation, determination of absolute configuration, and cytotoxicity.

A new complex iridoid, prismatomerin (1), has been isolated from the leaves of Prismatomeris tetrandra, together with the known glucoside gaertneroside (4). The structures of 1 and 4 were determined by spectroscopic analysis, notably 2D NMR techniques. The (1R,5S,8S,9S,10S)-(-) absolute configuration of prismatomerin (1) was determined by comparison of the vibrational circular dichroism (VCD) spectrum calculated using density functional theory and the experimental VCD spectrum of the O-acetyl derivative 3. Prismatomerin (1) showed remarkable antitumor activity and also interfered with mitotic spindle formation.

Determination of the absolute configurations of natural products via density functional theory calculations of vibrational circular dichroism, electronic circular dichroism and optical rotation: the schizozygane alkaloid schizozygine.

The development of density functional theory (DFT) methods for the calculation of vibrational circular dichroism (VCD), electronic circular dichroism (ECD), and transparent spectral region optical rotation (OR) has revolutionized the determination of the absolute configurations (ACs) of chiral molecules using these chiroptical properties. We report the first concerted application of DFT calculations of VCD, ECD, and OR to the determination of the AC of a natural product whose AC was previously undetermined. The natural product is the alkaloid schizozygine, isolated from Schizozygia caffaeoides. Comparison of DFT calculations of the VCD, ECD, and OR of schizozygine to experimental data leads, for each chiroptical technique, to the AC 2R,7S,20S,21S for the naturally occurring (+)-schizozygine. Three other alkaloids, schizogaline, schizogamine, and 6,7-dehydro-19beta-hydroxyschizozygine, have also been isolated from S. caffaeoides and shown to have structures closely related to schizozygine. Assuming a common biosynthetic pathway, their ACs are defined by that of schizozygine.

Insight into the activation mechanism of Escherichia coli octaprenyl pyrophosphate synthase derived from pre-steady-state kinetic analysis.

Octaprenyl pyrophosphate synthase (OPPs) catalyzes the sequential condensation of five molecules of isopentenyl pyrophosphate with farnesyl pyrophosphate to generate all-trans C40-octaprenyl pyrophosphate, which constitutes the side chain of ubiquinone. Due to the slow product release, a long-chain polyprenyl pyrophosphate synthase often requires detergent or another factor for optimal activity. Our previous studies in examining the activity enhancement of Escherichia coli undecaprenyl pyrophosphate synthase have demonstrated a switch of the rate-determining step from product release to isopentenyl pyrophosphate (IPP) condensation reaction in the presence of Triton [12]. In order to understand the mechanism of enzyme activation for E. coli OPPs, a single-turnover reaction was performed and the measured IPP condensation rate (2 s(-1)) was 100 times larger than the steady-state rate (0.02 s(-1)). The high molecular weight fractions and Triton could accelerate the steady-state rate by 3-fold (0.06 s(-1)) but insufficient to cause full activation (100-fold). A burst product formation was observed in enzyme multiple turnovers indicating a slow product release.