Organic Process Research & Development
Article
To this was added 1.5 volumes of 4.9 M aqueous NaCN
solution dropwise (total reaction volume concentration = 100
g/L of 1). These studies evaluated the reaction rate at 21, 30,
and 37 °C by HPLC analysis at 1, 2.5, 4.5, 7, and 24 h.13 This
study revealed that while both reactions at 21 °C and at 30 °C
showed complete consumption of the starting materials over
24 h, the reaction at 37 °C did not reach completion,
presumably due to enzyme deactivation from both the high
temperature and substrate concentration. Interestingly, lower
reaction concentrations of 68 g/L of 1 (0.5 M) at this higher
37 °C temperature were tolerated and allowed the reaction to
be complete, further highlighting this subtle interplay between
substrate concentration and temperature for this process.
Based on this observation, this effect was further investigated
using more challenging lower whole cell enzyme loadings (2%
by weight) to identify the preferred reaction conditions that
enabled robust performance while maximizing throughput
(Table 3). These results showed incomplete reactions for the
enzyme.16 This approach was preferred since it often can
explain the in vitro results while enabling exploration of the
enzyme’s multidimensional fitness landscapes, thereby allowing
evolutionary studies to focus on smaller data sets of mutants at
a fraction of the cost versus typical directed evolution
approaches.
We began our rational engineering of BCJ2315 by
constructing a homology model of the enzyme and submitting
it for extensive molecular dynamics (MD) simulations to better
sample the positions of all atoms in explicit solvent (Figure 1).
These studies found the active center is predominately
comprised of the anticipated catalytic triad (E48, K130 and
C164), with the following residues that comprise the binding
pocket: Y54, L56, W59, T134, E137, W165, W188, P189,
S190, F191, S192, A197, A198, G202, P203, N206, L246,
Q247, A248, G249, and G250.
The nitrilase reaction mechanism can be divided into three
steps:17 (1) The thiol group of the catalytic C164 which
undergoes nucleophilic attack on the electrophilic carbon of
the substrate nitrile group, forming a thioimidate intermediate;
(2) the thioimidate intermediate is hydrolyzed forming a
second acyl enzyme intermediate, with release of ammonia;
(3) and finally the acyl enzyme intermediate is hydrolyzed to
produce the carboxylic acid. In this overall process, the first
step of the reaction mechanism should be rate limiting and was
therefore chosen as the initial enzyme state for docking of the
(R)- and (S)-substrates.
Homology model validation was conducted by docking (R)-
2-chloromandelonitrile and (S)-2-chloromandelonitrile to WT
BCJ2315 and comparing computation findings with previously
reported experimental findings.7 This analysis found the
binding of (R)-2-chloromandelonitrile was better than that of
its enantiomer, as the coordination between the substrate
nitrile carbon and the C164 thiol group is shorter for (R)-2-
chloromandelonitrile.13 Docking of our substrates ((R)-2 and
(S)-2) in this validated enzyme model revealed that the
reaction coordinate for the first step is shorter and the angle of
the nucleophilic attack is better for (R)-2 and therefore
potentially closer to the Michaelis−Menten complex (Figure
2). Moreover, (R)-2 is able to stabilize stronger intermolecular
interactions within the active site via hydrogen bonding of the
S192 side chain with the oxygen of the methoxy group and also
by hydrogen bonding of W165 with the hydroxyl group of (R)-
2. Both interactions are not present in the docking pose of the
(S)-2. These findings again correlate to the experimentally
observed selective formation of 3 from (R)-2 by BCJ2315.
In silico model analysis focused predominately on identifying
favorable mutations of two general types: first, active site
mutations that could improve selectivity by (1) increasing the
binding affinity for the (R)-2 and (2) decrease the step 1
reaction coordinate of the (R)-2 while subsequently increasing
it for the (S)-2, and second, through B-factor analysis using
four replicate 50 ns MD simulations of BCJ2315 to identify
flexible residues/regions of the enzyme that could be modified
to improve structural stability.13 Additional mutants were also
selected from previous literature reports and through
coevolution analysis of nitrilase sequences.7a In total, 96
mutants were assessed experimentally in the first round of in
vitro screening which revealed that a substantial number of
mutants were able to improve selectivity and/or activity/
stability.13 A list of mutants which gave either increased
selectivity and conversion or an increase in both, relative to
WT, are detailed in Table 4. All mutants were expressed at the
Table 3. Effect of Aldehyde and Temperature on Reaction
a
Profile using Lyophilized Whole Cell BCJ2315
Aldehyde
Loading (g/L)
Temp % AN 3 + % AN ent-3 after
b
Entry
(°C)
24 h (%)
% ee
1
2
3
4
5
6
7
8
68
68
45
37
37
30
21
30
21
19−24
64
95
84
95
96
17
94
93
94
94
94
93
93
N/A
100
100
100
150
150
200
48−89
c
−
a
The reaction mixtures were prepared by mixing 7.5 volumesof
aqueous aldehyde-bisulfite adduct solution, prepared via reaction of
appropriate amount of 1 with 1.1 equiv of NaHSO3 in water, with 1
volume of ethyl acetate, followed by dropwise addition of 1 equiv of
NaCN in 1.5 volumes of water. BCJ2315 (2% by weight relative to 1)
was added, and the reaction mixture was stirred at the appropriate
b
temperature. Determined by HPLC. Full conversion was deemed
achieved when no 2 was detected at 283 nm. Conversion was
determined by [% AN 3 + % AN ent-3]/[% AN 1 + % AN 2 + % AN
c
3 + % AN ent-3].15 No 3 or ent-3 was detected.
process using ethyl acetate cosolvent at higher reaction
temperature and/or substrate concentrations. A further
complicating factor was that at ambient temperatures and
100 g/L of 1 loading, the aqueous solution of the bisulfite-
aldehyde adduct 4 used to prepare cyanohydrin 2 has a
concentration of 0.979 M, which is near its saturation point of
1.02 M.14 This was believed to be a contributing reason for the
observed inconsistent results when increasing the aldehyde 1
loading to 150 g/L (Table 3, entry 7). For this reason, the
target temperature and substrate concentration chosen were 21
°C and 100 g/L of 1 for the process, as it was expected to give
the most consistent results regarding reaction rate, total
turnover, and selectivity. As shown in Table 3, the ee of the
product was not affected by changes to either substrate loading
or temperature under these reaction conditions.
With suitable reaction conditions in hand, we began enzyme
engineering nitrilase BCJ23156 to fit the process with the goal
of improving substrate loading, enzyme robustness, and the
enantioselectivity of 3. A rational in silico approach was used to
identify key enzyme residues for modification based on our
previous success in efficiently re-engineering a transaminase
D
Org. Process Res. Dev. XXXX, XXX, XXX−XXX