PchB is an isochorismate pyruvate lyase (IPL) involved in siderophore biosynthesis in Pseudomonas aeruginosa. Besides catalyzing the [1,5]-sigmatropic rearrangement of isochorismate, PchB also has weak chorismate mutase (CM) activity, promoting the [3,3]-sigmatropic transformation of chorismate. It has been suggested that the secondary metabolism enzyme PchB evolved from a primary metabolism CM precursor. Here, we employed directed evolution to convert PchB (back) into an efficient CM. A total of seven residues around the active site differing between PchB and a prototypical CM from Escherichia coli were randomized, and the resulting gene library was subjected to selection for CM activity. After growth selection in an auxotrophic strain, a catalyst with 10-fold increased CM activity emerged. The improved enzyme was again randomized at three active site positions and subjected to selection, leading to a PchB variant with a kcat/Km of 96,000 M-1 s-1, which is 40 times higher than that of the parent enzyme and well within the range of dedicated natural CMs. The facile conversion of an IPL into a CM by directed evolution coincides with the fact that both reactions proceed through mechanistically interesting pericyclic processes, reaction types otherwise rarely used by enzymes. When probing typical established CMs for catalytic promiscuity, we discovered spurious IPL activity for the secreted CM from Mycobacterium tuberculosis. Our results hint at active site features, particularly a Val at the bottom of the substrate-binding pocket that may have served as a steppingstone for the evolution of IPL activity in a primordial CM.

NMR spectroscopy of biomolecules provides atomic level information into their structure, dynamics and interactions with their binding partners. However, signal attenuation from line broadening caused by fast relaxation and signal overlap often limits the application of NMR to large macromolecular systems. Here we leverage the slow relaxation properties of 13C nuclei attached to 19F in aromatic 19F-13C spin pairs as well as the spin-spin coupling between the fluorinated 13C nucleus and the hydrogen atom at the meta-position to record two-dimensional 1H-13CF correlation spectra with transverse relaxation-optimized spectroscopy selection on 13CF. To accomplish this, we synthesized [4-19F13Cζ; 3,5-2H2ε] Phe, engineered for optimal relaxation properties, and adapted a residue-specific route to incorporate this residue globally into proteins and a site-specific 4-19F Phe encoding strategy. This approach resulted in narrow linewidths for proteins ranging from 30 kDa to 180 kDa, enabling interaction studies with small-molecule ligands without requiring specialized 19F-compatible probes.

Chorismate mutase (CM) and cyclohexadienyl dehydratase (CDT) catalyze two subsequent reactions in the intracellular biosynthesis of phenylalanine. Surprisingly, exported CMs and CDTs exist in bacterial pathogens. Here, we report the discovery of novel and extremely rare exported bifunctional fusion enzymes, consisting of fused CM and CDT domains. Such enzymes were found in only nine bacterial species belonging to non-pathogenic γ- or β-proteobacteria. In γ-proteobacterial fusion enzymes, the CM domain is N-terminal to the CDT domain, whereas in β-proteobacteria the order is inversed. The CM domains share 15-20% sequence identity with the AroQ γ class CM holotype of Mycobacterium tuberculosis (*MtCM), and the CDT domains 40-60% identity with the exported monofunctional enzyme of Pseudomonas aeruginosa (PheC). In vitro kinetics revealed a K m <7 µM, much lower than for *MtCM, whereas kinetic parameters are similar for CDT domains and PheC. There is no feedback inhibition of CM or CDT by the pathway’s end product Phe, and no catalytic benefit of the domain fusion compared to engineered single-domain constructs. The fusion enzymes of Aequoribacter fuscus , Janthinobacterium sp. HH01, and Duganella sacchari were crystallized and their structures refined to 1.6, 1.7, and 2.4 Å resolution, respectively. Neither the crystal structures nor size-exclusion chromatography show evidence for substrate channeling or higher oligomeric structure that could account for cooperation of CM and CDT active sites. The genetic neighborhood with genes encoding transporter and substrate binding proteins suggests that these exported bifunctional fusion enzymes may participate in signaling systems rather than in the biosynthesis of Phe.

Unlike typical chorismate mutases, the enzyme from Mycobacterium tuberculosis (MtCM) has only low activity on its own. Remarkably, its catalytic efficiency kcat/Km can be boosted more than 100-fold by complex formation with a partner enzyme. Recently, an autonomously fully active MtCM variant was generated using directed evolution, and its structure was solved by X-ray crystallography. However, key residues were involved in crystal contacts, challenging the functional interpretation of the structural changes. Here, we address these challenges by microsecond molecular dynamics simulations, followed up by additional kinetic and structural analyses of selected sets of specifically engineered enzyme variants. A comparison of wild-type MtCM with naturally and artificially activated MtCMs revealed the overall dynamic profiles of these enzymes as well as key interactions between the C-terminus and the active site loop. In the artificially evolved variant of this model enzyme, this loop is preorganized and stabilized by Pro52 and Asp55, two highly conserved residues in typical, highly active chorismate mutases. Asp55 stretches across the active site and helps to appropriately position active site residues Arg18 and Arg46 for catalysis. The role of Asp55 can be taken over by another acidic residue, if introduced at position 88 close to the C-terminus of MtCM, as suggested by molecular dynamics simulations and confirmed by kinetic investigations of engineered variants.

Chorismate mutase (CM), an essential enzyme at the branch-point of the shikimate pathway, is required for the biosynthesis of phenylalanine and tyrosine in bacteria, archaea, plants, and fungi. MtCM, the CM from Mycobacterium tuberculosis, has less than 1% of the catalytic efficiency of a typical natural CM and requires complex formation with 3-deoxy-d-arabino-heptulosonate 7-phosphate synthase for high activity. To explore the full potential of MtCM for catalyzing its native reaction, we applied diverse iterative cycles of mutagenesis and selection, thereby raising kcat/Km 270-fold to 5 × 105m−1s−1, which is even higher than for the complex. Moreover, the evolutionarily optimized autonomous MtCM, which had 11 of its 90 amino acids exchanged, was stabilized compared with its progenitor, as indicated by a 9 °C increase in melting temperature. The 1.5 Å crystal structure of the top-evolved MtCM variant reveals the molecular underpinnings of this activity boost. Some acquired residues (e.g. Pro52 and Asp55) are conserved in naturally efficient CMs, but most of them lie beyond the active site. Our evolutionary trajectories reached a plateau at the level of the best natural enzymes, suggesting that we have exhausted the potential of MtCM. Taken together, these findings show that the scaffold of MtCM, which naturally evolved for mediocrity to enable inter-enzyme allosteric regulation of the shikimate pathway, is inherently capable of high activity.

The rational design of enzymes is an important goal for both fundamental and practical reasons. Here, we describe a process to learn the constraints for specifying proteins purely from evolutionary sequence data, design and build libraries of synthetic genes, and test them for activity in vivo using a quantitative complementation assay. For chorismate mutase, a key enzyme in the biosynthesis of aromatic amino acids, we demonstrate the design of natural-like catalytic function with substantial sequence diversity. Further optimization focuses the generative model toward function in a specific genomic context. The data show that sequence-based statistical models suffice to specify proteins and provide access to an enormous space of functional sequences. This result provides a foundation for a general process for evolution-based design of artificial proteins.

Corynebacterium glutamicum is widely used for the industrial production of amino acids, nucleotides, and vitamins. The shikimate pathway enzymes DAHP synthase (CgDS, Cg2391) and chorismate mutase (CgCM, Cgl0853) play a key role in the biosynthesis of aromatic compounds. Here we show that CgCM requires the formation of a complex with CgDS to achieve full activity, and that both CgCM and CgDS are feedback regulated by aromatic amino acids binding to CgDS. Kinetic analysis showed that Phe and Tyr inhibit CgCM activity by inter-enzyme allostery, whereas binding of Trp to CgDS strongly activates CgCM. Mechanistic insights were gained from crystal structures of the CgCM homodimer, tetrameric CgDS, and the heterooctameric CgCM-CgDS complex, refined to 1.1, 2.5, and 2.2 Å resolution, respectively. Structural details from the allosteric binding sites reveal that DAHP synthase is recruited as the dominant regulatory platform to control the shikimate pathway, similar to the corresponding enzyme complex from Mycobacterium tuberculosis.

DAHP synthase and chorismate mutase catalyze key steps in the shikimate biosynthetic pathway en route to aromatic amino acids. In Mycobacterium tuberculosis, chorismate mutase (MtCM; Rv0948c), located at the branch point toward phenylalanine and tyrosine, has poor activity on its own. However, it is efficiently activated by the first enzyme of the pathway, DAHP synthase (MtDS; Rv2178c), through formation of a non-covalent MtCM–MtDS complex. Here, we show how MtDS serves as an allosteric platform for feedback regulation of both enzymes, using X-ray crystallography, small-angle X-ray scattering, size-exclusion chromatography, and multi-angle light scattering. Crystal structures of the fully inhibited MtDS and the allosterically down-regulated MtCM–MtDS complex, solved at 2.8 and 2.7Å, respectively, reveal how effector binding at the internal MtDS subunit interfaces regulates the activity of MtDS and MtCM. While binding of all three metabolic end products to MtDS shuts down the entire pathway, the binding of phenylalanine jointly with tyrosine releases MtCM from the MtCM–MtDS complex, hence suppressing MtCM activation by 'inter-enzyme allostery'. This elegant regulatory principle, invoking a transient allosteric enzyme interaction, seems to be driven by dynamics and is likely a general strategy used by nature.

Computational methods have been recently applied to the design of protein–protein interfaces. Using this approach, a 61 amino acid long protein called Spider Roll was engineered to recognize the kinase domain of the human p21-activated kinase 1 (PAK1) with good specificity but modest affinity (KD=100μM). Here we show that this artificial protein can be optimized by yeast surface display and fluorescence-activated cell sorting. After three rounds of mutagenesis and screening, a diverse set of tighter binding variants was obtained. A representative binder, MSR7, has a >102-fold higher affinity for PAK1 when displayed on yeast and a 6 to 11-fold advantage when produced free in solution. In contrast to the starting Spider Roll protein, however, MSR7 unexpectedly exhibits characteristics typical of partially disordered proteins, including lower α-helical content, non-cooperative thermal denaturation, and NMR data showing peak broadening and poor signal dispersion. Although conformational disorder is increasingly recognized as an important property of proteins involved in cellular signaling and regulation, it is poorly modeled by current computational methods. Explicit consideration of structural flexibility may improve future protein designs and provide deeper insight into molecular events at protein–protein interfaces.

The shikimate pathway enzyme chorismate mutase converts chorismate into prephenate, a precursor of Tyr and Phe. The intracellular chorismate mutase (MtCM) of Mycobacterium tuberculosis is poorly active on its own, but becomes >100-fold more efficient upon formation of a complex with the first enzyme of the shikimate pathway, 3-deoxy-d-arabino-heptulosonate-7-phosphate synthase (MtDS). The crystal structure of the enzyme complex revealed involvement of C-terminal MtCM residues with the MtDS interface. Here we employed evolutionary strategies to probe the tolerance to substitution of the C-terminal MtCM residues from positions 84-90. Variants with randomized positions were subjected to stringent selection in vivo requiring productive interactions with MtDS for survival. Sequence patterns identified in active library members coincide with residue conservation in natural chorismate mutases of the AroQδ subclass to which MtCM belongs. An Arg-Gly dyad at positions 85 and 86, invariant in AroQδ sequences, was intolerant to mutation, whereas Leu88 and Gly89 exhibited a preference for small and hydrophobic residues in functional MtCM-MtDS complexes. In the absence of MtDS, selection under relaxed conditions identifies positions 84-86 as MtCM integrity determinants, suggesting that the more C-terminal residues function in the activation by MtDS. Several MtCM variants, purified using a novel plasmid-based T7 RNA polymerase gene expression system, showed that a diminished ability to physically interact with MtDS correlates with reduced activatability and feedback regulatory control by Tyr and Phe. Mapping critical protein-protein interaction sites by evolutionary strategies may pinpoint promising targets for drugs that interfere with the activity of protein complexes.

For more than half a century, transition state theory has provided a useful framework for understanding the origins of enzyme catalysis. As proposed by Pauling, enzymes accelerate chemical reactions by binding transition states tighter than substrates, thereby lowering the activation energy compared with that of the corresponding uncatalyzed process. This paradigm has been challenged for chorismate mutase (CM), a well-characterized metabolic enzyme that catalyzes the rearrangement of chorismate to prephenate. Calculations have predicted the decisive factor in CM catalysis to be ground state destabilization rather than transition state stabilization. Using X-ray crystallography, we show, in contrast, that a sluggish variant of Bacillus subtilis CM, in which a cationic active-site arginine was replaced by a neutral citrulline, is a poor catalyst even though it effectively preorganizes chorismate for the reaction. A series of high-resolution molecular snapshots of the reaction coordinate, including the apo enzyme, and complexes with substrate, transition state analog and product, demonstrate that an active site, which is only complementary in shape to a reactive substrate conformer, is insufficient for effective catalysis. Instead, as with other enzymes, electrostatic stabilization of the CM transition state appears to be crucial for achieving high reaction rates.

Evolutionary advances are often fueled by unanticipated innovation. Directed evolution of a computationally designed enzyme suggests that pronounced molecular changes can also drive the optimization of primitive protein active sites. The specific activity of an artificial retro-aldolase was boosted >4,400-fold by random mutagenesis and screening, affording catalytic efficiencies approaching those of natural enzymes. However, structural and mechanistic studies reveal that the engineered catalytic apparatus, consisting of a reactive lysine and an ordered water molecule, was unexpectedly abandoned in favor of a new lysine residue in a substrate-binding pocket created during the optimization process. Structures of the initial in silico design, a mechanistically promiscuous intermediate and one of the most evolved variants highlight the importance of loop mobility and supporting functional groups in the emergence of the new catalytic center. Such internal competition between alternative reactive sites may have characterized the early evolution of many natural enzymes.

The contemporary proteinogenic repertoire contains 20 amino acids with diverse functional groups and side chain geometries. Primordial proteins, in contrast, were presumably constructed from a subset of these building blocks. Subsequent expansion of the proteinogenic alphabet would have enhanced their capabilities, fostering the metabolic prowess and organismal fitness of early living systems. While the addition of amino acids bearing innovative functional groups directly enhances the chemical repertoire of proteomes, the inclusion of chemically redundant monomers is difficult to rationalize. Here, we studied how a simplified chorismate mutase evolves upon expanding its amino acid alphabet from nine to potentially 20 letters. Continuous evolution provided an enhanced enzyme variant that has only two point mutations, both of which extend the alphabet and jointly improve protein stability by >4 kcal/mol and catalytic activity tenfold. The same, seemingly innocuous substitutions (Ile→Thr, Leu→Val) occurred in several independent evolutionary trajectories. The increase in fitness they confer indicates that building blocks with very similar side chain structures are highly beneficial for fine-tuning protein structure and function.