Directed Evolution Mimics Allosteric Activation by Stepwise Tuning of the Conformational Ensemble

A naturally occurring standalone TrpB enzyme provides insights into allosteric communication within tryptophan synthase

Allosteric regulation of catalytic activity is a widespread property of multi-enzyme complexes. The tryptophan synthase is a prototypical allosteric enzyme where the constituting α (TrpA) and β (TrpB) subunits mutually activate each other in a manner that is incompletely understood. Experimental and computational studies have shown that LBCA-TrpB from the last bacterial common ancestor contains six residues (Res6) distal from the active site that allow for high stand-alone catalytic activity in the absence of a TrpA subunit. In the present study, a database search revealed that Res6 is also present in the extant plTrpB from Pelodictyon luteolum. The plTrpB enzyme showed a high stand-alone activity and only a moderate activation by plTrpA. The replacement of LBCA-Res6 in plTrpB with the consensus residues from a multiple sequence alignment yielded plTrpB-con, which showed a dramatically decreased stand-alone activity but was strongly stimulated by plTrpA. These findings suggest that the effect of these six key allosteric residues is largely independent of the protein context within a specific TrpB enzyme. Analysis of the conformational landscapes of plTrpB and plTrpB-con revealed that plTrpB in isolation displays efficient closure of both the active site and the communication (COMM) domain. In contrast, these catalytically competent states are destabilized in plTrpB-con but can be recovered by the addition of plTrpA. A correlation-based shortest path map (SPM) analysis reveals that the catalytically and allosterically relevant domains-specifically, the COMM domain in TrpB and loops 2 and 6 in TrpA-are tightly interconnected exclusively in plTrpA:plTrpB-con.

Identification and understanding of allostery hotspots in proteins: Integration of deep mutational scanning and multi-faceted computational analyses

Motivated by recent deep mutational scanning (DMS) experiments, we have carried out a diverse set of computations to better understand the distribution and contributions of allostery hotspot residues in a transcription factor, TetR. These include extensive atomistic simulations and free energy computations for different functional states of TetR, machine learning analysis of the DMS data and a statistical thermodynamic model for the experimental induction data for the WT protein and a handful of hotspot mutants. Collectively, these computations provided insights into the structural and energetic basis of allostery in TetR, and the distinct contributions of allostery hotspots. The results highlight that the allostery function (i.e., the induction activity) of TetR can be modulated by perturbing both inter-domain coupling and intra-domain properties, such as the population of the binding-competent conformation of each domain. This mechanistic degeneracy qualitatively explains the broad distribution of allostery hotspots across the protein structure observed in the DMS experiments, and also informs the design of strategies aimed at identifying allostery hotspots. The mechanistic framework and the multi-faceted computational approaches are expected to be applicable to the analysis of other allostery systems, especially those sharing the similar two-domain structural topology, and to the design of allostery modulators.

Advances in biosynthesis of 7-Dehydrocholesterol through de novo cell factory strategies

7-Dehydrocholesterol (7-DHC) is an important sterol for maintaining human health and is present in the skin. After sun exposure, 7-DHC in the skin is converted to vitamin D3 to strengthen the immune system. In recent years, synthetic biology has gained importance due to the effective and efficient production of various important compounds using microorganisms. Despite the understanding of the mechanisms and pathways of 7-DHC biosynthesis, achieving higher production yields remains a significant challenge. This review aims to provide a comprehensive overview of the current state of 7-DHC biosynthesis. Various synthetic strategies including optimization of rate-limiting enzymes, metabolic fluxes, redox balance, and subcellular localization are discussed. Moreover, the role of omics technology in designing important proteins and gene editing techniques for strain modification to efficiently synthesize 7-DHC will also be discussed.

Allostery

Allostery describes the ability of biological macromolecules to transmit signals spatially through the molecule from an allosteric site – a site that is distinct from orthosteric binding sites of primary, endogenous ligands – to the functional or active site. This review starts with a historical overview and a description of the classical example of allostery – hemoglobin – and other well-known examples (aspartate transcarbamoylase, Lac repressor, kinases, G-protein-coupled receptors, adenosine triphosphate synthase, and chaperonin). We then discuss fringe examples of allostery, including intrinsically disordered proteins and inter-enzyme allostery, and the influence of dynamics, entropy, and conformational ensembles and landscapes on allosteric mechanisms, to capture the essence of the field. Thereafter, we give an overview over central methods for investigating molecular mechanisms, covering experimental techniques as well as simulations and artificial intelligence (AI)-based methods. We conclude with a review of allostery-based drug discovery, with its challenges and opportunities: with the recent advent of AI-based methods, allosteric compounds are set to revolutionize drug discovery and medical treatments.

Allostery and Evolution: A Molecular Journey Through the Structural and Dynamical Landscape of an Enzyme Super Family

Allosteric regulation is a powerful mechanism for controlling the efficiency of enzymes. Deciphering the evolutionary mechanisms by which allosteric properties have been acquired in enzymes is of fundamental importance. We used the malate (MalDH) and lactate deydrogenases (LDHs) superfamily as model to elucidate this phenomenon. By introducing a few of mutations associated to the emergence of allosteric LDHs into the non-allosteric MalDH from Methanopyrus kandleri, we have gradually shifted its enzymatic profile toward that typical of allosteric LDHs. We first investigated the process triggering homotropic activation. The structures of the resulting mutants show the typical compact organization of the R-active state of LDHs, but a distorted (T-like) catalytic site demonstrating that they corresponds to hybrid states. Molecular dynamics simulations and free energy calculations confirmed the capability of these mutants to sample the T-inactive state. By adding a final single mutation to fine-tune the flexibility of the catalytic site, we obtained an enzyme with both sigmoid (homotropic) and hyperbolic (heterotropic) substrate activation profiles. Its structure shows a typical extended T-state as in LDHs, whereas its catalytic state has as a restored configuration favorable for catalysis. Free energy calculations indicate that the T and R catalytic site configurations are in an equilibrium that depends on solvent conditions. We observed long-range communication between monomers as required for allosteric activation. Our work links the evolution of allosteric regulation in the LDH/MDH superfamily to the ensemble model of allostery at molecular level, and highlights the important role of the underlying protein dynamics.

Understanding the Directed Evolution of a Natural-like Efficient Artificial Metalloenzyme

The artificial metalloenzyme containing iridium in place of iron along with four directed evolution mutations C317G, T213G, L69V, and V254L in a natural cytochrome P450 presents an important milestone in merging the extraordinary efficiency of biocatalysts with the versatility of small molecule chemical catalysts in catalyzing a new-to-nature carbene insertion reaction. This is a show-stopper enzyme, as it exhibits a catalytic efficiency similar to that of natural enzymes. Despite this remarkable discovery, there is no mechanistic and structural understanding as to why it displays extraordinary efficiency after the incorporation of the four active site mutations by directed evolution methods, which so far has been intractable to any experimental methods. In this study, we have deciphered how directed evolution mutations gradually alter the protein conformational ensemble to populate a catalytically active conformation to boost a multistep catalysis in a natural-like artificial metalloenzyme using large-scale molecular dynamics simulations, rigorous quantum chemical (QM), and multiscale quantum chemical/molecular mechanics (QM/MM) calculations. It reveals how evolution precisely positions the cofactor-substrate in an unusual but effective orientation within a reshaped active site in the catalytically active conformation stabilized by C-H···π interactions from more ordered mutated L69V and V254L residues to achieve preferential transition state stabilization compared to the ground state. This work essentially tracks down in atomistic detail the shift in the conformational ensemble of the highly active conformation from the less efficient single mutant to the most efficient quadruple mutant and offers valuable insights for designing better enzymes. The active conformation correctly reproduces the experimental barrier height and also accounts for the catalytic effect, which is in good agreement with experimental observations. Moreover, this conformation features an unusual bonding interaction in a metal-carbene species that preferentially stabilizes the rate-determining formation of an iridium porphyrin carbene intermediate to render the observed high catalytic rate acceleration. Our study provides crucial insights into the underlying rationale for directed evolution, reports the major catalytic role of nonelectrostatic interactions in enzyme catalysis different from the electrostatic model, and suggests a crucial principle toward designing enzymes with natural efficiency.

Altering Active-Site Loop Dynamics Enhances Standalone Activity of the Tryptophan Synthase Alpha Subunit

The α-subunit (TrpA) of the allosterically regulated bifunctional tryptophan synthase αββα enzyme catalyzes the retro-aldol cleavage of indole-glycerol phosphate (IGP) to d-glyceraldehyde 3-phosphate (G3P) and indole. The activity of the enzyme is highly dependent on the β-subunit (TrpB), which allosterically regulates and activates TrpA for enhanced function. This contrasts with the homologous BX1 enzyme from Zea mays that can catalyze the same reaction as TrpA without requiring the presence of any additional binding partner. In this study, we computationally evaluated and compared the conformational landscapes of the homologous ZmBX1 and ZmTrpA enzymes. Our results indicate that enhanced TrpA standalone activity requires the modulation of the conformational dynamics of two relevant active-site loops, loop 6 and 2, that need to be synchronized for accessing the catalytically activated closed state for IGP cleavage, as well as open states for favoring indole/G3P release. Taking as inspiration the evolutionary blueprint ZmBX1 and using our developed correlation-based tool shortest path map focused on the rate-determining conformational transition leading to the catalytically activated closed state, we computationally designed a variant named ZmTrpASPM4-L6BX1, which displays a 163-fold improvement in catalytic efficiency for the retro-aldol cleavage of IGP. This study showcases the importance of fine-tuning the conformational dynamics of active-site loops for altering and improving function, especially in those cases in which a conformational change is rate determining.

Continuous evolution of user-defined genes at 1 million times the genomic mutation rate

When nature evolves a gene over eons at scale, it produces a diversity of homologous sequences with patterns of conservation and change that contain rich structural, functional, and historical information about the gene. However, natural gene diversity accumulates slowly and likely excludes large regions of functional sequence space, limiting the information that is encoded and extractable. We introduce upgraded orthogonal DNA replication (OrthoRep) systems that radically accelerate the evolution of chosen genes under selection in yeast. When applied to a maladapted biosynthetic enzyme, we obtained collections of extensively diverged sequences with patterns that revealed structural and environmental constraints shaping the enzyme's activity. Our upgraded OrthoRep systems should support the discovery of factors influencing gene evolution, uncover previously unknown regions of fitness landscapes, and find broad applications in biomolecular engineering.

Harnessing conformational dynamics in enzyme catalysis to achieve nature-like catalytic efficiencies: the shortest path map tool for computational enzyme redesign

Enzymes exhibit diverse conformations, as represented in the free energy landscape (FEL). Such conformational diversity provides enzymes with the ability to evolve towards novel functions. The challenge lies in identifying mutations that enhance specific conformational changes, especially if located in distal sites from the active site cavity. The shortest path map (SPM) method, which we developed to address this challenge, constructs a graph based on the distances and correlated motions of residues observed in nanosecond timescale molecular dynamics (MD) simulations. We recently introduced a template based AlphaFold2 (tAF2) approach coupled with 10 nanosecond MD simulations to quickly estimate the conformational landscape of enzymes and assess how the FEL is shifted after mutation. In this study, we evaluate the potential of SPM when coupled with tAF2-MD in estimating conformational heterogeneity and identifying key conformationally-relevant positions. The selected model system is the beta subunit of tryptophan synthase (TrpB). We compare how the SPM pathways differ when integrating tAF2 with different MD simulation lengths from as short as 10 ns until 50 ns and considering two distinct Amber forcefield and water models (ff14SB/TIP3P versus ff19SB/OPC). The new methodology can more effectively capture the distal mutations found in laboratory evolution, thus showcasing the efficacy of tAF2-MD-SPM in rapidly estimating enzyme dynamics and identifying the key conformationally relevant hotspots for computational enzyme engineering.

Elucidation of the stereochemical mechanism of cystathionine γ-lyase reveals how substrate specificity constrains catalysis

Pyridoxal phosphate (PLP)-dependent enzymes play essential roles in metabolism and have found applications for organic synthesis and as enzyme therapeutics. The vinylglycine ketimine (VGK) subfamily hosts a growing set of enzymes that play diverse roles in primary and secondary metabolism. However, the molecular determinates of substrate specificity and the complex acid-base chemistry that enables VGK catalysis remain enigmatic. We use a recently discovered amino acid γ-lyase as a model system to probe catalysis in this enzyme family. We discovered that two stereochemically distinct proton transfer pathways occur. Combined kinetic and spectroscopic analysis revealed that progression through the catalytic cycle is correlated with the presence of an H-bond donor after Cγ of an amino acid substrate, suggesting substrate binding is kinetically coupled to a conformational change. High-resolution X-ray crystallography shows that cystathionine-γ-lyases generate an s-trans intermediate and that this geometry is likely conserved throughout the VGK family. An H-bond acceptor in the active site templates substrate binding but does so by pre-organizing substrates away from catalytically productive orientations. Mutagenesis eliminates this pre-organization, such that there is a relaxation of the substrate specificity, but an increase in k cat for diverse substrates. We exploit this information to perform preparative scale α,β,β-tri-deuteration of polar amino acids. Together, these data untangle a complex mode of substrate specificity and provide a foundation for the future study and applications of VGK enzymes.

The β-subunit of tryptophan synthase is a latent tyrosine synthase

Aromatic amino acids and their derivatives are diverse primary and secondary metabolites with critical roles in protein synthesis, cell structure and integrity, defense and signaling. All de novo aromatic amino acid production relies on a set of ancient and highly conserved chemistries. Here we introduce a new enzymatic transformation for L-tyrosine synthesis by demonstrating that the β-subunit of tryptophan synthase-which natively couples indole and L-serine to form L-tryptophan-can act as a latent 'tyrosine synthase'. A single substitution of a near-universally conserved catalytic residue unlocks activity toward simple phenol analogs and yields exclusive para carbon-carbon bond formation to furnish L-tyrosines. Structural and mechanistic studies show how a new active-site water molecule orients phenols for a nonnative mechanism of alkylation, with additional directed evolution resulting in a net >30,000-fold rate enhancement. This new biocatalyst can be used to efficiently prepare valuable L-tyrosine analogs at gram scales and provides the missing chemistry for a conceptually different pathway to L-tyrosine.

Ultrahigh Throughput Evolution of Tryptophan Synthase in Droplets via an Aptamer Sensor

Tryptophan synthase catalyzes the synthesis of a wide array of noncanonical amino acids and is an attractive target for directed evolution. Droplet microfluidics offers an ultrahigh throughput approach to directed evolution (up to 107 experiments per day), enabling the search for biocatalysts in wider regions of sequence space with reagent consumption minimized to the picoliter volume (per library member). While the majority of screening campaigns in this format on record relied on an optically active reaction product, a new assay is needed for tryptophan synthase. Tryptophan is not fluorogenic in the visible light spectrum and thus falls outside the scope of conventional droplet microfluidic readouts, which are incompatible with UV light detection at high throughput. Here, we engineer a tryptophan DNA aptamer into a sensor to quantitatively report on tryptophan production in droplets. The utility of the sensor was validated by identifying five-fold improved tryptophan synthases from ∼100,000 protein variants. More generally, this work establishes the use of DNA-aptamer sensors with a fluorogenic read-out in widening the scope of droplet microfluidic evolution.

Temporal Resolution of Activity-Related Solvation Dynamics in the TIM Barrel Enzyme Murine Adenosine Deaminase

Murine adenosine deaminase (mADA) is a prototypic system for studying the thermal activation of active site chemistry within the TIM barrel family of enzyme reactions. Previous temperature-dependent hydrogen deuterium exchange studies under various conditions have identified interconnected thermal networks for heat transfer from opposing protein-solvent interfaces to active site residues in mADA. One of these interfaces contains a solvent exposed helix-loop-helix moiety that presents the hydrophobic face of its long α-helix to the backside of bound substrate. Herein we pursue the time and temperature dependence of solvation dynamics at the surface of mADA, for comparison to established kinetic parameters that represent active site chemistry. We first created a modified protein devoid of native tryptophans with close to native kinetic behavior. Single site-specific tryptophan mutants were back inserted into each of the four positions where native tryptophans reside. Measurements of nanosecond fluorescence relaxation lifetimes and Stokes shift decays, that reflect time dependent environmental reoroganization around the photo-excited state of Trp*, display minimal temperature dependences. These regions serve as controls for the behavior of a new single tryptophan inserted into a solvent exposed region near the helix-loop-helix moiety located behind the bound substrate, Lys54Trp. This installed Trp displays a significantly elevated value for Ea ( k Stokes shift ) ; further, when Phe61 within the long helix positioned behind bound substrate is replaced by a series of aliphatic hydrophobic side chains, the trends in Ea ( k Stokes shift ) mirror the earlier reported impact of the same series of function-altering hydrophobic side chains on the activation energy of catalysis, Ea ( k cat ) .The reported experimental findings implicate a solvent initiated and rapid (>ns) protein restructuring that contributes to the enthalpic activation barrier to catalysis in mADA.

Flexible active-site loops fine-tune substrate specificity of hyperthermophilic metallo-oxidases

Hyperthermophilic ('superheat-loving') archaea found in high-temperature environments such as Pyrobaculum aerophilum contain multicopper oxidases (MCOs) with remarkable efficiency for oxidizing cuprous and ferrous ions. In this work, directed evolution was used to expand the substrate specificity of P. aerophilum McoP for organic substrates. Six rounds of error-prone PCR and DNA shuffling followed by high-throughput screening lead to the identification of a hit variant with a 220-fold increased efficiency (kcat/Km) than the wild-type for 2,2'-azino-bis(3-ethylbenzothiazoline-6-sulphonic acid) (ABTS) without compromising its intrinsic activity for metal ions. The analysis of the X-ray crystal structure reveals four proximal mutations close to the T1Cu active site. One of these mutations is within the 23-residues loop that occludes this site, a distinctive feature of prokaryotic MCOs. The increased flexibility of this loop results in an enlarged tunnel and one additional pocket that facilitates bulky substrate-enzyme interactions. These findings underscore the synergy between mutations that modulate the dynamics of the active-site loop enabling enhanced catalytic function. This study highlights the potential of targeting loops close to the T1Cu for engineering improvements suitable for biotechnological applications.

Epistasis arises from shifting the rate-limiting step during enzyme evolution of a β-lactamase

Epistasis, the non-additive effect of mutations, can provide combinatorial improvements to enzyme activity that substantially exceed the gains from individual mutations. Yet the molecular mechanisms of epistasis remain elusive, undermining our ability to predict pathogen evolution and engineer biocatalysts. Here we reveal how directed evolution of a β-lactamase yielded highly epistatic activity enhancements. Evolution selected four mutations that increase antibiotic resistance 40-fold, despite their marginal individual effects (≤2-fold). Synergistic improvements coincided with the introduction of super-stochiometric burst kinetics, indicating that epistasis is rooted in the enzyme's conformational dynamics. Our analysis reveals that epistasis stemmed from distinct effects of each mutation on the catalytic cycle. The initial mutation increased protein flexibility and accelerated substrate binding, which is rate-limiting in the wild-type enzyme. Subsequent mutations predominantly boosted the chemical steps by fine-tuning substrate interactions. Our work identifies an overlooked cause for epistasis: changing the rate-limiting step can result in substantial synergy that boosts enzyme activity.

Continuous evolution of user-defined genes at 1-million-times the genomic mutation rate

When nature maintains or evolves a gene's function over millions of years at scale, it produces a diversity of homologous sequences whose patterns of conservation and change contain rich structural, functional, and historical information about the gene. However, natural gene diversity likely excludes vast regions of functional sequence space and includes phylogenetic and evolutionary eccentricities, limiting what information we can extract. We introduce an accessible experimental approach for compressing long-term gene evolution to laboratory timescales, allowing for the direct observation of extensive adaptation and divergence followed by inference of structural, functional, and environmental constraints for any selectable gene. To enable this approach, we developed a new orthogonal DNA replication (OrthoRep) system that durably hypermutates chosen genes at a rate of >10 -4 substitutions per base in vivo . When OrthoRep was used to evolve a conditionally essential maladapted enzyme, we obtained thousands of unique multi-mutation sequences with many pairs >60 amino acids apart (>15% divergence), revealing known and new factors influencing enzyme adaptation. The fitness of evolved sequences was not predictable by advanced machine learning models trained on natural variation. We suggest that OrthoRep supports the prospective and systematic discovery of constraints shaping gene evolution, uncovering of new regions in fitness landscapes, and general applications in biomolecular engineering.

Improving the Activity of Tryptophan Synthetase via a Nucleic Acid Scaffold

Tryptophan synthetase (TSase), which functions as a tetramer, is a typical enzyme with a substrate channel effect, and shows excellent performance in the production of non-standard amino acids, histamine, and other biological derivatives. Based on previous work, we fused a mutant CE protein (colistin of E. coli, a polypeptide with antibacterial activity) sequence with the sequence of TSase to explore whether its catalytic activity could be enhanced, and we also analyzed whether the addition of a DNA scaffold was a feasible strategy. Here, dCE (CE protein without DNase activity) protein tags were constructed and fused to the TrapA and TrapB subunits of TSase, and the whole cell was used for the catalytic reaction. The results showed that after the dCE protein tag was fused to the TrapB subunit, its whole cell catalytic activity increased by 50%. Next, the two subunits were expressed separately, and the proteins were bound in vitro to ensure equimolar combination between the two subunits. After the dCE label was fused to TrapB, the activity of TSase assembled with TrapA also improved. A series of experiments revealed that the enzyme fused with dCE9 showed higher activity than the wild-type protein. In general, the activity of assembly TSase was optimal when the temperature was 50 °C and the pH was about 9.0. After a long temperature treatment, the enzyme maintained good activity. With the addition of exogenous nucleic acid, the activity of the enzyme increased. The maximum yield was 0.58 g/L, which was almost three times that of the wild-type TSase (0.21 g/L). The recombinant TSase constructed in this study with dCE fusion had the advantages of higher heat resistance and higher activity, and confirmed the feasibility of adding a nucleic acid scaffold, providing a new idea for the improvement of structurally similar enzymes.

Engineered Biocatalytic Synthesis of β-N-Substituted-α-Amino Acids

Non-canonical amino acids (ncAAs) are useful synthons for the development of new medicines, materials, and probes for bioactivity. Recently, enzyme engineering has been leveraged to produce a suite of highly active enzymes for the synthesis of β-substituted amino acids. However, there are few examples of biocatalytic N-substitution reactions to make α,β-diamino acids. In this study, we used directed evolution to engineer the β-subunit of tryptophan synthase, TrpB, for improved activity with diverse amine nucleophiles. Mechanistic analysis shows that high yields are hindered by product re-entry into the catalytic cycle and subsequent decomposition. Additional equivalents of l-serine can inhibit product reentry through kinetic competition, facilitating preparative scale synthesis. We show β-substitution with a dozen aryl amine nucleophiles, including demonstration on a gram scale. These transformations yield an underexplored class of amino acids that can serve as unique building blocks for chemical biology and medicinal chemistry.

Allostery, engineering and inhibition of tryptophan synthase

The final two steps of tryptophan biosynthesis are catalyzed by the enzyme tryptophan synthase (TS), composed of alpha (αTS) and beta (βTS) subunits. Recently, experimental and computational methods have mapped "allosteric networks" that connect the αTS and βTS active sites. In αTS, allosteric networks change across the catalytic cycle, which might help drive the conformational changes associated with its function. Directed evolution studies to increase catalytic function and expand the substrate profile of stand-alone βTS have also revealed the importance of αTS in modulating the conformational changes in βTS. These studies also serve as a foundation for the development of TS inhibitors, which can find utility against Mycobacterium tuberculosis and other bacterial pathogens.

Computational remodeling of an enzyme conformational landscape for altered substrate selectivity

Structural plasticity of enzymes dictates their function. Yet, our ability to rationally remodel enzyme conformational landscapes to tailor catalytic properties remains limited. Here, we report a computational procedure for tuning conformational landscapes that is based on multistate design of hinge-mediated domain motions. Using this method, we redesign the conformational landscape of a natural aminotransferase to preferentially stabilize a less populated but reactive conformation and thereby increase catalytic efficiency with a non-native substrate, resulting in altered substrate selectivity. Steady-state kinetics of designed variants reveals activity increases with the non-native substrate of approximately 100-fold and selectivity switches of up to 1900-fold. Structural analyses by room-temperature X-ray crystallography and multitemperature nuclear magnetic resonance spectroscopy confirm that conformational equilibria favor the target conformation. Our computational approach opens the door to targeted alterations of conformational states and equilibria, which should facilitate the design of biocatalysts with customized activity and selectivity.

Stereoselective amino acid synthesis by synergistic photoredox-pyridoxal radical biocatalysis

Developing synthetically useful enzymatic reactions that are not known in biochemistry and organic chemistry is an important challenge in biocatalysis. Through the synergistic merger of photoredox catalysis and pyridoxal 5'-phosphate (PLP) biocatalysis, we developed a pyridoxal radical biocatalysis approach to prepare valuable noncanonical amino acids, including those bearing a stereochemical dyad or triad, without the need for protecting groups. Using engineered PLP enzymes, either enantiomeric product could be produced in a biocatalyst-controlled fashion. Synergistic photoredox-pyridoxal radical biocatalysis represents a powerful platform with which to discover previously unknown catalytic reactions and to tame radical intermediates for asymmetric catalysis.

Allosteric Regulation of Glycogen Phosphorylase by Order/Disorder Transition of the 250' and 280s Loops

Allostery is a fundamental mechanism of protein activation, yet the precise dynamic changes that underlie functional regulation of allosteric enzymes, such as glycogen phosphorylase (GlyP), remain poorly understood. Despite being the first allosteric enzyme described, its structural regulation is still a challenging problem: the key regulatory loops of the GlyP active site (250' and 280s) are weakly stable and often missing density or have large b-factors in structural models. This led to the longstanding hypothesis that GlyP regulation is achieved through gating of the active site by (dis)order transitions, as first proposed by Barford and Johnson. However, testing this requires a quantitative measurement of weakly stable local structure which, to date, has been technically challenging in such a large protein. Hydrogen-deuterium-exchange mass spectrometry (HDX-MS) is a powerful tool for studying protein dynamics, and millisecond HDX-MS has the ability to measure site-localized stability differences in weakly stable structures, making it particularly valuable for investigating allosteric regulation in GlyP. Here, we used millisecond HDX-MS to measure the local structural perturbations of glycogen phosphorylase b (GlyPb), the phosphorylated active form (GlyPa), and the inhibited glucose-6 phosphate complex (GlyPb:G6P) at near-amino acid resolution. Our results support the Barford and Johnson hypothesis for GlyP regulation by providing insight into the dynamic changes of the key regulatory loops.

Magnetic-Field Dependence of LC-Photo-CIDNP in the Presence of Target Molecules Carrying a Quasi-Isolated Spin Pair

NMR spectroscopy is well known for its superb resolution, especially at high applied magnetic field. However, the sensitivity of this technique is very low. Liquid-state low-concentration photo-chemically-induced dynamic nuclear polarization (LC-photo-CIDNP) is a promising emerging methodology capable of enhancing NMR sensitivity in solution. LC-photo-CIDNP works well on solvent-exposed Trp and Tyr residues, either in isolation or within proteins. This study explores the magnetic-field dependence of the LC-photo-CIDNP experienced by two tryptophan isotopologs in solution upon in situ LED-mediated optical irradiation. Out of the two uniformly 13C,15N-labeled Trp (Trp-U-13C,15N) and Trp-α-13C-β,β,2,4,5,6,7-d7 species employed here, only the latter bears a quasi-isolated 1Hα-13Cα spin pair. Computer simulations of the predicted polarization due to geminate recombination of both species display a roughly bell-shaped field dependence. However, while Trp-U-13C,15N is predicted to show a maximum at ca. 500 MHz (11.7 T) and a fairly weak field dependence, Trp-α-13C-β,β,2,4,5,6,7-d7 is expected to display a much sharper field dependence accompanied by a dramatic polarization increase at lower field (ca. 200 MHz, 4.7 T). Experimental LC-photo-CIDNP studies on both Trp isotopologs at 1μM concentration, performed at selected fields, are consistent with the theoretical predictions. In summary, this study highlights the prominent field-dependence of LC-photo-CIDNP enhancements (ε ) experienced by Trp isotopologs bearing a quasi-isolated spin pair.

Bypassing evolutionary dead ends and switching the rate-limiting step of a human immunotherapeutic enzyme

The Trp metabolite kynurenine (KYN) accumulates in numerous solid tumours and mediates potent immunosuppression. Bacterial kynureninases (KYNases), which preferentially degrade kynurenine, can relieve immunosuppression in multiple cancer models, but immunogenicity concerns preclude their clinical use, while the human enzyme (HsKYNase) has very low activity for kynurenine and shows no therapeutic effect. Using fitness selections, we evolved a HsKYNase variant with 27-fold higher activity, beyond which exploration of >30 evolutionary trajectories involving the interrogation of >109 variants led to no further improvements. Introduction of two amino acid substitutions conserved in bacterial KYNases reduced enzyme fitness but potentiated rapid evolution of variants with ~500-fold improved activity and reversed substrate specificity, resulting in an enzyme capable of mediating strong anti-tumour effects in mice. Pre-steady-state kinetics revealed a switch in rate-determining step attributable to changes in both enzyme structure and conformational dynamics. Apart from its clinical significance, our work highlights how rationally designed substitutions can potentiate trajectories that overcome barriers in protein evolution.

Estimating conformational heterogeneity of tryptophan synthase with a template-based Alphafold2 approach

The three-dimensional structure of the enzymes provides very relevant information on the arrangement of the catalytic machinery and structural elements gating the active site pocket. The recent success of the neural network Alphafold2 in predicting the folded structure of proteins from the primary sequence with high levels of accuracy has revolutionized the protein design field. However, the application of Alphafold2 for understanding and engineering function directly from the obtained single static picture is not straightforward. Indeed, understanding enzymatic function requires the exploration of the ensemble of thermally accessible conformations that enzymes adopt in solution. In the present study, we evaluate the potential of Alphafold2 in assessing the effect of the mutations on the conformational landscape of the beta subunit of tryptophan synthase (TrpB). Specifically, we develop a template-based Alphafold2 approach for estimating the conformational heterogeneity of several TrpB enzymes, which is needed for enhanced stand-alone activity. Our results show the potential of Alphafold2, especially if combined with molecular dynamics simulations, for elucidating the changes induced by mutation in the conformational landscapes at a rather reduced computational cost, thus revealing its plausible application in computational enzyme design.

Substrate multiplexed protein engineering facilitates promiscuous biocatalytic synthesis

Enzymes with high activity are readily produced through protein engineering, but intentionally and efficiently engineering enzymes for an expanded substrate scope is a contemporary challenge. One approach to address this challenge is Substrate Multiplexed Screening (SUMS), where enzyme activity is measured on competing substrates. SUMS has long been used to rigorously quantitate native enzyme specificity, primarily for in vivo settings. SUMS has more recently found sporadic use as a protein engineering approach but has not been widely adopted by the field, despite its potential utility. Here, we develop principles of how to design and interpret SUMS assays to guide protein engineering. This rich information enables improving activity with multiple substrates simultaneously, identifies enzyme variants with altered scope, and indicates potential mutational hot-spots as sites for further engineering. These advances leverage common laboratory equipment and represent a highly accessible and customizable method for enzyme engineering.

Efficient engineering of enzymes for expanded substrate scope is currently challenging. Here, the authors develop simple principles of how to design and interpret Substrate Multiplexed Screening assays to guide protein engineering to enable activity improvements with simultaneously with multiple substrates.

Functional roles of enzyme dynamics in accelerating active site chemistry: Emerging techniques and changing concepts

With the growing acceptance of the contribution of protein conformational ensembles to enzyme catalysis and regulation, research in the field of protein dynamics has shifted toward an understanding of the atomistic properties of protein dynamical networks and the mechanisms and time scales that control such behavior. A full description of an enzymatic reaction coordinate is expected to extend beyond the active site and include site-specific networks that communicate with the protein/water interface. Advances in experimental tools for the spatial resolution of thermal activation pathways are being complemented by biophysical methods for visualizing dynamics in real time. An emerging multidimensional model integrates the impacts of bound substrate/effector on the distribution of protein substates that are in rapid equilibration near room temperature with reaction-specific protein embedded heat transfer conduits.

C-N Bond Forming Radical Rebound Is the Enantioselectivity-Determining Step in P411-Catalyzed Enantioselective C(sp3)-H Amination: A Combined Computational and Experimental Investigation

Engineered metalloenzymes represent promising catalysts for stereoselective C-H functionalization reactions. Recently, P450 enzymes have been evolved to allow for new-to-nature intramolecular C(sp3)-H amination reactions via a nitrene transfer mechanism, giving rise to diamine derivatives with excellent enantiocontrol. To shed light on the origin of enantioselectivity, a combined computational and experimental study was carried out. Hybrid quantum mechanics/molecular mechanics calculations were performed to investigate the activation energies and enantioselectivities of both the hydrogen atom transfer (HAT) and the subsequent C-N bond forming radical rebound steps. Contrary to previously hypothesized enantioinduction mechanisms, our calculations show that the radical rebound step is enantioselectivity-determining, whereas the preceding HAT step is only moderately stereoselective. Furthermore, the selectivity in the initial HAT is ablated by rapid conformational change of the radical intermediate prior to C-N bond formation. This finding is corroborated by our experimental study using a set of enantiomerically pure, monodeuterated substrates. Furthermore, classical and ab initio molecular dynamics simulations were carried out to investigate the conformational flexibility of the carbon-centered radical intermediate. This key radical species undergoes a facile conformational change in the enzyme active site from the pro-(R) to the pro-(S) configuration, whereas the radical rebound is slower due to the spin-state change and ring strain of the cyclization process, thereby allowing stereoablative C-N bond formation. Together, these studies revealed an underappreciated enantioinduction mechanism in biocatalytic C(sp3)-H functionalizations involving radical intermediates, opening up new avenues for the development of other challenging asymmetric C(sp3)-H functionalizations.

Rational Design of P450 aMOx for Improving Anti-Markovnikov Selectivity Based on the “Butterfly” Model

Aromatic aldehydes are important industrial raw materials mainly synthesized by anti-Markovnikov (AM) oxidation of corresponding aromatic olefins. The AM product selectivity remains a big challenge. P450 aMOx is the first reported enzyme that could catalyze AM oxidation of aromatic olefins. Here, we reported a rational design strategy based on the “butterfly” model of the active site of P450 aMOx. Constrained molecular dynamic simulations and a binding energy analysis of key residuals combined with an experimental alanine scan were applied. As a result, the mutant A275G showed high AM selectivity of >99%. The results also proved that the “butterfly” model is an effective design strategy for enzymes.

l-Serine Biosensor-Controlled Fermentative Production of l-Tryptophan Derivatives by Corynebacterium glutamicum

Simple Summary

l-tryptophan is an amino acid found in proteins. Its derivatives, such as hydroxylated or halogenated l-tryptophans, find applications in the chemical and pharmaceutical industries, for example, in therapeutic peptides. Biotechnology provides a sustainable way for the production of l-tryptophan and its derivatives. In the final reaction of l-tryptophan biosynthesis in bacteria, such as Corynebacterium glutamicum, another amino acid, l-serine, is incorporated. Here, we show that C. glutamicum TrpB is able to convert indole derivatives, which were added to cells synthesizing l-serine, to the corresponding l-tryptophan derivatives. The gene trpB was expressed under the control of the l-serine-responsive transcriptional activator SerR in the C. glutamicum cells engineered for this fermentation process.

Abstract

l-Tryptophan derivatives, such as hydroxylated or halogenated l-tryptophans, are used in therapeutic peptides and agrochemicals and as precursors of bioactive compounds, such as serotonin. l-Tryptophan biosynthesis depends on another proteinogenic amino acid, l-serine, which is condensed with indole-3-glycerophosphate by tryptophan synthase. This enzyme is composed of the α-subunit TrpA, which catalyzes the retro-aldol cleavage of indole-3-glycerol phosphate, yielding glyceraldehyde-3-phosphate and indole, and the β-subunit TrpB that catalyzes the β-substitution reaction between indole and l-serine to water and l-tryptophan. TrpA is reported as an allosteric actuator, and its absence severely attenuates TrpB activity. In this study, however, we showed that Corynebacterium glutamicum TrpB is catalytically active in the absence of TrpA. Overexpression of C. glutamicumtrpB in a trpBA double deletion mutant supported growth in minimal medium only when exogenously added indole was taken up into the cell and condensed with intracellularly synthesized l-serine. The fluorescence reporter gene of an l-serine biosensor, which was based on the endogenous transcriptional activator SerR and its target promoter P_serE, was replaced by trpB. This allowed for l-serine-dependent expression of trpB in an l-serine-producing strain lacking TrpA. Upon feeding of the respective indole derivatives, this strain produced the l-tryptophan derivatives 5-hydroxytryptophan, 7-bromotryptophan, and 5-fluorotryptophan.

Distal Mutations Shape Substrate-Binding Sites during Evolution of a Metallo-Oxidase into a Laccase

Laccases are in increasing demand as innovative solutions in the biorefinery fields. Here, we combine mutagenesis with structural, kinetic, and in silico analyses to characterize the molecular features that cause the evolution of a hyperthermostable metallo-oxidase from the multicopper oxidase family into a laccase (k _cat 273 s^-1 for a bulky aromatic substrate). We show that six mutations scattered across the enzyme collectively modulate dynamics to improve the binding and catalysis of a bulky aromatic substrate. The replacement of residues during the early stages of evolution is a stepping stone for altering the shape and size of substrate-binding sites. Binding sites are then fine-tuned through high-order epistasis interactions by inserting distal mutations during later stages of evolution. Allosterically coupled, long-range dynamic networks favor catalytically competent conformational states that are more suitable for recognizing and stabilizing the aromatic substrate. This work provides mechanistic insight into enzymatic and evolutionary molecular mechanisms and spots the importance of iterative experimental and computational analyses to understand local-to-global changes.

Directed evolution driving the generation of an efficient keratinase variant to facilitate the feather degradation

Keratinases can specifically degrade keratins, which widely exist in hair, horns, claws and human skin. There is a great interest in developing keratinase to manage keratin waste generated by the poultry industry and reusing keratin products in agriculture, medical treatment and feed industries. Degradation of keratin waste by keratinase is more environmentally friendly and more sustainable compared with chemical and physical methods. However, the wild-type keratinase-producing strains usually cannot meet the requirements of industrial production, and some are pathogenic, limiting their development and utilization. The main purpose of this study is to improve the catalytic performance of keratinase via directed evolution technology for the degradation of feathers. We first constructed a mutant library through error-prone PCR and screened variants with enhanced enzyme activity. The keratinase activity was further improved through fermentation conditions optimization and fed-batch strategies in a 7-L bioreactor. As a result, nine mutants with enhanced activity were identified and the highest enzyme activity was improved from 1150 to 8448 U/mL finally. The mutant achieved efficient biodegradation of feathers, increasing the degradation rate from 49 to 88%. Moreover, a large number of amino acids and soluble peptides were obtained as degradation products, which were excellent protein resources to feed. Therefore, the study provided a keratinase mutant with application potential in the management of feather waste and preparation of protein feed additive.

Engineered Tryptophan Synthase Balances Equilibrium Effects and Fast Dynamic Effects

Creating efficient and stable enzymes for catalysis in pharmaceutical and industrial laboratories is an important research goal. Arnold et al. used directed evolution to engineer a natural tryptophan synthase to create a mutant that is operable under laboratory conditions without the need for a natural allosteric effector. The use of directed evolution allows researchers to improve enzymes without understanding the structure-activity relationship. Here, we present a transition path sampling study of a key chemical transformation in the tryptophan synthase catalytic cycle. We observed that while directed evolution does mimic the natural allosteric effect from a stability perspective, fast protein dynamics associated with chemistry has been dramatically altered. This work provides further evidence of the role of protein dynamics in catalysis and clearly demonstrates the multifaceted complexity of mutations associated with protein engineering. This study also demonstrates a fascinating contrast between allosteric and stand-alone functions at the femtosecond time scale.

Imaging active site chemistry and protonation states: NMR crystallography of the tryptophan synthase α-aminoacrylate intermediate

NMR-assisted crystallography-the integrated application of solid-state NMR, X-ray crystallography, and first-principles computational chemistry-holds significant promise for mechanistic enzymology: by providing atomic-resolution characterization of stable intermediates in enzyme active sites, including hydrogen atom locations and tautomeric equilibria, NMR crystallography offers insight into both structure and chemical dynamics. Here, this integrated approach is used to characterize the tryptophan synthase α-aminoacrylate intermediate, a defining species for pyridoxal-5'-phosphate-dependent enzymes that catalyze β-elimination and replacement reactions. For this intermediate, NMR-assisted crystallography is able to identify the protonation states of the ionizable sites on the cofactor, substrate, and catalytic side chains as well as the location and orientation of crystallographic waters within the active site. Most notable is the water molecule immediately adjacent to the substrate β-carbon, which serves as a hydrogen bond donor to the ε-amino group of the acid-base catalytic residue βLys87. From this analysis, a detailed three-dimensional picture of structure and reactivity emerges, highlighting the fate of the L-serine hydroxyl leaving group and the reaction pathway back to the preceding transition state. Reaction of the α-aminoacrylate intermediate with benzimidazole, an isostere of the natural substrate indole, shows benzimidazole bound in the active site and poised for, but unable to initiate, the subsequent bond formation step. When modeled into the benzimidazole position, indole is positioned with C3 in contact with the α-aminoacrylate Cβ and aligned for nucleophilic attack. Here, the chemically detailed, three-dimensional structure from NMR-assisted crystallography is key to understanding why benzimidazole does not react, while indole does.

Biocatalytic Alkylation Chemistry: Building Molecular Complexity with High Selectivity

Biocatalysis has traditionally been viewed as a field that primarily enables access to chiral centers. This includes the synthesis of chiral alcohols, amines and carbonyl compounds, often through functional group interconversion via hydrolytic or oxidation-reduction reactions. This limitation is partly being overcome by the design and evolution of new enzymes. Here, we provide an overview of a recently thriving research field that we summarize as biocatalytic alkylation chemistry. In the past 3-4 years, numerous new enzymes have been developed that catalyze sp3 C-C/N/O/S bond formations. These enzymes utilize different mechanisms to generate molecular complexity by coupling simple fragments with high activity and selectivity. In many cases, the engineered enzymes perform reactions that are difficult or impossible to achieve with current small-molecule catalysts such as organocatalysts and transition-metal complexes. This review further highlights that the design of new enzyme function is particularly successful when off-the-shelf synthetic reagents are utilized to access non-natural reactive intermediates. This underscores how biocatalysis is gradually moving to a field that build molecules through selective bond forming reactions.

In Silico Identification and Experimental Validation of Distal Activity-Enhancing Mutations in Tryptophan Synthase

Allostery is a central mechanism for the regulation of multi-enzyme complexes. The mechanistic basis that drives allosteric regulation is poorly understood but harbors key information for enzyme engineering. In the present study, we focus on the tryptophan synthase complex that is composed of TrpA and TrpB subunits, which allosterically activate each other. Specifically, we develop a rational approach for identifying key amino acid residues of TrpB distal from the active site. Those residues are predicted to be crucial for shifting the inefficient conformational ensemble of the isolated TrpB to a productive ensemble through intra-subunit allosteric effects. The experimental validation of the conformationally driven TrpB design demonstrates its superior stand-alone activity in the absence of TrpA, comparable to those enhancements obtained after multiple rounds of experimental laboratory evolution. Our work evidences that the current challenge of distal active site prediction for enhanced function in computational enzyme design has become within reach.

Directed Evolution Improves the Enzymatic Synthesis of L-5-Hydroxytryptophan by an Engineered Tryptophan Synthase

L-5-Hydroxytryptophan is an important amino acid that is widely used in food and medicine. In this study, L-5-hydroxytryptophan was synthesized by a modified tryptophan synthase. A direct evolution strategy was applied to engineer tryptophan synthase from Escherichia coli to improve the efficiency of L-5-hydroxytryptophan synthesis. Tryptophan synthase was modified by error-prone PCR. A high-activity mutant enzyme (V231A/K382G) was obtained by a high-throughput screening method. The activity of mutant enzyme (V231A/K382G) is 3.79 times higher than that of its parent, and kcat/Km of the mutant enzyme (V231A/K382G) is 4.36 mM-1∙s-1. The mutant enzyme (V231A/K382G) reaction conditions for the production of L-5-hydroxytryptophan were 100 mmol/L L-serine at pH 8.5 and 35°C for 15 h, reaching a yield of L-5-hydroxytryptophan of 86.7%. Directed evolution is an effective strategy to increase the activity of tryptophan synthase.

A sustainable approach towards the three-component synthesis of unsubstituted 1H-imidazoles in the water at ambient conditions

A green protocol for the synthesis of unsubstituted imidazoles has been demonstrated herein. The reaction is realized using commercially available lipase enzyme, porcine pancreas lipase (PPL) in water. The reaction conditions are selective and mild which helped to tolerate a wide variety of functional groups to give the desired products in good chemical yields.[Formula: see text].

Substitution of a Surface-Exposed Residue Involved in an Allosteric Network Enhances Tryptophan Synthase Function in Cells

Networks of noncovalent amino acid interactions propagate allosteric signals throughout proteins. Tryptophan synthase (TS) is an allosterically controlled bienzyme in which the indole product of the alpha subunit (αTS) is transferred through a 25 Å hydrophobic tunnel to the active site of the beta subunit (βTS). Previous nuclear magnetic resonance and molecular dynamics simulations identified allosteric networks in αTS important for its function. We show here that substitution of a distant, surface-exposed network residue in αTS enhances tryptophan production, not by activating αTS function, but through dynamically controlling the opening of the indole channel and stimulating βTS activity. While stimulation is modest, the substitution also enhances cell growth in a tryptophan-auxotrophic strain of Escherichia coli compared to complementation with wild-type αTS, emphasizing the biological importance of the network. Surface-exposed networks provide new opportunities in allosteric drug design and protein engineering, and hint at potential information conduits through which the functions of a metabolon or even larger proteome might be coordinated and regulated.

Stability of Ligand-induced Protein Conformation Influences Affinity in Maltose-binding Protein

Our understanding of what determines ligand affinity of proteins is poor, even with high-resolution structures available. Both the non-covalent ligand-protein interactions and the relative free energies of available conformations contribute to the affinity of a protein for a ligand. Distant, non-binding site residues can influence the ligand affinity by altering the free energy difference between a ligand-free and ligand-bound conformation. Our hypothesis is that when different ligands induce distinct ligand-bound conformations, it should be possible to tweak their affinities by changing the free energies of the available conformations. We tested this idea for the maltose-binding protein (MBP) from Escherichia coli. We used single-molecule Förster resonance energy transfer (smFRET) to distinguish several unique ligand-bound conformations of MBP. We engineered mutations, distant from the binding site, to affect the stabilities of different ligand-bound conformations. We show that ligand affinity can indeed be altered in a conformation-dependent manner. Our studies provide a framework for the tuning of ligand affinity, apart from modifying binding site residues.

Exploring the Allosteric Territory of Protein Function

While the pervasiveness of allostery in proteins is commonly accepted, we further show the generic nature of allosteric mechanisms by analyzing here transmembrane ion-channel viroporin 3a and RNA-dependent RNA polymerase (RdRp) from SARS-CoV-2 along with metabolic enzymes isocitrate dehydrogenase 1 (IDH1) and fumarate hydratase (FH) implicated in cancers. Using the previously developed structure-based statistical mechanical model of allostery (SBSMMA), we share our experience in analyzing the allosteric signaling, predicting latent allosteric sites, inducing and tuning targeted allosteric response, and exploring the allosteric effects of mutations. This, yet incomplete list of phenomenology, forms a complex and unique allosteric territory of protein function, which should be thoroughly explored. We propose a generic computational framework, which not only allows one to obtain a comprehensive allosteric control over proteins but also provides an opportunity to approach the fragment-based design of allosteric effectors and drug candidates. The advantages of allosteric drugs over traditional orthosteric compounds, complemented by the emerging role of the allosteric effects of mutations in the expansion of the cancer mutational landscape and in the increased mutability of viral proteins, leave no choice besides further extensive studies of allosteric mechanisms and their biomedical implications.

The evolution and engineering of enzyme activity through tuning conformational landscapes

Proteins are dynamic molecules whose structures consist of an ensemble of conformational states. Dynamics contribute to protein function and a link to protein evolution has begun to emerge. This increased appreciation for the evolutionary impact of conformational sampling has grown from our developing structural biology capabilities and the exploration of directed evolution approaches, which have allowed evolutionary trajectories to be mapped. Recent studies have provided empirical examples of how proteins can evolve via conformational landscape alterations. Moreover, minor conformational substates have been shown to be involved in the emergence of new enzyme functions as they can become enriched through evolution. The role of remote mutations in stabilizing new active site geometries has also granted insight into the molecular basis underpinning poorly understood epistatic effects that guide protein evolution. Finally, we discuss how the growth of our understanding of remote mutations is beginning to refine our approach to engineering enzymes.

Computational investigations of allostery in aromatic amino acid biosynthetic enzymes

Allostery, in which binding of ligands to remote sites causes a functional change in the active sites, is a fascinating phenomenon observed in enzymes. Allostery can occur either with or without significant conformational changes in the enzymes, and the molecular basis of its mechanism can be difficult to decipher using only experimental techniques. Computational tools for analyzing enzyme sequences, structures, and dynamics can provide insights into the allosteric mechanism at the atomic level. Combining computational and experimental methods offers a powerful strategy for the study of enzyme allostery. The aromatic amino acid biosynthesis pathway is essential in microorganisms and plants. Multiple enzymes involved in this pathway are sensitive to feedback regulation by pathway end products and are known to use allostery to control their activities. To date, four enzymes in the aromatic amino acid biosynthesis pathway have been computationally investigated for their allosteric mechanisms, including 3-deoxy-d-arabino-heptulosonate 7-phosphate synthase, anthranilate synthase, chorismate mutase, and tryptophan synthase. Here we review the computational studies and findings on the allosteric mechanisms of these four enzymes. Results from these studies demonstrate the capability of computational tools and encourage future computational investigations of allostery in other enzymes of this pathway.

Towards Photochromic Azobenzene-Based Inhibitors for Tryptophan Synthase

Light regulation of drug molecules has gained growing interest in biochemical and pharmacological research in recent years. In addition, a serious need for novel molecular targets of antibiotics has emerged presently. Herein, the development of a photocontrollable, azobenzene-based antibiotic precursor towards tryptophan synthase (TS), an essential metabolic multienzyme complex in bacteria, is presented. The compound exhibited moderately strong inhibition of TS in its E configuration and five times lower inhibition strength in its Z configuration. A combination of biochemical, crystallographic, and computational analyses was used to characterize the inhibition mode of this compound. Remarkably, binding of the inhibitor to a hitherto-unconsidered cavity results in an unproductive conformation of TS leading to noncompetitive inhibition of tryptophan production. In conclusion, we created a promising lead compound for combatting bacterial diseases, which targets an essential metabolic enzyme, and whose inhibition strength can be controlled with light.

Biomolecular QM/MM Simulations: What Are Some of the "Burning Issues"?

QM/MM simulations have become an indispensable tool in many chemical and biochemical investigations. Considering the tremendous degree of success, including recognition by a 2013 Nobel Prize in Chemistry, are there still "burning challenges" in QM/MM methods, especially for biomolecular systems? In this short Perspective, we discuss several issues that we believe greatly impact the robustness and quantitative applicability of QM/MM simulations to many, if not all, biomolecules. We highlight these issues with observations and relevant advances from recent studies in our group and others in the field. Despite such limited scope, we hope the discussions are of general interest and will stimulate additional developments that help push the field forward in meaningful directions.

l-Threonine Transaldolase Activity Is Enabled by a Persistent Catalytic Intermediate

l-Threonine transaldolases (lTTAs) are a poorly characterized class of pyridoxal-5'-phosphate (PLP) dependent enzymes responsible for the biosynthesis of diverse β-hydroxy amino acids. Here, we study the catalytic mechanism of ObiH, an lTTA essential for biosynthesis of the β-lactone natural product obafluorin. Heterologously expressed ObiH purifies as a mixture of chemical states including a catalytically inactive form of the PLP cofactor. Photoexcitation of ObiH promotes the conversion of the inactive state of the enzyme to the active form. UV-vis spectroscopic analysis reveals that ObiH catalyzes the retro-aldol cleavage of l-threonine to form a remarkably persistent glycyl quinonoid intermediate, with a half-life of ∼3 h. Protonation of this intermediate is kinetically disfavored, enabling on-cycle reactivity with aldehydes to form β-hydroxy amino acids. We demonstrate the synthetic potential of ObiH via the single step synthesis of (2S,3R)-β-hydroxyleucine. To further understand the structural features underpinning this desirable reactivity, we determined the crystal structure of ObiH bound to PLP as the Schiff's base at 1.66 Å resolution. This high-resolution model revealed a unique active site configuration wherein the evolutionarily conserved Asp that traditionally H-bonds to the cofactor is swapped for a neighboring Glu. Molecular dynamics simulations combined with mutagenesis studies indicate that a structural rearrangement is associated with l-threonine entry into the catalytic cycle. Together, these data explain the basis for the unique reactivity of lTTA enzymes and provide a foundation for future engineering and mechanistic analysis.

Tryptophan Synthase: Biocatalyst Extraordinaire

Tryptophan synthase (TrpS) has emerged as a paragon of noncanonical amino acid (ncAA) synthesis and is an ideal biocatalyst for synthetic and biological applications. TrpS catalyzes an irreversible, C-C bond-forming reaction between indole and serine to make l-tryptophan; native TrpS complexes possess fairly broad specificity for indole analogues, but are difficult to engineer to extend substrate scope or to confer other useful properties due to allosteric constraints and their heterodimeric structure. Directed evolution freed the catalytically relevant TrpS β-subunit (TrpB) from allosteric regulation by its TrpA partner and has enabled dramatic expansion of the enzyme's substrate scope. This review examines the long and storied career of TrpS from the perspective of its application in ncAA synthesis and biocatalytic cascades.

Scalable continuous evolution for the generation of diverse enzyme variants encompassing promiscuous activities

Enzyme orthologs sharing identical primary functions can have different promiscuous activities. While it is possible to mine this natural diversity to obtain useful biocatalysts, generating comparably rich ortholog diversity is difficult, as it is the product of deep evolutionary processes occurring in a multitude of separate species and populations. Here, we take a first step in recapitulating the depth and scale of natural ortholog evolution on laboratory timescales. Using a continuous directed evolution platform called OrthoRep, we rapidly evolve the Thermotoga maritima tryptophan synthase β-subunit (TmTrpB) through multi-mutation pathways in many independent replicates, selecting only on TmTrpB's primary activity of synthesizing L-tryptophan from indole and L-serine. We find that the resulting sequence-diverse TmTrpB variants span a range of substrate profiles useful in industrial biocatalysis and suggest that the depth and scale of evolution that OrthoRep affords will be generally valuable in enzyme engineering and the evolution of biomolecular functions.