Computational studies on the sterol-like cyclization of a monodomain class II terpene cyclase.

Recently the first example of a class II terpene cyclase comprised of only a single domain was reported. Class II synthases are a diverse group of enzymes that catalyze exceptionally complex reactions, including the remarkable cyclization of steroids. This discovery of a single-domain enzyme being able to catalyze a steroid-like product contradicted the long-held tenet that complex class II cyclizations required double-domain enzymes. The proposed mechanism for the sterol-like cyclization of a monodomain class II terpene cyclase was studied computationally by using density functional theory (DFT). The complete pathway for the conversion of 5-geranyl-3,4-dihydroxybenzoate to the steroid-like pentacyclic product merosterolic acid A was elucidated. The formation of a tricyclic carbocation intermediate with three cyclohexane rings was found to be a concerted, but asynchronous, cyclization. The formation of the fourth ring proceeds with a low energy activation Friedel-Crafts reaction. Subsequent deprotonation of this pentacyclic system gave as the final product merosterolic acid. The overall conversion was found to be highly exothermic due to the conversion of three C-C double bonds to C-C single bonds.

Biosynthetic Mechanism of Lanosterol: Cyclization.

The remarkable cyclization mechanism of the formation of the 6-6-6-5 tetracyclic lanosterol (a key triterpenoid intermediate in the biosynthesis of cholesterol) from the acyclic 2,3-oxidosqualene catalyzed by oxidosqualene cyclase (OSC) has stimulated the interest of chemists and biologists for over a half century. Herein, the elaborate, state-of-the-art two-dimensional (2D) QM/MM MD simulations have clearly shown that the cyclization of the A-C rings involves a nearly concerted, but highly asynchronous cyclization, to yield a stable intermediate with "6-6-5" rings followed by the ring expansion of the C-ring concomitant with the formation of the D-ring to yield the "6-6-6-5" protosterol cation. The calculated reaction barrier of the rate-limiting step (≈22 kcal mol(-1)) is comparable to the experimental kinetic results. Furthermore all previous experimental mutagenic evidence is highly consistent with the identified reaction mechanism.

The concerted nature of the cyclization of squalene oxide to the protosterol cation.

Concerted A-C ring formation: A concerted, but highly asynchronous, pathway was identified for the formation of rings A-C in the biosynthetic conversion of squalene oxide to the prosterol cation, with ring B being formed in the required boat conformation.

Concerted, highly asynchronous, enzyme-catalyzed [4 + 2] cycloaddition in the biosynthesis of spinosyn A; computational evidence.

A theoretical study has been carried out on model systems to study a recently reported, (Nature, 2011, 473, 109) biosynthetic, [4 + 2] cycloaddition catalyzed by a stand-alone enzyme (the cyclase SpnF). It was suggested in this paper that SpnF is the first known example of a Diels-Alderase (DA). In the present study, for a model system of the substrate a transition structure was found with density functional calculations (DFT). In addition, the intrinsic reaction coordinate calculations indicated that the transition structure is that of a concerted, but highly asynchronous, DA reaction. Based on the DFT and Møller-Plesset second order calculations the activation energy was estimated to be about 15 kcal mol(-1). The results of a natural population analysis indicated that there is significant charge transfer in the transition state, and it is proposed that possibly the enzyme plays a dual role of not only folding the substrate into the proper conformation for the DA reaction to occur, but also lowering its activation energy by stabilization of the highly polarized transition structure.

Physical constraints on sesquiterpene diversity arising from cyclization of the eudesm-5-yl carbocation.

The biogenic origins of complex cyclic terpenes derive from the interplay of enzymes and the intrinsic reactivity of carbocation species at major branch-points along intramolecular cyclization pathways to ultimately determine the distribution of terpene skeletal types in nature. Solanaceous plants biosynthesize chemical defense compounds, largely derived from the eremophilane and spirovetivane-type sesquiterpenes. These hydrocarbon skeletons share a common biogenic origin, stemming from alternative Wagner-Meerwein rearrangements of the eudesm-5-yl carbocation during the cyclization of farnesyl pyrophosphate (FPP) catalyzed by sesquiterpene synthases. While the spirojatamane skeleton shares the same carbocation intermediate, this class of sesquiterpenes has not been reported in the Solanaceae and is exceedingly rare in nature. To investigate the physical basis for alternative rearrangements of the eudesm-5-yl carbocation, we carried out quantum mechanics (QM) analyses to calculate the allowable conformations, energies, and transition states linking conformers of the eudesm-5-yl carbocation to the eremophilene, spirovetivane, and spirojatamane skeletons. Additionally, we conducted parallel investigations on simplified decalin carbocation systems to examine the contribution of ring substituents to allowable conformations and rearrangement pathways. Our study reveals that ring substituents expand the conformational space accessible to the eudesm-5-yl carbocation while sterically blocking rearrangements in certain contexts. From our analysis, we define a conformational threshold for each possible rearrangement based on dihedral angles describing transition state geometry. Further, our calculations indicate that methylene migration rearrangements leading to spiro compounds are thermodynamically dominant in the eudesm-5-yl and simpler decalin carabocation systems. Interestingly, the theoretical abundance of sesquiterpene skeletal types arising from the intrinsic reactivity of the eudesm-5-yl carbocation stands in sharp contrast to their currently known natural abundance. The implications of these results for the catalytic tragectories catalyzed by sesquiterpene synthases are discussed.

Lanthanides caged by the organic chelates; structural properties.

The structure, in particular symmetry, geometry and morphology of organic chelates coordinated with the lanthanide ions are analyzed in the present review. This is the first part of a complete presentation of a theoretical description of the properties of systems, which are widely used in technology, but most of all, in molecular biology and medicine. The discussion is focused on the symmetry and geometry of the cages, since these features play a dominant role in the spectroscopic activity of the lanthanides caged by organic chelates. At the same time, the spectroscopic properties require more formal presentation in the language of Racah algebra, and deserve a separate analysis. In addition to the parent systems of DOTA, DOTP, EDTMP and CDTMP presented here, their modifications by various antennas are analyzed. The conclusions that have a strong impact upon the theory of the energy transfer and the sensitized luminescence of these systems are based on the results of numerical density functional theory calculations.

Compelling computational evidence for the concerted cyclization of the ABC rings of hopene from protonated squalene.

The long-standing question of what is the nature of the cyclization of squalene to form tetracyclic and pentacyclic triterpenes has been addressed computationally. Using the DFT method with an intrinsic reaction coordinate calculation, we find that the first three rings of protonated squalene were formed without the intermediacy of mono- or bicyclic carbocations. The cyclization, calculated in the gas phase, proceeds in a highly asynchronous, concerted reaction to yield the tricyclic, tertiary carbocation with a 5-membered C ring. The fourth double bond of squalene is not properly oriented for the ring expansion of the C ring in concert with the formation of the 5-membered ring.

Structural elucidation of cisoid and transoid cyclization pathways of a sesquiterpene synthase using 2-fluorofarnesyl diphosphates.

Sesquiterpene skeletal complexity in nature originates from the enzyme-catalyzed ionization of (trans,trans)-farnesyl diphosphate (FPP) (1a) and subsequent cyclization along either 2,3-transoid or 2,3-cisoid farnesyl cation pathways. Tobacco 5-epi-aristolochene synthase (TEAS), a transoid synthase, produces cisoid products as a component of its minor product spectrum. To investigate the cryptic cisoid cyclization pathway in TEAS, we employed (cis,trans)-FPP (1b) as an alternative substrate. Strikingly, TEAS was catalytically robust in the enzymatic conversion of (cis,trans)-FPP (1b) to exclusively (>/=99.5%) cisoid products. Further, crystallographic characterization of wild-type TEAS and a catalytically promiscuous mutant (M4 TEAS) with 2-fluoro analogues of both all-trans FPP (1a) and (cis,trans)-FPP (1b) revealed binding modes consistent with preorganization of the farnesyl chain. These results provide a structural glimpse into both cisoid and transoid cyclization pathways efficiently templated by a single enzyme active site, consistent with the recently elucidated stereochemistry of the cisoid products. Further, computational studies using density functional theory calculations reveal concerted, highly asynchronous cyclization pathways leading to the major cisoid cyclization products. The implications of these discoveries for expanded sesquiterpene diversity in nature are discussed.

Quantitative exploration of the catalytic landscape separating divergent plant sesquiterpene synthases.

Throughout molecular evolution, organisms create assorted chemicals in response to varying ecological niches. Catalytic landscapes underlie metabolic evolution, wherein mutational steps alter the biosynthetic properties of enzymes. Here we report the first systematic quantitative characterization of the catalytic landscape underlying the evolution of sesquiterpene chemical diversity. On the basis of our previous discovery of a set of nine naturally occurring amino acid substitutions that functionally interconverted orthologous sesquiterpene synthases from Nicotiana tabacum and Hyoscyamus muticus, we created a library of all possible residue combinations (2(9) = 512) in the N. tabacum enzyme. The product spectra of 418 active enzymes revealed a rugged landscape where several minimal combinations of the nine mutations encode convergent solutions to the interconversions of parental activities. Quantitative comparisons indicated context dependence for mutational effects--epistasis--in product specificity and promiscuity. These results provide a measure of the mutational accessibility of phenotypic variability in a diverging lineage of terpene synthases.

Lewis acid-catalyzed one-pot, three-component route to chiral 3,3'-bipyrroles.

3,3'-Bipyrroles 3 could be synthesized using a double Michael addition reaction involving diaroyl acetylene 1 and the appropriate 1,3-dicarbonyls 2 using ammonium acetate as a nitrogen source. The axial chirality of bipyrrole was anticipated from the X-ray crystal structure and DFT calculations and confirmed by separating the racemates on a chiral column and subsequent CD spectra of the enantiomers. The absolute configuration of the enantiomers was achieved by theoretical CD spectra calculation using the ZINDO method.

Thermal [1,5] hydrogen sigmatropic shifts in cis,cis-1,3-cyclononadienes probed by gas-phase kinetic studies and density functional theory calculations.

The kinetics of gas-phase thermal [1,5] hydrogen shifts interconverting the five isomeric mono-deuterium-labeled cis,cis-1,3-cyclononadienes have been followed at four temperatures from 240 to 287 degrees C. The activation parameters found were Ea = 37.1 +/- 0.8 kcal/mol, log A = 11.6 +/- 0.3, DeltaH++ = 36.0 +/- 0.8 kcal/mol, and DeltaS++ = -9.0 +/- 0.3 eu. Density functional theory based calculations have provided geometries and energies for the ground-state cyclononadiene conformational isomers, for the transition states linking one to another, and for the transition states for [1,5] hydrogen shifts responsible for isomerizations among the five labeled dienes. A generalized formulation of the Winstein-Holness equation is presented and applied to the complex system, one that involves 11 ground-state conformers, 10 transition states separating them, and five transition states for [1,5] hydrogen shifts. The value for the empirical Ea derived from calculated mole fractions of ground-state conformers and calculated energies for specific ground-state conformers and [1,5] hydrogen shift transition structures was 37.5 kcal/mol, in excellent agreement with the experimentally obtained activation energy. The significance of conformational options in various ground states and transition structures for the [1,5] hydrogen shifts is considerable, an inference that may well have general applicability.

Density-functional theory structures of 1,4,7,10-tetraazacyclododecane-1,4,7,10-tetraacetic acid complexes for ions across the lanthanide series.

The use of organically chelated lanthanides in diagnosis and treatment is a rapidly growing field in medicine. In order to gain a deeper understanding into the properties of these chelates, particularly spectroscopic, density-functional calculations have been performed on a series of lanthanide ions chelated with 1,4,7,10-tetraazacyclododecane-1,4,7,10-tetraacetic acid. Based on the results of these calculations, it has been concluded that the local symmetry experienced by the chelated lanthanide ion may be treated as being axial, which will make the interpretation of their spectroscopic properties greatly simplified. It has also been suggested that the so-called "capping" water molecule at the ninth coordination position of the lanthanide is hydrogen bonded to the acetate oxygens of the sidearms, rather than coordinated as the ninth ligand of the lanthanide.

Concerted nature of AB ring formation in the enzymatic cyclization of squalene to hopenes.

The cyclization of the A-B rings of squalene to hopene is studied computationally (DFT). A transition structure is found for a concerted, asynchronous pathway for the formation of chair-chair decalin carbocation. The computationally derived conformer leading to this asynchronous transition structure is remarkably similar to the analogous region of 2-azasqualene encapsulated by squalene-hopene cyclase recently reported by Schulz. A concerted A-B ring closure is likely to occur in the cyclization of squalene to hopene. [structure--see text]