Computational Evidence for the Enzymatic Transformation of 2

Publication Date (Web): July 12, 2018. Copyright © 2018 American Chemical Society. Cite this:ACS Earth Space Chem. XXXX, XXX, XXX-XXX ...
0 downloads 0 Views 1MB Size
Subscriber access provided by STEPHEN F AUSTIN STATE UNIV

Article

Computational Evidence for the Enzymatic Transformation of 2-hydroxypropylphosphonate to Methylphosphonate Yanwei Li, Xiaodan Wang, Ruiming Zhang, Junjie Wang, Zhongyue Yang, Likai Du, Xiaowen Tang, Qingzhu Zhang, and Wenxing Wang ACS Earth Space Chem., Just Accepted Manuscript • DOI: 10.1021/ acsearthspacechem.8b00070 • Publication Date (Web): 12 Jul 2018 Downloaded from http://pubs.acs.org on July 16, 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 26 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 Earth and Space Chemistry

1 2

Computational Evidence for the Enzymatic

3

Transformation of 2-hydroxypropylphosphonate to

4

Methylphosphonate

5 6 7

Yanwei Li†, Xiaodan Wang†, Ruiming Zhang†, Junjie Wang†, Zhongyue

8

Yang‡, Likai Du§*, Xiaowen Tangǁ, Qingzhu Zhang†*, Wenxing Wang†

9 †

10

Environment Research Institute, Shandong University, Jinan, 250100, P. R. China ‡

11

Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States

12 §

13

Hubei Key Laboratory of Agricultural Bioinformatics, College of Informatics, Huazhong Agricultural University, Wuhan, 430070, P.R. China.

14 15

ǁ

School of Pharmaceutical Sciences, Sun Yat-sen University, Guangzhou 510006, P. R. China

16 17

Keywords

18 19

Quantum

mechanics/molecular

20

Greenhouse

21

efficiency

22

___________________________________________________________

23

*

24

Fax: 86-531-8836 1990

gas

methane,

mechanics,

Enzymatic

Methylphosphonate,

transformation,

Enzymatic

Corresponding authors. E-mail: [email protected][email protected]

25 1   

ACS Paragon Plus Environment

ACS Earth and Space Chemistry 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 2 of 26

26 27

Abstract

28

Understanding the origins of the greenhouse gas methane in the ocean is

29

of great environmental importance, especially for global climate change and the flow

30

of

31

2-hydroxyethylphosphonate dioxygenase (HEPD) has been reported to catalyze the

32

transformation of 2-hydroxypropylphosphonate (2-HEP) to methylphosphonate

33

(MPn), a compound that can be easily transformed to methane by C–P lyase in marine

34

microbe. Here, HEPD E176H-catalyzed transformation of 2-HEP to MPn was

35

investigated at the molecular level by using QM/MM method. The results evidenced

36

the feasibility of the transformation of 2-HEP to MPn and highlighted that the

37

transformation contains five elementary steps: H-abstraction, O-O bond cleavage,

38

H-transfer, C-C bond cleavage, and MPn formation. H-abstraction was found to be the

39

rate-determining step with an energy barrier of 17.8 kcal/mol, which is in reasonable

40

accordance with the experimentally determined rate constant (0.38 s-1, correspond to

41

18.0 kcal/mol). Three intersystem crossing events were involved in H-abstraction,

42

H-transfer, and MPn formation steps. Residue electrostatic analysis on the rate

43

determining step suggests that proper mutation of Tyr174 may improve the enzymatic

44

efficiency.

carbon

within

the

earth

surface

system.

45 46 47 48 49 50 51 52

2   

ACS Paragon Plus Environment

A

mutant

(E176H)

of

Page 3 of 26 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 Earth and Space Chemistry

53

1. Introduction

54

Methane, a greenhouse gas with more than 20-fold more warming

55

potential than carbon dioxide, is supersaturated in the ocean relative to the atmosphere

56

1, 2

57

Understanding the origin of the excess methane in the ocean is of great environmental

58

importance, especially for global climate change and the flow of carbon within the

59

earth surface system 4. It has been hypothesized that methane comes from the

60

transformation of methylphosphonate (MPn) by carbon-phosphorus lyase

61

hypothesis was once challenged by the fact that no MPn was detected in the marine

62

ecosystem. Recently a breakthrough finding has supported that MPn can be vastly

63

produced

64

exopolysaccharides in ubiquitous marine archaeon Nitrosopumilus maritime, adding a

65

strong evidence to support the hypothesis that the excess methane in the ocean comes

66

from MPn 9, 10.

. The ocean contributes up to 4 % of the annual global methane emission 3.

by

methylphosphonate

synthase

(MPnS)

to

build

3-8

the

. This

polar

67

MPnS is a non-heme iron-dependent oxygenase that converts

68

2-hydroxyethylphosphonate (2-HEP) to MPn at a moderate rate constant of 0.18 s-1 11.

69

Lowering its catalytic efficiency may help control the methane emission from the

70

ocean and regulate the global climate while improving its efficiency may be helpful in

71

rapidly harnessing methane for industrial and daily use. It is thus of great significance

72

to understand the enzymatic processes both kinetically and structurally, however,

73

efforts to crystalize MPnS have not been successful

3   

ACS Paragon Plus Environment

12

. Lacking crystal structure

ACS Earth and Space Chemistry 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 26

74

hinders further molecular-level investigation of the catalytic mechanism of MPnS. It

75

has been reported that the active site of MPnS is highly conserved with an enzyme

76

known as 2-hydroxyethylphosphonate dioxygenase (HEPD), whose crystal structure

77

has been successfully deposited

78

2-HEP to MPn, just like MPnS does

79

alternative of MPnS to investigate the enzymatic transformation of 2-HEP to MPn.

12-16

. Actually, HEPD E176H mutant can convert 12

. HEPD E176H mutant can thus serve as an

80

The transformation of 2-HEP to hydroxymethylphosphonate (HMP) by

81

wild type HEPD (rather than HEPD E176H mutant) has been extensively investigated

82

through many state-of-the-art experimental techniques

83

theoretical investigations which used DFT and QM/MM methods to elucidate the

84

detailed mechanism at molecular level

85

such as a ferric superoxide was found to active 2-HEP and start the reaction, the

86

H-abstraction from 2-HEP was demonstrated to be the rate-determining step, and

87

valence state change [Fe(III)→Fe(IV)→Fe(III)] was spotted during the transformation

88

from 2-HEP to HMP. Based on these results, possible reaction mechanism of HEPD

89

E176H-catalyzed transformation of 2-HEP to MPn has been proposed

90

However, an in-depth understanding at molecular level is still not available.

91

Understanding the detailed mechanism is essential for regulating the transformation

92

efficiency of 2-HEP to MPn. For instance, with a deep understanding of the

93

mechanism, people can perform site-directed mutation studies to accelerate or

94

decelerate the reaction rate. This can be particularly useful for regulating the methane

95

emission from the ocean to the atmosphere or harnessing methane for industrial or

22-25

12-21

. There are also several

. Many valuable results were obtained,

4   

ACS Paragon Plus Environment

11, 12, 20

.

Page 5 of 26 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 Earth and Space Chemistry

96

daily use.

97

Here the detailed mechanism of HEPD E176H-catalyzed transformation

98

of 2-HEP to MPn was investigated by using QM/MM method. QM/MM is a

99

multiscale computational method that balances the accuracy and efficiency in

100

modelling enzymatic reactions. It has become an increasingly powerful tool to

101

complement experimental enzyme chemistry

102

feasibility of the transformation of 2-HEP to MPn and highlighted that the

103

transformation contains five elementary steps: H-abstraction, O-O bond cleavage,

104

H-transfer, C-C bond cleavage, and MPn formation. H-abstraction was found to be the

105

rate-determining step. Interestingly, three intersystem crossing (ISC) events were

106

observed along the transformation processes.

107

2. Methods

26-30

. The results evidenced the

The selected model of HEPD E176H was obtained based on the QM/MM

108

24, 25

109

calculations previously as described

110

work was similar to our previous studies and only important details of the methods are

111

briefly summarized here. The QM/MM calculations were performed by using

112

ChemShell 31 software, which interfaces QM program Turbomole 32 to MM program

113

DL-POLY 33. The charge shift model

114

used during the QM/MM calculations. The geometries of the intermediates and

115

transition states were optimized by using hybrid delocalized internal coordinates

116

optimizer and microiterative TS optimizer under the B3LYP/def-SV(P)//CHARMM22

117

level 36. Frequency calculations were performed to verify the one imaginary frequency

118

character of transition state structures, and the suitability of the transition vector was

34

. The QM/MM protocol used in the present

and electrostatic embedding method

5   

ACS Paragon Plus Environment

35

were

ACS Earth and Space Chemistry 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

119

also confirmed. Additional single point energy calculations were carried out at

120

B3LYP/def2-TZVPP//CHARMM22 level for a better description of the energy profiles.

121

Similar choice of method has been successfully applied in different systems by many

122

individual groups 37, 38. For HEPD E176H catalyzed reaction, the QM region contains

123

Fe, His129, His176, His182, 2-HEP, O2, and two water molecules (Scheme 1). This

124

resulted in 57 QM atoms. Atoms within 20 Å of Fe were allowed to move (~3700

125

MM atoms) while the rest of the system (~29000 MM atoms) was fixed during

126

QM/MM calculations. The total charge of the whole system is 0. Five different spin

127

multiplicities (singlet, triplet, quintet, septet, and nonet) were considered during

128

QM/MM calculations to determine the most feasible reaction pathway. Additional

129

computational details were provided in the Supporting Information.

130

3. Results and Discussion

131

By using QM/MM method, we have previously identified the detailed

132

reaction mechanism of wild type 2-hydroxyethylphosphonate dioxygenase (HEPD)

133

catalyzed transformation of 2-HEP to HMP

134

mechanism of HEPD E176H-catalyzed transformation of 2-HEP to MPn was fully

135

investigated. Five elementary steps were found: H-abstraction, O-O bond cleavage,

136

H-transfer, C-C bond cleavage, and MPn formation, as indicated in Scheme 1. Each

137

individual elementary step will be discussed in detail in the following paragraphs.

138

3.1 Reaction mechanism and energy profiles

24, 25

. Here the detailed reaction

139

It is critical to determine the most stable spin state in the reactant (R-1).

140

Literally, triplet, quintet or septet state was found to be more stable for non-heme 6   

ACS Paragon Plus Environment

Page 6 of 26

Page 7 of 26 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 Earth and Space Chemistry

22-25, 37, 38

141

iron-dependent systems

. Here the singlet, triplet, quintet, septet, and nonet

142

state have been systematically analyzed, and the order of relative energies is

143

determined as: septet