Computational design of a thermostable mutant of cocaine esterase via molecular dynamics simulations

Critical assessment of structure-based approaches to improve protein resistance in aqueous ionic liquids by enzyme-wide saturation mutagenesis

Ionic liquids (IL) and aqueous ionic liquids (aIL) are attractive (co-)solvents for green industrial processes involving biocatalysts, but often reduce enzyme activity. Experimental and computational methods are applied to predict favorable substitution sites and, most often, subsequent site-directed surface charge modifications are introduced to enhance enzyme resistance towards aIL. However, almost no studies evaluate the prediction precision with random mutagenesis or the application of simple data-driven filtering processes. Here, we systematically and rigorously evaluated the performance of 22 previously described structure-based approaches to increase enzyme resistance to aIL based on an experimental complete site-saturation mutagenesis library of Bacillus subtilis Lipase A (BsLipA) screened against four aIL. We show that, surprisingly, most of the approaches yield low gain-in-precision (GiP) values, particularly for predicting relevant positions: 14 approaches perform worse than random mutagenesis. Encouragingly, exploiting experimental information on the thermostability of BsLipA or structural weak spots of BsLipA predicted by rigidity theory yields GiP = 3.03 and 2.39 for relevant variants and GiP = 1.61 and 1.41 for relevant positions. Combining five simple-to-compute physicochemical and evolutionary properties substantially increases the precision of predicting relevant variants and positions, yielding GiP = 3.35 and 1.29. Finally, combining these properties with predictions of structural weak spots identified by rigidity theory additionally improves GiP for relevant variants up to 4-fold to ∼10 and sustains or increases GiP for relevant positions, resulting in a prediction precision of ∼90% compared to ∼9% in random mutagenesis. This combination should be applicable to other enzyme systems for guiding protein engineering approaches towards improved aIL resistance.

Systematically Scrutinizing the Impact of Substitution Sites on Thermostability and Detergent Tolerance for Bacillus subtilis Lipase A

Improving an enzyme's (thermo-)stability or tolerance against solvents and detergents is highly relevant in protein engineering and biotechnology. Recent developments have tended toward data-driven approaches, where available knowledge about the protein is used to identify substitution sites with high potential to yield protein variants with improved stability, and subsequently, substitutions are engineered by site-directed or site-saturation (SSM) mutagenesis. However, the development and validation of algorithms for data-driven approaches have been hampered by the lack of availability of large-scale data measured in a uniform way and being unbiased with respect to substitution types and locations. Here, we extend our knowledge on guidelines for protein engineering following a data-driven approach by scrutinizing the impact of substitution sites on thermostability or/and detergent tolerance for Bacillus subtilis lipase A (BsLipA) at very large scale. We systematically analyze a complete experimental SSM library of BsLipA containing all 3439 possible single variants, which was evaluated as to thermostability and tolerances against four detergents under respectively uniform conditions. Our results provide systematic and unbiased reference data at unprecedented scale for a biotechnologically important protein, identify consistently defined hot spot types for evaluating the performance of data-driven protein-engineering approaches, and show that the rigidity theory and ensemble-based approach Constraint Network Analysis yields hot spot predictions with an up to ninefold gain in precision over random classification.

Protein thermostability engineering

Using structure and sequence based analysis we can engineer proteins to increase their thermal stability.

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.

Lysine-Based Site-Directed Mutagenesis Increased Rigid β-Sheet Structure and Thermostability of Mesophilic 1,3-1,4-β-Glucanase

1,3-1,4-β-Glucanase is widely applied in the food industry, while its low thermostability often reduces its performance. In a previous study, chemical modification of surface lysine residues was proved to increase the thermostability of β-glucanase. To improve the thermostability, the mesophilic β-glucanase from Bacillus terquilensis was rationally engineered through site-directed mutagenesis of the 12 lysines into serines. The results showed that the K20S, K117S, and K165S mutants could both enhance the specific activities and thermostability of β-glucanase. The triple mutant (K20S/K117S/K165S) could increase the optimal temperature and T50 value by 15 and 14 °C, respectively. Five percent more structured residues were observed in the mutant, which formed new β-sheet structures in the concave side. Molecular dynamics simulation analysis showed that the flexibility in the mutation regions was decreased, which resulted in the overall rigidity of the β-glucanase. Therefore, the lysine-based site-directed mutagenesis is a simple and effective method for improving the thermostability of β-glucanase.

A computational strategy for altering an enzyme in its cofactor preference to NAD(H) and/or NADP(H)

Coenzyme engineering, especially for altered coenzyme specificity, has been a research hotspot for more than a decade. In the present study, a novel computational strategy that enhances the hydrogen-bond interaction between an enzyme and a coenzyme was developed and utilized to alter the coenzyme preference. This novel computational strategy only required the structure of the target enzyme. No other homologous enzymes were needed to achieve alteration in the coenzyme preference of a certain enzyme. Using our novel strategy, Gox2181 was reconstructed from exhibiting complete NADPH preference to exhibiting dual cofactor specificity for NADH and NADPH. Structure-guided Gox2181 mutants were designed in silico and molecular dynamics simulations were performed to evaluate the strength of hydrogen-bond interactions between the enzyme and the coenzyme NADPH. Three Gox2181 mutants displaying high structure stability and structural compatibility to NADH/NADPH were chosen for experimental confirmation. Among the three Gox2181 mutants, Gox2181-Q20R&D43S showed the highest enzymatic activity by utilizing NADPH as its coenzyme, which was even better than the wild-type enzyme. In addition, isothermal titration calorimetry analysis further verified that Gox2181-Q20R&D43S was able to interact with NADPH but the wild-type enzyme could not. This novel computational strategy represents an insightful approach for altering the cofactor preference of target enzymes.

Enhancement of proteolytic activity of a thermostable papain-like protease by structure-based rational design

Ervatamins (A, B and C) are papain-like cysteine proteases from the plant Ervatamia coronaria. Among Ervatamins, Ervatamin-C is a thermostable protease, but it shows lower catalytic efficiency. In contrast, Ervatamin-A which has a high amino acid sequence identity (∼90%) and structural homology (Cα rmsd 0.4 Å) with Ervatamin-C, has much higher catalytic efficiency (∼57 times). From the structural comparison of Ervatamin-A and -C, two residues Thr32 and Tyr67 in the catalytic cleft of Ervatamin-A have been identified whose contributions for higher activity of Ervatamin-A are established in our earlier studies. In this study, these two residues have been introduced in Ervatamin-C by site directed mutagenesis to enhance the catalytic efficiency of the thermostable protease. Two single mutants (S32T and A67Y) and one double mutant (S32T/A67Y) of Ervatamin-C have been generated and characterized. All the three mutants show ∼ 8 times higher catalytic efficiency (k cat/K m) than the wild-type. The thermostability of all the three mutant enzymes remained unchanged. The double mutant does not achieve the catalytic efficiency of the template enzyme Ervatamin-A. By modeling the structure of the double mutant and probing the role of active site residues by docking a substrate, the mechanistic insights of higher activity of the mutant protease have been addressed. The in-silico study demonstrates that the residues beyond the catalytic cleft also influence the substrate binding and positioning of the substrate at the catalytic centre, thus controlling the catalytic efficiency of an enzyme.

Plants as a source of butyrylcholinesterase variants designed for enhanced cocaine hydrolase activity

Cocaine addiction affects millions of people with disastrous personal and social consequences. Cocaine is one of the most reinforcing of all drugs of abuse, and even those who undergo rehabilitation and experience long periods of abstinence have more than 80% chance of relapse. Yet there is no FDA-approved treatment to decrease the likelihood of relapse in rehabilitated addicts. Recent studies, however, have demonstrated a promising potential treatment option with the help of the serum enzyme butyrylcholinesterase (BChE), which is capable of breaking down naturally occurring (-)-cocaine before the drug can influence the reward centers of the brain or affect other areas of the body. This activity of wild-type (WT) BChE, however, is relatively low. This prompted the design of variants of BChE which exhibit significantly improved catalytic activity against (-)-cocaine. Plants are a promising means to produce large amounts of these cocaine hydrolase variants of BChE, cheaply, safely with no concerns regarding human pathogens and functionally equivalent to enzymes derived from other sources. Here, in expressing cocaine-hydrolyzing mutants of BChE in Nicotiana benthamiana using the MagnICON virus-assisted transient expression system, and in reporting their initial biochemical analysis, we provide proof-of-principle that plants can express engineered BChE proteins with desired properties.

Reaction Pathway and Free Energy Profile for Cocaine Hydrolase-Catalyzed Hydrolysis of (-)-Cocaine

Reaction pathway of (-)-cocaine hydrolysis catalyzed by our recently discovered most efficient cocaine hydrolase, which is the A199S/F227A/S287G/A328W/Y332G mutant of human butyrylcholinesterase (BChE), and the corresponding free energy profile have been studied by performing first-principles pseudobond quantum mechanical/molecular mechanical (QM/MM)-free energy (FE) calculations. Based on the QM/MM-FE results, the catalytic hydrolysis process consists of four major reaction steps, including the nucleophilic attack on carbonyl carbon of (-)-cocaine benzoyl ester by hydroxyl group of S198, dissociation of (-)-cocaine benzoyl ester, nucleophilic attack on carbonyl carbon of (-)-cocaine benzoyl ester by water, and finally the dissociation between (-)-cocaine benzoyl group and S198 of the enzyme. The second reaction step is rate-determining. The calculated free energy barrier associated with the transition state for the rate-determining step is ~15.0 kcal/mol, which is in excellent agreement with the experimentally-derived activation free energy of ~14.7 kcal/mol. The mechanistic insights obtained from the present study will be valuable for rational design of more active cocaine hydrolase against (-)-cocaine. In particular, future efforts aiming at further increasing the catalytic activity of the enzyme against (-)-cocaine should focus on stabilization of the transition state for the second reaction step in which the benzoyl ester of (-)-cocaine dissociates.

Cocaine esterase-cocaine binding process and the free energy profiles by molecular dynamics and potential of mean force simulations

The combined molecular dynamics (MD) and potential of mean force (PMF) simulations have been performed to determine the free energy profiles for the binding process of (-)-cocaine interacting with wild-type cocaine esterase (CocE) and its mutants (T172R/G173Q and L119A/L169K/G173Q). According to the MD simulations, the general protein-(-)-cocaine binding mode is not affected by the mutations; e.g.. the benzoyl group of (-)-cocaine is always bound in a subsite composed of aromatic residues W151, W166, F261, and F408 and hydrophobic residue L407, while the carbonyl oxygen on the benzoyl group of (-)-cocaine is hydrogen-bonded with the oxyanion-hole residues Y44 and Y118. According to the PMF-calculated free energy profiles for the binding process, the binding free energies for (-)-cocaine with the wild-type, T172R/G173Q, and L119A/L169K/G173Q CocEs are predicted to be -6.4, -6.2, and -5.0 kcal/mol, respectively. The computational predictions are supported by experimental kinetic data, as the calculated binding free energies are in good agreement with the experimentally derived binding free energies, i.e., -7.2, -6.7, and -4.8 kcal/mol for the wild-type, T172R/G173Q, and L119A/L169K/G173Q, respectively. The reasonable agreement between the computational and experimental data suggests that the PMF simulations may be used as a valuable tool in new CocE mutant design that aims to decrease the Michaelis-Menten constant of the enzyme for (-)-cocaine.

Pharmacokinetic strategies for treatment of drug overdose and addiction

The pharmacokinetic treatment strategy targets the drug molecule itself, aiming to reduce drug concentration at the site of action, thereby minimizing any pharmacodynamic effect. This approach might be useful in the treatment of acute drug toxicity/overdose and in the long-term treatment of addiction. Phase IIa controlled clinical trials with anticocaine and antinicotine vaccines have shown good tolerability and some efficacy, but Phase IIb and III trials have been disappointing because of the failure to generate adequate antibody titers in most participants. Monoclonal antibodies against cocaine, methamphetamine and phencyclidine have shown promise in animal studies, as has enhancing cocaine metabolism with genetic variants of human butyrylcholinesterase, with a bacterial esterase, and with catalytic monoclonal antibodies. Pharmacokinetic treatments offer potential advantages in terms of patient adherence, absence of medication interactions and benefit for patients who cannot take standard medications.