Hein J. Wijma

Structure- and sequence-analysis inspired engineering of proteins for enhanced thermostability

Protein engineering strategies for increasing stability can be improved by replacing random mutagenesis and high-throughput screening by approaches that include bioinformatics and computational design. Mutations can be focused on regions in the structure that are most flexible and involved in the early steps of thermal unfolding. Sequence analysis can often predict the position and nature of stabilizing mutations, and may allow the reconstruction of thermostable ancestral sequences. Various computational tools make it possible to design stabilizing features, such as hydrophobic clusters and surface charges. Different methods for designing chimeric enzymes can also support the engineering of more stable proteins without the need of high-throughput screening.

Computationally designed libraries for rapid enzyme stabilization

The ability to engineer enzymes and other proteins to any desired stability would have wide-ranging applications. Here, we demonstrate that computational design of a library with chemically diverse stabilizing mutations allows the engineering of drastically stabilized and fully functional variants of the mesostable enzyme limonene epoxide hydrolase. First, point mutations were selected if they significantly improved the predicted free energy of protein folding. Disulfide bonds were designed using sampling of backbone conformational space, which tripled the number of experimentally stabilizing disulfide bridges. Next, orthogonal in silico screening steps were used to remove chemically unreasonable mutations and mutations that are predicted to increase protein flexibility. The resulting library of 64 variants was experimentally screened, which revealed 21 (pairs of) stabilizing mutations located both in relatively rigid and in flexible areas of the enzyme. Finally, combining 10–12 of these confirmed mutations resulted in multi-site mutants with an increase in apparent melting temperature from 50 to 85°C, enhanced catalytic activity, preserved regioselectivity and a >250-fold longer half-life. The developed Framework for Rapid Enzyme Stabilization by Computational libraries (FRESCO) requires far less screening than conventional directed evolution.

A Random-sequential Mechanism for Nitrite Binding and Active Site Reduction in Copper-containing Nitrite Reductase

The homotrimeric copper-containing nitrite reductase (NiR) contains one type-1 and one type-2 copper center per monomer. Electrons enter through the type-1 site and are shuttled to the type-2 site where nitrite is reduced to nitric oxide. To investigate the catalytic mechanism of NiR the effects of pH and nitrite on the turnover rate in the presence of three different electron donors at saturating concentrations were measured. The activity of NiR was also measured electrochemically by exploiting direct electron transfer to the enzyme immobilized on a graphite rotating disk electrode. In all cases, the steady-state kinetics fitted excellently to a random-sequential mechanism in which electron transfer from the type-1 to the type-2 site is rate-limiting. At low [NO–2] reduction of the type-2 site precedes nitrite binding, at high [NO–2] the reverse occurs. Below pH 6.5, the catalytic activity diminished at higher nitrite concentrations, in agreement with electron transfer being slower to the nitrite-bound type-2 site than to the water-bound type-2 site. Above pH 6.5, substrate activation is observed, in agreement with electron transfer to the nitrite-bound type-2 site being faster than electron transfer to the hydroxyl-bound type-2 site. To study the effect of slower electron transfer between the type-1 and type-2 site, NiR M150T was used. It has a type-1 site with a 125-mV higher midpoint potential and a 0.3-eV higher reorganization energy leading to an ∼50-fold slower intramolecular electron transfer to the type-2 site. The results confirm that NiR employs a random-sequential mechanism.

Computational design gains momentum in enzyme catalysis engineering

Computational protein design is becoming a powerful tool for tailoring enzymes for specific biotechnological applications. When applied to existing enzymes, computational re‐design makes it possible to obtain orders of magnitude improvement in catalytic activity towards a new target substrate. Computational methods also allow the design of completely new active sites that catalyze reactions that are not known to occur in biological systems. If initial designs display modest catalytic activity, which is often the case, this may be improved by iterative cycles of computational design or by follow‐up engineering through directed evolution. Compared to established protein engineering methods such as directed evolution and structure‐based mutagenesis, computational design allows for much larger jumps in sequence space; for example, by introducing more than a dozen mutations in a single step or by introducing loops that provide new functional interactions. Recent advances in the computational design toolbox, which include new backbone re‐design methods and the use of molecular dynamics simulations to better predict the catalytic activity of designed variants, will further enhance the use of computational tools in enzyme engineering.

Directed Evolution Strategies for Enantiocomplementary Haloalkane Dehalogenases: From Chemical Waste to Enantiopure Building Blocks

We used directed evolution to obtain enantiocomplementary haloalkane dehalogenase variants that convert the toxic waste compound 1,2,3‐trichloropropane (TCP) into highly enantioenriched (R)‐ or (S)‐2,3‐dichloropropan‐1‐ol, which can easily be converted into optically active epichlorohydrins—attractive intermediates for the synthesis of enantiopure fine chemicals. A dehalogenase with improved catalytic activity but very low enantioselectivity was used as the starting point. A strategy that made optimal use of the limited capacity of the screening assay, which was based on chiral gas chromatography, was developed. We used pair‐wise site‐saturation mutagenesis (SSM) of all 16 noncatalytic active‐site residues during the initial two rounds of evolution. The resulting best R‐ and S‐enantioselective variants were further improved in two rounds of site‐restricted mutagenesis (SRM), with incorporation of carefully selected sets of amino acids at a larger number of positions, including sites that are more distant from the active site. Finally, the most promising mutations and positions were promoted to a combinatorial library by using a multi‐site mutagenesis protocol with restricted codon sets. To guide the design of partly undefined (ambiguous) codon sets for these restricted libraries we employed structural information, the results of multiple sequence alignments, and knowledge from earlier rounds. After five rounds of evolution with screening of only 5500 clones, we obtained two strongly diverged haloalkane dehalogenase variants that give access to (R)‐epichlorohydrin with 90 % ee and to (S)‐epichlorohydrin with 97 % ee, containing 13 and 17 mutations, respectively, around their active sites.

Enantioselective enzymes by computational design and in silico screening.

Computational enzyme design holds great promise for providing new biocatalysts for synthetic chemistry. A strategy to design small mutant libraries of complementary enantioselective epoxide hydrolase variants for the production of highly enantioenriched (S,S)-diols and (R,R)-diols is developed. Key features of this strategy (CASCO, catalytic selectivity by computational design) are the design of mutations that favor binding of the substrate in a predefined orientation, the introduction of steric hindrance to prevent unwanted substrate binding modes, and ranking of designs by high-throughput molecular dynamics simulations. Using this strategy we obtained highly stereoselective mutants of limonene epoxide hydrolase after experimental screening of only 37 variants. The results indicate that computational methods can replace a substantial amount of laboratory work when developing enantioselective enzymes.

Biocatalytic and structural properties of a highly engineered halohydrin dehalogenase.

Two highly engineered halohydrin dehalogenase variants were characterized in terms of their performance in dehalogenation and epoxide cyanolysis reactions. Both enzyme variants outperformed the wild‐type enzyme in the cyanolysis of ethyl (S)‐3,4‐epoxybutyrate, a conversion yielding ethyl (R)‐4‐cyano‐3‐hydroxybutyrate, an important chiral building block for statin synthesis. One of the enzyme variants, HheC2360, displayed catalytic rates for this cyanolysis reaction enhanced up to tenfold. Furthermore, the enantioselectivity of this variant was the opposite of that of the wild‐type enzyme, both for dehalogenation and for cyanolysis reactions. The 37‐fold mutant HheC2360 showed an increase in thermal stability of 8 °C relative to the wild‐type enzyme. Crystal structures of this enzyme were elucidated with chloride and ethyl (S)‐3,4‐epoxybutyrate or with ethyl (R)‐4‐cyano‐3‐hydroxybutyrate bound in the active site. The observed increase in temperature stability was explained in terms of a substantial increase in buried surface area relative to the wild‐type HheC, together with enhanced interfacial interactions between the subunits that form the tetramer. The structures also revealed that the substrate binding pocket was modified both by substitutions and by backbone movements in loops surrounding the active site. The observed changes in the mutant structures are partly governed by coupled mutations, some of which are necessary to remove steric clashes or to allow backbone movements to occur. The importance of interactions between substitutions suggests that efficient directed evolution strategies should allow for compensating and synergistic mutations during library design.

Stabilization of cyclohexanone monooxygenase by a computationally designed disulfide bond spanning only one residue

Enzyme stability is an important parameter in biocatalytic applications, and there is a strong need for efficient methods to generate robust enzymes. We investigated whether stabilizing disulfide bonds can be computationally designed based on a model structure. In our approach, unlike in previous disulfide engineering studies, short bonds spanning only a few residues were included. We used cyclohexanone monooxygenase (CHMO), a Baeyer–Villiger monooxygenase (BVMO) from Acinetobacter sp. NCIMB9871 as the target enzyme. This enzyme has been the prototype BVMO for many biocatalytic studies even though it is notoriously labile. After creating a small library of mutant enzymes with introduced cysteine pairs and subsequent screening for improved thermostability, three stabilizing disulfide bonds were identified. The introduced disulfide bonds are all within 12 Å of each other, suggesting this particular region is critical for unfolding. This study shows that stabilizing disulfide bonds do not have to span many residues, as the most stabilizing disulfide bond, L323C–A325C, spans only one residue while it stabilizes the enzyme, as shown by a 6 °C increase in its apparent melting temperature.

Computational Library Design for Increasing Haloalkane Dehalogenase Stability

We explored the use of a computational design framework for the stabilization of the haloalkane dehalogenase LinB. Energy calculations, disulfide bond design, molecular dynamics simulations, and rational inspection of mutant structures predicted many stabilizing mutations. Screening of these in small mutant libraries led to the discovery of seventeen point mutations and one disulfide bond that enhanced thermostability. Mutations located in or contacting flexible regions of the protein had a larger stabilizing effect than mutations outside such regions. The combined introduction of twelve stabilizing mutations resulted in a LinB mutant with a 23 °C increase in apparent melting temperature (Tm,app, 72.5 °C) and an over 200‐fold longer half‐life at 60 °C. The most stable LinB variants also displayed increased compatibility with co‐solvents, thus allowing substrate conversion and kinetic resolution at much higher concentrations than with the wild‐type enzyme.

Reconstitution of the type-1 active site of the H145G/A variants of nitrite reductase by ligand insertion

Variants of the copper-containing nitrite reductase (NiR) of Alcaligenes faecalis S6 were constructed by site-directed mutagenesis, by which the C-terminal histidine ligand (His145) of the Cu in the type-1 site was replaced by an alanine or a glycine. The type-1 sites in the NiR variants as isolated, are in the reduced form, but can be oxidized in the presence of external ligands, like (substituted) imidazoles and chloride. The reduction potential of the type-1 site of NiR-H145A reconstituted with imidazole amounts to 505 mV vs NHE (20 degrees C, pH 7, 10 mM imidazole), while for the native type-1 site it amounts to 260 mV. XRD data on crystals of the reduced and oxidized NiR-H145A variant show that in the reduced type-1 site the metal is 3-coordinated, but in the oxidized form takes up a ligand from the solution. With the fourth (exogenous) ligand in place the type-1 site is able to accept electrons at about the same rate as the wt NiR, but it is unable to pass the electron onto the type-2 site, leading to loss of enzymatic activity. It is argued that the uptake of an electron by the mutated type-1 site is accompanied by a loss of the exogenous ligand and a concomitant rise of the redox potential. This rise effectively traps the electron in the type-1 site.