EnzymeML: seamless data flow and modeling of enzymatic data.

The design of biocatalytic reaction systems is highly complex owing to the dependency of the estimated kinetic parameters on the enzyme, the reaction conditions, and the modeling method. Consequently, reproducibility of enzymatic experiments and reusability of enzymatic data are challenging. We developed the XML-based markup language EnzymeML to enable storage and exchange of enzymatic data such as reaction conditions, the time course of the substrate and the product, kinetic parameters and the kinetic model, thus making enzymatic data findable, accessible, interoperable and reusable (FAIR). The feasibility and usefulness of the EnzymeML toolbox is demonstrated in six scenarios, for which data and metadata of different enzymatic reactions are collected and analyzed. EnzymeML serves as a seamless communication channel between experimental platforms, electronic lab notebooks, tools for modeling of enzyme kinetics, publication platforms and enzymatic reaction databases. EnzymeML is open and transparent, and invites the community to contribute. All documents and codes are freely available at https://enzymeml.org .

Selectivity and kinetic modeling of penicillin G acylase variants for the synthesis of cephalexin under a broad range of substrate concentrations.

The kinetics of cephalexin synthesis and hydrolysis of the activated acyl-donor precursor phenylglycine methyl ester (PGME) were characterized under a broad range of substrate concentrations. A previously developed model by Youshko-Svedas involving the formation of the acyl-enzyme complex followed by binding of the nucleophilic ß-lactam donor does not fully estimate the maximum reaction yields for cephalexin synthesis at different concentrations using initial-rate data. 7-aminodesacetoxycephalosporanic acid (7-ADCA) was discovered to be a potent inhibitor of cephalexin hydrolysis, which may account for the deviation from model predictions. Three kinetic models were compared for cephalexin synthesis, with the model incorporating competitive inhibition due to 7-ADCA yielding the best fit. Additionally, the ßF24A variant and Assemblase® did not exhibit significantly different kinetics for the synthesis of cephalexin compared to the wild-type, for the concentration range evaluated and for both initial-rate experiments and time-course synthesis experiments. Lastly, a continuous stirred-tank reactor for cephalexin synthesis was simulated using the model incorporating competitive inhibition by 7-ADCA, with clear tradeoffs observed between productivity, fractional yield, and PGME conversion.

Protein-inorganic calcium-phosphate supraparticles as a robust platform for enzyme co-immobilization.

Immobilization of enzymes provides many benefits, including facile separation and recovery of enzymes from reaction mixtures, enhanced stability, and co-localization of multiple enzymes. Calcium-phosphate-protein supraparticles imbued with a leucine zipper binding domain (ZR ) serve as a modular immobilization platform for enzymes fused to the complementary leucine zipper domain (ZE ). The zippers provide high-affinity, specific binding, separating enzymatic activity from the binding event. Using fluorescent model proteins (mCherryZE and eGFPZE ), an amine dehydrogenase (AmDHZE ), and a formate dehydrogenase (FDHZE ), the efficacy of supraparticles as a biocatalytic solid support was assessed. Supraparticles demonstrated several benefits as an immobilization support, including predictable loading of multiple proteins, structural integrity in a panel of solvents, and the ability to elute and reload proteins without damaging the support. The dual-enzyme reaction successfully converted ketone to amine on supraparticles, highlighting the efficacy of this system.

Tuning the Morphology of Protein-Inorganic Calcium-Phosphate Supraparticles via Directed Assembly.

Understanding the phenomena that govern complex interfacial and directed assemblies is essential for both control and scale-up of particle syntheses. The present work describes an effort to understand, control, and tune the formation of protein-inorganic calcium-phosphate supraparticles that are produced at an oscillating air-water interface created by end-over-end rotation of the synthesis solution. Supraparticles were synthesized under an array of different conditions that varied reagent concentration, the presence of additives, tube size, and rotational speed. Paired with a fluid mechanics model of the end-over-end rotation and dimensional analysis, the sensitivity of the synthesis to physicochemical and mechanical parameters was determined. Surface tension and bubble formation were found to be important criteria for changing the size distribution of supraparticles. Thresholds for the values of the Froude, Iribarren, and rotational Reynolds numbers were identified for narrowing particle size distribution. These results both guide the specific protein-inorganic supraparticle synthesis described here and inform future manipulation and scale-up of other complex interfacial colloidal assemblies.

Biomolecular Assemblies: Moving from Observation to Predictive Design.

Biomolecular assembly is a key driving force in nearly all life processes, providing structure, information storage, and communication within cells and at the whole organism level. These assembly processes rely on precise interactions between functional groups on nucleic acids, proteins, carbohydrates, and small molecules, and can be fine-tuned to span a range of time, length, and complexity scales. Recognizing the power of these motifs, researchers have sought to emulate and engineer biomolecular assemblies in the laboratory, with goals ranging from modulating cellular function to the creation of new polymeric materials. In most cases, engineering efforts are inspired or informed by understanding the structure and properties of naturally occurring assemblies, which has in turn fueled the development of predictive models that enable computational design of novel assemblies. This Review will focus on selected examples of protein assemblies, highlighting the story arc from initial discovery of an assembly, through initial engineering attempts, toward the ultimate goal of predictive design. The aim of this Review is to highlight areas where significant progress has been made, as well as to outline remaining challenges, as solving these challenges will be the key that unlocks the full power of biomolecules for advances in technology and medicine.

The Short-chain Dehydrogenase/Reductase Engineering Database (SDRED): A classification and analysis system for a highly diverse enzyme family.

The Short-chain Dehydrogenases/Reductases Engineering Database (SDRED) covers one of the largest known protein families (168 150 proteins). Assignment to the superfamilies of Classical and Extended SDRs was achieved by global sequence similarity and by identification of family-specific sequence motifs. Two standard numbering schemes were established for Classical and Extended SDRs that allow for the determination of conserved amino acid residues, such as cofactor specificity determining positions or superfamily specific sequence motifs. The comprehensive sequence dataset of the SDRED facilitates the refinement of family-specific sequence motifs. The glycine-rich motifs for Classical and Extended SDRs were refined to improve the precision of superfamily classification. In each superfamily, the majority of sequences formed a tightly connected sequence network and belonged to a large homologous family. Despite their different sequence motifs and their different sequence length, the two sequence networks of Classical and Extended SDRs are not separate, but connected by edges at a threshold of 40% sequence similarity, indicating that all SDRs belong to a large, connected network. The SDRED is accessible at https://sdred.biocatnet.de/.

Photoirradiation Generates an Ultrastable 8-Formyl FAD Semiquinone Radical with Unusual Properties in Formate Oxidase.

Formate oxidase (FOX) was previously shown to contain a noncovalently bound 8-formyl FAD (8-fFAD) cofactor. However, both the absorption spectra and the kinetic parameters previously reported for FOX are inconsistent with more recent reports. The ultraviolet-visible (UV-vis) absorption spectrum reported in early studies closely resembles the spectra observed for protein-bound 8-formyl flavin semiquinone species, thus suggesting FOX may be photosensitive. Therefore, the properties of dark and light-exposed FOX were investigated using steady-state kinetics and site-directed mutagenesis analysis along with inductively coupled plasma optical emission spectroscopy, UV-vis absorption spectroscopy, circular dichroism spectroscopy, liquid chromatography and mass spectrometry, and electron paramagnetic resonance (EPR) spectroscopy. Surprisingly, these experimental results demonstrate that FOX is deactivated in the presence of light through generation of an oxygen stable, anionic (red) 8-fFAD semiquinone radical capable of persisting either in an aerobic environment for multiple weeks or in the presence of a strong reducing agent like sodium dithionite. Herein, we study the photoinduced formation of the 8-fFAD semiquinone radical in FOX and report the first EPR spectrum of this radical species. The stability of the 8-fFAD semiquinone radical suggests FOX to be a model enzyme for probing the structural and mechanistic features involved in stabilizing flavin semiquinone radicals. It is likely that the photoinduced formation of a stable 8-fFAD semiquinone radical is a defining characteristic of 8-formyl flavin-dependent enzymes. Additionally, a better understanding of the radical stabilization process may yield a FOX enzyme with more robust activity and broader industrial usefulness.

Mechanistic studies of formate oxidase from Aspergillus oryzae: A novel member of the glucose-Methanol-choline oxidoreductase enzyme superfamily that oxidizes carbon acids.

Formate oxidase (FOX) from Aspergillus oryzae is the only GMC member that oxidizes a carbon acid rather than alcohols; thus, its catalytic mechanism may be different from that of other GMC members. We have used pH, solvent viscosity, and deuterium kinetic isotope effects, to investigate the catalytic mechanism of FOX. The enzyme followed a Bi-Bi sequential steady-state kinetic mechanism. The kcat value was pH-independent between pH 2.8 and 6.8, suggesting a lack of ionizable groups in kinetic step(s) that limit the overall turnover of the enzyme. The kcat/Kformate value decreased from a value of 10,000 M-1s-1 at low pH with a pKa value of 4.4, consistent with the requirement of a protonated group for substrate binding. An inverse viscosity dependence on the kcat/Kformate value indicated an isomerization of the Michaelis complex. The kcat/Koxygen value was 340,000 M-1s-1 and pH independent up to pH 6.0. The Dkcat and D(kcat/Kformate) values were 2.5 and 1.9, respectively, indicating that substrate CH bond cleavage is rate-limiting for FOX catalysis. Analytical ultracentrifugation indicated a concentration dependence of the oligomeric state of FOX. The appkred,H value was ∼75% that of kcat,H, indicating that the anaerobic reduction of FOX was dependent on the oligomeric state of FOX.

A high-throughput pH-based colorimetric assay: application focus on alpha/beta hydrolases.

Research involving α/ß hydrolases, including α-amino acid ester hydrolase and cocaine esterase, has been limited by the lack of an online high throughput screening assay. The development of a high throughput screening assay capable of detecting α/ß hydrolase activity toward specific substrates and/or chemical reactions (e.g., hydrolysis in lieu of amidase activity and/or synthesis instead of thioesterase activity) is of interest in a broad set of scientific questions and applications. Here we present a general framework for pH-based colorimetric assays, as well as the mathematical considerations necessary to estimate de novo the experimental response required to assign a 'hit' or a 'miss,' in the absence of experimental standard curves. This combination is valuable for screening the hydrolysis and synthesis activity of α/ß hydrolases on a variety of substrates, and produces data comparable to the current standard technique involving High Performance Liquid Chromatography (HPLC). In contrast to HPLC, this assay enables screening experiments to be performed with greater efficiency.

Informing Efforts to Develop Nitroreductase for Amine Production.

Nitroreductases (NRs) hold promise for converting nitroaromatics to aromatic amines. Nitroaromatic reduction rate increases with Hammett substituent constant for NRs from two different subgroups, confirming substrate identity as a key determinant of reactivity. Amine yields were low, but compounds yielding amines tend to have a large π system and electron withdrawing substituents. Therefore, we also assessed the prospects of varying the enzyme. Several different subgroups of NRs include members able to produce aromatic amines. Comparison of four NR subgroups shows that they provide contrasting substrate binding cavities with distinct constraints on substrate position relative to the flavin. The unique architecture of the NR dimer produces an enormous contact area which we propose provides the stabilization needed to offset the costs of insertion of the active sites between the monomers. Thus, we propose that the functional diversity included in the NR superfamily stems from the chemical versatility of the flavin cofactor in conjunction with a structure that permits tremendous active site variability. These complementary properties make NRs exceptionally promising enzymes for development for biocatalysis in prodrug activation and conversion of nitroaromatics to valuable aromatic amines. We provide a framework for identifying NRs and substrates with the greatest potential to advance.

Enzyme-Mediated Conversion of Flavin Adenine Dinucleotide (FAD) to 8-Formyl FAD in Formate Oxidase Results in a Modified Cofactor with Enhanced Catalytic Properties.

Flavins, including flavin adenine dinucleotide (FAD), are fundamental catalytic cofactors that are responsible for the redox functionality of a diverse set of proteins. Alternatively, modified flavin analogues are rarely found in nature as their incorporation typically results in inactivation of flavoproteins, thus leading to the disruption of important cellular pathways. Here, we report that the fungal flavoenzyme formate oxidase (FOX) catalyzes the slow conversion of noncovalently bound FAD to 8-formyl FAD and that this conversion results in a nearly 10-fold increase in formate oxidase activity. Although the presence of an enzyme-bound 8-formyl FMN has been reported previously as a result of site-directed mutagenesis studies of lactate oxidase, FOX is the first reported case of 8-formyl FAD in a wild-type enzyme. Therefore, the formation of the 8-formyl FAD cofactor in formate oxidase was investigated using steady-state kinetics, site-directed mutagenesis, ultraviolet-visible, circular dichroism, and fluorescence spectroscopy, liquid chromatography with mass spectrometry, and computational analysis. Surprisingly, the results from these studies indicate not only that 8-formyl FAD forms spontaneously and results in the active form of FOX but also that its autocatalytic formation is dependent on a nearby arginine residue, R87. Thus, this work describes a new enzyme cofactor and provides insight into the little-understood mechanism of enzyme-mediated 8α-flavin modifications.

Contributions of the Prion Protein Sequence, Strain, and Environment to the Species Barrier.

Amyloid propagation requires high levels of sequence specificity so that only molecules with very high sequence identity can form cross-ß-sheet structures of sufficient stringency for incorporation into the amyloid fibril. This sequence specificity presents a barrier to the transmission of prions between two species with divergent sequences, termed a species barrier. Here we study the relative effects of protein sequence, seed conformation, and environment on the species barrier strength and specificity for the yeast prion protein Sup35p from three closely related species of the Saccharomyces sensu stricto group; namely, Saccharomyces cerevisiae, Saccharomyces bayanus, and Saccharomyces paradoxus. Through in vivo plasmid shuffle experiments, we show that the major characteristics of the transmission barrier and conformational fidelity are determined by the protein sequence rather than by the cellular environment. In vitro data confirm that the kinetics and structural preferences of aggregation of the S. paradoxus and S. bayanus proteins are influenced by anions in accordance with their positions in the Hofmeister series, as observed previously for S. cerevisiae. However, the specificity of the species barrier is primarily affected by the sequence and the type of anion present during the formation of the initial seed, whereas anions present during the seeded aggregation process typically influence kinetics rather than the specificity of prion conversion. Therefore, our work shows that the protein sequence and the conformation variant (strain) of the prion seed are the primary determinants of cross-species prion specificity both in vivo and in vitro.

Green fluorescent protein with anionic tryptophan-based chromophore and long fluorescence lifetime.

Spectral diversity of fluorescent proteins, crucial for multiparameter imaging, is based mainly on chemical diversity of their chromophores. Recently we have reported, to our knowledge, a new green fluorescent protein WasCFP-the first fluorescent protein with a tryptophan-based chromophore in the anionic state. However, only a small portion of WasCFP molecules exists in the anionic state at physiological conditions. In this study we report on an improved variant of WasCFP, named NowGFP, with the anionic form dominating at 37°C and neutral pH. It is 30% brighter than enhanced green fluorescent protein (EGFP) and exhibits a fluorescence lifetime of 5.1 ns. We demonstrated that signals of NowGFP and EGFP can be clearly distinguished by fluorescence lifetime in various models, including mammalian cells, mouse tumor xenograft, and Drosophila larvae. NowGFP thus provides an additional channel for multiparameter fluorescence lifetime imaging microscopy of green fluorescent proteins.

Imidazolium Catalysts Formed by an Iterative Synthetic Process as a Model System for Chemical Evolution.

Processes exhibiting diversity and selection would have been necessary to promote chemical evolution on early Earth. In this work, a model process was developed using non-kinetic selection to synthesize and isolate small molecule imidazolium catalysts. These catalysts were purified by affinity chromatography and recycled back into the process, forming a product feedback loop. In dimethylformamide, the catalysts activated the coupling of formaldehyde to short chain sugars. This sugar mixture was reacted with aniline, acetic acid, and paraformaldehyde to generate new catalysts. Thus chemical diversity was produced through non-selective, multi-component synthesis. Applying sequential dilution-reaction-purification cycles it was demonstrated that this process can function independently of starting catalyst. Over three process cycles, the initiator catalyst is effectively diluted out as a new catalyst population emerges to take its place. This system offers an alternative viewpoint for chemical evolution via the generation of small molecule organocatalysts.

Engineering towards nitroreductase functionality in ene-reductase scaffolds.

Nitroreductases (NRs) and ene-reductases (ERs) both utilize flavin mononucleotide cofactors but catalyze distinct reactions. NRs reduce nitroaromatics, whereas ERs reduce unsaturated C=C double bonds, and these functionalities are known to somewhat overlap. Recent studies on the ER xenobiotic reductase A (XenA) from Pseudomonas putida demonstrated the possibility of increasing NR activity with active site modifications. Structural comparison between NRs and ERs led us to hypothesize that active site cavity size plays an important role in determining enzyme functionality. Residues of ER KYE1 from Kluyveromyces lactis were selected to increase the binding pocket size, compensate for hydrogen bonding pattern changes, and eliminate ER activity. Single variants were screened, and promising mutations were combined. Variant F296A/Y275A showed a 100-fold improvement in NR specific activity over wild-type, and variant H191A/F296A/Y375A exhibited complete conversion to a NR.

Hidden photoinduced reactivity of the blue fluorescent protein mKalama1.

Understanding the photoinduced dynamics of fluorescent proteins is essential for their applications in bioimaging. Despite numerous studies on the ultrafast dynamics, the delayed response of these proteins, which often results in population of kinetically trapped dark states of various origins, is largely unexplored. Here, by using transient absorption spectroscopy spanning the time scale from picoseconds to seconds, we reveal a hidden reactivity of the bright blue-light emitting protein mKalama1 previously thought to be inert. This protein shows no excited-state proton transfer during its nanosecond excited-state lifetime; however, its tyrosine-based chromophore undergoes deprotonation coupled to non-radiative electronic relaxation. Such deprotonation causes distinct optical absorption changes in the broad UV-to-NIR spectral range (ca. 300-800 nm); the disappearance of the transient absorption signal has a complex nature and spans the whole microsecond-to-second time scale. The mechanisms underlying the relaxation kinetics are disclosed based on the X-ray structural analysis of mKalama1 and the high-level electronic structure calculations of proposed intermediates in the photocycle. We conclude that the non-radiative excited-state decay includes two major branches: internal conversion coupled to intraprotein proton transfer, where a conserved residue E222 serves as the proton acceptor; and ionization induced by two consecutive resonant absorption events, followed by deprotonation of the chromophore radical cation to bulk solvent through a novel water-mediated proton-wire pathway. Our findings open up new perspectives on the dynamics of fluorescent proteins as tracked by its optical transient absorption in the time domain extending up to seconds.

Ion-specific effects on prion nucleation and strain formation.

Ordered, fibrous, self-seeding aggregates of misfolded proteins known as amyloids are associated with important diseases in mammals and control phenotypic traits in fungi. A given protein may adopt multiple amyloid conformations, known as variants or strains, each of which leads to a distinct disease pattern or phenotype. Here, we study the effect of Hofmeister ions on amyloid nucleation and strain generation by the prion domain-containing fragment (Sup35NM) of a yeast protein Sup35p. Strongly hydrated anions (kosmotropes) initiate nucleation quickly and cause rapid fiber elongation, whereas poorly hydrated anions (chaotropes) delay nucleation and mildly affect the elongation rate. For the first time, we demonstrate that kosmotropes favor formation of amyloid strains that are characterized by lower thermostability and higher frangibility in vitro and stronger phenotypic and proliferation patterns effectively in vivo as compared with amyloids formed in chaotropes. These phenomena point to inherent differences in the biochemistry of Hofmeister ions. Our work shows that the ionic composition of a solution not only influences the kinetics of amyloid nucleation but also determines the amyloid strain that is preferentially formed.

In situ mixed donor synthesis of ampicillin with ethylene glycol co-solvent.

Ethylene glycol has the ability to act as a reactant when employed as a co-solvent for the enzymatic synthesis of ampicillin from (R)-phenylglycine methyl ester and 6-aminopenicillanic acid. The side reaction positively affects yield because its product, (R)-phenylglycine hydroxyethyl ester, is an intermediate for ampicillin synthesis.

Stabilizing biocatalysts.

The area of biocatalysis itself is in rapid development, fueled by both an enhanced repertoire of protein engineering tools and an increasing list of solved problems. Biocatalysts, however, are delicate materials that hover close to the thermodynamic limit of stability. In many cases, they need to be stabilized to survive a range of challenges regarding temperature, pH value, salt type and concentration, co-solvents, as well as shear and surface forces. Biocatalysts may be delicate proteins, however, once stabilized, they are efficiently active enzymes. Kinetic stability must be achieved to a level satisfactory for large-scale process application. Kinetic stability evokes resistance to degradation and maintained or increased catalytic efficiency of the enzyme in which the desired reaction is accomplished at an increased rate. However, beyond these limitations, stable biocatalysts can be operated at higher temperatures or co-solvent concentrations, with ensuing reduction in microbial contamination, better solubility, as well as in many cases more favorable equilibrium, and can serve as more effective templates for combinatorial and data-driven protein engineering. To increase thermodynamic and kinetic stability, immobilization, protein engineering, and medium engineering of biocatalysts are available, the main focus of this work. In the case of protein engineering, there are three main approaches to enhancing the stability of protein biocatalysts: (i) rational design, based on knowledge of the 3D-structure and the catalytic mechanism, (ii) combinatorial design, requiring a protocol to generate diversity at the genetic level, a large, often high throughput, screening capacity to distinguish 'hits' from 'misses', and (iii) data-driven design, fueled by the increased availability of nucleotide and amino acid sequences of equivalent functionality.

Optically modulatable blue fluorescent proteins.

Blue fluorescent proteins (BFPs) offer visualization of protein location and behavior, but often suffer from high autofluorescent background and poor signal discrimination. Through dual-laser excitation of bright and photoinduced dark states, mutations to the residues surrounding the BFP chromophore enable long-wavelength optical modulation of BFP emission. Such dark state engineering enables violet-excited blue emission to be increased upon lower energy, green coillumination. Turning this green coillumination on and off at a specific frequency dynamically modulates collected blue fluorescence without generating additional background. Interpreted as transient photoconversion between neutral cis and anionic trans chromophoric forms, mutations tune photoisomerization and ground state tautomerizations to enable long-wavelength depopulation of the millisecond-lived, spectrally shifted dark states. Single mutations to the tyrosine-based blue fluorescent protein T203V/S205V exhibit enhanced modulation depth and varied frequency. Importantly, analogous single point mutations in the nonmodulatable BFP, mKalama1, creates a modulatable variant. Building modulatable BFPs offers opportunities for improved BFP signal discrimination vs background, greatly enhancing their utility.