Ketosynthase mutants enable short-chain fatty acid biosynthesis in E. coli.

Cells build fatty acids in tightly regulated assembly lines, or fatty acid synthases (FASs), in which ß-ketoacyl-acyl carrier protein (ACP) synthases (KSs) catalyze sequential carbon-carbon bond forming reactions that generate acyl-ACPs of varying lengths-precursors for a diverse set of lipids and oleochemicals. To date, most efforts to control fatty acid synthesis in engineered microbes have focused on modifying termination enzymes such as acyl-ACP thioesterases, which release free fatty acids from acyl-ACPs. Changes to the substrate specificity of KSs provide an alternative-and, perhaps, more generalizable-approach that focuses on controlling the acyl-ACPs available for downstream products. This study combines mutants of FabF and FabB, the two elongating KSs of the E. coli FAS, with in vitro and in vivo analyses to explore the use of KS mutants to control fatty acid synthesis. In vitro, single amino acid substitutions in the gating loop and acyl binding pocket of FabF shifted the product profiles of reconstituted FASs toward short chains and showed that KS mutants, alone, can cause large shifts in average length (i.e., 6.5-13.5). FabB, which is essential for unsaturated fatty acid synthesis, blunted this effect in vivo, but exogenously added cis-vaccenic acid (C18:1) enabled sufficient transcriptional repression of FabB to restore it. Strikingly, a single mutant of FabB afforded titers of octanoic acid as high as those generated by an engineered thioesterase. Findings indicate that fatty acid synthesis must be decoupled from microbial growth to resolve the influence of KS mutants on fatty acid profiles but show that these mutants offer a versatile approach for tuning FAS outputs.

A kinetic rationale for functional redundancy in fatty acid biosynthesis.

Cells build fatty acids with biocatalytic assembly lines in which a subset of enzymes often exhibit overlapping activities (e.g., two enzymes catalyze one or more identical reactions). Although the discrete enzymes that make up fatty acid pathways are well characterized, the importance of catalytic overlap between them is poorly understood. We developed a detailed kinetic model of the fatty acid synthase (FAS) of Escherichia coli and paired that model with a fully reconstituted in vitro system to examine the capabilities afforded by functional redundancy in fatty acid synthesis. The model captures-and helps explain-the effects of experimental perturbations to FAS systems and provides a powerful tool for guiding experimental investigations of fatty acid assembly. Compositional analyses carried out in silico and in vitro indicate that FASs with multiple partially redundant enzymes enable tighter (i.e., more independent and/or broader range) control of distinct biochemical objectives-the total production, unsaturated fraction, and average length of fatty acids-than FASs with only a single multifunctional version of each enzyme (i.e., one enzyme with the catalytic capabilities of two partially redundant enzymes). Maximal production of unsaturated fatty acids, for example, requires a second dehydratase that is not essential for their synthesis. This work provides a kinetic, control-theoretic rationale for the inclusion of partially redundant enzymes in fatty acid pathways and supplies a valuable framework for carrying out detailed studies of FAS kinetics.

Kinetically guided, ratiometric tuning of fatty acid biosynthesis.

Cellular metabolism is a nonlinear reaction network in which dynamic shifts in enzyme concentration help regulate the flux of carbon to different products. Despite the apparent simplicity of these biochemical adjustments, their influence on metabolite biosynthesis tends to be context-dependent, difficult to predict, and challenging to exploit in metabolic engineering. This study combines a detailed kinetic model with a systematic set of in vitro and in vivo analyses to explore the use of enzyme concentration as a control parameter in fatty acid synthesis, an essential metabolic process with important applications in oleochemical production. Compositional analyses of a modeled and experimentally reconstituted fatty acid synthase (FAS) from Escherichia coli indicate that the concentration ratio of two native enzymes-a promiscuous thioesterase and a ketoacyl synthase-can tune the average length of fatty acids, an important design objective of engineered pathways. The influence of this ratio is sensitive to the concentrations of other FAS components, which can narrow or expand the range of accessible chain lengths. Inside the cell, simple changes in enzyme concentration can enhance product-specific titers by as much as 125-fold and elicit shifts in overall product profiles that rival those of thioesterase mutants. This work develops a kinetically guided approach for using ratiometric adjustments in enzyme concentration to control the product profiles of FAS systems and, broadly, provides a detailed framework for understanding how coordinated shifts in enzyme concentration can afford tight control over the outputs of nonlinear metabolic pathways.

A kinetic framework for modeling oleochemical biosynthesis in Escherichia coli.

Microorganisms build fatty acids with biocatalytic assembly lines, or fatty acid synthases (FASs), that can be repurposed to produce a broad set of fuels and chemicals. Despite their versatility, the product profiles of FAS-based pathways are challenging to adjust without experimental iteration, and off-target products are common. This study uses a detailed kinetic model of the Escherichia coli FAS as a foundation to model nine oleochemical pathways. These models provide good fits to experimental data and help explain unexpected results from in vivo studies. An analysis of pathways for alkanes and fatty acid ethyl esters (FAEEs), for example, suggests that reductions in titer caused by enzyme overexpression-an experimentally consistent phenomenon-can result from shifts in metabolite pools that are incompatible with the substrate specificities of downstream enzymes, and a focused examination of multiple alcohol pathways indicates that coordinated shifts in enzyme concentrations provide a general means of tuning the product profiles of pathways with promiscuous components. The study concludes by integrating all models into a graphical user interface. The models supplied by this work provide a versatile kinetic framework for studying oleochemical pathways in different biochemical contexts.

Autocatalytic, bistable, oscillatory networks of biologically relevant organic reactions.

Networks of organic chemical reactions are important in life and probably played a central part in its origin. Network dynamics regulate cell division, circadian rhythms, nerve impulses and chemotaxis, and guide the development of organisms. Although out-of-equilibrium networks of chemical reactions have the potential to display emergent network dynamics such as spontaneous pattern formation, bistability and periodic oscillations, the principles that enable networks of organic reactions to develop complex behaviours are incompletely understood. Here we describe a network of biologically relevant organic reactions (amide formation, thiolate-thioester exchange, thiolate-disulfide interchange and conjugate addition) that displays bistability and oscillations in the concentrations of organic thiols and amides. Oscillations arise from the interaction between three subcomponents of the network: an autocatalytic cycle that generates thiols and amides from thioesters and dialkyl disulfides; a trigger that controls autocatalytic growth; and inhibitory processes that remove activating thiol species that are produced during the autocatalytic cycle. In contrast to previous studies that have demonstrated oscillations and bistability using highly evolved biomolecules (enzymes and DNA) or inorganic molecules of questionable biochemical relevance (for example, those used in Belousov-Zhabotinskii-type reactions), the organic molecules we use are relevant to metabolism and similar to those that might have existed on the early Earth. By using small organic molecules to build a network of organic reactions with autocatalytic, bistable and oscillatory behaviour, we identify principles that explain the ways in which dynamic networks relevant to life could have developed. Modifications of this network will clarify the influence of molecular structure on the dynamics of reaction networks, and may enable the design of biomimetic networks and of synthetic self-regulating and evolving chemical systems.

Optogenetic Analysis of Allosteric Control in Protein Tyrosine Phosphatases.

Allosteric regulation enables dynamic adjustments to protein function that permit tight control over cellular biochemistry. Discrepancies in the allosteric systems of related proteins can thus reveal important differences in their susceptibilities to influential stimuli (e.g., allosteric ligands, mutations, or post-translational modifications). This study uses an optogenetic actuator as a tool to compare the allosteric systems of two structurally related regulatory proteins: protein tyrosine phosphatase 1B (PTP1B) and T-cell protein tyrosine phosphatase (TCPTP). It begins with an interesting observation: The fusion of a protein light switch to the allosterically influential α7 helix of PTP1B permits optical modulation of its catalytic activity, but a similar fusion to TCPTP does not. A subsequent analysis of different PTP chimeras shows that replacing regions of TCPTP with homologous regions from PTP1B can enhance photocontrol; as TCPTP becomes more "PTP1B-like", its photosensitivity increases. Interestingly, the structural changes required for photocontrol also enhance the sensitivity of TCPTP to other allosteric inputs, notably, an allosteric inhibitor and a newly reported activating mutation. Our findings indicate that the allosteric functionality of the α7 helix of PTP1B is not conserved across the PTP family and highlight residues necessary to transfer this functionality to other PTPs. More broadly, our results suggest that simple gene fusion events can strengthen allosteric communication within individual protein domains and describe an intriguing application for optogenetic actuators as structural probes-a sort of physically disruptive "ratchet"-for studying protein allostery.

Evolutionarily Conserved Allosteric Communication in Protein Tyrosine Phosphatases.

Protein tyrosine phosphatases (PTPs) are an important class of regulatory enzymes that exhibit aberrant activities in a wide range of diseases. A detailed mapping of allosteric communication in these enzymes could, thus, reveal the structural basis of physiologically relevant-and, perhaps, therapeutically informative-perturbations (i.e., mutations, post-translational modifications, or binding events) that influence their catalytic states. This study combines detailed biophysical studies of protein tyrosine phosphatase 1B (PTP1B) with bioinformatic analyses of the PTP family to examine allosteric communication in this class of enzymes. Results of X-ray crystallography, molecular dynamics simulations, and sequence-based statistical analyses indicate that PTP1B possesses a broadly distributed allosteric network that is evolutionarily conserved across the PTP family, and findings from both kinetic studies and mutational analyses show that this network is functionally intact in sequence-diverse PTPs. The allosteric network resolved in this study reveals new sites for targeting allosteric inhibitors of PTPs and helps explain the functional influence of a diverse set of disease-associated mutations.

Abietane-Type Diterpenoids Inhibit Protein Tyrosine Phosphatases by Stabilizing an Inactive Enzyme Conformation.

Protein tyrosine phosphatases (PTPs) contribute to a striking variety of human diseases, yet they remain vexingly difficult to inhibit with uncharged, cell-permeable molecules; no inhibitors of PTPs have been approved for clinical use. This study uses a broad set of biophysical analyses to evaluate the use of abietane-type diterpenoids, a biologically active class of phytometabolites with largely nonpolar structures, for the development of pharmaceutically relevant PTP inhibitors. Results of nuclear magnetic resonance analyses, mutational studies, and molecular dynamics simulations indicate that abietic acid can inhibit protein tyrosine phosphatase 1B, a negative regulator of insulin signaling and an elusive drug target, by binding to its active site in a non-substrate-like manner that stabilizes the catalytically essential WPD loop in an inactive conformation; detailed kinetic studies, in turn, show that minor changes in the structures of abietane-type diterpenoids (e.g., the addition of hydrogens) can improve potency (i.e., lower IC50) by 7-fold. These findings elucidate a previously uncharacterized mechanism of diterpenoid-mediated inhibition and suggest, more broadly, that abietane-type diterpenoids are a promising source of structurally diverse-and, intriguingly, microbially synthesizable-molecules on which to base the design of new PTP-inhibiting therapeutics.

Warning signals for eruptive events in spreading fires.

Spreading fires are noisy (and potentially chaotic) systems in which transitions in dynamics are notoriously difficult to predict. As flames move through spatially heterogeneous environments, sudden shifts in temperature, wind, or topography can generate combustion instabilities, or trigger self-stabilizing feedback loops, that dramatically amplify the intensities and rates with which fires propagate. Such transitions are rarely captured by predictive models of fire behavior and, thus, complicate efforts in fire suppression. This paper describes a simple, remarkably instructive physical model for examining the eruption of small flames into intense, rapidly moving flames stabilized by feedback between wind and fire (i.e., "wind-fire coupling"-a mechanism of feedback particularly relevant to forest fires), and it presents evidence that characteristic patterns in the dynamics of spreading flames indicate when such transitions are likely to occur. In this model system, flames propagate along strips of nitrocellulose with one of two possible modes of propagation: a slow, structured mode, and a fast, unstructured mode sustained by wind-fire coupling. Experimental examination of patterns in dynamics that emerge near bifurcation points suggests that symptoms of critical slowing down (i.e., the slowed recovery of the system from perturbations as it approaches tipping points) warn of impending transitions to the unstructured mode. Findings suggest that slowing responses of spreading flames to sudden changes in environment (e.g., wind, terrain, temperature) may anticipate the onset of intense, feedback-stabilized modes of propagation (e.g., "blowup fires" in forests).

Water-Restructuring Mutations Can Reverse the Thermodynamic Signature of Ligand Binding to Human Carbonic Anhydrase.

This study uses mutants of human carbonic anhydrase (HCAII) to examine how changes in the organization of water within a binding pocket can alter the thermodynamics of protein-ligand association. Results from calorimetric, crystallographic, and theoretical analyses suggest that most mutations strengthen networks of water-mediated hydrogen bonds and reduce binding affinity by increasing the enthalpic cost and, to a lesser extent, the entropic benefit of rearranging those networks during binding. The organization of water within a binding pocket can thus determine whether the hydrophobic interactions in which it engages are enthalpy-driven or entropy-driven. Our findings highlight a possible asymmetry in protein-ligand association by suggesting that, within the confines of the binding pocket of HCAII, binding events associated with enthalpically favorable rearrangements of water are stronger than those associated with entropically favorable ones.

Interactions between Hofmeister anions and the binding pocket of a protein.

This paper uses the binding pocket of human carbonic anhydrase II (HCAII, EC 4.2.1.1) as a tool to examine the properties of Hofmeister anions that determine (i) where, and how strongly, they associate with concavities on the surfaces of proteins and (ii) how, upon binding, they alter the structure of water within those concavities. Results from X-ray crystallography and isothermal titration calorimetry show that most anions associate with the binding pocket of HCAII by forming inner-sphere ion pairs with the Zn(2+) cofactor. In these ion pairs, the free energy of anion-Zn(2+) association is inversely proportional to the free energetic cost of anion dehydration; this relationship is consistent with the mechanism of ion pair formation suggested by the "law of matching water affinities". Iodide and bromide anions also associate with a hydrophobic declivity in the wall of the binding pocket. Molecular dynamics simulations suggest that anions, upon associating with Zn(2+), trigger rearrangements of water that extend up to 8 Å away from their surfaces. These findings expand the range of interactions previously thought to occur between ions and proteins by suggesting that (i) weakly hydrated anions can bind complementarily shaped hydrophobic declivities, and that (ii) ion-induced rearrangements of water within protein concavities can (in contrast with similar rearrangements in bulk water) extend well beyond the first hydration shells of the ions that trigger them. This study paints a picture of Hofmeister anions as a set of structurally varied ligands that differ in size, shape, and affinity for water and, thus, in their ability to bind to­and to alter the charge and hydration structure of­polar, nonpolar, and topographically complex concavities on the surfaces of proteins.

A single-molecule analysis reveals morphological targets for cellulase synergy.

The mechanisms of enzyme activity on solid substrates are not well understood. Unlike enzyme catalysis in aqueous solutions, enzyme activity on surfaces is complicated by adsorption steps and structural heterogeneities that make enzyme-substrate interactions difficult to characterize. Cellulase enzymes, which catalyze the depolymerization of cellulose, show binding specificities for different cellulose surface morphologies, but the influence of these specificities on the activity of multienzyme mixtures has remained unclear. We developed a metric to quantify binding-target arrangements determined by photoactivated localization microscopy, and we used that metric to show that combinations of cellulases designed to bind within similar but nonidentical morphologies can have synergistic activity. This phenomenon cannot be explained with the binary crystalline or amorphous classifications commonly used to characterize cellulase-binding targets. Our results reveal a strategy for improving the activity of cellulolytic mixtures and demonstrate a versatile method for investigating protein organization on heterogeneous surfaces.

Analysis of neutral mutational drift in an allosteric enzyme.

Neutral mutational drift is an important source of biological diversity that remains underexploited in fundamental studies of protein biophysics. This study uses a synthetic transcriptional circuit to study neutral drift in protein tyrosine phosphatase 1B (PTP1B), a mammalian signaling enzyme for which conformational changes are rate limiting. Kinetic assays of purified mutants indicate that catalytic activity, rather than thermodynamic stability, guides enrichment under neutral drift, where neutral or mildly activating mutations can mitigate the effects of deleterious ones. In general, mutants show a moderate activity-stability tradeoff, an indication that minor improvements in the activity of PTP1B do not require concomitant losses in its stability. Multiplexed sequencing of large mutant pools suggests that substitutions at allosterically influential sites are purged under biological selection, which enriches for mutations located outside of the active site. Findings indicate that the positional dependence of neutral mutations within drifting populations can reveal the presence of allosteric networks and illustrate an approach for using synthetic transcriptional systems to explore these mutations in regulatory enzymes.

Genetically Encoded Detection of Biosynthetic Protease Inhibitors.

Proteases are an important class of drug targets that continue to drive inhibitor discovery. These enzymes are prone to resistance mutations, yet their promise for treating viral diseases and other disorders continues to grow. This study develops a general approach for detecting microbially synthesized protease inhibitors and uses it to screen terpenoid pathways for inhibitory compounds. The detection scheme relies on a bacterial two-hybrid (B2H) system that links protease inactivation to the transcription of a swappable reporter gene. This system, which can accomodate multiple biochemical outputs (i.e., luminescence and antibiotic resistance), permitted the facile incorporation of four disease-relevant proteases. A B2H designed to detect the inactivation of the main protease of severe acute respiratory syndrome coronavirus 2 enabled the identification of a terpenoid inhibitor of modest potency. An analysis of multiple pathways that make this terpenoid, however, suggested that its production was necessary but not sufficient to confer a survival advantage in growth-coupled assays. This finding highlights an important challenge associated with the use of genetic selection to search for inhibitors─notably, the influence of pathway toxicity─and underlines the value of including multiple pathways with overlapping product profiles in pathway screens. This study provides a detailed experimental framework for using microbes to screen libraries of biosynthetic pathways for targeted protease inhibitors.

Biophysical Rationale for the Selective Inhibition of PTP1B over TCPTP by Nonpolar Terpenoids.

Protein tyrosine phosphatases (PTPs) are emerging drug targets for many diseases, including cancer, autoimmunity, and neurological disorders. A high degree of structural similarity between their catalytic domains, however, has hindered the development of selective pharmacological agents. Our previous research uncovered two unfunctionalized terpenoid inhibitors that selectively inhibit PTP1B over T-cell PTP (TCPTP), two PTPs with high sequence conservation. Here, we use molecular modeling, with supporting experimental validation, to study the molecular basis of this unusual selectivity. Molecular dynamics (MD) simulations suggest that PTP1B and TCPTP share a h-bond network that connects the active site to a distal allosteric pocket; this network stabilizes the closed conformation of the catalytically essential WPD loop, which it links to the L-11 loop and neighboring α3 and α7 helices on the other side of the catalytic domain. Terpenoid binding to either of two proximal C-terminal sites─an α site and a ß site─can disrupt the allosteric network; however, binding to the α site forms a stable complex only in PTP1B. In TCPTP, two charged residues disfavor binding at the α site in favor of binding at the ß site, which is conserved between the two proteins. Our findings thus indicate that minor amino acid differences at the poorly conserved α site enable selective binding, a property that might be enhanced with chemical elaboration, and illustrate more broadly how minor differences in the conservation of neighboring─yet functionally similar─allosteric sites can affect the selectivity of inhibitory scaffolds (e.g., fragments).

A biophysical rationale for the selective inhibition of PTP1B over TCPTP by nonpolar terpenoids.

Protein tyrosine phosphatases (PTPs) are emerging drug targets for many diseases, including type 2 diabetes, obesity, and cancer. However, a high degree of structural similarity between the catalytic domains of these enzymes has made the development of selective pharmacological inhibitors an enormous challenge. Our previous research uncovered two unfunctionalized terpenoid inhibitors that selectively inhibit PTP1B over TCPTP, two PTPs with high sequence conservation. Here, we use molecular modeling with experimental validation to study the molecular basis of this unusual selectivity. Molecular dynamics (MD) simulations indicate that PTP1B and TCPTP contain a conserved h-bond network that connects the active site to a distal allosteric pocket; this network stabilizes the closed conformation of the catalytically influential WPD loop, which it links to the L-11 loop and α 3 and α 7 helices-the C-terminal side of the catalytic domain. Terpenoid binding to either of two proximal allosteric sites-an α site and a ß site-can disrupt the allosteric network. Interestingly, binding to the α site forms a stable complex with only PTP1B; in TCPTP, where two charged residues disfavor binding at the α site, the terpenoids bind to the ß site, which is conserved between the two proteins. Our findings indicate that minor amino acid differences at the poorly conserved α site enable selective binding, a property that might be enhanced with chemical elaboration, and illustrate, more broadly, how minor differences in the conservation of neighboring-yet functionally similar-allosteric sites can have very different implications for inhibitor selectivity.

Initial- and processive-cut products reveal cellobiohydrolase rate limitations and the role of companion enzymes.

Efforts to improve the activity of cellulases, which catalyze the hydrolysis of insoluble cellulose, have been hindered by uncertainty surrounding the mechanistic origins of rate-limiting phenomena and by an incomplete understanding of complementary enzyme function. In particular, direct kinetic measurements of individual steps occurring after enzymes adsorb to the cellulose surface have proven to be experimentally elusive. This work describes an experimental and analytical approach, derived from a detailed mechanistic model of cellobiohydrolase action, for determining rates of initial- and processive-cut product generation by Trichoderma longibrachiatum cellobiohydrolase I (TlCel7A) as it catalyzes the hydrolysis of bacterial microcrystalline cellulose (BMCC) alone and in the presence of Talaromyces emersonii endoglucanase II (TemGH5). This analysis revealed that the rate of TlCel7A-catalyzed hydrolysis of crystalline cellulose is limited by the rate of enzyme complexation with glycan chains, which is shown to be equivalent to the rate of initial-cut product generation. This rate is enhanced in the presence of endoglucanase enzymes. The results confirm recent reports about the role of morphological obstacles in enzyme processivity and also provide the first direct evidence that processive length may be increased by the presence of companion enzymes, including small amounts of TemGH5. The findings of this work indicate that efforts to improve cellobiohydrolase activity should focus on enhancing the enzyme's ability to complex with cellulose chains, and the analysis employed provides a new technique for investigating the mechanism by which companion enzymes influence cellobiohydrolase activity.

Analysis of Three Architectures for Controlling PTP1B with Light.

Photosensory domains are powerful tools for placing proteins under optical control, but their integration into light-sensitive chimeras is often challenging. Many designs require structural iterations, and direct comparisons of alternative approaches are rare. This study uses protein tyrosine phosphatase 1B (PTP1B), an influential regulatory enzyme, to compare three architectures for controlling PTPs with light: a protein fusion, an insertion chimera, and a split construct. All three designs permitted optical control of PTP1B activity in vitro (i.e., kinetic assays of purified enzyme) and in mammalian cells; photoresponses measured under both conditions, while different in magnitude, were linearly correlated. The fusion- and insertion-based architectures exhibited the highest dynamic range and maintained native localization patterns in mammalian cells. A single insertion architecture enabled optical control of both PTP1B and TCPTP, but not SHP2, where the analogous chimera was active but not photoswitchable. Findings suggest that PTPs are highly tolerant of domain insertions and support the use of in vitro screens to evaluate different optogenetic designs.