Search

Biocatalysis is becoming a data science. High-throughput experimentation generates a rapidly increasing stream of biocatalytic data, which is the raw material for mechanistic and novel data-driven modeling approaches for the predictive design of improved biocatalysts and novel bioprocesses. The holistic and molecular understanding of enzymatic reaction systems will enable us to identify and overcome kinetic bottlenecks and shift the thermodynamics of a reaction. The full characterization and modeling of reaction systems is a community effort; therefore, published methods and results should be findable, accessible, interoperable, and reusable (FAIR), which is achieved by developing standardized data exchange formats, by a complete and reproducible documentation of experimentation, by collaborative platforms for developing sustainable software and for analyzing data, and by repositories for publishing results together with raw data. The FAIRification of biocatalysis is a prerequisite to developing highly automated laboratory infrastructures that improve the reproducibility of scientific results and reduce the time and costs required to develop novel synthesis routes.

Carbohydrate-active enzymes (CAZymes) can be found in all domains of life and play a crucial role in metabolic and physiological processes. CAZymes often possess a modular structure, comprising not only catalytic domains but also associated domains such as carbohydrate-binding modules (CBMs) and linker domains. By exploring the modular diversity of CAZy families, catalysts with novel properties can be discovered and further insight in their biological functions and evolutionary relationships can be obtained. Here we present the carbohydrate-active enzyme domain analysis tool (CANDy), an assembly of several novel scripts, tools and databases that allows users to analyze the domain architecture of all protein sequences in a given CAZy family. CANDy's usability is shown on glycoside hydrolase family 48, a small yet underexplored family containing multi-domain enzymes. Our analysis reveals the existence of 35 distinct domain assemblies, including eight known architectures, with the remaining assemblies awaiting characterization. Moreover, we substantiate the occurrence of horizontal gene transfer from prokaryotes to insect orthologs and provide evidence for the subsequent removal of auxiliary domains, likely through a gene fission event. CANDy is available at https://github.com/PyEED/CANDy., Competing Interests: The authors declare no potential conflict of interest., (Copyright: © 2024 Windels et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Protein lysine methyltransferases (PKMTs) transfer up to three methyl groups to the side chains of lysine residues in proteins and fulfill important regulatory functions by controlling protein stability, localization and protein/protein interactions. The methylation reactions are highly regulated, and aberrant methylation of proteins is associated with several types of diseases including neurologic disorders, cardiovascular diseases, and various types of cancer. This review describes novel insights into the catalytic machinery of various PKMTs achieved by the combined application of biochemical experiments and simulation approaches during the last years, focusing on clinically relevant and well-studied enzymes of this group like DOT1L, SMYD1-3, SET7/9, G9a/GLP, SETD2, SUV420H2, NSD1/2, different MLLs and EZH2. Biochemical experiments have unraveled many mechanistic features of PKMTs concerning their substrate and product specificity, processivity and the effects of somatic mutations observed in PKMTs in cancer cells. Structural data additionally provided information about the substrate recognition, enzyme-substrate complex formation, and allowed for simulations of the substrate peptide interaction and mechanism of PKMTs with atomistic resolution by molecular dynamics and hybrid quantum mechanics/molecular mechanics methods. These simulation technologies uncovered important mechanistic details of the PKMT reaction mechanism including the processes responsible for the deprotonation of the target lysine residue, essential conformational changes of the PKMT upon substrate binding, but also rationalized regulatory principles like PKMT autoinhibition. Further developments are discussed that could bring us closer to a mechanistic understanding of catalysis of this important class of enzymes in the near future. The results described here illustrate the power of the investigation of enzyme mechanisms by the combined application of biochemical experiments and simulation technologies.

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 ., (© 2023. The Author(s), under exclusive licence to Springer Nature America, Inc.)

The Systems Biology community has taken numerous actions to develop data and modeling standards towards FAIR data and model handling. Nevertheless, the debate about incentives and rewards for individual researchers to make their results reproducible is ongoing. Here, we pose the specific question of whether reproducible models have a higher impact in terms of citations. Therefore, we statistically analyze 328 published models recently classified by Tiwari et al. based on their reproducibility. For hypothesis testing, we use a flexible Bayesian approach that provides complete distributional information for all quantities of interest and can handle outliers. The results show that in the period from 2013, i.e., 10 years after the introduction of SBML, to 2020, the group of reproducible models is significantly more cited than the non-reproducible group. We show that differences in journal impact factors do not explain this effect and that this effect increases with additional standardization of data and error model integration via PEtab. Overall, our statistical analysis demonstrates the long-term merits of reproducible modeling for the individual researcher in terms of citations. Moreover, it provides evidence for the increased use of reproducible models in the scientific community., (© 2023. The Author(s).)

Petroleum-based plastics are durable and accumulate in all ecological niches. Knowledge on enzymatic degradation is sparse. Today, less than 50 verified plastics-active enzymes are known. First examples of enzymes acting on the polymers polyethylene terephthalate (PET) and polyurethane (PUR) have been reported together with a detailed biochemical and structural description. Furthermore, very few polyamide (PA) oligomer active enzymes are known. In this article, the current known enzymes acting on the synthetic polymers PET and PUR are briefly summarized, their published activity data were collected and integrated into a comprehensive open access database. The Plastics-Active Enzymes Database (PAZy) represents an inventory of known and experimentally verified enzymes that act on synthetic fossil fuel-based polymers. Almost 3000 homologs of PET-active enzymes were identified by profile hidden Markov models. Over 2000 homologs of PUR-active enzymes were identified by BLAST. Based on multiple sequence alignments, conservation analysis identified the most conserved amino acids, and sequence motifs for PET- and PUR-active enzymes were derived., (© 2022 The Authors. Proteins: Structure, Function, and Bioinformatics published by Wiley Periodicals LLC.)

Lipases are widely used enzymes that catalyze hydrolysis and alcoholysis of fatty acid esters. At high concentrations of small alcohols such as methanol or ethanol, many lipases are inhibited by the substrate. The molecular basis of the inhibition of Candida antarctica lipase B (CALB) by methanol was investigated by unbiased molecular dynamics (MD) simulations, and the substrate binding kinetics was analyzed by Markov state models (MSMs). The modeled fluxes of productive methanol binding at concentrations between 50 mM and 5.5 M were in good agreement with the experimental activity profile of CALB, with a peak at 300 mM. The kinetic and structural analysis uncovered the molecular basis of CALB inhibition. Beyond 300 mM, the kinetic bottleneck results from crowding of methanol in the substrate access channel, which is caused by the gradual formation of methanol patches close to Leu140 (helix α5), Leu278, and Ile285 (helix α10) at a distance of 4-5 Å from the active site. Our findings demonstrate the usefulness of unbiased MD simulations to study enzyme-substrate interactions at realistic substrate concentrations and the feasibility of scale-bridging by an MSM analysis to derive kinetic information.

Expansins have the remarkable ability to loosen plant cell walls and cellulose material without showing catalytic activity and therefore have potential applications in biomass degradation. To support the study of sequence-structure-function relationships and the search for novel expansins, the Expansin Engineering Database (ExED, https://exed.biocatnet.de) collected sequence and structure data on expansins from Bacteria, Fungi, and Viridiplantae, and expansin-like homologues such as carbohydrate binding modules, glycoside hydrolases, loosenins, swollenins, cerato-platanins, and EXPNs. Based on global sequence alignment and protein sequence network analysis, the sequences are highly diverse. However, many similarities were found between the expansin domains. Newly created profile hidden Markov models of the two expansin domains enable standard numbering schemes, comprehensive conservation analyses, and genome annotation. Conserved key amino acids in the expansin domains were identified, a refined classification of expansins and carbohydrate binding modules was proposed, and new sequence motifs facilitate the search of novel candidate genes and the engineering of expansins., (© 2020 The Authors. Proteins: Structure, Function, and Bioinformatics published by Wiley Periodicals LLC.)

The β-hydroxyacid dehydrogenase from Thermocrinus albus (Ta-βHAD), which catalyzes the NADP + -dependent oxidation of β-hydroxyacids, was engineered to accept imines as substrates. The catalytic activity of the proton-donor variant K189D was further increased by the introduction of two nonpolar flanking residues (N192 L, N193 L). Engineering the putative alternative proton donor (D258S) and the gate-keeping residue (F250 A) led to a switched substrate specificity as compared to the single and triple variants. The two most active Ta-βHAD variants were applied to biocatalytic asymmetric reductions of imines at elevated temperatures and enabled enhanced product formation at a reaction temperature of 50 °C., (© 2020 The Authors. Published by Wiley-VCH GmbH.)

Multicopper oxidases (MCOs) use copper ions as cofactors to oxidize a variety of substrates while reducing oxygen to water. MCOs have been identified in various taxa, with notable occurrences in fungi. The role of these fungal MCOs in lignin degradation sparked an interest due to their potential for application in biofuel production and various other industries. MCOs consist of different protein domains, which led to their classification into two-, three-, and six-domain MCOs. The previously established Laccase and Multicopper Oxidase Engineering Database (https://lcced.biocatnet.de) was updated and now includes 51 058 sequences and 229 structures of MCOs. Sequences and structures of all MCOs were systematically compared. All MCOs consist of cupredoxin-like domains. Two-domain MCOs are formed by the N- and C-terminal domain (domain N and C), while three-domain MCOs have an additional domain (M) in between, homologous to domain C. The six-domain MCOs consist of alternating domains N and C, each three times. Two standard numbering schemes were developed for the copper-binding domains N and C, which facilitated the identification of conserved positions and a comparison to previously reported results from mutagenesis studies. Two sequence motifs for the copper binding sites were identified per domain. Their modularity, depending on the placement of the T1-copper binding site, was demonstrated. Protein sequence networks showed relationships between two- and three-domain MCOs, allowing for family-specific annotation and inference of evolutionary relationships., (© 2020 The Authors. Proteins: Structure, Function, and Bioinformatics published by Wiley Periodicals, Inc.)

The family of NAD(P)H-dependent short-chain dehydrogenases/reductases (SDRs) comprises numerous biocatalysts capable of C=O or C=C reduction. The highly homologous noroxomaritidine reductase (NR) from Narcissus sp. aff. pseudonarcissus and Zt_SDR from Zephyranthes treatiae, however, are SDRs with an extended imine substrate scope. Comparison with a similar SDR from Asparagus officinalis (Ao_SDR) exhibiting keto-reducing activity, yet negligible imine-reducing capability, and mining the Short-Chain Dehydrogenase/Reductase Engineering Database indicated that NR and Zt_SDR possess a unique active-site composition among SDRs. Adapting the active site of Ao_SDR accordingly improved its imine-reducing capability. By applying the same strategy, an unrelated SDR from Methylobacterium sp. 77 (M77_SDR) with distinct keto-reducing activity was engineered into a promiscuous enzyme with imine-reducing activity, thereby confirming that the ability to reduce imines can be rationally introduced into members of the "classical" SDR enzyme family. Thus, members of the SDR family could be a promising starting point for protein approaches to generate new imine-reducing enzymes., (© 2020 The Authors. Published by Wiley-VCH Verlag GmbH & Co. KGaA.)

The enzymatic, asymmetric reduction of imines is catalyzed by imine reductases (IREDs), members of the short-chain dehydrogenase/reductase (SDR) family, and β-hydroxy acid dehydrogenase (βHAD) variants. Systematic evaluation of the structures and substrate-binding sites of the three enzyme families has revealed four common principles for imine reduction: structurally conserved cofactor-binding domains; tyrosine, aspartate, or glutamate as proton donor; at least four characteristic flanking residues that adapt the donor's pK a and polarize the substrate; and a negative electrostatic potential in the substrate-binding site to stabilize the transition state. As additional catalytically relevant positions, we propose alternative proton donors in IREDs and βHADs as well as proton relays in IREDs, βHADs, and SDRs. The functional role of flanking residues was experimentally confirmed by alanine scanning of the imine-reducing SDR from Zephyranthes treatiae. Mutating the "gatekeeping" phenylalanine at standard position 200 resulted in a tenfold increase in imine-reducing activity., (© 2020 The Authors. Published by Wiley-VCH Verlag GmbH & Co. KGaA.)

Calculating free energies of binding (Δ G bind ) between ligands and their target protein is of major interest to drug discovery and safety, yet it is still associated with several challenges and difficulties. Linear interaction energy (LIE) is an efficient in silico method for Δ G bind computation. LIE models can be trained and used to directly calculate binding affinities from interaction energies involving ligands in the bound and unbound states only, and LIE can be combined with statistical weighting to calculate Δ G bind for flexible proteins that may bind their ligands in multiple orientations. Here, we investigate if LIE predictions can be effectively improved by explicitly including the entropy of (de)solvation into our free-energy calculations. For that purpose, we combine LIE calculations for the protein-ligand-bound state with explicit free-energy perturbation to rigorously compute the unbound ligand's solvation free energy. We show that for 28 Cytochrome P450 2A6 (CYP2A6) ligands, coupling LIE with alchemical solvation free-energy calculation helps to improve obtained correlation between computed and reference (experimental) binding data.

Thiamine diphosphate-dependent decarboxylases catalyze both cleavage and formation of CC bonds in various reactions, which have been assigned to different homologous sequence families. This work compares 53 ThDP-dependent decarboxylases with known crystal structures. Both sequence and structural information were analyzed synergistically and data were analyzed for global and local properties by means of statistical approaches (principle component analysis and principal coordinate analysis) enabling complexity reduction. The different results obtained both locally and globally, that is, individual positions compared with the overall protein sequence or structure, revealed challenges in the assignment of separated homologous families. The methods applied herein support the comparison of enzyme families and the identification of functionally relevant positions. The findings for the family of ThDP-dependent decarboxylases underline that global sequence identity alone is not sufficient to distinguish enzyme function. Instead, local sequence similarity, defined by comparisons of structurally equivalent positions, allows for a better navigation within several groups of homologous enzymes. The differentiation between homologous sequences is further enhanced by taking structural information into account, such as BioGPS analysis of the active site properties or pairwise structural superimpositions. The methods applied herein are expected to be transferrable to other enzyme families, to facilitate family assignments for homologous protein sequences., (© 2019 Wiley Periodicals, Inc.)

Chorismate and isochorismate constitute branch-point intermediates in the biosynthesis of many aromatic metabolites in microorganisms and plants. To obtain unnatural compounds, we modified the route to menaquinone in Escherichia coli. We propose a model for the binding of isochorismate to the active site of MenD ((1R,2S, 5S,6S)-2-succinyl-5-enolpyruvyl-6-hydroxycyclohex-3-ene-1-carboxylate (SEPHCHC) synthase) that explains the outcome of the native reaction with α-ketoglutarate. We have rationally designed variants of MenD for the conversion of several isochorismate analogues. The double-variant Asn117Arg-Leu478Thr preferentially converts (5S,6S)-5,6-dihydroxycyclohexa-1,3-diene-1-carboxylate (2,3-trans-CHD), the hydrolysis product of isochorismate, with a >70-fold higher ratio than that for the wild type. The single-variant Arg107Ile uses (5S,6S)-6-amino-5-hydroxycyclohexa-1,3-diene-1-carboxylate (2,3-trans-CHA) as substrate with >6-fold conversion compared to wild-type MenD. The novel compounds have been made accessible in vivo (up to 5.3 g L -1 ). Unexpectedly, as the identified residues such as Arg107 are highly conserved (>94 %), some of the designed variations can be found in wild-type SEPHCHC synthases from other bacteria (Arg107Lys, 0.3 %). This raises the question for the possible natural occurrence of as yet unexplored branches of the shikimate pathway., (© 2019 The Authors. Published by Wiley-VCH Verlag GmbH & Co. KGaA.)

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/., (© 2019 Wiley Periodicals, Inc.)

Aqueous solutions of Candida antarctica lipase B (CALB) were simulated considering three different water models (SPC/E, TIP3P, TIP4P) by a series of molecular dynamics (MD) simulations of three different box sizes (L = 9, 14, and 19 nm) to determine the diffusion coefficient, the water viscosity and the protein density. The protein-water systems were equilibrated for 500 ns, followed by 100 ns production runs which were analysed. The diffusional properties of CALB were characterized by the Stokes radius (R S ), which was derived from the diffusion coefficient and the viscosity. R S was compared to the geometric radius (R G ) of CALB, which was derived from the protein density. R S and R G differed by 0.27 nm for SPC/E and by 0.40 and 0.39 nm for TIP3P and TIP4P, respectively, which characterizes the thickness of the diffusive hydration layer on the protein surface. The simulated hydration layer of CALB resulted in agreement with those experimentally determined for other seven different proteins of comparable size. By avoiding the most common pitfalls, protein diffusion can be reliably simulated: simulating different box sizes to account for the finite size effect, equilibrating the protein-water system sufficiently, and using the complete production run for the determination of the diffusion coefficient.

The effect of water content on the static and dynamic properties of the deep eutectic solvent glyceline is studied using molecular dynamics (MD) simulations. Static properties are additionally calculated using the PC-SAFT equation of state. Force fields calibrated on water-free glyceline show predictive power for density and water activity over the entire composition range. In contrast, the PC-SAFT approach using pseudo one-component or two-component modelling strategies performed better for the density or the water activity, respectively. The MD simulations show that at low water content, the hydrogen-bond network between glycerol molecules as well as between glycerol and the cholinium cation is hardly affected by the water molecules while at higher water content, glycerol-glycerol hydrogen bonds are replaced by glycerol-water hydrogen bonds indicating the formation of an aqueous solution accompanied by a strong decrease of the shear viscosity. At the same time, the thermodynamic activity of water increases such that the MD simulations are able to guide the optimal composition with respect to requirements in biocatalytic applications such as low viscosity and low water activity. The combined application of PC-SAFT to efficiently predict static properties and molecular dynamics simulations to predict static and dynamic properties offers a powerful framework in solvent design applications.

Chorismatases catalyse the cleavage of chorismate, yielding (dihydroxy-)benzoate derivatives, which often constitute starter units for pharmaceutically relevant secondary metabolites. Depending on their products, chorismatases have been classified into three different subfamilies. These can be assigned using a set of amino acid residues in the active site. Here, we describe five new chorismatases, two of them members of a new subfamily, which has been discovered through correlation analysis of homologous protein sequences. The enzymes from the new subfamily produce exclusively 4-hydroxybenzoate, the same compound as produced by the structurally unrelated chorismate lyases. This showcase of convergent evolution is an example of the existence of more than one pathway to central building blocks. In contrast to chorismate lyases, however, chorismatases do not suffer from product inhibition (up to 2 mM 4-HBA), while the remaining kinetic parameters are in the same range; this makes them an interesting alternative for biocatalytic applications.

The ω-Transaminase Engineering Database (oTAED) was established as a publicly accessible resource on sequences and structures of the biotechnologically relevant ω-transaminases (ω-TAs) from Fold types I and IV. The oTAED integrates sequence and structure data, provides a classification based on fold type and sequence similarity, and applies a standard numbering scheme to identify equivalent positions in homologous proteins. The oTAED includes 67 210 proteins (114 655 sequences) which are divided into 169 homologous families based on global sequence similarity. The 44 and 39 highly conserved positions which were identified in Fold type I and IV, respectively, include the known catalytic residues and a large fraction of glycines and prolines in loop regions, which might have a role in protein folding and stability. However, for most of the conserved positions the function is still unknown. Literature information on positions that mediate substrate specificity and stereoselectivity was systematically examined. The standard numbering schemes revealed that many positions which have been described in different enzymes are structurally equivalent. For some positions, multiple functional roles have been suggested based on experimental data in different enzymes. The proposed standard numbering schemes for Fold type I and IV ω-TAs assist with analysis of literature data, facilitate annotation of ω-TAs, support prediction of promising mutation sites, and enable navigation in ω-TA sequence space. Thus, it is a useful tool for enzyme engineering and the selection of novel ω-TA candidates with desired biochemical properties., (© 2018 Wiley Periodicals, Inc.)

Lipase-catalyzed transesterification of triglycerides and alcohols to obtain biodiesel is an environmentally friendly and sustainable route for fuels production since, besides proceeding in mild reaction conditions, it allows for the use of low-cost feedstocks that contain water and free fatty acids, for example non-edible oils and waste oils. This review article reports recent advances in the field and focus in particular on a major issue in the enzymatic process, the inactivation of most lipases caused by methanol, the preferred acyl acceptor used for alcoholysis. The recent results about immobilization of enzymes on nano-materials and the use of whole-cell biocatalysts, as well as the use of cell-surface display technologies and metabolic engineering strategies for microbial production of biodiesel are described. It is discussed also insight into the effects of methanol on lipases obtained by modeling approaches and report on studies aimed at mining novel alcohol stable enzymes or at improving robustness in existing ones by protein engineering., (© 2018 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Macrokinetic Michaelis-Menten models based on thermodynamic activity provide insights into enzyme kinetics because they separate substrate-enzyme from substrate-solvent interactions. Kinetic parameters are estimated from experimental progress curves of enzyme-catalyzed reactions. Three pitfalls are discussed: deviations between thermodynamic and concentration-based models, product effects on the substrate activity coefficient, and product inhibition., (Copyright © 2017 Elsevier Ltd. All rights reserved.)

Esterases receive special attention because of their wide distribution in biological systems and environments and their importance for physiology and chemical synthesis. The prediction of esterases' substrate promiscuity level from sequence data and the molecular reasons why certain such enzymes are more promiscuous than others remain to be elucidated. This limits the surveillance of the sequence space for esterases potentially leading to new versatile biocatalysts and new insights into their role in cellular function. Here, we performed an extensive analysis of the substrate spectra of 145 phylogenetically and environmentally diverse microbial esterases, when tested with 96 diverse esters. We determined the primary factors shaping their substrate range by analyzing substrate range patterns in combination with structural analysis and protein-ligand simulations. We found a structural parameter that helps rank (classify) the promiscuity level of esterases from sequence data at 94% accuracy. This parameter, the active site effective volume, exemplifies the topology of the catalytic environment by measuring the active site cavity volume corrected by the relative solvent accessible surface area (SASA) of the catalytic triad. Sequences encoding esterases with active site effective volumes (cavity volume/SASA) above a threshold show greater substrate spectra, which can be further extended in combination with phylogenetic data. This measure provides also a valuable tool for interrogating substrates capable of being converted. This measure, found to be transferred to phosphatases of the haloalkanoic acid dehalogenase superfamily and possibly other enzymatic systems, represents a powerful tool for low-cost bioprospecting for esterases with broad substrate ranges, in large scale sequence data sets.

The tail-to-head terpene cyclization is arguably one of the most complex reactions found in nature. The hydrogen-bond-based resorcinarene capsule represents the first man-made enzyme-like catalyst that is capable of catalyzing this reaction. Based on noncovalent interactions between the capsule and the substrate, the product selectivity can be tuned by using different leaving groups. A detailed mechanistic investigation was performed to elucidate the reaction mechanism. For the cyclization of geranyl acetate, it was found that the cleavage of the leaving group is the rate-determining step. Furthermore, the studies revealed that trace amounts of acid are required as cocatalyst. A series of control experiments demonstrate that a synergistic interplay between the supramolecular capsule and the acid traces is required for catalytic activity.

The experimentally determined Michaelis constant K m results from a combination of two effects: the recognition of the substrate by the enzyme and the interactions between substrate and solvent. For a solvent-independent analysis of substrate specificity, the thermodynamic activity of the substrate, rather than its concentration, must be considered., (Copyright © 2017 Elsevier Ltd. All rights reserved.)

Candida antarctica lipase B (CalB) acts as a lipase when adsorbed to an acylglyceride interface and as an esterase when exposed to an aqueous environment. The effect of the molecular self-assembly nanostructure of triglyceride-water interfaces on structural conformations of adsorbed CalB and the implications to its catalytic function were studied by molecular dynamics simulations. Systems of CalB adsorbed to interfaces and solvated in water were compared. The two environments induced relative motions of helices α5 and α10 that resulted in open and closed conformations. The open conformation was stabilized by interactions between the polar and nonpolar amino acids of α5 and α10 and the nanostructure of triglyceride aggregates, which self-assembled into crystalline-like patterns of alternating polar and nonpolar lamellae. Thus, the structure of CalB has been adapted by evolution to the geometric constraints imposed by the interface nanostructure for optimized catalytic activity. Helices α5 and α10 have two functions. As mobile elements, they ensure access of bulky substrates to the active site in the open conformation. As a part of the active site pocket, they ensure binding of substrate molecules in a productive orientation near the active site. In water, access to the binding site is limited, and the smaller substrate binding site is beneficial for the binding of small, water-soluble substrates. The CalB crystal structure commonly used for protein engineering studies represents an intermediate state between open and closed, and may thus not be adequate to assess the function of CalB, neither as lipase nor as esterase.

The development of novel enzymes for biocatalytic processes requires knowledge on substrate profile and selectivity; this can be derived from databases and from publications. Often, these sources lack time-course data for the substrate or product, and an unambiguous link between experiment and enzyme sequence. The lack of integrated, original data hampers the comprehensive analysis of enzyme kinetics and the evaluation of sequence-function relationships. In order to accelerate enzyme engineering, BioCatNet integrates protein sequence, protein structure, and experimental data for a given enzyme family. BioCatNet explicitly assigns the enzyme sequence to the experimental data, which consists of information on reaction conditions and time-course data. BioCatNet facilitates the consistent documentation of reaction conditions, the archiving of time-course data, and the efficient exchange of experimental data among collaborators. Data integration is demonstrated for three case studies by using the TEED (Thiamine diphosphate-dependent Enzymes Engineering Database)., (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim.)

The binding of solvent molecules to a protein surface was modeled by molecular dynamics simulations of of Candida antarctica (C. antarctica) lipase B in binary mixtures of water, methanol, and toluene. Two models were analyzed: a competitive Langmuir model which assumes identical solvent binding sites with a different affinity toward water (KWat), methanol (KMet), and toluene (KTol) and a competitive Langmuir model with an additional interaction between free water and already bound water (KWatWat). The numbers of protein-bound molecules of both components of a binary mixture were determined for different compositions as a function of their thermodynamic activities in the bulk phase, and the binding constants were simultaneously fitted to the six binding curves (two components of three different mixtures). For both Langmuir models, the values of KWat, KMet, and KTol were highly correlated. The highest binding affinity was found for methanol, which was almost 4-fold higher than the binding affinities of water and toluene (KMet ≫ KWat ≈ KTol). Binding of water was dominated by the water-water interaction (KWatWat). Even for the three protein surface patches of highest water affinity, the binding affinity of methanol was 2-fold higher than water and 8-fold higher than toluene (KMet > KWat > KTol). The Langmuir model provides insights into the protein destabilizing mechanism of methanol which has a high binding affinity toward the protein surface. Thus, destabilizing solvents compete with intraprotein interactions and disrupt the tertiary structure. In contrast, benign solvents such as water or toluene have a low affinity toward the protein surface. Water is a special solvent: only few water molecules bind directly to the protein; most water molecules bind to already bound water molecules thus forming water patches. A quantitative mechanistic model of protein-solvent interactions that includes competition and miscibility of the components contributes a robust basis for solvent and protein engineering.

The regioselective terminal hydroxylation of alkanes and fatty acids is of great interest in a variety of industrial applications, such as in cosmetics, in fine chemicals, and in the fragrance industry. The chemically challenging activation and oxidation of non-activated C-H bonds can be achieved with cytochrome P450 enzymes. CYP153AM.aq. -CPRBM3 is an artificial fusion construct consisting of the heme domain from Marinobacter aquaeolei and the reductase domain of CYP102A1 from Bacillus megaterium. It has the ability to hydroxylate medium- and long-chain fatty acids selectively at their terminal positions. However, the activity of this interesting P450 construct needs to be improved for applications in industrial processes. For this purpose, the design of mutant libraries including two consecutive steps of mutagenesis is demonstrated. Targeted positions and residues chosen for substitution were based on semi-rational protein design after creation of a homology model of the heme domain of CYP153AM.aq. , sequence alignments, and docking studies. Site-directed mutagenesis was the preferred method employed to address positions within the binding pocket, whereas diversity was created with the aid of a degenerate codon for amino acids located at the substrate entrance channel. Combining the successful variants led to the identification of a double variant-G307A/S233G-that showed alterations of one position within the binding pocket and one position located in the substrate access channel. This double variant showed twofold increased activity relative to the wild type for the terminal hydroxylation of medium-chain-length fatty acids. This variant furthermore showed improved activity towards short- and long-chain fatty acids and enhanced stability in the presence of higher concentrations of fatty acids., (© 2016 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Chiral amines are valuable building blocks for the production of a variety of pharmaceuticals, agrochemicals and other specialty chemicals. Only recently, imine reductases (IREDs) were discovered which catalyze the stereoselective reduction of imines to chiral amines. Although several IREDs were biochemically characterized in the last few years, knowledge of the reaction mechanism and the molecular basis of substrate specificity and stereoselectivity is limited. To gain further insights into the sequence-function relationships, the Imine Reductase Engineering Database (www.IRED.BioCatNet.de) was established and a systematic analysis of 530 putative IREDs was performed. A standard numbering scheme based on R-IRED-Sk was introduced to facilitate the identification and communication of structurally equivalent positions in different proteins. A conservation analysis revealed a highly conserved cofactor binding region and a predominantly hydrophobic substrate binding cleft. Two IRED-specific motifs were identified, the cofactor binding motif GLGxMGx(5 )[ATS]x(4) Gx(4) [VIL]WNR[TS]x(2) [KR] and the active site motif Gx[DE]x[GDA]x[APS]x(3){K}x[ASL]x[LMVIAG]. Our results indicate a preference toward NADPH for all IREDs and explain why, despite their sequence similarity to β-hydroxyacid dehydrogenases (β-HADs), no conversion of β-hydroxyacids has been observed. Superfamily-specific conservations were investigated to explore the molecular basis of their stereopreference. Based on our analysis and previous experimental results on IRED mutants, an exclusive role of standard position 187 for stereoselectivity is excluded. Alternatively, two standard positions 139 and 194 were identified which are superfamily-specifically conserved and differ in R- and S-selective enzymes., (© 2016 Wiley Periodicals, Inc.)

The Lactamase Engineering Database (www.LacED.uni-stuttgart.de) was developed to facilitate the classification and analysis of TEM β-lactamases. The current version contains 474 TEM variants. Two hundred fifty-nine variants form a large scale-free network of highly connected point mutants. The network was divided into three subnetworks which were enriched by single phenotypes: one network with predominantly 2be and two networks with 2br phenotypes. Fifteen positions were found to be highly variable, contributing to the majority of the observed variants. Since it is expected that a considerable fraction of the theoretical sequence space is functional, the currently sequenced 474 variants represent only the tip of the iceberg of functional TEM β-lactamase variants which form a huge natural reservoir of highly interconnected variants. Almost 50% of the variants are part of a quartet. Thus, two single mutations that result in functional enzymes can be combined into a functional protein. Most of these quartets consist of the same phenotype, or the mutations are additive with respect to the phenotype. By predicting quartets from triplets, 3,916 unknown variants were constructed. Eighty-seven variants complement multiple quartets and therefore have a high probability of being functional. The construction of a TEM β-lactamase network and subsequent analyses by clustering and quartet prediction are valuable tools to gain new insights into the viable sequence space of TEM β-lactamases and to predict their phenotype. The highly connected sequence space of TEM β-lactamases is ideally suited to network analysis and demonstrates the strengths of network analysis over tree reconstruction methods., (Copyright © 2016, American Society for Microbiology. All Rights Reserved.)

Cytochrome P450 monooxygenases can catalyse the stereoselective C-H activation of a very broad range of substrates. Prediction and control of enantioselectivity of this enzyme class is of great interest for the synthesis of high-value chiral molecules. Here we have used a combination of molecular dynamics simulations and experimental screening to study the enantioselectivity of a library of active-site mutants of chimeric P450cam-RhFRed towards the benzylic hydroxylation of structurally related regioisomers of ethylmethylbenzene. Small variations either in substrate structure or in enzyme active site architecture were shown to lead to dramatic changes in enantioselectivity; this was broadly in agreement with computational predictions. In addition to validating computational approaches, these studies have provided us with a deeper understanding of effects that might control stereoselectivity in these biooxidation reactions., (© 2016 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Cytochrome P450 monooxygenases (CYPs) are a large, highly diverse protein family with a common fold. The sequences, structures, and functions of CYPs have been extensively studied resulting in more than 53,000 scientific articles. A sequence-based literature mining algorithm was designed to systematically analyze this wealth of information on SNPs, designed mutations, structural interactions, or functional roles of individual residues. Structurally corresponding positions in different CYPs were compared and universal selectivity-determining positions were identified. Based on the Cytochrome P450 Engineering Database (www.CYPED.BioCatNet.de) and a standard numbering scheme for all CYPs, 4000 residues in 168 CYPs mentioned in 2400 articles could be assigned to 440 structurally corresponding standard positions of the CYP fold, covering 96% of all standard positions. Seventeen individual standard positions were mentioned in the context of more than 32 different CYPs. The majority of these most frequently mentioned positions are located on the six substrate recognition sites and are involved in control of selectivity, such as the well-studied position 87 in CYP102A1 (P450(BM-3)) which was mentioned in the articles on 63 different CYPs. The recurrent citation of the 17 frequently mentioned positions for different CYPs suggests their universal functional relevance., (© 2015 Wiley Periodicals, Inc.)

The underlying molecular mechanisms of macroscopic excess properties were studied by molecular dynamics simulations for different compositions of methanol-water mixtures. Structural data (nearest neighbor relationships, clustering analysis) and dynamic data (hydrogen bond lifetimes, rotational autocorrelation, translational diffusion) were evaluated. Nearest neighbor relationships provide quantitative evidence and a pictorial description of incomplete mixing at the molecular level as a source for mixture anomalies, while a comparative study of water surrounding methyl moieties versus water in the bulk-like environment provides evidence against the hydrophobicity model of clathrate-like hydration. Furthermore, the formation or breakdown of the system-wide hydrogen bonding network at a critical threshold of approximately equimolar mixture is perceived to separate the mixture system into two hydrogen bonding regimes: hydrogen-bonded water clusters embedded in methanol for mixtures with low water content and methanol molecules within a system-wide hydrogen-bonded water network for mixtures with high water content.

Interactions of lipases with hydrophobic substrate-water interfaces are of great interest to design improved lipase variants and engineer reaction conditions. This chapter describes the necessary steps to carry out molecular dynamics simulations of Candida antarctica lipase B at tributyrin-water interface using the GROMACS simulation software. Special attention is drawn to the preparation of the protein and the substrate-water interface and to the analysis of the obtained trajectory.

The Omega-loop of TEM beta-lactamase is involved in substrate recognition and catalysis. Its dynamical properties and interaction with water molecules were investigated by performing multiple molecular dynamics simulations of up to 50 ns. Protein flexibility was assessed by calculating the root mean-square fluctuations and the generalized order parameter, S(2). The residues in secondary structure elements are highly ordered, whereas loop regions are more flexible, which is in agreement with previous experimental observations. Interestingly, the Omega-loop (residues 161-179) is rigid with order parameters similar to secondary structure elements, with the exception of the tip of the loop (residues 173-177) that has a considerably higher flexibility and performs an opening and closing motion on the 50-ns timescale. The rigidity of the main part of the Omega-loop is mediated by stabilizing and highly conserved water bridges inside a cavity lined by the Omega-loop and residues 65-69 of the protein core. In contrast, the flexible tip of the Omega-loop lacks these interactions. Hydration of the cavity and exchange of the water molecules with the bulk solvent occurs via two pathways: the flexible tip that serves as a door to the cavity, and a temporary water channel involving the side chain of Arg(164).

The Medium-Chain Dehydrogenase/Reductase Engineering Database (MDRED, http://www.mdred.uni-stuttgart.de) has been established to serve as an analysis tool for a systematic investigation of sequence-structure-function relationships. It includes sequence and structure information of 2684 and 42 medium-chain dehydrogenases/reductases (MDRs), respectively. Although MDRs are very diverse in sequence, they have a conserved tertiary structure. MDRs are assigned to 199 homologous families and 29 superfamilies. For each family, annotated multiple sequence alignments are provided, and functionally relevant residues are annotated. Twenty-five superfamilies were classified as zinc-containing MDRs, four as non-zinc-containing MDRs. For the zinc-containing MDRs, three subclasses were identified by systematic analysis of a variable loop region, the quaternary structure determining loop (QSDL): the class of short, medium, and long QSDL, which include 11, 3, and 5 superfamilies, respectively. The length of the QSDL is predictive for tetramer (short QSDL) and dimer (long QSDL) formation. The class of medium QSDL includes both tetrameric and dimeric MDRs. The shape of the substrate-binding site is highly conserved in all zinc-containing MDRs with the exception of two variable regions, the substrate recognition sites (SRS): two residues located on the QSDL (SRS1) and, for the class of long QSDL, one residue located in the catalytic domain (SRS2). The MDRED is the first online-accessible resource of MDRs that integrates information on sequence, structure, and function. Annotation of functionally relevant residues assist the understanding of sequence-structure-function relationships. Thus, the MDRED serves as a valuable tool to identify potential hotspots for engineering properties such as substrate specificity.

A set of 49 high-resolution (

Immunoliposomes generated by coupling of antibodies to the liposomal surface allow for an active targeting of entrapped compounds to diseased areas. Single-chain Fv fragments (scFv) represent the smallest part of an antibody containing the entire antigen-binding site. They can be coupled in a defined and site-directed manner through genetically engineered cysteine residues, for example, those added at the C-terminus. Here, we have performed a comparative analysis of various scFv' variants with cysteine residues present at the end of a C-terminal extension of varying length and composition (HC variants) or introduced in the linker sequence connecting the variable heavy and light chain domain (LC variants). Using a scFv fragment directed against fibroblast activation protein (FAP) as a model antibody, we could show that all variants can be employed for the generation of active immunoliposomes, although the presence of three additional cysteine residues in one scFv' molecule resulted in decreased binding of immunoliposomes compared to that of immunoliposomes generated with scFv' molecules containing only one additional cysteine residue. In order to further improve the scFv' format by reducing the number of additional amino acid residues, we also generated molecules with the hexahistidyl-tag incorporated into the linker sequence together with a cysteine residue either at position 1 or 3 of the linker sequence (LCH variants). These newly designed scFv' molecules may be particularly suitable for the generation of immunoliposomes and other antibody conjugates, limiting the number of additional residues in these antibody molecules to a minimum.

The mouse strain H-2(aw18) shows typical characteristics of 21-hydroxylase deficiency (21-OHD). A deletion of the active Cyp21a1 gene has been postulated; however, the changes on the nucleotide level are still unknown. To investigate whether this animal model, the only one available, is suitable for studying congenital adrenal hyperplasia in man, a detailed analysis of the Cyp21 locus has been performed to ascertain the genetic cause of 21-OHD in H-2(aw18) mice. We demonstrate that 21-OHD is caused by unequal crossing over between the active Cyp21a1 gene and the pseudogene resulting in a hybrid Cyp21a1-Cyp21a2-p gene including a partial deletion of Cyp21a1. Next to several pseudogene-specific point mutations, various novel missense mutations and a nonsense mutation are present. Enzyme activity for each point mutation has been determined in vitro and the structure-function relationship has been studied by sequence conservation analysis and a three-dimensional murine 21-hydroxylase protein (Cyp21) structure model. The mutations are classified in three classes: I, no or minor decrease in enzyme activity: R238Q, P465L, R361K, A362V, P458L; II, loss of enzyme activity caused by inefficient electron flux: R346H, R400C; III, loss of activity due to deficient substrate binding: I462F, L464F. The combination of in vitro protein expression and three-dimensional structure modeling provides a valuable tool to understand the role of the different mutations and polymorphisms on the resulting enzyme activity. The underlying genetic mechanisms are also known to be responsible for 21-OHD in humans, so rodent 21-OHD turns out to be an excellent genetic model for studying the human disease.

Family GH16 glycoside hydrolases can be assigned to five subgroups according to their substrate specificities, including xyloglucan transglucosylases/hydrolases (XTHs), (1,3)-beta-galactanases, (1,4)-beta-galactanases/kappa-carrageenases, "nonspecific" (1,3/1,3;1,4)-beta-D-glucan endohydrolases, and (1,3;1,4)-beta-D-glucan endohydrolases. A structured family GH16 glycoside hydrolase database has been constructed (http://www.ghdb.uni-stuttgart.de) and provides multiple sequence alignments with functionally annotated amino acid residues and phylogenetic trees. The database has been used for homology modeling of seven glycoside hydrolases from the GH16 family with various substrate specificities, based on structural coordinates for (1,3;1,4)-beta-D-glucan endohydrolases and a kappa-carrageenase. In combination with multiple sequence alignments, the models predict the three-dimensional (3D) dispositions of amino acid residues in the substrate-binding and catalytic sites of XTHs and (1,3/1,3;1,4)-beta-d-glucan endohydrolases; there is no structural information available in the databases for the latter group of enzymes. Models of the XTHs, compared with the recently determined structure of a Populus tremulos x tremuloides XTH, reveal similarities with the active sites of family GH11 (1,4)-beta-D-xylan endohydrolases. From a biological viewpoint, the classification, molecular modeling and a new 3D structure of the P. tremulos x tremuloides XTH establish structural and evolutionary connections between XTHs, (1,3;1,4)-beta-D-glucan endohydrolases and xylan endohydrolases. These findings raise the possibility that XTHs from higher plants could be active not only on cell wall xyloglucans, but also on (1,3;1,4)-beta-D-glucans and arabinoxylans, which are major components of walls in grasses. A role for XTHs in (1,3;1,4)-beta-D-glucan and arabinoxylan modification would be consistent with the apparent overrepresentation of XTH sequences in cereal expressed sequence tags databases.

The thienopyridine derivatives ticlopidine and clopidogrel are inhibitors of ADP-induced platelet aggregation. Pharmacological activity of these prodrugs depends on cytochrome P450 (P450)-dependent oxidation to the active antithrombotic agent. In this study, we investigated the interaction potential of clopidogrel and ticlopidine by using human liver microsomes and recombinantly expressed P450 isoforms. Both clopidogrel and ticlopidine inhibited CYP2B6 with highest potency and CYP2C19 with lower potency. Clopidogrel also inhibited CYP2C9, and ticlopidine also inhibited CYP1A2, with lower potency. Inhibition of CYP2B6 was time- and concentration-dependent, and as shown by dialysis experiments, it was irreversible and dependent on NADPH, suggesting a mechanism-based mode of action. Inactivation was of nonpseudo-firstorder type with maximal rates of inactivation (K(inact)) for clopidogrel and ticlopidine in microsomes (recombinant CYP2B6) of 0.35 (1.5 min(-1)) and 0.5 min(-1) (0.8 min(-1)), respectively, and half-maximal inactivator concentrations (KI) were 0.5 microM (1.1 microM) for clopidogrel and 0.2 microM (0.8 microM) for ticlopidine. Inhibition was attenuated by the presence of alternative active site ligands but not by nucleophilic trapping agents or reactive oxygen scavengers, further supporting mechanism-based action. A chemical mechanism is discussed based on the known metabolic activation of clopidogrel and on the finding that hemoprotein integrity of recombinant CYP2B6 was not affected by irreversible inhibition. These results suggest the possibility of drug interactions between thienopyridine derivates and drug substrates of CYP2B6 and CYP2C19.

Two models, a purely nonbonded model and a cationic dummy atom approach, were examined for the modeling of the binuclear zinc-containing IMP-1 metallo-beta-lactamase in complex with a mercaptocarboxylate inhibitor. The cationic dummy atom approach had substantial advantages as it maintained the initial, experimentally determined geometry of the metal-containing active site during molecular dynamics simulations in water. The method was extended to the modeling of the free enzyme and the enzyme in complex with a cephalosporin substrate docked in an intermediate structure. For all three systems, the modeled complexes and the tetrahedral coordination of the zinc ions were stable. The average zinc-zinc distance increased by approximately 1 A in the substrate complex compared with the inhibitor complex and the free enzyme in which a hydroxide ion acts as a bridging ligand. Thus, the zinc ions are predicted to undergo a back and forth movement upon the cycle of hydrolysis. In contrast to previous assumptions, no interaction of the Asn167 side chain with the bound cephalosporin substrate was observed. Our observations are in agreement with quantum-mechanical calculations and experimental data and indicate that the cationic dummy atom approach is useful to model zinc-containing metallo-beta-lactamases as free proteins, in complex with inhibitors and in complex with substrates.