A conserved lysine in beta-lactam synthetase assists ring cyclization: Implications for clavam and carbapenem biosynthesis

beta-Lactam synthetase (beta-LS) is the paradigm of a growing class of enzymes that form the critical beta-lactam ring in the clavam and carbapenem antibiotics. beta-LS catalyzes a two-stage reaction in which N(2)-(2-carboxyethyl)-L-arginine is first adenylated, and then undergoes intramolecular ring closure. It was previously shown that the forward kinetic commitment to beta-lactam formation is high, and that the overall rate of reaction is partially limited to a protein conformational change rather than to the chemical step alone of closing the strained ring. beta-Lactam formation was evaluated on the basis of X-ray crystal structures, site-specific mutation, and kinetic and computational studies. The combined evidence clearly points to a reaction coordinate involving the formation of a tetrahedral transition state/intermediate stabilized by a conserved Lys. The combination of substrate preorganization, a well-stabilized transition state and an excellent leaving group facilitates this acyl substitution to account for the strong forward commitment to catalysis and to lower the barrier of four-membered ring formation to the magnitude of a protein conformational change.

Absence of the aflatoxin biosynthesis gene, norA, allows accumulation of deoxyaflatoxin B1 in Aspergillus flavus cultures

Biosynthesis of the highly toxic and carcinogenic aflatoxins in select Aspergillus species from the common intermediate O-methylsterigmatocystin has been postulated to require only the cytochrome P450 monooxygenase, OrdA (AflQ). We now provide evidence that the aryl alcohol dehydrogenase NorA (AflE) encoded by the aflatoxin biosynthetic gene cluster in Aspergillus flavus affects the accumulation of aflatoxins in the final steps of aflatoxin biosynthesis. Mutants with inactive norA produced reduced quantities of aflatoxin B(1) (AFB(1)), but elevated quantities of a new metabolite, deoxyAFB(1). To explain this result, we suggest that, in the absence of NorA, the AFB(1) reduction product, aflatoxicol, is produced and is readily dehydrated to deoxyAFB(1) in the acidic medium, enabling us to observe this otherwise minor toxin produced in wild-type A. flavus.

Non-heme iron oxygenases generate natural structural diversity in carbapenem antibiotics

Carbapenems are a clinically important antibiotic family. More than 50 naturally occurring carbapenam/ems are known and are distinguished primarily by their C-2/C-6 side chains where many are only differentiated by the oxidation states of these substituents. With a limited palette of variations the carbapenem family comprises a natural combinatorial library, and C-2/C-6 oxidation is associated with increased efficacy. We demonstrate that ThnG and ThnQ encoded by the thienamycin gene cluster in Streptomyces cattleya oxidize the C-2 and C-6 moieties of carbapenems, respectively. ThnQ stereospecifically hydroxylates PS-5 (5) giving N-acetyl thienamycin (2). ThnG catalyzes sequential desaturation and sulfoxidation of PS-5 (5), giving PS-7 (7) and its sulfoxide (9). The enzymes are relatively substrate selective but are proposed to give rise to the oxidative diversity of carbapenems produced by S. cattleya, and orthologues likely function similarly in allied streptomyces. Elucidating the roles of ThnG and ThnQ will focus further investigations of carbapenem antibiotic biosynthesis.

Production of octaketide polyenes by the calicheamicin polyketide synthase CalE8: implications for the biosynthesis of enediyne core structures

Enediyne antibiotics are categorized according to the presence of either a 9- or 10-membered ring within their polyketide-derived core structures. Recent literature reports have favored the notion that biosynthetic divergence of the two structural families is determined by the enediyne polyketide synthases (PKSs) alone. We now disclose the simultaneous in vitro production of three octaketide polyenes by biosynthetic enzymes for the 10-membered enediyne calicheamicin gamma(1)(I), including the elusive beta-keto acid precursor to a previously described C15 methyl hexaenone. Alongside these two polyene products, we have additionally detected a hydrocarbon heptaene previously isolated only from 9-membered enediyne systems. The discovery of the heptaene in the calicheamicin system promotes a more convergent model for the early steps of enediyne biosynthesis. Furthermore, the synthesis of this set of octaketides by the enediyne PKS CalE8 and thioesterase CalE7 suggests, in contrast to recent biosynthetic proposals, that accessory enzymes may be necessary to initiate differentiation to 9- or 10-membered enediyne precursors, either by modulation of enediyne PKS activity or by interception and modification of polyketide chain-extension intermediates.

Structural basis for biosynthetic programming of fungal aromatic polyketide cyclization

Polyketides are a class of natural products with diverse structures and biological activities. The structural variability of aromatic products of fungal nonreducing, multidomain iterative polyketide synthases (NR-PKS group of IPKSs) results from regiospecific cyclizations of reactive poly-beta-keto intermediates. How poly-beta-keto species are synthesized and stabilized, how their chain lengths are determined, and, in particular, how specific cyclization patterns are controlled have been largely inaccessible and functionally unknown until recently. A product template (PT) domain is responsible for controlling specific aldol cyclization and aromatization of these mature polyketide precursors, but the mechanistic basis is unknown. Here we present the 1.8 A crystal structure and mutational studies of a dissected PT monodomain from PksA, the NR-PKS that initiates the biosynthesis of the potent hepatocarcinogen aflatoxin B(1) in Aspergillus parasiticus. Despite having minimal sequence similarity to known enzymes, the structure displays a distinct 'double hot dog' (DHD) fold. Co-crystal structures with palmitate or a bicyclic substrate mimic illustrate that PT can bind both linear and bicyclic polyketides. Docking and mutagenesis studies reveal residues important for substrate binding and catalysis, and identify a phosphopantetheine localization channel and a deep two-part interior binding pocket and reaction chamber. Sequence similarity and extensive conservation of active site residues in PT domains suggest that the mechanistic insights gleaned from these studies will prove general for this class of IPKSs, and lay a foundation for defining the molecular rules controlling NR-PKS cyclization specificity.

A catalytic asymmetric route to carbapenems

Efficient syntheses of N-acetyl thienamycin and epithienamycin A in their readily deprotected form are reported where three contiguous stereocenters are established in a single catalytic asymmetric azetidinone-forming reaction. These examples are a template for synthesizing C-5/C-6 cis or trans carbapenems with independent control of the C-8 stereocenter. A library of oxidatively and sterochemically defined azetidinone precursors to a variety of naturally occurring carbapenems and potential biosynthetic intermediates has been prepared to facilitate studies of carbapenem antibiotic biosynthesis.

A conserved tyrosyl-glutamyl catalytic dyad in evolutionarily linked enzymes: carbapenam synthetase and beta-lactam synthetase

Beta-lactam-synthesizing enzymes carbapenam synthetase (CPS) and beta-lactam synthetase (beta-LS) are evolutionarily linked to a common ancestor, asparagine synthetase B (AS-B). These three relatives catalyze substrate acyl-adenylation and nucleophilic acyl substitution by either an external (AS-B) or internal (CPS, beta-LS) nitrogen source. Unlike AS-B, crystal structures of CPS and beta-LS revealed a putative Tyr-Glu dyad (CPS, Y345/E380; beta-LS, Y348/E382) proposed to deprotonate the respective internal nucleophile. CPS and beta-LS site-directed mutagenesis (Y345/8A, Y345/8F, E380/2D, E380/2Q, E380A) resulted in the reduction of their catalytic efficiency, with Y345A, E380A, and E382Q producing undetectable amounts of beta-lactam product. However, [(32)P]PP(i)-ATP exchange assays demonstrated Y345A and E380A undergo the first half-reaction, with the remaining active mutants showing decreased forward commitment to beta-lactam cyclization. pH-rate profiles of CPS and beta-LS supported the importance of a Tyr-Glu dyad in beta-lactam formation and suggested its reverse protonation in beta-LS. The kinetics of CPS double-site mutants reinforced the synergism of Tyr-Glu in catalysis. Furthermore, significant solvent isotope effects on k(cat) ((D)k(cat)) for Y345F (1.9) and Y348F (1.7) maintained the assignment of Y345/8 in proton transfer. A proton inventory on Y348F determined its (D)(k(cat)/K(m)) = 0.2 to arise from multiple reactant-state fractionation factors, presumably from water molecule(s) replacing the missing Tyr hydroxyl. The role of a CPS and beta-LS Tyr-Glu catalytic dyad was solidified by a significant decrease in mutant k(cat) viscosity dependence with respect to the wild-type enzymes. The evolutionary relation and potential for engineered biosynthesis were demonstrated by beta-LS acting as a carbapenam synthetase.

Design and synthesis of small molecule glycerol 3-phosphate acyltransferase inhibitors

The incidence of obesity and other diseases associated with an increased triacylglycerol mass is growing rapidly, particularly in the United States. Glycerol 3-phosphate acyltransferase (GPAT) catalyzes the rate-limiting step of glycerolipid biosynthesis, the acylation of glycerol 3-phosphate with saturated long-chain acyl-CoAs. In an effort to produce small molecule inhibitors of this enzyme, a series of benzoic and phosphonic acids was designed and synthesized. In vitro testing of this series has led to the identification of several compounds, in particular 2-(nonylsulfonamido)benzoic acid (15g), possessing moderate GPAT inhibitory activity in an intact mitochondrial assay.

Design and synthesis of a beta-lactamase activated 5-fluorouracil prodrug

An efficient synthesis of a 5-fluorouracil-cephalosporin prodrug is described for use against colorectal and other cancers in antibody and gene-directed therapies. The compound shows stability in aqueous media until specifically activated by beta-lactamase (betaL). The kinetic parameters of the 5-fluorouracil-cephalosporin conjugate were determined in the presence of Enterobacter cloacae P99 betaL (ECl betaL) revealing a K(m)=95.4 microM and V(max)=3.21 microMol min(-1) mg(-1). The data compare favorably to related systems that have been reported and enable testing of this prodrug against cancer cell lines in vitro and in vivo.

Dissection of the stepwise mechanism to beta-lactam formation and elucidation of a rate-determining conformational change in beta-lactam synthetase

Clavulanic acid is a widely used beta-lactamase inhibitor whose key beta-lactam core is formed by beta-lactam synthetase. beta-Lactam synthetase exhibits a Bi-Ter mechanism consisting of two chemical steps, acyl-adenylation followed by beta-lactam formation. 32PPi-ATP exchange assays showed the first irreversible step of catalysis is acyl-adenylation. From a small, normal solvent isotope effect (1.38 +/- 0.04), it was concluded that beta-lactam synthesis contributes at least partially to kcat. Site-specific mutation of Lys-443 identified this residue as the ionizable group at pKa approximately 8.1 apparent in the pH-kcat profile that stabilizes the beta-lactam-forming step. Viscosity studies demonstrated that a protein conformational change was also partially rate-limiting on kcat attenuating the observed solvent isotope effect on beta-lactam formation. Adherence to Kramers' theory gave a slope of 1.66 +/- 0.08 from a plot of log(o kcat/kcat) versus log(eta/eta(o)) consistent with opening of a structured loop visible in x-ray data preceding product release. Internal "friction" within the enzyme contributes to a slope of > 1 in this analysis. Correspondingly, earlier in the catalytic cycle ordering of a mobile active site loop upon substrate binding was manifested by an inverse solvent isotope effect (0.67 +/- 0.15) on kcat/Km. The increased second-order rate constant in heavy water was expected from ordering of this loop over the active site imposing torsional strain. Finally, an Eyring plot displayed a large enthalpic change accompanying loop movement (DeltaH approximately 20 kcal/mol) comparable to the chemical barrier of beta-lactam formation.

Deconstruction of iterative multidomain polyketide synthase function

PksA, which initiates biosynthesis of the environmental carcinogen aflatoxin B1, is one of the multidomain iterative polyketide synthases (IPKSs), a large, poorly understood family of biosynthetic enzymes. We found that dissection of PksA and its reconstitution from selected sets of domains allows the accumulation and characterization of advanced octaketide intermediates bound to the enzyme, permitting the reactions controlled by individual catalytic domains to be identified. A product template (PT) domain unites with the ketosynthase and thioesterase in this IPKS system to assemble precisely seven malonyl-derived building blocks to a hexanoyl starter unit and mediate a specific cyclization cascade. Because the PT domain is common among nonreducing IPKSs, these mechanistic features should prove to be general for IPKS-catalyzed production of aromatic polyketides.

Starter unit specificity directs genome mining of polyketide synthase pathways in fungi

Search of the protein database with the aflatoxin pathway polyketide synthase (PKS) revealed putative PKSs in the pathogenic fungi Coccidioides immitis and Coccidioides posadasii that could require partnerships with a pair of fatty acid synthase (FAS) subunits for the biosynthesis of fatty acid-polyketide hybrid metabolites. A starter unit:acyl-carrier protein transacylase (SAT) domain was discovered in the nonreducing PKS. This domain is thought to accept the fatty acid product from the FAS to initiate polyketide synthesis. We expressed the C. immitis SAT domain in Escherichia coli and showed that this domain, unlike that from the aflatoxin pathway PKS, transferred octanoyl-CoA four times faster than hexanoyl-CoA. The SAT domain also formed a covalent octanoyl intermediate and transferred this group to a free-standing ACP domain. Our results suggest that C. immitis/posadasii, both human fungal pathogens, contain a FAS/PKS cluster with functional similarity to the aflatoxin cluster found in Aspergillus species. Dissection of the PKS and determination of in vitro SAT domain specificity provides a tool to uncover the growing number of similar sequenced pathways in fungi, and to guide elucidation of the fatty acid-polyketide hybrid metabolites that they produce.

Observation of an acryloyl-thiamin diphosphate adduct in the first step of clavulanic acid biosynthesis

The first committed biosynthetic step toward clavulanic acid, the clinically important beta-lactamase inhibitor, is catalyzed by the thiamin diphosphate (ThDP)-dependent enzyme N2-(2-carboxyethyl)arginine synthase (CEAS). This protein carries out a unique reaction among ThDP-dependent processes in which a C-N bond is formed, and an electrophilic acryloyl-thiazolium intermediate of ThDP is proposed to be involved, unlike the nucleophilic enamine species typically generated by this class of enzymes. Here we present evidence for the existence of the putative acryloyl adduct and report the unexpected observation of a long-wavelength chromophore (lambda = 433 nm), which we attribute to this enzyme-bound species. Chemical models were synthesized that both confirm its expected absorption (lambda = 310-320 nm) and exclude self-condensation and intramolecular imine formation with the cofactor as its cause. Circular dichroism experiments and others discount charge transfer as a likely explanation for the approximately 120 nm red shift of the chromophore ( approximately 25 kcal). Examples are well-known of charged molecules that exhibit significantly red-shifted UV-visible spectra compared to their neutral forms as, for example, polyene cations and dyes such as indigo and the cyanines. Rhodopsin is the classic biochemical example where the protein (opsin)-bound protonated Schiff base of retinal displays a remarkable range of red-shifted absorptions modulated by the protein environment. Similar tuning of the chromophoric behavior of the enzyme-bound CEAS acryloyl.ThDP species may be occurring.

Pharmacological stimulation of brain carnitine palmitoyl-transferase-1 decreases food intake and body weight

Inhibition of brain carnitine palmitoyl-transferase-1 (CPT-1) is reported to decrease food intake and body weight in rats. Yet, the fatty acid synthase (FAS) inhibitor and CPT-1 stimulator C75 produces hypophagia and weight loss when given to rodents intracerebroventricularly (icv). Thus roles and relative contributions of altered brain CPT-1 activity and fatty acid oxidation in these phenomena remain unclarified. We administered compounds that target FAS or CPT-1 to mice by single icv bolus and examined acute and prolonged effects on feeding and body weight. C75 decreased food intake rapidly and potently at all doses (1-56 nmol) and dose dependently inhibited intake on day 1. Dose-dependent weight loss on day 1 persisted through 4 days of postinjection monitoring. The FAS inhibitor cerulenin produced dose-dependent (560 nmol) hypophagia for 1 day, weight loss for 2 days, and weight regain to vehicle control by day 3. The CPT-1 inhibitor etomoxir (32, 320 nmol) did not alter overall day 1 feeding. However, etomoxir attenuated the hypophagia produced by C75, indicating that CPT-1 stimulation is important for C75's effect. A novel compound, C89b, was characterized in vitro as a selective stimulator of CPT-1 that does not affect fatty acid synthesis. C89b (100, 320 nmol) decreased feeding in mice for 3 days and produced persistent weight loss for 6 days without producing conditioned taste aversion. Similarly, intraperitoneal administration decreased feeding and body weight without producing conditioned taste aversion. These results suggest a role for brain CPT-1 in the regulation of energy balance and implicate CPT-1 stimulation as a pharmacological approach to weight loss.

Rate-limiting steps and role of active site Lys443 in the mechanism of carbapenam synthetase

Carbapenam synthetase (hereafter named CPS) catalyzes the formation of the beta-lactam ring in the biosynthetic pathway to (5R)-carbapen-2-em-3-carboxylate, the simplest of the carbapenem antibiotics. Kinetic studies showed remarkable tolerance to substrate stereochemistry in the turnover rate but did not distinguish between chemistry and a nonchemical step such as product release or conformational change as being rate-determining. Also, X-ray structural studies and modest sequence homology to beta-lactam synthetase, an enzyme that catalyzes the formation of a monocyclic beta-lactam ring in a similar ATP/Mg2+-dependent reaction, implicate K443 as an essential residue for substrate binding and intermediate stabilization. In these experiments, we use pH-rate profiles, deuterium solvent isotope effects, and solvent viscosity measurements to examine the rate-limiting step in this complex overall process of substrate adenylation and intramolecular ring formation. Mutagenesis and chemical rescue demonstrate that K443 is the general acid visible in the pH-rate profile of the wild-type CPS-catalyzed reaction. On the basis of these results, we propose a mechanism in which the rate-limiting step is beta-lactam ring formation coupled to a protein conformational change and underscore the role of K443 throughout the reaction.

Fatty acid synthase inhibition activates AMP-activated protein kinase in SKOV3 human ovarian cancer cells

Fatty acid synthase (FAS), the enzyme responsible for the de novo synthesis of fatty acids, is highly expressed in ovarian cancers and most common human carcinomas. Inhibition of FAS and activation of AMP-activated protein kinase (AMPK) have been shown to be cytotoxic to human cancer cells in vitro and in vivo. In this report, we explore the cytotoxic mechanism of action of FAS inhibition and show that C93, a synthetic FAS inhibitor, increases the AMP/ATP ratio, activating AMPK in SKOV3 human ovarian cancer cells, which leads to cytotoxicity. As a physiologic consequence of AMPK activation, acetyl-CoA carboxylase (ACC), the rate-limiting enzyme of fatty acid synthesis, was phosphorylated and inhibited whereas glucose oxidation was increased. Despite these attempts to conserve energy, the AMP/ATP ratio increased with worsening cellular redox status. Pretreatment of SKOV3 cells with compound C, an AMPK inhibitor, substantially rescued the cells from C93 cytotoxicity, indicating its dependence on AMPK activation. 5-(Tetradecyloxy)-2-furoic acid, an ACC inhibitor, did not activate AMPK despite inhibiting fatty acid synthesis pathway activity and was not significantly cytotoxic to SKOV3 cells. This indicates that substrate accumulation from FAS inhibition triggering AMPK activation, not end-product depletion of fatty acids, is likely responsible for AMPK activation. C93 also exhibited significant antitumor activity and apoptosis against SKOV3 xenografts in athymic mice without significant weight loss or cytotoxicity to proliferating cellular compartments such as bone marrow, gastrointestinal tract, or skin. Thus, pharmacologic FAS inhibition selectively activates AMPK in ovarian cancer cells, inducing cytotoxicity while sparing most normal human tissues from the pleiotropic effects of AMPK activation.

Identification of a starter unit acyl-carrier protein transacylase domain in an iterative type I polyketide synthase

Polyketides are a class of natural products that exhibit a wide range of functional and structural diversity. They include antibiotics, immunosuppressants, antifungals, antihypercholesterolemics, and cytotoxins. Polyketide synthases (PKSs) use chemistry similar to fatty acid synthases (FASs), although building block variation and differing extents of reduction of the growing polyketide chain underlie their biosynthetic versatility. In contrast to the well studied sequential modular type I PKSs, less is known about how the iterative type I PKSs carry out and control chain initiation, elongation, folding, and cyclization during polyketide processing. Domain structure analysis of a group of related fungal, nonreducing PKSs has revealed well defined N-terminal domains longer than commonly seen for FASs and modular PKSs. Predicted structure of this domain disclosed a region similar to malonyl-CoA:acyl-carrier protein (ACP) transacylases (MATs). MATs play a key role transferring precursor CoA thioesters from solution onto FASs and PKSs for chain elongation. On the basis of site-directed mutagenesis, radiolabeling, and kinetics experiments carried out with individual domains of the norsolorinic acid PKS, we propose that the N-terminal domain is a starter unit:ACP transacylase (SAT domain) that selects a C(6) fatty acid from a dedicated yeast-like FAS and transfers this unit onto the PKS ACP, leading to the production of the aflatoxin precursor, norsolorinic acid. These findings could indicate a much broader role for SAT domains in starter unit selection among nonreducing iterative, fungal PKSs, and they provide a biochemical rationale for the classical acetyl "starter unit effect."

Metabolic engineering of the E. coli l-phenylalanine pathway for the production of d-phenylglycine (d-Phg)

D-phenylglycine (D-Phg) is an important side chain building block for semi-synthetic penicillins and cephalosporins such as ampicillin and cephalexin. To produce d-Phg ultimately from glucose, metabolic engineering was applied. Starting from phenylpyruvate, which is the direct precursor of L-phenylalanine, an artificial D-Phg biosynthesis pathway was created. This three-step route is composed of the enzymes hydroxymandelate synthase (HmaS), hydroxymandelate oxidase (Hmo), and the stereoinverting hydroxyphenylglycine aminotransferase (HpgAT). Together they catalyse the conversion of phenylpyruvate via mandelate and phenylglyoxylate to D-Phg. The corresponding genes were obtained from Amycolatopsis orientalis, Streptomyces coelicolor, and Pseudomonas putida. Combined expression of these activities in E. coli strains optimized for the production of L-phenylalanine resulted in the first completely fermentative production of D-Phg.

New Images Evoke FAScinating Questions

Two recent papers in Science reported the X-ray structures of the large, organizationally distinct animal and fungal fatty acid synthases at 5 A. These new structural insights have unexpected implications for enzyme function for the other "iterative" and "assembly line" megasynthases.

Rational strain improvement for enhanced clavulanic acid production by genetic engineering of the glycolytic pathway in Streptomyces clavuligerus

Clavulanic acid is a potent beta-lactamase inhibitor used to combat resistance to penicillin and cephalosporin antibiotics. There is a demand for high-yielding fermentation strains for industrial production of this valuable product. Clavulanic acid biosynthesis is initiated by the condensation of L-arginine and D-glyceraldehyde-3-phosphate (G3P). To overcome the limited G3P pool and improve clavulanic acid production, we genetically engineered the glycolytic pathway in Streptomyces clavuligerus. Two genes (gap1 and gap2) whose protein products are distinct glyceraldehyde-3-phosphate dehydrogenases (GAPDHs) were inactivated in S. clavuligerus by targeted gene disruption. A doubled production of clavulanic acid was consistently obtained when gap1 was disrupted, and reversed by complementation. Addition of arginine to the cultured mutant further improved clavulanic acid production giving a greater than 2-fold increase over wild type, suggesting that arginine became limiting for biosynthesis. This is the first reported application of genetic engineering to channel precursor flux to improve clavulanic acid production.

Hexanoate synthase, a specialized type I fatty acid synthase in aflatoxin B1 biosynthesis

In fungi, fatty acids are biosynthesized by large multifunctional enzyme complexes, the fatty acid synthases (FASs), which catalyze chain assembly in an iterative manner. Many fungal secondary metabolites contain fatty acid moieties, and it is often unclear whether they are recruited from primary metabolism or are biosynthesized de novo by secondary metabolic FASs. The most convincing evidence of such a dedicated FAS comes from the biosyntheses of aflatoxin (AF) and sterigmatocystin (ST) in certain species of the filamentous fungus Aspergillus. Incorporation studies in AF and genetic analyses of ST and AF biosynthesis strongly suggest that their biosyntheses begin with the production of a C6 fatty acid by a specialized FAS. The genes encoding the alpha (hexA) and beta (hexB) subunits of this hexanoate synthase (HexS) from the AF pathway in Aspergillus parsiticus SU-1 were cloned and both their gDNAs and cDNAs were sequenced and their transcriptional ends analyzed. Translated amino acid sequences are predicted to result in proteins of 181.3 and 210.5 kDa, for HexA and HexB, respectively. Comparison of the HexA and HexB sequences with those of the ST FAS subunits and primary metabolic FASs indicated that the secondary metabolic enzymes are members of a well-defined subclass of the FAS family. Phylogenetic predictions and an analysis of GC-bias in AF and ST pathway genes compared with primary metabolic Aspergillus genes were used as a basis to propose a route for the evolution of the AF and ST clusters.