Prediction of the solvent affecting site and the computational design of stable Candida antarctica lipase B in a hydrophilic organic solvent

Enhanced Activity and Stability of an Acetyl Xylan Esterase in Hydrophilic Alcohols through Site-Directed Mutagenesis

Current demands for the development of suitable biocatalysts showing high process performance is stimulated by the need to replace current chemical synthesis with cleaner alternatives. A drawback to the use of biocatalysts for unique applications is their low performance in industrial conditions. Hence, enzymes with improved performance are needed to achieve innovative and sustainable biocatalysis. In this study, we report the improved performance of an engineered acetyl xylan esterase (BaAXE) in a hydrophilic organic solvent. The structure of BaAXE was partitioned into a substrate-binding region and a solvent-affecting region. Using a rational design approach, charged residues were introduced at protein surfaces in the solvent-affecting region. Two sites present in the solvent-affecting region, A12D and Q143E, were selected for site-directed mutagenesis, which generated the mutants MUT12, MUT143 and MUT12-143. The mutants MUT12 and MUT143 reported lower Km (0.29 mM and 0.27 mM, respectively) compared to the wildtype (0.41 mM). The performance of the mutants in organic solvents was assessed after enzyme incubation in various strengths of alcohols. The mutants showed improved activity and stability compared to the wild type in low strengths of ethanol and methanol. However, the activity of MUT143 was lost in 40% methanol while MUT12 and MUT12-143 retained over 70% residual activity in this environment. Computational analysis links the improved performance of MUT12 and MUT12-143 to novel intermolecular interactions that are absent in MUT143. This work supports the rationale for protein engineering to augment the characteristics of wild-type proteins and provides more insight into the role of charged residues in conferring stability.

Solvent tolerant enzymes in extremophiles: Adaptations and applications

Non-aqueous enzymology has always drawn attention due to the wide range of unique possibilities in biocatalysis. In general, the enzymes do not or insignificantly catalyze substrate in the presence of solvents. This is due to the interfering interactions of the solvents between enzyme and water molecules at the interface. Therefore, information about solvent-stable enzymes is scarce. Yet, solvent-stable enzymes prove quite valuable in the present day biotechnology. The enzymatic hydrolysis of the substrates in solvents synthesizes commercially valuable products, such as peptides, esters, and other transesterification products. Extremophiles, the most valuable yet not extensively explored candidates, can be an excellent source to investigate this avenue. Due to inherent structural attributes, many extremozymes can catalyze and maintain stability in organic solvents. In the present review, we aim to consolidate information about the solvent-stable enzymes from various extremophilic microorganisms. Further, it would be interesting to learn about the mechanism adapted by these microorganisms to sustain solvent stress. Various approaches to protein engineering are used to enhance catalytic flexibility and stability and broaden biocatalysis's prospects under non-aqueous conditions. It also describes strategies to achieve optimal immobilization with minimum inhibition of the catalysis. The proposed review would significantly aid our understanding of non-aqueous enzymology.

Improving ethanol tolerance of ethyl carbamate hydrolase by diphasic high pressure molecular dynamic simulations

Ethyl carbamate (EC) is mainly found in fermented foods and fermented alcoholic beverages, which could cause carcinogenic potential to humans. Reducing EC is one of the key research priorities to address security of fermented foods. Enzymatic degradation of EC with EC hydrolase in food is the most reliable and efficient method. However, poor tolerance to ethanol severely hinders application of EC hydrolase. In this study, the mutants of EC hydrolase were screened by diphasic high pressure molecular dynamic simulations (dHP-MD). The best variant with remarkable improvement in specific activity and was H68A/K70R/S325N, whose specific activity was approximately 3.42-fold higher than WT, and relative enzyme activity under 20% (v/v) was 5.02-fold higher than WT. Moreover, the triple mutant increased its stability by acquiring more hydration shell and forming extra hydrogen bonds. Furthermore, the ability of degrading EC of the immobilized triple mutant was both detected in mock wine and under certain reaction conditions. The stability of immobilized triple mutant and WT were both improved, and immobilized triple mutant degraded nearly twice as much EC as that of immobilized WT. Overall, dHP-MD was proved to effectively improve enzyme activity and ethanol tolerance for extent application at industrial scale.

Engineering a Bacillus subtilis esterase for selective hydrolysis of d, l-menthyl acetate in an organic solvent-free system†

Esterase/lipase-catalyzed selective hydrolysis of d, l-menthyl esters has become one of the promising approaches for producing l-menthol, one of the most important flavoring chemicals with extensive uses. However, the activity and l-enantioselectivity of the biocatalyst are not sufficient for meeting the industrial requirements. Herein, a highly active para-nitrobenzyl esterase from Bacillus subtilis 168 (pnbA-BS) was cloned and then engineered to enhance its l-enantioselectivity. On the basis of the strategy tailoring the steric exclusion effect and structural flexibility of the region adjacent to the substrate, the substitution of Ala400 to Pro caused a remarkable improvement in the E value from 1.0 to 466.6. The variant A400P was purified and further confirmed with strict l-enantioselectivity in the selective hydrolysis of d, l-menthyl acetate, whereas the improved l-enantioselectivity caused decreased activity. To develop an efficient, easy-to-use, and green methodology, organic solvent was omitted and substrate constant feeding was integrated into the whole-cell catalyzed system. During the catalytic process, the selective hydrolysis of 1.0 M d, l-menthyl acetate in 14 h offered a conversion of 48.9%, e.e._p value of >99%, and space-time yield of 160.52 g (l d)⁻¹.

Esterase/lipase-catalyzed selective hydrolysis of d, l-menthyl esters has become one of the promising approaches for producing l-menthol, one of the most important flavoring chemicals with extensive uses.

Single-sequence protein structure prediction using a language model and deep learning

AlphaFold2 and related computational systems predict protein structure using deep learning and co-evolutionary relationships encoded in multiple sequence alignments (MSAs). Despite high prediction accuracy achieved by these systems, challenges remain in (1) prediction of orphan and rapidly evolving proteins for which an MSA cannot be generated; (2) rapid exploration of designed structures; and (3) understanding the rules governing spontaneous polypeptide folding in solution. Here we report development of an end-to-end differentiable recurrent geometric network (RGN) that uses a protein language model (AminoBERT) to learn latent structural information from unaligned proteins. A linked geometric module compactly represents Cα backbone geometry in a translationally and rotationally invariant way. On average, RGN2 outperforms AlphaFold2 and RoseTTAFold on orphan proteins and classes of designed proteins while achieving up to a 106-fold reduction in compute time. These findings demonstrate the practical and theoretical strengths of protein language models relative to MSAs in structure prediction.

Using Molecular Simulation to Guide Protein Engineering for Biocatalysis in Organic Solvents

Biocatalysis in organic solvents (OSs) is very appealing for the industry in producing bulk and/or fine chemicals, such as pharmaceuticals, biodiesel, and fragrances. The poor performance of enzymes in OSs (e.g., reduced activity, insufficient stability, and deactivation) negates OSs' excellent solvent properties. Molecular dynamics (MD) simulations provide a complementary method to study the relationship between enzymes dynamics and the stability in OSs. Here we describe computational procedure for MD simulation of enzymes in OSs with an example of Bacillus subtilis lipase A (BSLA) in dimethyl sulfoxide (DMSO) cosolvent with software GROMACS. We discuss main essential practical issues considered (such as choice of force field, parameterization, simulation setup, and trajectory analysis). The core part of this protocol (enzyme-OS system setup, analysis of structural-based and solvation-based observables) is transferable to other enzymes and any OS systems. Combining with experimental studies, the obtained molecular knowledge is most likely to guide researchers to access rational protein engineering approaches to tailor OS resistant enzymes and expand the scope of biocatalysis in OS media. Finally, we discuss potential solutions to overcome the remaining challenges of computational biocatalysis in OSs and briefly draw future directions for further improvement in this field.

Improving the organic solvent resistance of lipase a from Bacillus subtilis in water-ethanol solvent through rational surface engineering

Given that lipase is an enzyme applicable in various industrial fields and water-miscible organic solvents are important reaction media for developing industrial-scale biocatalysis, a structure-based strategy was explored to stabilize lipase A from Bacillus subtilis in a water-ethanol cosolvent. Site-directed mutagenesis of ethanol-interacting sites resulted in 4 mutants, i.e., Ser16Gly, Ala38Gly, Ala38Thr, and Leu108Asn, which were stable in 50% ethanol and had up to 1.8-fold higher stability than the wild-type. In addition, Leu108Asn was more thermostable at 45 °C than the wild type. The results discussed in this study not only provide insights into strategies for enzyme engineering to improve organic solvent resistance but also suggest perspectives on pioneering routes for constructing enzyme-based biorefineries to produce value-added fuels and chemicals.

CompassR Yields Highly Organic-Solvent-Tolerant Enzymes through Recombination of Compatible Substitutions

The CompassR (computer-assisted recombination) rule enables, among beneficial substitutions, the identification of those that can be recombined in directed evolution. Herein, a recombination strategy is systematically investigated to minimize experimental efforts and maximize possible improvements. In total, 15 beneficial substitutions from Bacillus subtilis lipase A (BSLA), which improves resistance to the organic cosolvent 1,4-dioxane (DOX), were studied to compare two recombination strategies, the two-gene recombination process (2GenReP) and the in silico guided recombination process (InSiReP), employing CompassR. Remarkably, both strategies yielded a highly DOX-resistant variant, M4 (I12R/Y49R/E65H/N98R/K122E/L124K), with up to 14.6-fold improvement after screening of about 270 clones. M4 has a remarkably enhanced resistance in 60 % (v/v) acetone (6.0-fold), 30 % (v/v) ethanol (2.1-fold), and 60 % (v/v) methanol (2.4-fold) compared with wild-type BSLA. Molecular dynamics simulations revealed that attracting water molecules by charged surface substitutions is the main driver for increasing the DOX resistance of BSLA M4. Both strategies and obtained molecular knowledge can likely be used to improve the properties of other enzymes with a similar α/β-hydrolase fold.

A novel heptapeptide derived from Crassostrea gigas shows anticoagulant activity by targeting for thrombin active domain

A novel food-derived anticoagulant heptapeptides (P-3-CG) was isolated and characterized from oyster (Crassostrea gigas) pepsin hydrolysate. P-3-CG competed with fibrinogen against thrombin active domain by a spontaneous and exothermic reaction which was entropically driven. The residue Lys7 of P-3-CG anchored thrombin S₁ pocket strongly, which inhibited fibrinogen binding to the thrombin, then blocked the conversion of fibrinogen to fibrin. The fibrinogen clotting time was prolonged to 27.55 s, and the reciprocally authenticated results of dynamic light scattering and scanning electron microscope further explained for fibrinogen clotting time extension. Inhibition of amidolytic activity of thrombin was affected significantly by reaction time and P-3-CG concentration. Furthermore, P-3-CG prolonged activated partial thromboplastin time significantly in vitro/vivo, and decreased the mortality which was confirmed by pulmonary pathological slide results. The obtained results demonstrated that P-3-CG may potentially serve as an alternative food-derived anticoagulant peptide that could be utilized for thrombosis prevention.

Characterization of a Novel Moderately Thermophilic Solvent-Tolerant Esterase Isolated From a Compost Metagenome Library

A novel esterase, EstCS1, was isolated from a compost metagenomics library. The EstCS1 protein, which consists of 309 amino acid residues with an anticipated molecular mass of 34 kDa, showed high amino acid sequence identities to predicted esterases and alpha/beta hydrolases (59%) from some cultured bacteria and to predicted lipases/esterases from uncultured bacteria. The phylogenetic analysis suggested that the EstCS1 belongs to the hormone-sensitive lipase family of lipolytic enzyme classification and contains a catalytic triad including Ser155–Asp255–His285. The Ser155 residue of the catalytic triad in the EstCS1 was located in the consensus active-site motif, GXSXG. Besides, a conserved HGGG motif placed in an oxyanion hole of the hormone-sensitive lipase family was discovered, too. The EstCS1 demonstrated the highest activity toward p-nitrophenyl propionate (C3) and caproate (C6) and was normally stable up to 60°C with optimal activity at 50°C. In addition, an optimal activity was observed at pH 8, and the EstCS1 possessed its stability within the pH range between 5 and 10. Interestingly, EstCS1 had an outstanding stability in up to 30% (v/v) organic solvents and activity over 50% in the presence of 50% (v/v) acetone, ethanol, dimethyl sulfoxide (DMSO), and N,N-dimethylformamide. The EstCS1 hydrolyzed sterically hindered tertiary alcohol esters of t-butyl acetate and linalyl acetate. Considering the properties, such as the moderate thermostability, stability against organic solvents, and activity toward esters of tertiary alcohols, the EstCS1 will be worthwhile to be used for organic synthesis and related industrial applications.

Algorithm-based coevolution network identification reveals key functional residues of the α/β hydrolase subfamilies

Covariant residues identified by computational algorithms have provided new insights into enzyme evolutionary routes. However, the reliability and accuracy of routine statistical coupling analysis (SCA) are unable to satisfy the needs of protein engineering because SCA depends only on sequence information. Here, we set up a new SCA algorithm, SCA.SIM, by integrating structure information and MD simulation data. The more reliable covariant residues with high-quality scores are obtained from sequence alignment weighted by residual movement for eight related subfamilies, belonging to α/β hydrolase family, with Candida antarctica lipase B (CALB). The 38 predicted covariant residues are tested for function by high-throughput quantitative evaluation in combination with activity and thermostability assays of a mutant library and deep sequencing. Based on the landscapes of both activity and thermostability, most mutants play key roles in catalysis, and some mutants gain 2.4- to 6-fold increase in half-life at 50°C and 9- to 12-fold improvement in catalytic efficiency. The activity of double mutants for A225F/T103A is higher than those of A225F and T103A which means that SCA.SIM method might be useful for identifying the allosteric coupling. The SCA.SIM algorithm can be used for protein coevolution and enzyme engineering research.

Effect of organic solvent on enzymatic degradation of cyclic PBS-based polymers by lipase N435

Poly(butylene succinate-co-cyclohexane dimethanol succinate) (P(BS-co-CHDMS)) and poly(butylene succinate-co-butanediol cyclohexanedicarboxylic acid) (P(BS-co-BCHDA)) were catalytically degraded by Candida antarctica lipase Novozyme 435 (N435) in CHCl₃ and THF. The results indicated that the degradation rate was P(BS-co-CHDMS) > P(BS-co-BCHDA) > poly(butylene succinate) (PBS). The degradation rate of copolyesters was higher in CHCl₃ than in THF, the highest degradation rate of 67% being obtained for P(BS-co-CHDMS). Hence, the CHCl₃ solvent is more suitable for the enzyme-catalytic degradation of copolyesters, since the lipase can easier recognize the butylene succinate (BS-), (butanediol cyclohexanedicarboxylic acid) (BCA-), and (cyclohexane dimethanol succinate-type) (CMS-type) ester bonds in this solvent. Moreover, it can recognize the CMS-type ester bonds with a higher specificity than the (butanediol cyclohexanedicarboxylic acid type) (BCA-type) ester bonds. Molecular simulation results indicated that the structure of the lipase was stable in CHCl₃ and THF. However, CHCl₃ proved to be more suitable for a stable activity of the enzyme. The active pocket contains acyl-binding hydrophilic residues which are recognized by the substrate. The increase in the content of saturated cycles can increase the hydrophobicity of the substrate and thus, the amount of substrate bond to enzyme active site is increased, which facilitates the enzymatic degradation of copolyesters.

Changes of Thermostability, Organic Solvent, and pH Stability in Geobacillus zalihae HT1 and Its Mutant by Calcium Ion

Thermostable T1 lipase from Geobacillus zalihae has been crystallized using counter-diffusion method under space and Earth conditions. The comparison of the three-dimensional structures from both crystallized proteins show differences in the formation of hydrogen bond and ion interactions. Hydrogen bond and ion interaction are important in the stabilization of protein structure towards extreme temperature and organic solvents. In this study, the differences of hydrogen bond interactions at position Asp43, Thr118, Glu250, and Asn304 and ion interaction at position Glu226 was chosen to imitate space-grown crystal structure, and the impact of these combined interactions in T1 lipase-mutated structure was studied. Using space-grown T1 lipase structure as a reference, subsequent simultaneous mutation D43E, T118N, E226D, E250L, and N304E was performed on recombinant wild-type T1 lipase (wt-HT1) to generate a quintuple mutant term as 5M mutant lipase. This mutant lipase shared similar characteristics to its wild-type in terms of optimal pH and temperature. The stability of mutant 5M lipase improved significantly in acidic and alkaline pH as compared to wt-HT1. 5M lipase was highly stable in organic solvents such as dimethyl sulfoxide (DMSO), methanol, and n-hexane compared to wt-HT1. Both wild-type and mutant lipases were found highly activated in calcium as compared to other metal ions due to the presence of calcium-binding site for thermostability. The presence of calcium prolonged the half-life of mutant 5M and wt-HT1, and at the same time increased their melting temperature (T_m). The melting temperature of 5M and wt-HT1 lipases increased at 8.4 and 12.1 °C, respectively, in the presence of calcium as compared to those without. Calcium enhanced the stability of mutant 5M in 25% (v/v) DMSO, n-hexane, and n-heptane. The lipase activity of wt-HT1 also increased in 25% (v/v) ethanol, methanol, acetonitrile, n-hexane, and n-heptane in the presence of calcium. The current study showed that the accumulation of amino acid substitutions D43E, T118N, E226D, E250L, and N304E produced highly stable T1 mutant when hydrolyzing oil in selected organic solvents such as DMSO, n-hexane, and n-heptane. It is also believed that calcium ion plays important role in regulating lipase thermostability.

Theoretical and experimental studies on the conformational changes of organic solvent-stable protease from Bacillus sphaericus DS11 in methanol/water mixtures

If natural proteases are used in organic synthesis, they are often inactivated or give a low rate of reaction in non-aqueous or aqueous-organic media. Therefore, to reveal the molecular mechanism governing the stability of proteases in organic solvents and increase protease stability in those systems is of intriguing interest. In the present study, the activity and conformational changes of an organic solvent-stable protease (OSP) from Bacillus sphaericus DS11 in different concentrations of methanol were investigated by measuring fluorescence, UV-Vis spectra, circular dichroism (CD), and conducting molecular dynamics (MD) simulations. The OSP expanded with increasing methanol concentration. The methanol molecules were able to enter into the OSP, leading to microenvironmental changes around the aromatic amino acids. More hydrophobic groups were exposed to the solvents at high methanol concentrations, and the original hydrophobic interaction in the protein decreased, thus resulting in the secondary and tertiary structure change in the OSP. Our results provide helpful insight into the molecular mechanism of the OSP tolerance to organic solvent and indicate directions for future work to design and engineer proteases that are stable at high organic solvent concentrations.

Solvent stable microbial lipases: current understanding and biotechnological applications

OBJECTIVE

This review examines on our current understanding of microbial lipase solvent tolerance, with a specific focus on the molecular strategies employed to improve lipase stability in a non-aqueous environment.

RESULTS

It provides an overview of known solvent tolerant lipases and of approaches to improving solvent stability such as; enhancing stabilising interactions, modification of residue flexibility and surface charge alteration. It shows that judicious selection of lipase source supplemented by appropriate enzyme stabilisation, can lead to a wide application spectrum for lipases.

CONCLUSION

Organic solvent stable lipases are, and will continue to be, versatile and adaptable biocatalytic workhorses commonly employed for industrial applications in the food, pharmaceutical and green manufacturing industries.

pH-Induced Rotation of Lidless Lipase LipA from Bacillus subtilis at Lipase-Detergent Interface

Lipases exhibit a unique process during the catalysis of the hydrolysis of triglyceride substrates called interfacial activation. Surfactants are used as cosolvents with water not only to offer a less polar environment to the lipases needed for their interfacial activation but also to solvate the substrate which are poorly soluble in water. However, the presence of detergent in the medium can affect both the lipase and the substrate, making the construction of a microkinetic model for lipase activity in the presence of the detergent difficult. Herein, we study the interfacial activation of a lidless lipase LipA from Bacillus subtilis using extensive atomistic molecular dynamics simulations at different concentrations of the surfactant, Thesit (C12E8), at two pH values. Residues which bind to the monomeric detergent are found to be the same as the ones which have been reported earlier to bind to the substrate. Very importantly, a pH-induced rotation of the enzyme with respect to surfactant aggregate has been observed which not only explains the experimentally observed pH-dependent enzymatic activity of this lidless lipase, but also suggests its reorientation at an aqueous-lipodophilic interface.

Lid opening and conformational stability of T1 Lipase is mediated by increasing chain length polar solvents

The dynamics and conformational landscape of proteins in organic solvents are events of potential interest in nonaqueous process catalysis. Conformational changes, folding transitions, and stability often correspond to structural rearrangements that alter contacts between solvent molecules and amino acid residues. However, in nonaqueous enzymology, organic solvents limit stability and further application of proteins. In the present study, molecular dynamics (MD) of a thermostable Geobacillus zalihae T1 lipase was performed in different chain length polar organic solvents (methanol, ethanol, propanol, butanol, and pentanol) and water mixture systems to a concentration of 50%. On the basis of the MD results, the structural deviations of the backbone atoms elucidated the dynamic effects of water/organic solvent mixtures on the equilibrium state of the protein simulations in decreasing solvent polarity. The results show that the solvent mixture gives rise to deviations in enzyme structure from the native one simulated in water. The drop in the flexibility in H₂O, MtOH, EtOH and PrOH simulation mixtures shows that greater motions of residues were influenced in BtOH and PtOH simulation mixtures. Comparing the root mean square fluctuations value with the accessible solvent area (SASA) for every residue showed an almost correspondingly high SASA value of residues to high flexibility and low SASA value to low flexibility. The study further revealed that the organic solvents influenced the formation of more hydrogen bonds in MtOH, EtOH and PrOH and thus, it is assumed that increased intraprotein hydrogen bonding is ultimately correlated to the stability of the protein. However, the solvent accessibility analysis showed that in all solvent systems, hydrophobic residues were exposed and polar residues tended to be buried away from the solvent. Distance variation of the tetrahedral intermediate packing of the active pocket was not conserved in organic solvent systems, which could lead to weaknesses in the catalytic H-bond network and most likely a drop in catalytic activity. The conformational variation of the lid domain caused by the solvent molecules influenced its gradual opening. Formation of additional hydrogen bonds and hydrophobic interactions indicates that the contribution of the cooperative network of interactions could retain the stability of the protein in some solvent systems. Time-correlated atomic motions were used to characterize the correlations between the motions of the atoms from atomic coordinates. The resulting cross-correlation map revealed that the organic solvent mixtures performed functional, concerted, correlated motions in regions of residues of the lid domain to other residues. These observations suggest that varying lengths of polar organic solvents play a significant role in introducing dynamic conformational diversity in proteins in a decreasing order of polarity.

Modulating Mobility: a Paradigm for Protein Engineering?

Proteins are highly mobile structures. In addition to gross conformational changes occurring on, for example, ligand binding, they are also subject to constant thermal motion. The mobility of a protein varies through its structure and can be modulated by ligand binding and other events. It is becoming increasingly clear that this mobility plays an important role in key functions of proteins including catalysis, allostery, cooperativity, and regulation. Thus, in addition to an optimum structure, proteins most likely also require an optimal dynamic state. Alteration of this dynamic state through protein engineering will affect protein function. A dramatic example of this is seen in some inherited metabolic diseases where alternation of residues distant from the active site affects the mobility of the protein and impairs function. We postulate that using molecular dynamics simulations, experimental data or a combination of the two, it should be possible to engineer the mobility of active sites. This may be useful in, for example, increasing the promiscuity of enzymes. Thus, a paradigm for protein engineering is suggested in which the mobility of the active site is rationally modified. This might be combined with more "traditional" approaches such as altering functional groups in the active site.

Geometric Simulation Approach for Grading and Assessing the Thermostability of CALBs

Candida antarctica lipase B (CALB) is a known stable and highly active enzyme used widely in biodiesel synthesis. In this work, the stability of native (4K6G) and mutant (4K5Q) CALB was studied through various structural parameters using conformational sampling approach. The contours of polar surface area and surface area of mutant CALB were 11357.67 Å(2) and 30007.4 Å(2), respectively, showing an enhanced stability compared to native CALB with a statistically significant P value of < 0.0001. Moreover, simulated thermal denaturation of CALB, a process involving dilution of hydrogen bond, significantly shielded against different intervals of energy application in mutant CALB revealing its augmentation of structural rigidity against native CALB. Finally, computational docking analysis showed an increase in the binding affinity of CALB and its substrate (triglyceride) in mutant CALB with Atomic Contact Energy (ACE) of -91.23 kcal/mol compared to native CALB (ACE of -70.3 kcal/mol). The computational observations proposed that the use of mutant CALB (4K5Q) could serve as a best template for production of biodiesel in the future. Additionally, it can also be used as a template to identify efficient thermostable lipases through further mutations.

An investigation of nitrile transforming enzymes in the chemo-enzymatic synthesis of the taxol sidechain

Paclitaxel (taxol) is an antimicrotubule agent widely used in the treatment of cancer. Taxol is prepared in a semisynthetic route by coupling the N-benzoyl-(2R,3S)-3-phenylisoserine sidechain to the baccatin III core structure. Precursors of the taxol sidechain have previously been prepared in chemoenzymatic approaches using acylases, lipases, and reductases, mostly featuring the enantioselective, enzymatic step early in the reaction pathway. Here, nitrile hydrolysing enzymes, namely nitrile hydratases and nitrilases, are investigated for the enzymatic hydrolysis of two different sidechain precursors. Both sidechain precursors, an openchain α-hydroxy-β-amino nitrile and a cyanodihydrooxazole, are suitable for coupling to baccatin III directly after the enzymatic step. An extensive set of nitrilases and nitrile hydratases was screened towards their activity and selectivity in the hydrolysis of two taxol sidechain precursors and their epimers. A number of nitrilases and nitrile hydratases converted both sidechain precursors and their epimers.

Single-Point Mutation with a Rotamer Library Toolkit: Toward Protein Engineering

Protein engineers have long been hard at work to harness biocatalysts as a natural source of regio-, stereo-, and chemoselectivity in order to carry out chemistry (reactions and/or substrates) not previously achieved with these enzymes. The extreme labor demands and exponential number of mutation combinations have induced computational advances in this domain. The first step in our virtual approach is to predict the correct conformations upon mutation of residues (i.e., rebuilding side chains). For this purpose, we opted for a combination of molecular mechanics and statistical data. In this work, we have developed automated computational tools to extract protein structural information and created conformational libraries for each amino acid dependent on a variable number of parameters (e.g., resolution, flexibility, secondary structure). We have also developed the necessary tool to apply the mutation and optimize the conformation accordingly. For side-chain conformation prediction, we obtained overall average root-mean-square deviations (RMSDs) of 0.91 and 1.01 Å for the 18 flexible natural amino acids within two distinct sets of over 3000 and 1500 side-chain residues, respectively. The commonly used dihedral angle differences were also evaluated and performed worse than the state of the art. These two metrics are also compared. Furthermore, we generated a family-specific library for kinases that produced an average 2% lower RMSD upon side-chain reconstruction and a residue-specific library that yielded a 17% improvement. Ultimately, since our protein engineering outlook involves using our docking software, Fitted/Impacts, we applied our mutation protocol to a benchmarked data set for self- and cross-docking. Our side-chain reconstruction does not hinder our docking software, demonstrating differences in pose prediction accuracy of approximately 2% (RMSD cutoff metric) for a set of over 200 protein/ligand structures. Similarly, when docking to a set of over 100 kinases, side-chain reconstruction (using both general and biased conformation libraries) had minimal detriment to the docking accuracy.