Navigating within thiamine diphosphate-dependent decarboxylases: Sequences, structures, functional positions, and binding sites.

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.

Investigating the Role of Conformational Effects on Laccase Stability and Hyperactivation under Stress Conditions.

Fungal laccase from Steccherinum ochraceum 1833 displays remarkable stability under different harsh conditions: organic/buffer mixtures, thermal treatment, and microwave radiation. The behavior is particularly significant in the light of the sharp inactivation observed for two different fungal laccases. Laccase from S. ochraceum 1833 also displays hyperactivation under mild thermal treatment (60 °C). Molecular dynamics simulations at 80 °C explained how this laccase retains the geometry of the electron transfer pathway, thereby assuring electron transfer through the copper ions and thus maintaining its catalytic activity at high temperature. Spectroscopic studies revealed that the thermal activation corresponds to specific conformational changes in the protein. The results indicate that this laccase is potentially applicable under denaturing conditions that might be beneficial for the biotransformation of recalcitrant substrates.

Efficient immobilisation of industrial biocatalysts: criteria and constraints for the selection of organic polymeric carriers and immobilisation methods.

Efficient immobilisation protocols are the result of perfect matching of factors depending on the enzyme, the process and the support for immobilisation. Physical-chemical phenomena, such as partition, solvation and diffusion, strongly affect the efficiency of the biocatalyst in each specific reaction system. Therefore, tailored solutions must be developed for each specific process of interest. Indeed, direct investigation of what occurs at the molecular level in a reaction catalysed by an immobilised enzyme is a quite formidable task and observed differences in the performance of immobilised biocatalysts must be interpreted very carefully. In any study dealing with enzyme immobilisation the prerequisite is the rigorous planning and reporting of experiments, being aware of the complexity of these multi-phase systems.

Biocatalytic Heteroaromatic Amide Formation in Water Enabled by a Catalytic Tetrad and Two Access Tunnels.

The amide moiety belongs to the most common motives in pharmaceutical chemistry, present in many prescribed small-molecule pharmaceuticals. Methods for its manufacture are still in high demand, especially using water/buffer as a solvent and avoiding stoichiometric amounts of activation reagents. Herein, we identified from a library of lipases/esterases/acyltransferases and variants thereof a lipase originating from Sphingomonas sp. HXN-200 (SpL) able to form amides in aqueous solution starting from a broad scope of sterically demanding heteroaromatic ethyl esters as well as aliphatic amines, reaching isolated yields up to 99% on preparative scale and space time yields of up to 864 g L-1 d-1; thus, in selected cases, the amide was formed within minutes. The enzyme features an aspartate next to the canonical serine of the catalytic triad, which was essential for amide formation. Furthermore, the enzyme structure revealed two tunnels to the active site, presumably one for the ester and one for the amine, which permit the bringing together of the sterically demanding heteroaromatic esters and the amine in the active site. This work shows that biocatalytic amide formation starting from various five- and six-membered heteroaromatic ethyl esters in the buffer can serve as a platform for preparative amide synthesis.

Molecular Mechanism of Methanol Inhibition in CALB-Catalyzed Alcoholysis: Analyzing Molecular Dynamics Simulations by a Markov State Model.

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.

Combined Linear Interaction Energy and Alchemical Solvation Free-Energy Approach for Protein-Binding Affinity Computation.

Calculating free energies of binding (ΔGbind) 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 ΔGbind 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 ΔGbind 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.

Diverse effects of aqueous polar co-solvents on Candida antarctica lipase B.

Biocatalysis in mixtures of water and co-solvents represents an opportunity to expand the application of enzymes. However, in the presence of organic solvents, enzymes can undergo reversible inhibition, inactivation, or aggregation. In this work, we studied the effects of three co-solvents (methanol, acetone, and dimethyl sulfoxide - DMSO) on the function and structure of the recombinant Candida antarctica lipase B (rCALB), a widely used enzyme in biotechnological applications. The effects of co-solvents on rCALB were investigated by steady-state kinetics experiments, biophysical assays and by molecular dynamics simulations in the presence and upon incubation with the three co-solvents. Methanol and acetone were found to act as competitive inhibitors of rCALB and to promote its aggregation, whereas DMSO is a non-essential activator of rCALB.

Simulation of protein diffusion: a sensitive probe of protein-solvent interactions.

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 (RS), which was derived from the diffusion coefficient and the viscosity. RS was compared to the geometric radius (RG) of CALB, which was derived from the protein density. RS and RG 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.

Modelling of substrate access and substrate binding to cephalosporin acylases.

Semisynthetic cephalosporins are widely used antibiotics currently produced by different chemical steps under harsh conditions, which results in a considerable amount of toxic waste. Biocatalytic synthesis by the cephalosporin acylase from Pseudomonas sp. strain N176 is a promising alternative. Despite intensive engineering of the enzyme, the catalytic activity is still too low for a commercially viable process. To identify the bottlenecks which limit the success of protein engineering efforts, a series of MD simulations was performed to study for two acylase variants (WT, M6) the access of the substrate cephalosporin C from the bulk to the active site and the stability of the enzyme-substrate complex. In both variants, cephalosporin C was binding to a non-productive substrate binding site (E86α, S369ß, S460ß) at the entrance to the binding pocket, preventing substrate access. A second non-productive binding site (G372ß, W376ß, L457ß) was identified within the binding pocket, which competes with the active site for substrate binding. Noteworthy, substrate binding to the protein surface followed a Langmuir model resulting in binding constants K = 7.4 and 9.2 mM for WT and M6, respectively, which were similar to the experimentally determined Michaelis constants KM = 11.0 and 8.1 mM, respectively.

Enantioselective production of 3-hydroxy metabolites of tibolone by yeast reduction.

The enantioselective reduction of tibolone into the corresponding 3alpha-hydroxy or 3beta-hydroxy metabolite can be controlled by choosing suited strains of yeasts and biotransformation conditions. A restricted screening performed among 52 yeasts showed that the 3alpha-epimer was preferentially obtained with high epimeric purity with various strains (i.e. with Kluyveromyces lactis CBS 2359), while only Saccharomyces cerevisiae CBS 3093 gave the 3beta-epimer as major product. The reduction of tibolone with K. lactis CBS 2359 and S. cerevisiae CBS 3093 was optimised. S. cerevisiae CBS 3093 furnished a 96:4 ratio of 3beta/3alpha with complete molar conversion within 72h when the initial concentration of substrate was below 2.5g/L. K. lactis CBS 2359 gave a 99:1 ratio of 3alpha/3beta with complete conversion in 64h.

Interpretation of cytochrome P450 monooxygenase kinetics by modeling of thermodynamic activity.

The experimentally determined Michaelis constant Kmc results from a combination of two effects: the recognition of the substrate by the enzyme and the molecular interactions between substrate and solvent. By separating substrate recognition from solvent effects, the thermodynamic activity-based Michaelis constant Kma allows for an unambiguous comparison of how different substrates fit into the substrate binding site. Kma of a poorly water-soluble substrate is calculated from the experimentally determined concentration-based Michaelis constant Kmc and its activity coefficient at infinite dilution γ∞. Comparing the Kma of different substrates instead of the experimentally determined Kmc prevents misinterpretations of the molecular basis of enzyme-substrate interactions. While n-octane showed the lowest Kmc value of six P450BM-3 substrates, its Kma was 500 fold higher than aniline, indicating that the binding of n-octane is mainly driven by its low water solubility, while binding of aniline is driven by its shape complementarity. For three substrates (aniline, oct-1-yne, n-octane), γ∞ was reliably calculated by molecular dynamics simulations, either in binary substrate-water mixtures or in ternary mixtures including DMSO as cosolvent. It is demonstrated that the widely used DMSO has a considerable effect on the measured Kmc value. Depending on the substrate, addition of 10% v/v DMSO increases Kmc by up to a factor of 11. To make biocatalytic experiments reproducible, it is therefore of utmost importance to carefully report the reaction conditions. The reliable simulation of activity coefficients in complex mixtures allows for an unambiguous comparison of enzyme-substrate interactions and provides a predictive tool for the design of biocatalytic processes.

The Closure of the Cycle: Enzymatic Synthesis and Functionalization of Bio-Based Polyesters.

The polymer industry is under pressure to mitigate the environmental cost of petrol-based plastics. Biotechnologies contribute to the gradual replacement of petrol-based chemistry and the development of new renewable products, leading to the closure of carbon circle. An array of bio-based building blocks is already available on an industrial scale and is boosting the development of new generations of sustainable and functionally competitive polymers, such as polylactic acid (PLA). Biocatalysts add higher value to bio-based polymers by catalyzing not only their selective modification, but also their synthesis under mild and controlled conditions. The ultimate aim is the introduction of chemical functionalities on the surface of the polymer while retaining its bulk properties, thus enlarging the spectrum of advanced applications.

BioGPS descriptors for rational engineering of enzyme promiscuity and structure based bioinformatic analysis.

A new bioinformatic methodology was developed founded on the Unsupervised Pattern Cognition Analysis of GRID-based BioGPS descriptors (Global Positioning System in Biological Space). The procedure relies entirely on three-dimensional structure analysis of enzymes and does not stem from sequence or structure alignment. The BioGPS descriptors account for chemical, geometrical and physical-chemical features of enzymes and are able to describe comprehensively the active site of enzymes in terms of "pre-organized environment" able to stabilize the transition state of a given reaction. The efficiency of this new bioinformatic strategy was demonstrated by the consistent clustering of four different Ser hydrolases classes, which are characterized by the same active site organization but able to catalyze different reactions. The method was validated by considering, as a case study, the engineering of amidase activity into the scaffold of a lipase. The BioGPS tool predicted correctly the properties of lipase variants, as demonstrated by the projection of mutants inside the BioGPS "roadmap".

Lipases immobilization for effective synthesis of biodiesel starting from coffee waste oils.

Immobilized lipases were applied to the enzymatic conversion of oils from spent coffee ground into biodiesel. Two lipases were selected for the study because of their conformational behavior analysed by Molecular Dynamics (MD) simulations taking into account that immobilization conditions affect conformational behavior of the lipases and ultimately, their efficiency upon immobilization. The enzymatic synthesis of biodiesel was initially carried out on a model substrate (triolein) in order to select the most promising immobilized biocatalysts. The results indicate that oils can be converted quantitatively within hours. The role of the nature of the immobilization support emerged as a key factor affecting reaction rate, most probably because of partition and mass transfer barriers occurring with hydrophilic solid supports. Finally, oil from spent coffee ground was transformed into biodiesel with yields ranging from 55% to 72%. The synthesis is of particular interest in the perspective of developing sustainable processes for the production of bio-fuels from food wastes and renewable materials. The enzymatic synthesis of biodiesel is carried out under mild conditions, with stoichiometric amounts of substrates (oil and methanol) and the removal of free fatty acids is not required.

Mechanochemical activation of vincamine mediated by linear polymers: assessment of some "critical" steps.

The aim of the research was to investigate three "critical steps" that deserve particular attention during the mechanochemical activation of vincamine. The first step consisted in the selection of the best polymeric carrier/most affine stabiliser between linear PVP and NaCMC by using the GRID and the GRID based AutoDock software packages which permit to calculate their surface features and interactions. Moreover, the calculation of the partial and total solubility parameters supported the results obtained by GRID and AutoDock software. Then, after the selection of linear PVP-K30 as the suitable carrier, the influence of process and formulation variables on the amorphisation degree and solubility enhancement was studied, to select the most suitable process conditions and formulation parameters. Subsequently, the best performing samples were widely characterised using XRPD, TEM and SSNMR (including the proton relaxation ((1)H T1 NMR) time) techniques. These studies highlighted that all the coground samples were nanocrystalline solid dispersions indicating a dramatic difference between the amorphisation capacities of linear PVP-K30 and cross-linked PVP, used in previous analogous experiences. In particular, (13)C, (15)N and (1)H T1 NMR data point to a description of the system as a dispersion of nanocrystals in the polymer. In these dispersions vincamine is in a disordered crystalline state due to extensive interactions and contacts with PVP-K30 but the main hydrogen bonding motif characterising its packing remains. Again, differently from cross-linked PVP, dissolution studies revealed that linear PVP-K30 was able to promote a complete in vitro solubilisation of vincamine in some coground samples. What is more important, by using a linear polymer, drug-to-polymer and milling time variables appeared less influent on the solid state and in vitro properties of the composites. Finally, stability studies conducted for a period of 1year highlighted the high physical stability of the selected samples.

Endo- and exo-inulinases: enzyme-substrate interaction and rational immobilization.

Three-dimensional models of exoinulinase from Bacillus stearothermophilus and endoinulinase from Aspergillus niger were built up by means of homology modeling. The crystal structure of exoinulinase from Aspergillus awamori was used as a template, which is the sole structure of inulinase resolved so far. Docking and molecular dynamics simulations were performed to investigate the differences between the two inulinases in terms of substrate selectivity. The analysis of the structural differences between the two inulinases provided the basis for the explanation of their different regio-selectivity and for the understanding of enzyme-substrate interactions. Surface analysis was performed to point out structural features that can affect the efficiency of enzymes also after immobilization. The computational analysis of the three-dimensional models proved to be an effective tool for acquiring information and allowed to formulate an optimal immobilized biocatalyst even more active that the native one, thus enabling the full exploitation of the catalytic potential of these enzymes.