The Moderately Efficient Enzyme: Futile Encounters and Enzyme Floppiness

A finely balanced order-disorder equilibrium sculpts the folding-binding landscape of an antibiotic sequestering protein

TipA, a MerR family transcription factor from Streptomyces lividans, promotes antibiotic resistance by sequestering broad-spectrum thiopeptide-based antibiotics, thus counteracting their inhibitory effect on ribosomes. TipAS, a minimal binding motif which is expressed as an isoform of TipA, harbors a partially disordered N-terminal subdomain that folds upon binding multiple antibiotics. The extent and nature of the underlying molecular heterogeneity in TipAS that shapes its promiscuous folding-function landscape is an open question and is critical for understanding antibiotic-sequestration mechanisms. Here, combining equilibrium and time-resolved experiments, statistical modeling, and simulations, we show that the TipAS native ensemble exhibits a pre-equilibrium between binding-incompetent and binding-competent substates, with the fully folded state appearing only as an excited state under physiological conditions. The binding-competent state characterized by a partially structured N-terminal subdomain loses structure progressively in the physiological range of temperatures, swells on temperature increase, and displays slow conformational exchange across multiple conformations. Binding to the bactericidal antibiotic thiostrepton follows a combination of induced-fit and conformational-selection-like mechanisms, via partial binding and concomitant stabilization of the binding-competent substate. These ensemble features are evolutionarily conserved across orthologs from select bacteria that infect humans, underscoring the functional role of partial disorder in the native ensemble of antibiotic-sequestering proteins belonging to the MerR family.

Determinants of CRISPR Cas12a nuclease activation by DNA and RNA targets

The RNA-guided CRISPR-associated (Cas) enzyme Cas12a cleaves specific double-stranded (ds-) or single-stranded (ss-) DNA targets (in cis), unleashing non-specific ssDNA cleavage (in trans). Though this trans-activity is widely coopted for diagnostics, little is known about target determinants promoting optimal enzyme performance. Using quantitative kinetics, we show formation of activated nuclease proceeds via two steps whereby rapid binding of Cas12a ribonucleoprotein to target is followed by a slower allosteric transition. Activation does not require a canonical protospacer-adjacent motif (PAM), nor is utilization of such PAMs predictive of high trans-activity. We identify several target determinants that can profoundly impact activation times, including bases within the PAM (for ds- but not ssDNA targets) and sequences within and outside those complementary to the spacer, DNA topology, target length, presence of non-specific DNA, and ribose backbone itself, uncovering previously uncharacterized cleavage of and activation by RNA targets. The results provide insight into the mechanism of Cas12a activation, with direct implications on the role of Cas12a in bacterial immunity and for Cas-based diagnostics.

A map of the rubisco biochemical landscape

Rubisco is the primary CO2 fixing enzyme of the biosphere yet has slow kinetics. The roles of evolution and chemical mechanism in constraining the sequence landscape of rubisco remain debated. In order to map sequence to function, we developed a massively parallel assay for rubisco using an engineered E. coli where enzyme function is coupled to growth. By assaying >99% of single amino acid mutants across CO2 concentrations, we inferred enzyme velocity and CO2 affinity for thousands of substitutions. We identified many highly conserved positions that tolerate mutation and rare mutations that improve CO2 affinity. These data suggest that non-trivial kinetic improvements are readily accessible and provide a comprehensive sequence-to-function mapping for enzyme engineering efforts.

Protein structure-function continuum model: Emerging nexuses between specificity, evolution, and structure

The rationale for replacing the old binary of structure-function with the trinity of structure, disorder, and function has gained considerable ground in recent years. A continuum model based on the expanded form of the existing paradigm can now subsume importance of both conformational flexibility and intrinsic disorder in protein function. The disorder is actually critical for understanding the protein-protein interactions in many regulatory processes, formation of membrane-less organelles, and our revised notions of specificity as amply illustrated by moonlighting proteins. While its importance in formation of amyloids and function of prions is often discussed, the roles of intrinsic disorder in infectious diseases and protein function under extreme conditions are also becoming clear. This review is an attempt to discuss how our current understanding of protein function, specificity, and evolution fit better with the continuum model. This integration of structure and disorder under a single model may bring greater clarity in our continuing quest for understanding proteins and molecular mechanisms of their functionality.

Discretised Flux Balance Analysis for Reaction-Diffusion Simulation of Single-Cell Metabolism

Metabolites have to diffuse within the sub-cellular compartments they occupy to specific locations where enzymes are, so reactions could occur. Conventional flux balance analysis (FBA), a method based on linear programming that is commonly used to model metabolism, implicitly assumes that all enzymatic reactions are not diffusion-limited though that may not always be the case. In this work, we have developed a spatial method that implements FBA on a grid-based system, to enable the exploration of diffusion effects on metabolism. Specifically, the method discretises a living cell into a two-dimensional grid, represents the metabolic reactions in each grid element as well as the diffusion of metabolites to and from neighbouring elements, and simulates the system as a single linear programming problem. We varied the number of rows and columns in the grid to simulate different cell shapes, and the method was able to capture diffusion effects at different shapes. We then used the method to simulate heterogeneous enzyme distribution, which suggested a theoretical effect on variability at the population level. We propose the use of this method, and its future extensions, to explore how spatiotemporal organisation of sub-cellular compartments and the molecules within could affect cell behaviour.

Theoretical Improvements in Enzyme Efficiency Associated with Noisy Rate Constants and Increased Dissipation

Previous studies have revealed the extraordinarily large catalytic efficiency of some enzymes. High catalytic proficiency is an essential accomplishment of biological evolution. Natural selection led to the increased turnover number, kcat, and enzyme efficiency, kcat/KM, of uni-uni enzymes, which convert a single substrate into a single product. We added or multiplied random noise with chosen rate constants to explore the correlation between dissipation and catalytic efficiency for ten enzymes: beta-galactosidase, glucose isomerase, β-lactamases from three bacterial strains, ketosteroid isomerase, triosephosphate isomerase, and carbonic anhydrase I, II, and T200H. Our results highlight the role of biological evolution in accelerating thermodynamic evolution. The catalytic performance of these enzymes is proportional to overall entropy production-the main parameter from irreversible thermodynamics. That parameter is also proportional to the evolutionary distance of β-lactamases PC1, RTEM, and Lac-1 when natural or artificial evolution produces the optimal or maximal possible catalytic efficiency. De novo enzyme design and attempts to speed up the rate-limiting catalytic steps may profit from the described connection between kinetics and thermodynamics.

The cofactor challenge in synthetic methylotrophy: bioengineering and industrial applications

Methanol is a promising feedstock for industrial bioproduction: it can be produced renewably and has high solubility and limited microbial toxicity. One of the key challenges for its bio-industrial application is the first enzymatic oxidation step to formaldehyde. This reaction is catalysed by methanol dehydrogenases (MDH) that can use NAD⁺, O₂ or pyrroloquinoline quinone (PQQ) as an electron acceptor. While NAD-dependent MDH are simple to express and have the highest energetic efficiency, they exhibit mediocre kinetics and poor thermodynamics at ambient temperatures. O₂-dependent methanol oxidases require high oxygen concentrations, do not conserve energy and thus produce excessive heat as well as toxic H₂O₂. PQQ-dependent MDH provide a good compromise between energy efficiency and good kinetics that support fast growth rates without any drawbacks for process engineering. Therefore, we argue that this enzyme class represents a promising solution for industry and outline engineering strategies for the implementation of these complex systems in heterologous hosts.

Orthogonal glycolytic pathway enables directed evolution of noncanonical cofactor oxidase

Noncanonical cofactor biomimetics (NCBs) such as nicotinamide mononucleotide (NMN+) provide enhanced scalability for biomanufacturing. However, engineering enzymes to accept NCBs is difficult. Here, we establish a growth selection platform to evolve enzymes to utilize NMN+-based reducing power. This is based on an orthogonal, NMN+-dependent glycolytic pathway in Escherichia coli which can be coupled to any reciprocal enzyme to recycle the ensuing reduced NMN+. With a throughput of >106 variants per iteration, the growth selection discovers a Lactobacillus pentosus NADH oxidase variant with ~10-fold increase in NMNH catalytic efficiency and enhanced activity for other NCBs. Molecular modeling and experimental validation suggest that instead of directly contacting NCBs, the mutations optimize the enzyme's global conformational dynamics to resemble the WT with the native cofactor bound. Restoring the enzyme's access to catalytically competent conformation states via deep navigation of protein sequence space with high-throughput evolution provides a universal route to engineer NCB-dependent enzymes.

Formulation mitigations for particle formation induced by enzymatic hydrolysis of polysorbate 20 in protein-based drug products: insights from a full-factorial longitudinal study

Hydrolytic degradation of the polysorbate 20 (PS20) surfactant in protein-based liquid formulations releases free fatty acids (FFAs), which can accumulate to form particles in drug products during real-time (long-term) storage. To identify formulation conditions that mitigate the risk of particle formation, we conducted a longitudinal study using purified recombinant monoclonal antibody (mAb) formulated in 24 conditions. In this real-time stability study at 5 °C, three key formulation parameters—mAb concentration, initial PS20 concentration, and pH—were varied across representative ranges in a full-factorial design. A longitudinal regression analysis was used to evaluate the effects of these parameters and their interactions on PS20 degradation (via measurements of PS20, FFAs, and PS20 ester distribution) and on particle formation (via visible particle observations and subvisible particle counts). The time-dependent onset of visible particles trended with the rise in subvisible particle counts and FFA levels and fall in PS20 concentration. In the ranges studied here, lower mAb concentration and higher initial PS20 concentration delayed the onset of particles, whereas pH had a negligible effect. These observations were consistent with the general trends predicted by our previously published FFA solubility model. Taken together, these findings highlight the complex relationships between formulation parameters, PS20 degradation, and particle formation.

Evolution avoids a pathological stabilizing interaction in the immune protein S100A9

Stability constrains evolution. While much is known about constraints on destabilizing mutations, less is known about the constraints on stabilizing mutations. We recently identified a mutation in the innate immune protein S100A9 that provides insight into such constraints. When introduced into human S100A9, M63F simultaneously increases the stability of the protein and disrupts its natural ability to activate Toll-like receptor 4. Using chemical denaturation, we found that M63F stabilizes a calcium-bound conformation of hS100A9. We then used NMR to solve the structure of the mutant protein, revealing that the mutation distorts the hydrophobic binding surface of hS100A9, explaining its deleterious effect on function. Hydrogen-deuterium exchange (HDX) experiments revealed stabilization of the region around M63F in the structure, notably Phe37. In the structure of the M63F mutant, the Phe37 and Phe63 sidechains are in contact, plausibly forming an edge-face π-stack. Mutating Phe37 to Leu abolished the stabilizing effect of M63F as probed by both chemical denaturation and HDX. It also restored the biological activity of S100A9 disrupted by M63F. These findings reveal that Phe63 creates a molecular staple with Phe37 that stabilizes a nonfunctional conformation of the protein, thus disrupting function. Using a bioinformatic analysis, we found that S100A9 proteins from different organisms rarely have Phe at both positions 37 and 63, suggesting that avoiding a pathological stabilizing interaction indeed constrains S100A9 evolution. This work highlights an important evolutionary constraint on stabilizing mutations, namely, that they must avoid inappropriately stabilizing nonfunctional protein conformations.

Enzyme Characterization of Pro-virulent SntA, a Cell Wall-Anchored Protein of Streptococcus suis, With Phosphodiesterase Activity on cyclic-di-AMP at a Level Suited to Limit the Innate Immune System

Streptococcus suis and Streptococcus agalactiae evade the innate immune system of the infected host by mechanisms mediated by cell wall-anchored proteins: SntA and CdnP, respectively. The former has been reported to interfere with complement responses, and the latter dampens STING-dependent type-I interferon (IFN) response by hydrolysis of bacterial cyclic-di-AMP (c-di-AMP). Both proteins are homologous but, while CdnP has been studied as a phosphohydrolase, the enzyme activities of SntA have not been investigated. The core structure of SntA was expressed in Escherichia coli as a GST-tagged protein that, after affinity purification, was characterized as phosphohydrolase with a large series of substrates. This included 3′-nucleotides, 2′,3′-cyclic nucleotides, cyclic and linear dinucleotides, and a variety of phosphoanhydride or phosphodiester compounds, most of them previously considered as substrates of E. coli CpdB, a periplasmic protein homologous to SntA and CdnP. Catalytic efficiency was determined for each SntA substrate, either by dividing parameters k_cat/K_M obtained from saturation curves or directly from initial rates at low substrate concentrations when saturation curves could not be obtained. SntA is concluded to act as phosphohydrolase on two groups of substrates with efficiencies higher or lower than ≈ 10⁵ M^–1 s^–1 (average value of the enzyme universe). The group with k_cat/K_M ≥ 10⁵ M^–1 s^–1 (good substrates) includes 3′-nucleotides, 2′,3′-cyclic nucleotides, and linear and cyclic dinucleotides (notably c-di-AMP). Compounds showing efficiencies <10⁴ M^–1 s^–1 are considered poor substrates. Compared with CpdB, SntA is more efficient with its good substrates and less efficient with its poor substrates; therefore, the specificity of SntA is more restrictive. The efficiency of the SntA activity on c-di-AMP is comparable with the activity of CdnP that dampens type-I IFN response, suggesting that this virulence mechanism is also functional in S. suis. SntA modeling revealed that Y530 and Y633 form a sandwich with the nitrogen base of nucleotidic ligands in the substrate-binding site. Mutants Y530A-SntA, Y633A-SntA, and Y530A+Y633A-SntA were obtained and kinetically characterized. For orientation toward the catalytic site, one tyrosine is enough, although this may depend on the substrate being attacked. On the other hand, both tyrosines are required for the efficient binding of good SntA substrates.

Diffusion control in biochemical specificity

Biochemical specificity is critical in enzyme function, evolution, and engineering. Here we employ an established kinetic model to dissect the effects of reactant geometry and diffusion on product formation speed and accuracy in the presence of cognate (correct) and near-cognate (incorrect) substrates. Using this steady-state model for spherical geometries, we find that, for distinct kinetic regimes, the speed and accuracy of the reactions are optimized on different regions of the geometric landscape. From this model we deduce that accuracy can be strongly dependent on reactant geometric properties even for chemically limited reactions. Notably, substrates with a specific geometry and reactivity can be discriminated by the enzyme with higher efficacy than others through purely diffusive effects. For similar cognate and near-cognate substrate geometries (as is the case for polymerases or the ribosome), we observe that speed and accuracy are maximized in opposing regions of the geometric landscape. We also show that, in relevant environments, diffusive effects on accuracy can be substantial even far from extreme kinetic conditions. Finally, we find how reactant chemical discrimination and diffusion can be related to simultaneously optimize steady-state flux and accuracy. These results highlight how diffusion and geometry can be employed to enhance reaction speed and discrimination, and similarly how they impose fundamental restraints on these quantities.

Glycan processing in the Golgi: optimal information coding and constraints on cisternal number and enzyme specificity

Many proteins that undergo sequential enzymatic modification in the Golgi cisternae are displayed at the plasma membrane as cell identity markers. The modified proteins, called glycans, represent a molecular code. The fidelity of this glycan code is measured by how accurately the glycan synthesis machinery realises the desired target glycan distribution for a particular cell type and niche. In this paper, we construct a simplified chemical synthesis model to quantitatively analyse the tradeoffs between the number of cisternae, and the number and specificity of enzymes, required to synthesize a prescribed target glycan distribution of a certain complexity to within a given fidelity. We find that to synthesize complex distributions, such as those observed in real cells, one needs to have multiple cisternae and precise enzyme partitioning in the Golgi. Additionally, for fixed number of enzymes and cisternae, there is an optimal level of specificity (promiscuity) of enzymes that achieves the target distribution with high fidelity. The geometry of the fidelity landscape in the multidimensional space of the number and specificity of enzymes, inter-cisternal transfer rates, and number of cisternae, provides a measure for robustness and identifies stiff and sloppy directions. Our results show how the complexity of the target glycan distribution and number of glycosylation enzymes places functional constraints on the Golgi cisternal number and enzyme specificity.

Single-Molecule Sampling of Dihydrofolate Reductase Shows Kinetic Pauses and an Endosteric Effect Linked to Catalysis

The ability to sample multiple reactions on the same single enzyme is important to link rare intermediates with catalysis and to unravel the role of conformational changes. Despite decades of efforts, however, the single-molecule characterization of nonfluorogenic enzymes during multiple catalytic turnovers has been elusive. Here, we show that nanopore currents allow sampling the dynamic exchange between five structural intermediates during E. coli dihydrofolate reductase (DHFR) catalysis. We found that an endosteric effect promotes the binding of the substrate to the enzyme with a specific hierarchy. The chemical step then switched the enzyme from the closed to the occluded conformation, which in turn promotes the release of the reduced cofactor NADP⁺. Unexpectedly, only a few reactive complexes lead to catalysis. Furthermore, second-long catalytic pauses were observed, possibly reflecting an off-path conformation generated during the reaction. Finally, the free energy from multiple cofactor binding events were required to release the product and switch DHFR back to the reactive conformer. This catalytic fueled concerted mechanism is likely to have evolved to improve the catalytic efficiency of DHFR under the high concentrations of NADP⁺ in E. coli and might be a general feature for complex enzymatic reactions where the binding and release of the products must be tightly controlled.

Substrate Dynamics Contribute to Enzymatic Specificity in Human and Bacterial Methionine Adenosyltransferases

Protein conformational changes can facilitate the binding of noncognate substrates and underlying promiscuous activities. However, the contribution of substrate conformational dynamics to this process is comparatively poorly understood. Here, we analyze human (hMAT2A) and Escherichia coli (eMAT) methionine adenosyltransferases that have identical active sites but different substrate specificity. In the promiscuous hMAT2A, noncognate substrates bind in a stable conformation to allow catalysis. In contrast, noncognate substrates sample stable productive binding modes less frequently in eMAT owing to altered mobility in the enzyme active site. Different cellular concentrations of substrates likely drove the evolutionary divergence of substrate specificity in these orthologues. The observation of catalytic promiscuity in hMAT2A led to the detection of a new human metabolite, methyl thioguanosine, that is produced at elevated levels in a cancer cell line. This work establishes that identical active sites can result in different substrate specificity owing to the effects of substrate and enzyme dynamics.

Structural dynamics in the evolution of a bilobed protein scaffold

Proteins conduct numerous complex biological functions by use of tailored structural dynamics. The molecular details of how these emerged from ancestral peptides remains mysterious. How does nature utilize the same repertoire of folds to diversify function? To shed light on this, we analyzed bilobed proteins with a common structural core, which is spread throughout the tree of life and is involved in diverse biological functions such as transcription, enzymatic catalysis, membrane transport, and signaling. We show here that the structural dynamics of the structural core differentiate predominantly via terminal additions during a long-period evolution. This diversifies substrate specificity and, ultimately, biological function.

Novel biophysical tools allow the structural dynamics of proteins and the regulation of such dynamics by binding partners to be explored in unprecedented detail. Although this has provided critical insights into protein function, the means by which structural dynamics direct protein evolution remain poorly understood. Here, we investigated how proteins with a bilobed structure, composed of two related domains from the periplasmic-binding protein–like II domain family, have undergone divergent evolution, leading to adaptation of their structural dynamics. We performed a structural analysis on ∼600 bilobed proteins with a common primordial structural core, which we complemented with biophysical studies to explore the structural dynamics of selected examples by single-molecule Förster resonance energy transfer and Hydrogen–Deuterium exchange mass spectrometry. We show that evolutionary modifications of the structural core, largely at its termini, enable distinct structural dynamics, allowing the diversification of these proteins into transcription factors, enzymes, and extracytoplasmic transport-related proteins. Structural embellishments of the core created interdomain interactions that stabilized structural states, reshaping the active site geometry, and ultimately altered substrate specificity. Our findings reveal an as-yet-unrecognized mechanism for the emergence of functional promiscuity during long periods of evolution and are applicable to a large number of domain architectures.

Resource uptake and the evolution of moderately efficient enzymes

Enzymes speed up reactions that would otherwise be too slow to sustain the metabolism of selfreplicators. Yet, most enzymes seem only moderately efficient, exhibiting kinetic parameters orders of magnitude lower than their expected physically achievable maxima and spanning over surprisingly large ranges of values. Here, we question how these parameters evolve using a mechanistic model where enzyme efficiency is a key component of individual competition for resources. We show that kinetic parameters are under strong directional selection only up to a point, above which enzymes appear to evolve under near-neutrality, thereby confirming the qualitative observation of other modeling approaches. While the existence of a large fitness plateau could potentially explain the extensive variation in enzyme features reported, we show using a population genetics model that such a widespread distribution is an unlikely outcome of evolution on a common landscape, as mutation–selection–drift balance occupy a narrow area even when very moderate biases towards lower efficiency are considered. Instead, differences in the evolutionary context encountered by each enzyme should be involved, such that each evolves on an individual, unique landscape. Our results point to drift and effective population size playing an important role, along with the kinetics of nutrient transporters, the tolerance to high concentrations of intermediate metabolites, and the reversibility of reactions. Enzyme concentration also shapes selection on kinetic parameters, but we show that the joint evolution of concentration and efficiency does not yield extensive variance in evolutionary outcomes when documented costs to protein expression are applied.

The number of catalytic cycles in an enzyme's lifetime and why it matters to metabolic engineering

The continuous replacement of enzymes and other proteins appropriates up to half the maintenance energy budget in microorganisms and plants. High enzyme replacement rates therefore cut the productivity of biosystems ranging from microbial fermentations to crops. However, yardsticks to assess what drives enzyme protein replacement and guidelines on how to reduce it are lacking. Accordingly, we compared enzymes’ life spans across kingdoms using a new yardstick (catalytic cycles until replacement [CCR]) and related CCR to enzyme reaction chemistry. We concluded that 1) many enzymes fail due to collateral damage from the reaction they catalyze, and 2) such damage and its attendant enzyme replacement costs are mitigable by engineering and are therefore promising targets for synthetic biology.

Metabolic engineering uses enzymes as parts to build biosystems for specified tasks. Although a part’s working life and failure modes are key engineering performance indicators, this is not yet so in metabolic engineering because it is not known how long enzymes remain functional in vivo or whether cumulative deterioration (wear-out), sudden random failure, or other causes drive replacement. Consequently, enzymes cannot be engineered to extend life and cut the high energy costs of replacement. Guided by catalyst engineering, we adopted catalytic cycles until replacement (CCR) as a metric for enzyme functional life span in vivo. CCR is the number of catalytic cycles that an enzyme mediates in vivo before failure or replacement, i.e., metabolic flux rate/protein turnover rate. We used estimated fluxes and measured protein turnover rates to calculate CCRs for ∼100–200 enzymes each from Lactococcus lactis, yeast, and Arabidopsis. CCRs in these organisms had similar ranges (<10³ to >10⁷) but different median values (3–4 × 10⁴ in L. lactis and yeast versus 4 × 10⁵ in Arabidopsis). In all organisms, enzymes whose substrates, products, or mechanisms can attack reactive amino acid residues had significantly lower median CCR values than other enzymes. Taken with literature on mechanism-based inactivation, the latter finding supports the proposal that 1) random active-site damage by reaction chemistry is an important cause of enzyme failure, and 2) reactive noncatalytic residues in the active-site region are likely contributors to damage susceptibility. Enzyme engineering to raise CCRs and lower replacement costs may thus be both beneficial and feasible.

A computational workflow for the expansion of heterologous biosynthetic pathways to natural product derivatives

Plant natural products (PNPs) and their derivatives are important but underexplored sources of pharmaceutical molecules. To access this untapped potential, the reconstitution of heterologous PNP biosynthesis pathways in engineered microbes provides a valuable starting point to explore and produce novel PNP derivatives. Here, we introduce a computational workflow to systematically screen the biochemical vicinity of a biosynthetic pathway for pharmaceutical compounds that could be produced by derivatizing pathway intermediates. We apply our workflow to the biosynthetic pathway of noscapine, a benzylisoquinoline alkaloid (BIA) with a long history of medicinal use. Our workflow identifies pathways and enzyme candidates for the production of (S)-tetrahydropalmatine, a known analgesic and anxiolytic, and three additional derivatives. We then construct pathways for these compounds in yeast, resulting in platforms for de novo biosynthesis of BIA derivatives and demonstrating the value of cheminformatic tools to predict reactions, pathways, and enzymes in synthetic biology and metabolic engineering.

Mainly on the Plane: Deep Subsurface Bacterial Proteins Bind and Alter Clathrate Structure

Gas clathrates are both a resource and a hindrance. They store massive quantities of natural gas but also can clog natural gas pipelines, with disastrous consequences. Eco-friendly technologies for controlling and modulating gas clathrate growth are needed. Type I Antifreeze Proteins (AFPs) from cold-water fish have been shown to bind to gas clathrates via repeating motifs of threonine and alanine. We tested whether proteins encoded in the genomes of bacteria native to natural gas clathrates bind to and alter clathrate morphology. We identified putative clathrate-binding proteins (CBPs) with multiple threonine/alanine motifs in a putative operon (cbp) in metagenomes from natural clathrate deposits. We recombinantly expressed and purified five CbpA proteins, four of which were stable, and experimentally confirmed that CbpAs bound to tetrahydrofuran (THF) clathrate, a low-pressure analogue for structure II gas clathrate. When grown in the presence of CbpAs, the THF clathrate was polycrystalline and platelike instead of forming single, octahedral crystals. Two CbpAs yielded branching clathrate crystals, similar to the effect of Type I AFP, while the other two produced hexagonal crystals parallel to the [1 1 1] plane, suggesting two distinct binding modes. Bacterial CBPs may find future utility in industry, such as maintaining a platelike structure during gas clathrate transportation.

Words of advice: teaching enzyme kinetics

Enzymology is concerned with the study of enzyme structure, function, regulation and kinetics. It is an interdisciplinary subject that can be treated as an exclusive sphere of exhaustive inquiry within mathematical, physico-chemical and biological sciences. Hence, teaching of enzymology, in general, and enzyme kinetics, in particular, should be undertaken in an interdisciplinary manner for a holistic appreciation of this subject. Further, analogous examples from everyday life should form an integral component of the teaching for an intuitive grasp of the subject matter. Furthermore, simulation-based appreciation of enzyme kinetics should be preferred over simplifying assumptions and approximations of traditional enzyme kinetics teaching. In this Words of Advice, I outline the domain depth of enzymology across the various disciplines and provide initial ideas on how appropriate analogies can provide firm insights into the subject. Further, I demonstrate how an intuitive feel for the subject can help not only in grasping abstract concepts but also extending it in experimental design and subsequent interpretation. Use of simulations in grasping complex concepts is also advocated given the advantages this medium offers over traditional approaches involving images and molecular models. Furthermore, I discuss the merits of incorporating the historical backdrop of major discoveries in enzymological teaching. We, at AstraZeneca, have experimented with this approach with the desired outcome of generating interest in the subject from people practising diverse disciplines.

A mechanistic view of enzyme evolution

New enzyme functions often evolve through the recruitment and optimization of latent promiscuous activities. How do mutations alter the molecular architecture of enzymes to enhance their activities? Can we infer general mechanisms that are common to most enzymes, or does each enzyme require a unique optimization process? The ability to predict the location and type of mutations necessary to enhance an enzyme's activity is critical to protein engineering and rational design. In this review, via the detailed examination of recent studies that have shed new light on the molecular changes underlying the optimization of enzyme function, we provide a mechanistic perspective of enzyme evolution. We first present a global survey of the prevalence of activity-enhancing mutations and their distribution within protein structures. We then delve into the molecular solutions that mediate functional optimization, specifically highlighting several common mechanisms that have been observed across multiple examples. As distinct protein sequences encounter different evolutionary bottlenecks, different mechanisms are likely to emerge along evolutionary trajectories toward improved function. Identifying the specific mechanism(s) that need to be improved upon, and tailoring our engineering efforts to each sequence, may considerably improve our chances to succeed in generating highly efficient catalysts in the future.

Predicting nitroimidazole antibiotic resistance mutations in Mycobacterium tuberculosis with protein engineering

Our inability to predict which mutations could result in antibiotic resistance has made it difficult to rapidly identify the emergence of resistance, identify pre-existing resistant populations, and manage our use of antibiotics to effectively treat patients and prevent or slow the spread of resistance. Here we investigated the potential for resistance against the new antitubercular nitroimidazole prodrugs pretomanid and delamanid to emerge in Mycobacterium tuberculosis, the causative agent of tuberculosis (TB). Deazaflavin-dependent nitroreductase (Ddn) is the only identified enzyme within M. tuberculosis that activates these prodrugs, via an F420H2-dependent reaction. We show that the native menaquinone-reductase activity of Ddn is essential for emergence from hypoxia, which suggests that for resistance to spread and pose a threat to human health, the native activity of Ddn must be at least partially retained. We tested 75 unique mutations, including all known sequence polymorphisms identified among ~15,000 sequenced M. tuberculosis genomes. Several mutations abolished pretomanid and delamanid activation in vitro, without causing complete loss of the native activity. We confirmed that a transmissible M. tuberculosis isolate from the hypervirulent Beijing family already possesses one such mutation and is resistant to pretomanid, before being exposed to the drug. Notably, delamanid was still effective against this strain, which is consistent with structural analysis that indicates delamanid and pretomanid bind to Ddn differently. We suggest that the mutations identified in this work be monitored for informed use of delamanid and pretomanid treatment and to slow the emergence of resistance.

Defining optimal electron transfer partners for light-driven cytochrome P450 reactions

Plants and cyanobacteria are promising heterologous hosts for metabolic engineering, and particularly suited for expression of cytochrome P450 (P450s), enzymes that catalyse key steps in biosynthetic pathways leading to valuable natural products such as alkaloids, terpenoids and phenylpropanoids. P450s are often difficult to express and require a membrane-bound NADPH-dependent reductase, complicating their use in metabolic engineering and bio-production. We previously demonstrated targeting of heterologous P450s to thylakoid membranes both in N. benthamiana chloroplasts and cyanobacteria, and functional substitution of their native reductases with the photosynthetic apparatus via the endogenous soluble electron carrier ferredoxin. However, because ferredoxin acts as a sorting hub for photosynthetic reducing power, there is fierce competition for reducing equivalents, which limits photosynthesis-driven P450 output. This study compares the ability of four electron carriers to increase photosynthesis-driven P450 activity. These carriers, three plant ferredoxins and a flavodoxin-like engineered protein derived from cytochrome P450 reductase, show only modest differences in their electron transfer to our model P450, CYP79A1 in vitro. However, only the flavodoxin-like carrier supplies appreciable reducing power in the presence of competition for reduced ferredoxin, because it possesses a redox potential that renders delivery of reducing equivalents to endogenous processes inefficient. We further investigate the efficacy of these electron carrier proteins in vivo by expressing them transiently in N. benthamiana fused to CYP79A1. All but one of the fusion enzymes show improved sequestration of photosynthetic reducing power. Fusion with the flavodoxin-like carrier offers the greatest improvement in this comparison - nearly 25-fold on a per protein basis. Thus, this study demonstrates that synthetic electron transfer pathways with optimal redox potentials can alleviate the problem of endogenous competition for reduced ferredoxin and sets out a new metabolic engineering strategy useful for producing valuable natural products.

Increasing the Enzyme Loading Capacity of Porous Supports by a Layer-by-Layer Immobilization Strategy Using PEI as Glue

A new strategy to increase the enzyme-loading capacity of porous supports was investigated. Lipase from Pseudomonas fluorescens (PFL) was immobilized on octyl-agarose (OA) beads and treated with polyethyleneimine (PEI). Then, PFL was immobilized on the previous PFL layer. Next, the biocatalyst was coated with PEI and a third layer of PFL was added. Sodium dodecyl sulfate polyacrylamide electrophoresis showed that the amount of PFL proportionally increased with each enzyme layer; however, the effects on biocatalyst activity were not as clear. Hydrolyzing 50 mM of triacetin at 25 °C, the activity of the three-layer biocatalyst was even lower than that of the bi-layer one; on the contrary its activity was higher when the activity was measured at 4 °C in the presence of 30% acetonitrile (that reduced the activity and thus the relevance of the substrate diffusion limitations). That is, the advantage of the multilayer formation depends on the specific activity of the enzyme and on the diffusion limitations of the substrate. When octyl agarose (OA)-PFL-PEI-PFL preparation was treated with glutaraldehyde, the activity was reduced, although the enzyme stability increased and the immobilization of the last PFL layer offered results similar to the one obtained using the three-layer preparation without glutaraldehyde modification (90%).

Exploiting enzyme promiscuity to shape plant specialized metabolism

The amazing variability of plant metabolism and its rapid divergence during evolution pose fundamental questions as to the driving forces, mechanisms, and players in metabolic differentiation. This review examines concepts that help us understand adaptive pathway evolution, with a particular emphasis on plant specialized metabolism, previously often termed secondary metabolism. Following a general introduction to pathway and metabolite evolution, the focus is directed to enzyme promiscuity and its classification. Promiscuous enzymes (or substrates), 'silent' elements of the metabolome, and the 'underground metabolism' may be used and combined to evolve 'new' metabolic pathways. It appears that new pathways rarely appear from scratch, but instead emerge from 'floppy' enzymes and elements of a 'messy' metabolism, and in this way a range of metabolites is generated, some of which may provide benefits to the plant.

Local frustration around enzyme active sites

Conflicting biological goals often meet in the specification of protein sequences for structure and function. Overall, strong energetic conflicts are minimized in folded native states according to the principle of minimal frustration, so that a sequence can spontaneously fold, but local violations of this principle open up the possibility to encode the complex energy landscapes that are required for active biological functions. We survey the local energetic frustration patterns of all protein enzymes with known structures and experimentally annotated catalytic residues. In agreement with previous hypotheses, the catalytic sites themselves are often highly frustrated regardless of the protein oligomeric state, overall topology, and enzymatic class. At the same time a secondary shell of more weakly frustrated interactions surrounds the catalytic site itself. We evaluate the conservation of these energetic signatures in various family members of major enzyme classes, showing that local frustration is evolutionarily more conserved than the primary structure itself.

Carbon Monoxide Dehydrogenases

Carbon monoxide dehydrogenases (CODHs) catalyze the reversible oxidation of CO with water to CO₂, two electrons, and two protons. Two classes of CODHs exist, having evolved from different scaffolds featuring active sites built from different transition metals. The basic properties of both classes are described in this overview chapter.

Tools and strategies for constructing cell-free enzyme pathways

Single enzyme systems or engineered microbial hosts have been used for decades but the notion of assembling multiple enzymes into cell-free synthetic pathways is a relatively new development. The extensive possibilities that stem from this synthetic concept makes it a fast growing and potentially high impact field for biomanufacturing fine and platform chemicals, pharmaceuticals and biofuels. However, the translation of individual single enzymatic reactions into cell-free multi-enzyme pathways is not trivial. In reality, the kinetics of an enzyme pathway can be very inadequate and the production of multiple enzymes can impose a great burden on the economics of the process. We examine here strategies for designing synthetic pathways and draw attention to the requirements of substrates, enzymes and cofactor regeneration systems for improving the effectiveness and sustainability of cell-free biocatalysis. In addition, we comment on methods for the immobilisation of members of a multi-enzyme pathway to enhance the viability of the system. Finally, we focus on the recent development of integrative tools such as in silico pathway modelling and high throughput flux analysis with the aim of reinforcing their indispensable role in the future of cell-free biocatalytic pathways for biomanufacturing.

Synergistic effects of functionally distinct substitutions in β-lactamase variants shed light on the evolution of bacterial drug resistance

The CTX-M β-lactamases have emerged as the most widespread extended-spectrum β-lactamases (ESBLs) in Gram-negative bacteria. These enzymes rapidly hydrolyze cefotaxime, but not the related cephalosporin, ceftazidime. ESBL variants have evolved, however, that provide enhanced ceftazidime resistance. We show here that a natural variant at a nonactive site, i.e. second-shell residue N106S, enhances enzyme stability but reduces catalytic efficiency for cefotaxime and ceftazidime and decreases resistance levels. However, when the N106S variant was combined with an active-site variant, D240G, that enhances enzyme catalytic efficiency, but decreases stability, the resultant double mutant exhibited higher resistance levels than predicted on the basis of the phenotypes of each variant. We found that this epistasis is due to compensatory effects, whereby increased stability provided by N106S overrides its cost of decreased catalytic activity. X-ray structures of the variant enzymes in complex with cefotaxime revealed conformational changes in the active-site loop spanning residues 103-106 that were caused by the N106S substitution and relieve steric strain to stabilize the enzyme, but also alter contacts with cefotaxime and thereby reduce catalytic activity. We noted that the 103-106 loop conformation in the N106S-containing variants is different from that of WT CTX-M but nearly identical to that of the non-ESBL, TEM-1 β-lactamase, having a serine at the 106 position. Therefore, residue 106 may serve as a "switch" that toggles the conformations of the 103-106 loop. When it is serine, the loop is in the non-ESBL, TEM-like conformation, and when it is asparagine, the loop is in a CTX-M-like, cefotaximase-favorable conformation.

Conformational dynamics and enzyme evolution

Enzymes are dynamic entities, and their dynamic properties are clearly linked to their biological function. It follows that dynamics ought to play an essential role in enzyme evolution. Indeed, a link between conformational diversity and the emergence of new enzyme functionalities has been recognized for many years. However, it is only recently that state-of-the-art computational and experimental approaches are revealing the crucial molecular details of this link. Specifically, evolutionary trajectories leading to functional optimization for a given host environment or to the emergence of a new function typically involve enriching catalytically competent conformations and/or the freezing out of non-competent conformations of an enzyme. In some cases, these evolutionary changes are achieved through distant mutations that shift the protein ensemble towards productive conformations. Multifunctional intermediates in evolutionary trajectories are probably multi-conformational, i.e. able to switch between different overall conformations, each competent for a given function. Conformational diversity can assist the emergence of a completely new active site through a single mutation by facilitating transition-state binding. We propose that this mechanism may have played a role in the emergence of enzymes at the primordial, progenote stage, where it was plausibly promoted by high environmental temperatures and the possibility of additional phenotypic mutations.

Evolution of cyclohexadienyl dehydratase from an ancestral solute-binding protein

The emergence of enzymes through the neofunctionalization of noncatalytic proteins is ultimately responsible for the extraordinary range of biological catalysts observed in nature. Although the evolution of some enzymes from binding proteins can be inferred by homology, we have a limited understanding of the nature of the biochemical and biophysical adaptations along these evolutionary trajectories and the sequence in which they occurred. Here we reconstructed and characterized evolutionary intermediate states linking an ancestral solute-binding protein to the extant enzyme cyclohexadienyl dehydratase. We show how the intrinsic reactivity of a desolvated general acid was harnessed by a series of mutations radiating from the active site, which optimized enzyme-substrate complementarity and transition-state stabilization and minimized sampling of noncatalytic conformations. Our work reveals the molecular evolutionary processes that underlie the emergence of enzymes de novo, which are notably mirrored by recent examples of computational enzyme design and directed evolution.

Rubisco is not really so bad

Ribulose-1,5-bisphosphate carboxylase/oxygenase (Rubisco) is the most widespread carboxylating enzyme in autotrophic organisms. Its kinetic and structural properties have been intensively studied for more than half a century. Yet important aspects of the catalytic mechanism remain poorly understood, especially the oxygenase reaction. Because of its relatively modest turnover rate (a few catalytic events per second) and the competitive inhibition by oxygen, Rubisco is often viewed as an inefficient catalyst for CO₂ fixation. Considerable efforts have been devoted to improving its catalytic efficiency, so far without success. In this review, we re-examine Rubisco's catalytic performance by comparison with other chemically related enzymes. We find that Rubisco is not especially slow. Furthermore, considering both the nature and the complexity of the chemical reaction, its kinetic properties are unremarkable. Although not unique to Rubisco, oxygenation is not systematically observed in enolate and enamine forming enzymes and cannot be considered as an inevitable consequence of the mechanism. It is more likely the result of a compromise between chemical and metabolic imperatives. We argue that a better description of Rubisco mechanism is still required to better understand the link between CO₂ and O₂ reactivity and the rationale of Rubisco diversification and evolution.

Structural and Mechanistic Basis for Extended-Spectrum Drug-Resistance Mutations in Altering the Specificity of TEM, CTX-M, and KPC β-lactamases

The most common mechanism of resistance to β-lactam antibiotics in Gram-negative bacteria is the production of β-lactamases that hydrolyze the drugs. Class A β-lactamases are serine active-site hydrolases that include the common TEM, CTX-M, and KPC enzymes. The TEM enzymes readily hydrolyze penicillins and older cephalosporins. Oxyimino-cephalosporins, such as cefotaxime and ceftazidime, however, are poor substrates for TEM-1 and were introduced, in part, to circumvent β-lactamase-mediated resistance. Nevertheless, the use of these antibiotics has lead to evolution of numerous variants of TEM with mutations that significantly increase the hydrolysis of the newer cephalosporins. The CTX-M enzymes emerged in the late 1980s and hydrolyze penicillins and older cephalosporins and derive their name from the ability to also hydrolyze cefotaxime. The CTX-M enzymes, however, do not efficiently hydrolyze ceftazidime. Variants of CTX-M enzymes, however, have evolved that exhibit increased hydrolysis of ceftazidime. Finally, the KPC enzyme emerged in the 1990s and is characterized by its broad specificity that includes penicillins, most cephalosporins, and carbapenems. The KPC enzyme, however, does not efficiently hydrolyze ceftazidime. As with the TEM and CTX-M enzymes, variants have recently evolved that extend the spectrum of KPC β-lactamase to include ceftazidime. This review discusses the structural and mechanistic basis for the expanded substrate specificity of each of these enzymes that result from natural mutations that confer oxyimino-cephalosporin resistance. For the TEM enzyme, extended-spectrum mutations act by establishing new interactions with the cephalosporin. These mutations increase the conformational heterogeneity of the active site to create sub-states that better accommodate the larger drugs. The mutations expanding the spectrum of CTX-M enzymes also affect the flexibility and conformation of the active site to accommodate ceftazidime. Although structural data are limited, extended-spectrum mutations in KPC may act by mediating new, direct interactions with substrate and/or altering conformations of the active site. In many cases, mutations that expand the substrate profile of these enzymes simultaneously decrease the thermodynamic stability. This leads to the emergence of additional global suppressor mutations that help correct the stability defects leading to increased protein expression and increased antibiotic resistance.

Directions for Optimization of Photosynthetic Carbon Fixation: RuBisCO's Efficiency May Not Be So Constrained After All

The ubiquitous enzyme Ribulose 1,5-bisphosphate carboxylase-oxygenase (RuBisCO) fixes atmospheric carbon dioxide within the Calvin-Benson cycle that is utilized by most photosynthetic organisms. Despite this central role, RuBisCO's efficiency surprisingly struggles, with both a very slow turnover rate to products and also impaired substrate specificity, features that have long been an enigma as it would be assumed that its efficiency was under strong evolutionary pressure. RuBisCO's substrate specificity is compromised as it catalyzes a side-fixation reaction with atmospheric oxygen; empirical kinetic results show a trend to tradeoff between relative specificity and low catalytic turnover rate. Although the dominant hypothesis has been that the active-site chemistry constrains the enzyme's evolution, a more recent study on RuBisCO stability and adaptability has implicated competing selection pressures. Elucidating these constraints is crucial for directing future research on improving photosynthesis, as the current literature casts doubt on the potential effectiveness of site-directed mutagenesis to improve RuBisCO's efficiency. Here we use regression analysis to quantify the relationships between kinetic parameters obtained from empirical data sets spanning a wide evolutionary range of RuBisCOs. Most significantly we found that the rate constant for dissociation of CO₂ from the enzyme complex was much higher than previous estimates and comparable with the corresponding catalytic rate constant. Observed trends between relative specificity and turnover rate can be expressed as the product of negative and positive correlation factors. This provides an explanation in simple kinetic terms of both the natural variation of relative specificity as well as that obtained by reported site-directed mutagenesis results. We demonstrate that the kinetic behaviour shows a lesser rather than more constrained RuBisCO, consistent with growing empirical evidence of higher variability in relative specificity. In summary our analysis supports an explanation for the origin of the tradeoff between specificity and turnover as due to competition between protein stability and activity, rather than constraints between rate constants imposed by the underlying chemistry. Our analysis suggests that simultaneous improvement in both specificity and turnover rate of RuBisCO is possible.

New insights into the molecular characteristics behind the function of Renilla luciferase

Renilla Luciferase (RLuc) is a blue light emitter protein which can be applied as a valuable tool in medical diagnosis. But due to lack of the crystal structure of RLuc-ligand complex, the functional motions and catalytic mechanism of this enzyme remain largely unknown. In the present study, the active site properties and the ligand-receptor interactions of the native RLuc and its red-shifted light emitting variant (Super RLuc 8) were investigated using molecular docking approach, molecular dynamics (MD) analysis, and MM-PBSA method. The detailed analysis of the main clusters led to identifying a lid-like structure and its functional motions. Furthermore, an induced-fit mechanism is proposed where ligand-binding induces conformational changes of the active site. Our findings give an insight into the deeper understanding of RLuc conformational changes during binding steps and ligand-receptor pattern. Moreover, our work broaden our understanding of how active site geometry is adjusted to support the catalytic activity and red-shifted light emission in Super RLuc 8.

Enzyme engineering: reaching the maximal catalytic efficiency peak

The practical need for highly efficient enzymes presents new challenges in enzyme engineering, in particular, the need to improve catalytic turnover (k_cat) or efficiency (k_cat/K_M) by several orders of magnitude. However, optimizing catalysis demands navigation through complex and rugged fitness landscapes, with optimization trajectories often leading to strong diminishing returns and dead-ends. When no further improvements are observed in library screens or selections, it remains unclear whether the maximal catalytic efficiency of the enzyme (the catalytic 'fitness peak') has been reached; or perhaps, an alternative combination of mutations exists that could yield additional improvements. Here, we discuss fundamental aspects of the process of catalytic optimization, and offer practical solutions with respect to overcoming optimization plateaus.

Structural Evidence for the Dopamine-First Mechanism of Norcoclaurine Synthase

Norcoclaurine synthase (NCS) is a Pictet-Spenglerase that catalyzes the first key step in plant benzylisoquinoline alkaloid metabolism, a compound family that includes bioactive natural products such as morphine. The enzyme has also shown great potential as a biocatalyst for the formation of chiral isoquinolines. Here we present new high-resolution X-ray crystallography data describing Thalictrum flavum NCS bound to a mechanism-inspired ligand. The structure supports two key features of the NCS “dopamine-first” mechanism: the binding of dopamine catechol to Lys-122 and the position of the carbonyl substrate binding site at the active site entrance. The catalytically vital residue Glu-110 occupies a previously unobserved ligand-bound conformation that may be catalytically significant. The potential roles of inhibitory binding and alternative amino acid conformations in the mechanism have also been revealed. This work significantly advances our understanding of the NCS mechanism and will aid future efforts to engineer the substrate scope and catalytic properties of this useful biocatalyst.