Analysis of Poly(ethylene terephthalate) degradation kinetics of evolved IsPETase variants using a surface crowding model.

Poly(ethylene terephthalate) (PET) is a major plastic polymer utilized in the single-use and textile industries. The discovery of PET-degrading enzymes (PETases) has led to an increased interest in the biological recycling of PET in addition to mechanical recycling. IsPETase from Ideonella sakaiensis is a candidate catalyst, but little is understood about its structure-function relationships with regards to PET degradation. To understand the effects of mutations on IsPETase productivity, we develop a directed evolution assay to identify mutations beneficial to PET film degradation at 30 °C. IsPETase also displays enzyme concentration-dependent inhibition effects, and surface crowding has been proposed as a causal phenomenon. Based on total internal reflectance fluorescence microscopy and adsorption experiments, IsPETase is likely experiencing crowded conditions on PET films. Molecular dynamics simulations of IsPETase variants reveal a decrease in active site flexibility in free enzymes and reduced probability of productive active site formation in substrate-bound enzymes under crowding. Hence, we develop a surface crowding model to analyze the biochemical effects of three hit mutations (T116P, S238N, S290P) that enhanced ambient temperature activity and/or thermostability. We find that T116P decreases susceptibility to crowding, resulting in higher PET degradation product accumulation despite no change in intrinsic catalytic rate. In conclusion, we show that a macromolecular crowding-based biochemical model can be used to analyze the effects of mutations on properties of PETases and that crowding behavior is a major property to be targeted for enzyme engineering for improved PET degradation.

Flexible active-site loops fine-tune substrate specificity of hyperthermophilic metallo-oxidases.

Hyperthermophilic ('superheat-loving') archaea found in high-temperature environments such as Pyrobaculum aerophilum contain multicopper oxidases (MCOs) with remarkable efficiency for oxidizing cuprous and ferrous ions. In this work, directed evolution was used to expand the substrate specificity of P. aerophilum McoP for organic substrates. Six rounds of error-prone PCR and DNA shuffling followed by high-throughput screening lead to the identification of a hit variant with a 220-fold increased efficiency (kcat/Km) than the wild-type for 2,2'-azino-bis(3-ethylbenzothiazoline-6-sulphonic acid) (ABTS) without compromising its intrinsic activity for metal ions. The analysis of the X-ray crystal structure reveals four proximal mutations close to the T1Cu active site. One of these mutations is within the 23-residues loop that occludes this site, a distinctive feature of prokaryotic MCOs. The increased flexibility of this loop results in an enlarged tunnel and one additional pocket that facilitates bulky substrate-enzyme interactions. These findings underscore the synergy between mutations that modulate the dynamics of the active-site loop enabling enhanced catalytic function. This study highlights the potential of targeting loops close to the T1Cu for engineering improvements suitable for biotechnological applications.

Rationally Guided Improvement of NOV1 Dioxygenase for the Conversion of Lignin-Derived Isoeugenol to Vanillin.

Biocatalysis is a key tool in both green chemistry and biorefinery fields. NOV1 is a dioxygenase that catalyzes the one-step, coenzyme-free oxidation of isoeugenol into vanillin and holds enormous biotechnological potential for the complete valorization of lignin as a sustainable starting material for biobased chemicals, polymers, and materials. This study integrates computational, kinetic, structural, and biophysical approaches to characterize a new NOV1 variant featuring improved activity and stability compared to those of the wild type. The S283F replacement results in a 2-fold increased turnover rate (kcat) for isoeugenol and a 4-fold higher catalytic efficiency (kcat/Km) for molecular oxygen compared to those of the wild type. Furthermore, the variant exhibits a half-life that is 20-fold higher than that of the wild type, which most likely relates to the enhanced stabilization of the iron cofactor in the active site. Molecular dynamics supports this view, revealing that the S283F replacement decreases the optimal pKa and favors conformations of the iron-coordinating histidines compatible with an increased level of binding to iron. Importantly, whole cells containing the S283F variant catalyze the conversion of ≤100 mM isoeugenol to vanillin, yielding >99% molar conversion yields within 24 h. This integrative strategy provided a new enzyme for biotechnological applications and mechanistic insights that will facilitate the future design of robust and efficient biocatalysts.

The ins and outs of vanillyl alcohol oxidase: Identification of ligand migration paths.

Vanillyl alcohol oxidase (VAO) is a homo-octameric flavoenzyme belonging to the VAO/PCMH family. Each VAO subunit consists of two domains, the FAD-binding and the cap domain. VAO catalyses, among other reactions, the two-step conversion of p-creosol (2-methoxy-4-methylphenol) to vanillin (4-hydroxy-3-methoxybenzaldehyde). To elucidate how different ligands enter and exit the secluded active site, Monte Carlo based simulations have been performed. One entry/exit path via the subunit interface and two additional exit paths have been identified for phenolic ligands, all leading to the si side of FAD. We argue that the entry/exit path is the most probable route for these ligands. A fourth path leading to the re side of FAD has been found for the co-ligands dioxygen and hydrogen peroxide. Based on binding energies and on the behaviour of ligands in these four paths, we propose a sequence of events for ligand and co-ligand migration during catalysis. We have also identified two residues, His466 and Tyr503, which could act as concierges of the active site for phenolic ligands, as well as two other residues, Tyr51 and Tyr408, which could act as a gateway to the re side of FAD for dioxygen. Most of the residues in the four paths are also present in VAO's closest relatives, eugenol oxidase and p-cresol methylhydroxylase. Key path residues show movements in our simulations that correspond well to conformations observed in crystal structures of these enzymes. Preservation of other path residues can be linked to the electron acceptor specificity and oligomerisation state of the three enzymes. This study is the first comprehensive overview of ligand and co-ligand migration in a member of the VAO/PCMH family, and provides a proof of concept for the use of an unbiased method to sample this process.

Distal mutations enhance efficiency of free and immobilized NOV1 dioxygenase for vanillin synthesis.

Protein engineering is crucial to improve enzymes' efficiency and robustness for industrial biocatalysis. NOV1 is a bacterial dioxygenase that holds biotechnological potential by catalyzing the one-step oxidation of the lignin-derived isoeugenol into vanillin, a popular flavoring agent used in food, cleaning products, cosmetics and pharmaceuticals. This study aims to enhance NOV1 activity and operational stability through the identification of distal hotspots, located at more than 9 Šfrom the active site using Zymspot, a tool that predicts advantageous distant mutations, streamlining protein engineering. A total of 41 variants were constructed using site-directed mutagenesis and the six most active enzyme variants were then recombined. Two variants, with two and three mutations, showed nearly a 10-fold increase in activity and up to 40-fold higher operational stability than the wild-type. Furthermore, these variants show 90-100 % immobilization efficiency in metal affinity resins, compared to approximately 60 % for the wild-type. In bioconversions where 50 mM of isoeugenol was added stepwise over 24-h cycles, the 1D2 variant produced approximately 144 mM of vanillin after six reaction cycles, corresponding to around 22 mg, indicating a 35 % molar conversion yield. This output was around 2.5 times higher than that obtained using the wild-type. Our findings highlight the efficacy of distal protein engineering in enhancing enzyme functions like activity, stability, and metal binding selectivity, thereby fulfilling the criteria for industrial biocatalysts. This study provides a novel approach to enzyme optimization that could have significant implications for various biotechnological applications.

In Silico Assessment of Biomolecule Reactivity with Leachables.

Leachables in pharmaceutical products may react with biomolecule active pharmaceutical ingredients (APIs), for example, monoclonal antibodies (mAb), peptides, and ribonucleic acids (RNA), potentially compromising product safety and efficacy or impacting quality attributes. This investigation explored a series of in silico models to screen extractables and leachables to assess their possible reactivity with biomolecules. These in silico models were applied to collections of known leachables to identify functional and structural chemical classes likely to be flagged by these in silico approaches. Flagged leachable functional classes included antimicrobials, colorants, and film-forming agents, whereas specific chemical classes included epoxides, acrylates, and quinones. In addition, a dataset of 22 leachables with experimental data indicating their interaction with insulin glargine was used to evaluate whether one or more in silico methods are fit-for-purpose as a preliminary screen for assessing this biomolecule reactivity. Analysis of the data showed that the sensitivity of an in silico screen using multiple methodologies was 80%-90% and the specificity was 58%-92%. A workflow supporting the use of in silico methods in this field is proposed based on both the results from this assessment and best practices in the field of computational modeling and quality risk management.

Electron transfer in proteins: theory, applications and future perspectives.

The study of electron transfer (ET) by means of computational techniques has experienced a great development in the last few decades. In particular, understanding the atomic details of its mechanism in complex biological systems is currently possible with a large range of different in silico modelling tools. We review here some theories and representative major contributions to this development. We also underline some of our group's main inputs, focusing on long range and protein-protein electron transfer, and analyse future perspectives. At the end of the article, we emphasize the importance of the basic electron transfer knowledge in the frame of medical and bioengineering applications: mitochondrial therapeutic targets, bioengineering for clean energy, and biosensors.

Enzyme immobilization studied through molecular dynamic simulations.

In recent years, simulations have been used to great advantage to understand the structural and dynamic aspects of distinct enzyme immobilization strategies, as experimental techniques have limitations in establishing their impact at the molecular level. In this review, we discuss how molecular dynamic simulations have been employed to characterize the surface phenomenon in the enzyme immobilization procedure, in an attempt to decipher its impact on the enzyme features, such as activity and stability. In particular, computational studies on the immobilization of enzymes using i) nanoparticles, ii) self-assembled monolayers, iii) graphene and carbon nanotubes, and iv) other surfaces are covered. Importantly, this thorough literature survey reveals that, while simulations have been primarily performed to rationalize the molecular aspects of the immobilization event, their use to predict adequate protocols that can control its impact on the enzyme properties is, up to date, mostly missing.

Corrigendum to "Unveiling molecular details behind improved activity at neutral to alkaline pH of an engineered DyP-type peroxidase" Comput Struct Biotechnol J, 20 (2022), 3899-3910.

[This corrects the article DOI: 10.1016/j.csbj.2022.07.032.].

Mechanistic insights into glycoside 3-oxidases involved in C-glycoside metabolism in soil microorganisms.

C-glycosides are natural products with important biological activities but are recalcitrant to degradation. Glycoside 3-oxidases (G3Oxs) are recently identified bacterial flavo-oxidases from the glucose-methanol-coline (GMC) superfamily that catalyze the oxidation of C-glycosides with the concomitant reduction of O2 to H2O2. This oxidation is followed by C-C acid/base-assisted bond cleavage in two-step C-deglycosylation pathways. Soil and gut microorganisms have different oxidative enzymes, but the details of their catalytic mechanisms are largely unknown. Here, we report that PsG3Ox oxidizes at 50,000-fold higher specificity (kcat/Km) the glucose moiety of mangiferin to 3-keto-mangiferin than free D-glucose to 2-keto-glucose. Analysis of PsG3Ox X-ray crystal structures and PsG3Ox in complex with glucose and mangiferin, combined with mutagenesis and molecular dynamics simulations, reveal distinctive features in the topology surrounding the active site that favor catalytically competent conformational states suitable for recognition, stabilization, and oxidation of the glucose moiety of mangiferin. Furthermore, their distinction to pyranose 2-oxidases (P2Oxs) involved in wood decay and recycling is discussed from an evolutionary, structural, and functional viewpoint.

Computational Enzyme Design at Zymvol.

Directed evolution is the most recognized methodology for enzyme engineering. The main drawback resides in its random nature and in the limited sequence exploration; both require screening of thousands (if not millions) of variants to achieve a target function. Computer-driven approaches can limit laboratorial screening to a few hundred candidates, enabling and accelerating the development of industrial enzymes. In this book chapter, the technology adopted at Zymvol is described. An overview of the current development and future directions in the company is also provided.

Distal Mutations Shape Substrate-Binding Sites during Evolution of a Metallo-Oxidase into a Laccase.

Laccases are in increasing demand as innovative solutions in the biorefinery fields. Here, we combine mutagenesis with structural, kinetic, and in silico analyses to characterize the molecular features that cause the evolution of a hyperthermostable metallo-oxidase from the multicopper oxidase family into a laccase (k cat 273 s-1 for a bulky aromatic substrate). We show that six mutations scattered across the enzyme collectively modulate dynamics to improve the binding and catalysis of a bulky aromatic substrate. The replacement of residues during the early stages of evolution is a stepping stone for altering the shape and size of substrate-binding sites. Binding sites are then fine-tuned through high-order epistasis interactions by inserting distal mutations during later stages of evolution. Allosterically coupled, long-range dynamic networks favor catalytically competent conformational states that are more suitable for recognizing and stabilizing the aromatic substrate. This work provides mechanistic insight into enzymatic and evolutionary molecular mechanisms and spots the importance of iterative experimental and computational analyses to understand local-to-global changes.

Unveiling molecular details behind improved activity at neutral to alkaline pH of an engineered DyP-type peroxidase.

DyP-type peroxidases (DyPs) are microbial enzymes that catalyze the oxidation of a wide range of substrates, including synthetic dyes, lignin-derived compounds, and metals, such as Mn2+ and Fe2+, and have enormous biotechnological potential in biorefineries. However, many questions on the molecular basis of enzyme function and stability remain unanswered. In this work, high-resolution structures of PpDyP wild-type and two engineered variants (6E10 and 29E4) generated by directed evolution were obtained. The X-ray crystal structures revealed the typical ferredoxin-like folds, with three heme access pathways, two tunnels, and one cavity, limited by three long loops including catalytic residues. Variant 6E10 displays significantly increased loops' flexibility that favors function over stability: despite the considerably higher catalytic efficiency, this variant shows poorer protein stability compared to wild-type and 29E4 variants. Constant-pH MD simulations revealed a more positively charged microenvironment near the heme pocket of variant 6E10, particularly in the neutral to alkaline pH range. This microenvironment affects enzyme activity by modulating the pK a of essential residues in the heme vicinity and should account for variant 6E10 improved activity at pH 7-8 compared to the wild-type and 29E4 that show optimal enzymatic activity close to pH 4. Our findings shed light on the structure-function relationships of DyPs at the molecular level, including their pH-dependent conformational plasticity. These are essential for understanding and engineering the catalytic properties of DyPs for future biotechnological applications.

Neutral and acidic hydrolysis reactions of the third generation anticancer drug oxaliplatin.

The hydrolysis of oxaliplatin, a third generation anticancer drug, is expected to play an important role in the activation of this compound before it reaches DNA. The first and second hydrolysis corresponding to the addition of the first water molecule concomitant with the ring-opening, followed by addition of a second water and loss of the monodentate oxalato ligand, respectively, were studied combining density functional theory (DFT) with the conductor-like dielectric continuum model (CPCM) approach. The reaction was studied in neutral and acidic conditions, and all stationary points have been identified. The computed potential energy surfaces show that, for the neutral hydrolysis, the ring-opening reaction is the rate-limiting process, with an activation barrier of about 28 kcal/mol. For the acid degradation in water, according to experimental data, the reaction is expected to proceed in a faster biphasic process, and the rate-limiting process is the ligand detachment that occurs with a barriers of about 22 kcal/mol. According to the calculated results, we expect that the reaction is favored in acidic conditions and that the monoaquated complex should be the species reacting with DNA.

Kauniolide synthase is a P450 with unusual hydroxylation and cyclization-elimination activity.

Guaianolides are an important class of sesquiterpene lactones with unique biological and pharmaceutical properties. They have been postulated to be derived from germacranolides, but for years no progress has been made in the elucidation of their biosynthesis that requires an unknown cyclization mechanism. Here we demonstrate the isolation and characterization of a cytochrome P450 from feverfew (Tanacetum parthenium), kauniolide synthase. Kauniolide synthase catalyses the formation of the guaianolide kauniolide from the germacranolide substrate costunolide. Unlike most cytochrome P450s, kauniolide synthase combines stereoselective hydroxylation of costunolide at the C3 position, with water elimination, cyclization and regioselective deprotonation. This unique mechanism of action is supported by in silico modelling and docking experiments. The full kauniolide biosynthesis pathway is reconstructed in the heterologous hosts Nicotiana benthamiana and yeast, paving the way for biotechnological production of guaianolide-type sesquiterpene lactones.

Simulating Substrate Recognition and Oxidation in Laccases: From Description to Design.

To meet the very specific requirements demanded by industry, proteins must be appropriately tailored. Engineering laccases, to improve the oxidation of small molecules, with applications in multiple fields, is, however, a difficult task. Most efforts have concentrated on increasing the redox potential of the enzyme, but in recent work, we have pursued an alternate strategy to engineering these biocatalysts. In particular, we have found that redesigning substrate binding at the T1 pocket, guided by in silico methodologies, to be a more consistent option. In this work, we evaluate the robustness of our computational approach to estimate activity, emphasizing the importance of the binding event in laccase reactivity. Strengths and weaknesses of the protocol are discussed along with its potential for scoring large numbers of protein sequences and thus its significance in protein engineering.

Oxidoreductases on their way to industrial biotransformations.

Fungi produce heme-containing peroxidases and peroxygenases, flavin-containing oxidases and dehydrogenases, and different copper-containing oxidoreductases involved in the biodegradation of lignin and other recalcitrant compounds. Heme peroxidases comprise the classical ligninolytic peroxidases and the new dye-decolorizing peroxidases, while heme peroxygenases belong to a still largely unexplored superfamily of heme-thiolate proteins. Nevertheless, basidiomycete unspecific peroxygenases have the highest biotechnological interest due to their ability to catalyze a variety of regio- and stereo-selective monooxygenation reactions with H2O2 as the source of oxygen and final electron acceptor. Flavo-oxidases are involved in both lignin and cellulose decay generating H2O2 that activates peroxidases and generates hydroxyl radical. The group of copper oxidoreductases also includes other H2O2 generating enzymes - copper-radical oxidases - together with classical laccases that are the oxidoreductases with the largest number of reported applications to date. However, the recently described lytic polysaccharide monooxygenases have attracted the highest attention among copper oxidoreductases, since they are capable of oxidatively breaking down crystalline cellulose, the disintegration of which is still a major bottleneck in lignocellulose biorefineries, along with lignin degradation. Interestingly, some flavin-containing dehydrogenases also play a key role in cellulose breakdown by directly/indirectly "fueling" electrons for polysaccharide monooxygenase activation. Many of the above oxidoreductases have been engineered, combining rational and computational design with directed evolution, to attain the selectivity, catalytic efficiency and stability properties required for their industrial utilization. Indeed, using ad hoc software and current computational capabilities, it is now possible to predict substrate access to the active site in biophysical simulations, and electron transfer efficiency in biochemical simulations, reducing in orders of magnitude the time of experimental work in oxidoreductase screening and engineering. What has been set out above is illustrated by a series of remarkable oxyfunctionalization and oxidation reactions developed in the frame of an intersectorial and multidisciplinary European RTD project. The optimized reactions include enzymatic synthesis of 1-naphthol, 25-hydroxyvitamin D3, drug metabolites, furandicarboxylic acid, indigo and other dyes, and conductive polyaniline, terminal oxygenation of alkanes, biomass delignification and lignin oxidation, among others. These successful case stories demonstrate the unexploited potential of oxidoreductases in medium and large-scale biotransformations.

New Monte Carlo Based Technique To Study DNA-Ligand Interactions.

We present a new all-atom Monte Carlo technique capable of performing quick and accurate DNA-ligand conformational sampling. In particular, and using the PELE software as a frame, we have introduced an additional force field, an implicit solvent, and an anisotropic network model to effectively map the DNA energy landscape. With these additions, we successfully generated DNA conformations for a test set composed of six DNA fragments of A-DNA and B-DNA. Moreover, trajectories generated for cisplatin and its hydrolysis products identified the best interacting compound and binding site, producing analogous results to microsecond molecular dynamics simulations. Furthermore, a combination of the Monte Carlo trajectories with Markov State Models produced noncovalent binding free energies in good agreement with the published molecular dynamics results, at a significantly lower computational cost. Overall our approach will allow a quick but accurate sampling of DNA-ligand interactions.

Atomic picture of ligand migration in toluene 4-monooxygenase.

Computational modeling combined with mutational and activity assays was used to underline the substrate migration pathways in toluene 4-monooxygenase, a member of the important family of bacterial multicomponent monooxygenases (BMMs). In all structurally defined BMM hydroxylases, several hydrophobic cavities in the α-subunit map a preserved path from the protein surface to the diiron active site. Our results confirm the presence of two pathways by which different aromatic molecules can enter/escape the active site. While the substrate is observed to enter from both channels, the more hydrophilic product is withdrawn mainly from the shorter channel ending at residues D285 and E214. The long channel ends in the vicinity of S395, whose variants have been seen to affect activity and specificity. These mutational effects are clearly reproduced and rationalized by the in silico studies. Furthermore, the combined computational and experimental results highlight the importance of residue F269, which is located at the intersection of the two channels.