Enzyme Engineering Based on X-ray Structures and Kinetic Profiling

State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic and Developmental Sciences, MOE-LSB & MOE-LSC, Scho...
0 downloads 4 Views 3MB Size
Subscriber access provided by - Access paid by the | UCSB Libraries

Enzyme engineering based on X-ray structures and kinetic profiling of substrate libraries: alcohol dehydrogenases for stereospecific synthesis of a broad range of chiral alcohols Yao Nie, Shanshan Wang, Yan Xu, Shenggan Luo, Yi-Lei Zhao, Rong Xiao, Gaetano Montelione, John F. Hunt, and Thomas Szyperski ACS Catal., Just Accepted Manuscript • DOI: 10.1021/acscatal.8b00364 • Publication Date (Web): 04 Apr 2018 Downloaded from http://pubs.acs.org on April 4, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Catalysis

1

Enzyme engineering based on X-ray structures and

2

kinetic profiling of substrate libraries: alcohol

3

dehydrogenases for stereospecific synthesis of a

4

broad range of chiral alcohols

5

Yao Nie,†# Shanshan Wang,†‡# Yan Xu,†§* Shenggan Luo,¶ Yi-Lei Zhao,¶* Rong Xiao,∥⫫ Gaetano

6

T. Montelione,∥ John F. Hunt,± and Thomas Szyperski⊥

7



8

Jiangnan University, Wuxi 214122, China.

9



School of Biotechnology, Key laboratory of Industrial Biotechnology, Ministry of Education,

School of Biological Science and Engineering, Shannxi University of Technology, Hanzhong

10

723001, China.

11

§

12

China.

13



14

Metabolic and Developmental Sciences, MOE-LSB & MOE-LSC, School of Life Sciences and

15

Biotechnology, Shanghai Jiao Tong University, Shanghai 200240, China

16

∥Center

17

Biochemistry, Rutgers, The State University of New Jersey, Piscataway, NJ 08854, USA

18

±

State Key Laboratory of Food Science and Technology, Jiangnan University, Wuxi 214122,

State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of

for Advanced Biotechnology and Medicine, Department of Molecular Biology and

Department of Biological Sciences, Columbia University, New York, NY 10027, USA

ACS Paragon Plus Environment

1

ACS Catalysis 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

19



20

USA

Page 2 of 36

Department of Chemistry, The State University of New York at Buffalo, Buffalo, NY 14260,

21 22

*E-mail: [email protected] (Y. Xu); [email protected] (Y.-L. Zhao)

23

24

ABSTRACT: The narrow substrate scope of naturally occurring alcohol dehydrogenases

25

(ADHs) greatly limits the enzymatic synthesis of important chiral alcohols. Based on X-ray

26

crystal structures and kinetic profiling of a substrate library, we engineered variants of the

27

stereospecific alcohol dehydrogenase from Candida parapsilopsis. This resulted in a set of four

28

mutant enzymes which enable the asymmetric reduction of a broad range of prochiral ketones,

29

including valuable pharmaceuticals and fine chemicals. The engineering strategy of this study

30

paves the way for creating additional ADHs tailored for production of complex chiral alcohols.

31

KEYWORDS: alcohol dehydrogenase, X-ray crystal structure, substrate library, kinetic

32

profiling, stereoselectivity, substrate specificity

33 34

ACS Paragon Plus Environment

2

Page 3 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

35

ACS Catalysis

INTRODUCTION

36

Chiral alcohols are versatile building blocks for synthesis of pharmaceuticals and fine

37

chemicals, such as antidepressants, anti-asthmatics, cholesterol-lowering agents, adrenergic

38

receptor agonists, α- or β-adrenergic drugs, etc.1-3 For example, (R)-ethyl-4,4,4-trifluoro-3-

39

hydroxybutanoate is the precursor for synthesizing befloxatone, a potent anti-depressant, and (S)-

40

ethyl-3-hydroxy-3-phenylpropionate can be used to produce fluoxetine, tomoxetine, and

41

nisoxetine.4,5 Alcohol dehydrogenases (ADH, EC 1.1.1.1-394) catalyze the reduction of prochiral

42

ketones with very high chemo-, regio-, and stereoselectivitiy, and offer a valuable synthetic route

43

for these alcohols.6-8 However, most naturally occurring ADHs exhibit a rather narrow substrate

44

scope and thus do not reduce most of the ketone intermediates required for industrial

45

applications.9

46

Protein engineering provides approaches to expand the substrate scope of enzymes.10 In

47

particular, structure-based redesign has been used to optimize the substrate scope for

48

biotechnological important lipases,11,12 epoxide hydrolases13,14 and a large variety of

49

oxidoreductases.15,16 In spite of an urgent demand, however, comparably few engineered ADHs

50

have been reported so far.17-21

51

Previously, we characterized the alcohol dehydrogenase/carbonyl reductase (EC 1.1.1.1)

52

from Candida parapsilosis (CpRCR) as a robust, highly stereoselective biocatalyst following

53

Prelog’s rule for ketone reduction,22-24 and determined the X-ray crystal structures in complex

54

with ligands.25 However, CpRCR likewise exhibits a rather narrow substrate scope.26

55

Here, we present the engineering of a set of CpRCR variants enabling stereospecific

56

synthesis of broad range of chiral alcohols. We selected a representative substrate library

ACS Paragon Plus Environment

3

ACS Catalysis 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 36

57

containing 12 aryl ketones, 8 β-ketoesters, and 4 aliphatic ketones (Scheme 1), and measured the

58

Michaelis-Menten kinetic parameters27 for wild-type (wt) CpRCR in order to accurately profile

59

substrate specificity. Furthermore, we measured the enantiomeric excess (e.e.) for the products to

60

profile stereoselectivity. Next, based on crystal structures and the profiling, we engineered

61

CpRCR variants by replacing bulky residues forming the active site. The resulting catalytic

62

efficiency and stereoselectivity was likewise profiled.

63 64

Scheme 1.Substrate library consisting of ketones (1a-24a), where substituent RL is larger than

65

RS. The CpRCR catalyzed reduction with NADH (shown at the top) yields the corresponding

66

chiral alcohols (1b-24b), illustrated assuming that stereoselective reduction follows Prelog’s

67

rule.22 The colors used for the substrates indicate how efficiently they are reduced by wt CpRCR

68

(see Figure 1).

69 70

ACS Paragon Plus Environment

4

Page 5 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Catalysis

71

RESULTS AND DISCUSSION

72

Kinetic profiling of wt CpRCR towards the substrate library

73 74

Figure 1. Catalytic profile of wt CpRCR towards 1a-24a of the substrate library (Scheme 1). A:

75

Michaelis-Menten parameters. From left (21a) to right (19a), substrates are sorted by increasing

76

molecular volume (for details see Table S1). For each substrate, Km (left axis), kcat (middle axis),

77

and kcat/Km (right axis) are represented, respectively, by open, black and grey bars. Standard

78

deviations for kcat and Km were obtained from triplicate measurements. B: Logarithmic graph

79

(log2-scale) of kcat with the increase of molecular volume.

80

Figure 1A provides the Michaelis-Menten kinetic parameters of wt CpRCR obtained for our

81

substrate library (Scheme 1). The kcat value is considered to be the most important parameter for

82

catalytic performance of enzymes28 under substrate saturation, as commonly encountered in

83

biotechnological processes. Based on the kcat values, the substrates were thus divided into three

84

classes (Figure 1B):

85

Class (1): “good substrates” with high kcat > 4 s-1, located in the orange region of Figure 1B,

86

includes the four substrates13a, 1a, 8a, and 14a. Compound 13a, the smallest β-ketoester, is

87

most efficiently reduced by wt CpRCR exhibiting the highest kcat (8.98 s-1) along with good

ACS Paragon Plus Environment

5

ACS Catalysis 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 36

88

affinity (Km 0.25 mM), which results in high catalytic efficiency (kcat/Km 35.92 mM-1 s-1).

89

Molecule 14a with a chlorinated methyl group as RS performs comparably. Substrates 1a and 8a,

90

aryl ketones, exhibit kcat > 7 s-1 and good affinities (Km 0.50 mM and Km 0.95 mM, respectively).

91

Hence, the substrates of class (1) are reduced very efficiently by wt CpRCR.

92

Class (2): “moderate substrates” with kcat in the range 2-4 s-1, located in the yellow region of

93

Figure 1B, includes 8 of the substrates. Catalytic activities are reduced by varying degrees

94

towards the 1a-derivatives 2a, 6a, 7a, and 9a bearing substituents on the benzene ring (‘RL‘) and

95

at the methyl group (‘RS‘). β-ketoester 17a with a phenyl group as RL exhibits good affinity (Km

96

0.18 mM) and high catalytic efficiency (kcat/Km 14.28 mM-1s-1). Unexpectedly, for aliphatic

97

ketones (21a, 23a, 24a), catalytic efficiency increases gradually along with the increase of the

98

length of carbon chain of RL. Substrate 24a bearing hexyl group is more suitable for binding to

99

the active site of wt CpRCR, as reflected by a lower Km (0.34 mM) and higher kcat (3.45 s-1) when

100

compared with the propyl, butyl, and pentyl congeners.

101

Class (3): “poor substrates” with kcat < 2 s-1, located in the blue region of Figure 1B,

102

includes half (12) of the substrates of our library (Scheme 1). Although 3a and 4a have the same

103

substituent at different C-position compared to 2a, the kcat values of 3a and 4a are reduced by

104

varying degrees. Moreover, molecule 5a bearing bulky bromo substituent has a very low kcat

105

value (0.83 s-1). The substrates 10a-12a, which have bulky aryl ketones, and β-ketoesters such as

106

trihalogenated compounds (15a and 16a) and ethylbenzoyl acetate derivatives (18a and 19a), are

107

also poor substrates for reduction with wt CpRCR.

108

The kinetic data (Figure 1A) effectively profile the substrate specificity of wt CpRCR.

109

Importantly, log(kcat) values are moderately anti-correlated with the molecular volume of the

ACS Paragon Plus Environment

6

Page 7 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Catalysis

110

substrate (Figure 1B). This finding suggests that steric hindrance affects catalytic efficiency

111

significantly. Since wt CpRCR exhibits kcat values which can be considered to be sufficiently

112

high for biotechnology applications only for the four class (1) substrates (Scheme 1 and Figure

113

1), engineering of CpRCR variants is required to improve catalytic efficiency towards class (2)

114

and (3) substrates.

115

Redesign of substrate-binding cavities

116

We recently described the high-resolution X-ray crystal structures25 of wt CpRCR bound to

117

NAD+ at 2.15 Å resolution (PDB ID: 3WLE), CpRCR mutant H49A bound to class 1 substrate

118

8a at 2.95 Å (PDB ID: 3WNQ), and wt CpRCR bound to its product (R)-8b at 2.30Å (PDB ID:

119

3WLF). These structures validated our detailed understanding of the catalytic mechanism of

120

ADHs.8 The active site contains both a large and a small hydrophobic cavity close to the

121

catalytic center, as well as the NADH binding site. The small cavity is formed mainly by the

122

bulky side chains of F285 and W286, together with a segment of NADH (Figure 2). We

123

anticipated that these two aromatic residues play an important role in determining substrate

124

specificity. In the large cavity, the bulky side chain of residue W116 exhibits a hydrophobic

125

interaction with the phenyl ring of substrate 8a in 3WNQ. Furthermore, significant spatial

126

rearrangements of the side chains of W116, F285 and W286, creating a more “closed”

127

conformation, were observed upon substrate binding (Figure S1).

128

Based on our kinetic profiling data (Figure 1) and high-resolution structural information

129

(Figure 2), we hypothesized that replacement of W116 or F285/W286 by Ala may (i) increase

130

the sizes of the small and large cavities, respectively, and (ii) provide increased conformational

131

flexibility of the active site to allow binding of a broader range of ketone substrates. In turn, we

ACS Paragon Plus Environment

7

ACS Catalysis 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 36

132

anticipated that such engineered enzymes may be suitable for catalyzing the reduction of class (2)

133

and (3) substrates with greatly improved efficiency.

134 135

Figure 2. Substrate-binding cavities in CpRCR (PDB ID: 3WLE) and key residues selected for

136

engineering. A: Surfaces of the large (magenta) and small cavity (yellow) of CpRCR. B: Stick

137

representation of the side chains of W116 (large cavity; magenta), F285, W286 and a segment of

138

NAD+ (small cavity; yellow).

139

Comparison of catalytic performance of the variants with wt CpRCR

140

To test our hypothesis, we prepared samples of CpRCR mutants W116A, F285A, W286A,

141

and F285A/W286A. Measurements of Michaelis-Menten parameters (Figures S2-S5) revealed

142

dramatic changes in substrate specificity profiles and kcat values relative to wt CpRCR (Figure 3).

143

For all four class (1) substrates (orange region in Figure 3), the ratios of kcat values (mutant/wt

144

CpRCR) were around or below 1.0 for all mutants. This finding validates our conclusion that the

145

active site of wt CpRCR is well tailored for small class (1) substrates. In contrast, for all class (2)

146

substrates except 21a and 24a (yellow region), at least one of the mutant enzymes exhibits higher

147

catalytic efficiency than wt CpRCR. Remarkably, kcat values for W116A, F285A, W286A, and

148

F285A/W286A are increased to 2.6, 1.6, 4.2, and 4.4-fold towards 2a, respectively. For all class

ACS Paragon Plus Environment

8

Page 9 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Catalysis

149

(3) substrates except for 12a and 15a (blue region), there is likewise at least one mutant which

150

exhibits higher catalytic efficiency than wt CpRCR, and for nine substrates (22a, 3a, 4a, 5a, 20a,

151

11a, 16a, 18a, 19a) the increase of kcat for one of the mutants is at least ~2 or higher. Specifically,

152

W116A shows 2.9 and 3.0-fold increases of kcat towards substrates 16a and 19a, F285A/W286A

153

exhibits 2.8, 2.9, 2.6, and 3.0-fold increases of kcat towards 22a, 4a, 11a, and 18a, and F285A

154

shows 5.5, 8.8, and 8.5-fold increases of kcat for 3a, 20a, and 18a, respectively.

155 156

Figure 3. Ratios of kcat values of mutants (W116A, F285A, W286A, and F285A/W286A) and wt

157

CpRCR towards substrates 1a-24a (Scheme 1) sorted by class (1) (orange), class (2) (yellow),

158

and class (3) (blue) (Figure 1).

159

Overall, for all but four of the 20 class (2) and (3) substrates (Figure 1), at least one of the

160

engineered mutants exhibits enhanced catalytic efficiency under substrate saturation (Figure 3).

161

The enhanced catalytic performance can likely be attributed to a reshaped active site of the

162

mutants: while wt CpRCR is suitable only for the smaller substrates, our engineered mutants can

ACS Paragon Plus Environment

9

ACS Catalysis 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 36

163

accommodate also more bulky substituents and exhibit a kinetic substrate profile which neatly

164

complements wt CpRCR (Figure 3).

165

Comparison of stereoselectivity of the variants with wt CpRCR

166

Since active site engineering may well affect stereoselectivity of ketone reduction and thus

167

enantiopurity of the alcohols,25 we also measured the optical purity of chiral alcohols resulting

168

from reduction of compounds constituting our substrate library (Scheme 1) by either wt or

169

mutant CpRCRs (Figure 4). The products of wt CpRCR-mediated reduction of 1a-24a follow

170

Prelog’s rule.22 Likewise, W116A exhibits Prelog stereopreference for all ketones, while

171

providing for some of them even higher stereoselectivity than wt CpRCR. This finding indicates

172

that the W116A mutation within the large cavity does not significantly affect stereoselectivity.

173

On the contrary, the W286A and F285A/W286A mutations dramatically affect stereoselectivity

174

(Figure 4). W286A exhibits for 14 ketones Prelog, and for 10 ketones anti-Prelog

175

stereoselectivity, with the e.e. values rarely exceeding 75%. F285A/W286A exhibits Prelog

176

specificity for only 6 ketones and anti-Prelog stereoselectivity for 18 ketones. Remarkably, for

177

eight of the anti-Prelog ketones (1a, 14a, 24a, 20a, 12a, 11a, 16a, and 17a), the e.e. exceeds

178

80%, so that chiral alcohols can be synthesized which cannot be efficiently obtained with most of

179

known naturally occurring ADHs. These findings show that reshaping of the small cavity of

180

CpRCR can alter stereoselectivity. Substrate 12a represents the most striking example of our

181

library: wt and W116A CpRCR yield the Prelog product with 100%e.e., while F285A, W286A,

182

and F285A/W286A yield the anti-Prelog product with 100%e.e.

ACS Paragon Plus Environment

10

Page 11 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Catalysis

183 184

Figure 4. Optical purity (%e.e.) of products from asymmetric reduction of 1a-24a, sorted by

185

increasing molecular volume (Table S1) from left to right, catalyzed by purified wt CpRCR and

186

mutants (for assay details see Table S2). Positive values of e.e. are expected following Prelog’s

187

rule, while negative values indicate an anti-Prelog stereopreference.

188

Overall, the combination of structure-based protein engineering with kinetic profiling of a

189

substrate library represents a powerful approach for creating a set of four ADH variants offering

190

broad substrate specificity and distinct stereoselectivity. Importantly, the kinetic profiling

191

enables experimental assessment of the impact of entropic effects and structural rearrangements

192

which are difficult to capture with computational methods.

193

Computational study on substrate specificity and stereoselectivity

194

To rationalize experimentally observed kinetic parameters (Figures 1 and 3, Figures S2-S5)

195

and stereoselectivity (Figure 4) we initiated a comprehensive computational modeling study.

196

First, to assess the impact of the mutations on the substrate-binding pocket, we performed 100 ns

197

molecular dynamics (MD) simulations for both the wt CpRCR-NADH complex and the mutant

ACS Paragon Plus Environment

11

ACS Catalysis 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 36

198

congeners. Pocket volume distribution analysis (Figure 5) revealed that: (i) W286A (small cavity

199

mutation) exhibits an average pocket volume and distribution which is quite similar to the wild

200

type; (ii) F285A (small cavity mutation) and F285A/W286A (two small cavity mutations) exhibit

201

larger average pocket volumes as well as broader distributions; and (iii) W116A (large cavity

202

mutation) exhibits a smaller average pocket volume and a narrower distribution. As expected,

203

these findings indicate that the mutations result in significant conformational rearrangements

204

coming along with a significant and varying degree of ‘collapse’ of the binding pocket, which

205

results from the replacement of the aromatic sidechains with a methyl group. Moreover, the

206

broadening of the volume distributions suggest that some of these mutations result in increased

207

conformational flexibility, which may enable the accommodation of a broader range of

208

substrates.

209 210

Figure 5. Distributions for pocket volumes of wt CpRCR and the mutants W116A, F285A,

211

W286A, and F285A/W286A derived from 100 ns MD simulations.

212

Second, we implemented a quantum mechanics/molecular mechanics (QM/MM)29 protocol

213

enabling us to eventually develop a rationale for substrate stereoselectivity (Figure 4). The

ACS Paragon Plus Environment

12

Page 13 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Catalysis

214

protocol was used to calculate the pro-S and pro-R transition states (TSs) for substrate 20a in wt

215

CpRCR which lead, respectively, to ‘Prelog’ and ‘anti-Prelog’ products. It was quite

216

straightforward to determine that, due to limited space in the active site pocket, the pro-R TS

217

could be modeled only by ‘flipping’ of the RS group by a rotation around the carbonyl-RL bond

218

(Figure 6). The similar result corresponding to the pro-S and pro-R TSs for substrate 20a in wt

219

CpRCR was also obtained to the double mutant F285A/W286A (Figures S6). Consistent with the

220

experimentally observed Prelog stereospecificity (Figure 4), the energy of the pro-S TS was

221

calculated to be about 11.3 kcal/mol lower than that of the pro-R. Notably, the ‘RS group flip’

222

might represent a distinct mechanism to explain inversion of stereospecificity in ADHs and may

223

guide future engineering studies. Our QM/MM protocol shall be applied for other substrates as

224

well as mutant enzymes to shed new light on the stereospecificity of ADHs in general.

225 226

Figure 6. The superposition of pro-S (yellow, expected according the Prelog’s rule) and pro-R

227

(green) transition states of substrate 20a in wt CpRCR.

228 229

ACS Paragon Plus Environment

13

ACS Catalysis 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 36

230 231

CONCLUSIONS B

A Ratio of kcat

0.5 0.8 1 1.5

2

3

4

F285A/W286A F285A/W286A W116A F285A W286A F285A/W286A F285A F285A/W286A F285A/W286A F285A F285A/W286A F285A/W286A W116A W116A W116A F285A/W286A F285A F285A F286A F285A/W286A W116A F285A/W286A F285A W116A

21a 13a 1a 22a 8a 2a 3a 4a 6a 15a 23a 5a 9a 7a 14a 24a 20a 10a 12a 11a 16a 17a 18a 19a

232

5

W116A F285A W286A

F285A/ W286A

e.e. (%)

100 90 75 50 25 0 -25 -50 -75 -90 -100

21a 13a 1a 22a 8a 2a 3a 4a 6a 15a 23a 5a 9a 7a 14a 24a 20a 10a 12a 11a 16a 17a 18a 19a

W116A WT W116A WT WT W116A WT WT W116A W116A W116A WT WT W116A WT WT W116A WT WT WT W116A WT WT WT WT

W116A F285A W286A

F285A/W286A F285A/W286A F285A/W286A W286A F285A/W286A F285A/W286A F285A/W286A F285A/W286A F285A/W286A F285A/W286A F285A/W286A F285A/W286A F285A/W286A F285A/W286A F285A/W286A F285A/W286A F285A/W286A F285A/W286A F285A/W286A F285A/W286A F285A/W286A F285A/W286A F285A/W286A F285A/W286A

F285A/ W286A

233

Figure 7. Comparison of catalytic performance of wt CpRCR with four engineered CpRCR

234

mutants. A: Ratio of mutant to wt CpRCR kcat values for each compound in the substrate library

235

(Scheme 1). The mutant with the largest ratio is indicated along with the corresponding color

236

code in the right-most column. B: Enantiomeric excess (e.e.) of wt CpRCR and the four mutants

237

for each compound. Positive e.e. values indicate synthesis following the Prelog’s rule, while the

238

negative values represent anti-Prelog synthesis (see Table S2 for details). The enzyme with the

239

highest positive or the lowest negative e.e. value is indicated along with the corresponding color

240

code in the right-most two columns.

241

We successfully altered the substrate specificity of CpRCR to allow reduction of a

242

broadened range of substrates, including significantly more bulky ketones (Figure 7). The

243

resulting set of four CpRCR mutants is now available for biotechnological applications (notably,

244

these four mutants were the only ones we expressed and tested). Furthermore, our strategy

ACS Paragon Plus Environment

14

Page 15 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Catalysis

245

allowed us to create mutants which either retain or reverse stereoselectivity of wt CpRCR, and

246

we initiated a computational study to rationalize our findings. In particular, the reversal of

247

stereoselectivity observed for selected variants and substrates enables the synthesis of chiral

248

alcohols which thus far can not be obtained with known naturally occurring ADHs. Figure 7

249

exemplifies the high complementarity of wt CpRCR together with four specific mutants for

250

synthesis of a broad range of alcohols with varying stereoselectivity. Depending on the desired

251

substrate and stereoselectivity, this analysis enables one to choose the best CpRCR variant. It is

252

evident that the engineering strategy described here can be employed to generate a suitably

253

expanded set of CpRCR mutant proteins.

254 255

EXPERIMENTAL SECTION

256

Materials

257

Enzymes, vectors, oligonucleotides and other reagents for DNA cloning and amplification

258

were obtained from Takara-Bio Co., Japan and Novagen Co., USA. Substrates (1a-24a) were

259

purchased from Sigma-Aldrich (USA) or TCI (Shanghai). NADH, NAD+, (R)-/(S)-1a,(R)-/(S)-8a

260

were purchased from Sigma-Aldrich (USA). All other racemic alcohol standards were prepared

261

by reduction of ketones with sodium borohydride. Hexane and isopropanol used for High

262

Performance Liquid Chromatography (HPLC) were of chromatographic grade (Honeywell Co.,

263

USA). All other chemicals used were of analytical grade and commercially available.

264

Expression and purification of wt CpRCR and variants

ACS Paragon Plus Environment

15

ACS Catalysis 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 36

265

Recombinant strains of wt CpRCR and the variants (F285A, W286A, and F285A/W286A)

266

were constructed previously using the pET32 Xa/LIC vector.25 The mutation W116A was

267

introduced by using the protocol of the Agilent QuikChange Lightning Site-Directed

268

Mutagenesis Kit. Importantly, the four mutants (W116A, F285A, W286A, and F285A/W286A)

269

were

270

GTAAAAACGCATTTGGTGATGCGTTCGGATTGGGGTACGATGGTG - 3’, W116A_R: 5’-

271

CACCATCGTACCCCAATCCGAACGCATCACCAAATGCGTTTTTAC - 3’. Plasmid DNA

272

was prepared and sequenced to validate positive clones. Recombinant plasmids were transformed

273

in Escherichia coli BL21 trxB (DE3). Cells were grown in LB medium containing 100 µg mL-1

274

ampicillin, 50 µg mL-1 kanamycin, and 0.2 mM zinc acetate. The protein was expressed without

275

induction of isopropyl β-D-1-thiogalactopyranoside (IPTG) at 28 °C for 12 hours.30 Cells were

276

harvested by centrifugation at 7,000 ×g and resuspended in the lysis buffer (pH 8.0) containing

277

25 mM Tris-HCl, 150 mM NaCl, and 20 mM imidazole. Cell lysates were prepared with a

278

French Press, and then centrifuged at 17,000 ×g to remove cell debris. Proteins were purified by

279

FPLC using a Ni-NTA column. The buffer and gradient for the Ni-NTA column were 25 mM

280

Tris-HCl, 150 mM NaCl, and 20-500 mM imidazole, pH 8.0. His-tagged CpRCR was eluted at

281

about 120 mM imidazole. The protein solution was dialyzed against the buffer containing 25

282

mM Tris-HCl, 150 mM NaCl, pH 8.0, and then subjected to Factor Xa digestion to remove the

283

His-tag. The mixture was then passed through Ni-NTA column again and the untagged CpRCR

284

was eluted with 5 mM imidazole-containing buffer. The purified protein was finally concentrated

285

to 3 mg mL-1 in the buffer (pH 8.0) containing 25 mM Tris-HCl and 150 mM NaCl.

286

Determination of kinetic parameters

the

only

ones

constructed

for

our

study.

Primers

were:

W116A_F:

5’-

ACS Paragon Plus Environment

16

Page 17 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Catalysis

287

The standard assay mixture in 100 µL for the enzyme activity of the purified wt CpRCR and

288

the variants comprised 0.1 M potassium phosphate buffer (pH 6.5), 5 mM ZnCl2, 1 mM NADH,

289

5 mM substrates, and an appropriate amount of enzyme. The decrease in NADH absorbance at

290

340 nm was recorded for 3 min at 30 °C. A molar extinction coefficient of 6,220 L mol-1 cm-1

291

was used for NADH. One unit (U) of enzyme activity was defined as 1 µmol of NADH

292

consumed per minute under the assay conditions. Kinetic parameters of wt CpRCR and the

293

variants towards various substrates were assayed by measuring initial velocity at various

294

concentrations of substrates and cofactor.25 To determine the apparent Km value, the

295

concentration of carbonyl compounds was varied from 0.1 to 8 mM with a fixed concentration of

296

NADH at 0.05 mM, 0.1 mM, and 0.25 mM, respectively. The reaction follows a Theorell-

297

Chance BiBi mechanism.27 The Michaelis-Menten kinetic parameters were obtained from a non-

298

linear least-squres fit of the Michaelis-Menten equation to the experimental data and further

299

validated from double reciprocal Lineweaver–Burk plots. All the data were averaged from three

300

replicates for each substrate and cofactor concentration and significant differences (p