Processive Degradation of Crystalline Cellulose by a Multimodular

Apr 4, 2018 - Processive hydrolysis of crystalline cellulose by cellulases is a critical step for lignocellulose deconstruction. The classic Trichoder...
0 downloads 0 Views 2MB Size
Subscriber access provided by Warwick University Library

Processive Degradation of Crystalline Cellulose by a Multimodular Endoglucanase via a Wirewalking Mode Kun-Di Zhang, Wen Li, Yefei Wang, Yan-Lin Zheng, Fang-Cheng Tan, Xiao-Qing Ma, Lishan Yao, Edward A. Bayer, Lushan Wang, and Fuli Li Biomacromolecules, Just Accepted Manuscript • DOI: 10.1021/acs.biomac.8b00340 • Publication Date (Web): 04 Apr 2018 Downloaded from http://pubs.acs.org on April 8, 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 38 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

Biomacromolecules

1

Processive Degradation of Crystalline Cellulose by a

2

Multimodular Endoglucanase via a Wirewalking

3

Mode

4

Kun-Di Zhanga,b, Wen Lib, Ye-Fei Wanga, Yan-Lin Zhengc, Fang-Cheng Tana, Xiao-Qing Mae,

5

Li-Shan Yaoa, Edward A. Bayerd, Lu-Shan Wangb,*, Fu-Li Lia,*

6

a

7

Institute of Bioenergy and Bioprocess Technology, Chinese Academy of Sciences, 266101

8

Qingdao, P. R. China; bState Key Laboratory of Microbial Technology, Shandong University,

9

Jinan, 250100, P. R. China; cCollege of Mathematics and Systems Science, Shandong University

Shandong Provincial Key Laboratory of Energy Genetics, Key Laboratory of Biofuels, Qingdao

10

of Science and Technology, Qingdao, 266590, P. R. China; dDepartment of Biomolecular

11

Sciences, The Weizmann Institute of Science, Rehovot, 76100, Israel; eShandong Provincial Key

12

Laboratory of Synthetic Biology

13

KEYWORDS: GH9, CBM3, Open cleft, Cellotetraose, Clostridium cellulosi

14

15

ABSTRACT: Processive hydrolysis of crystalline cellulose by cellulases is a critical step for

16

lignocellulose deconstruction. The classic Trichoderma reesei exoglucanase TrCel7A, which has

17

a closed active-site tunnel, starts each processive run by threading the tunnel with a cellulose

18

chain. Loop regions are necessary for tunnel conformation, resulting in weak thermostability of

ACS Paragon Plus Environment

1

Biomacromolecules 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 38

1

fungal exoglucanases. However, endoglucanase CcCel9A, from the thermophilic bacterium

2

Clostridium cellulosi, comprises a glycoside hydrolase (GH) family 9 module with an open cleft

3

and five carbohydrate-binding modules (CBMs) and hydrolyzes crystalline cellulose

4

processively. How CcCel9A and other similar GH9 enzymes bind to the smooth surface of

5

crystalline cellulose to achieve processivity is still unknown. Our results demonstrate that the C-

6

terminal CBM3b and three CBMX2s enhance productive adsorption to cellulose, while the

7

CBM3c adjacent to the GH9 is tightly bound to 11 glucosyl units, thereby extending the catalytic

8

cleft to 17 subsites, which facilitates decrystallization by forming a supramodular binding

9

surface. In the open cleft, the strong interaction forces between substrate-binding subsites and

10

glucosyl rings enable cleavage of the hydrogen bonds and extraction of a single cellulose chain.

11

In addition, subsite −4 is capable of drawing the chain to its favored location. Cellotetraose is

12

released from the open cleft as the initial product to achieve high processivity, which is further

13

hydrolyzed to cellotriose, cellobiose and glucose by the catalytic cleft of the endoglucanase. On

14

this basis, we propose a wirewalking mode for processive degradation of crystalline cellulose by

15

an endoglucanase, which provides insights for rational design of industrial cellulases.

16

Introduction

17

Cellulose is the earth’s most abundant biopolymer and is responsible for the structural

18

integrity of plants 1. Processive cellulases degrade crystalline cellulose and are defined

19

by their ability to remain bound to their substrates and continuously repeat cycles of

20

catalysis before dissociating

21

reacting with adjacent sites without dissociating from the cellulose chain

22

two types of processive cellulases: exoglucanases and processive endoglucanases.

23

Generally, their structures are similar and consist of a glycoside hydrolase (GH) module

2-3

. They do so by binding and sliding along substrates and 4-6

. There are

ACS Paragon Plus Environment

2

Page 3 of 38 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

Biomacromolecules

1

and one or more carbohydrate-binding modules (CBMs). The cellulases attach to the

2

hydrophobic surface of the crystalline substrate by the CBMs, thus forming an enzyme–

3

substrate complex that produces cello-oligosaccharides in a successive stepwise

4

manner as the enzyme slides along a single cellulose chain. However, there is a

5

distinguishing structural feature of active-site architecture between the two processive

6

enzyme types – exoglucanases have a closed tunnel while endoglucanases have an

7

open cleft.

8

One of the best studied exocellulases is Cel7A of the fungus, Trichoderma reesei

9

(TrCel7A). It is composed of a GH7 catalytic module and a CBM1. An early molecular

10

dynamics study suggested that the enzyme–substrate interactions produced forces that

11

enabled the enzyme to slide along the cellulose surface with a single substrate chain

12

channeled through its catalytic tunnel, which was referred to as threading mode

13

visual study using high-speed atomic force microscopy suggested that the movement

14

was accompanied by catalytic activity and that the CBM enhanced the concentration of

15

enzyme molecules on the substrate 10. Other studies have highlighted the importance of

16

hydrophobic interactions between amino acid residues in the tunnel and the cellulose

17

chain

18

processivity of GHs is directly related to oligosaccharide binding free energy, and a

19

mathematical model has been derived to quantify the binding free energy of the enzyme

20

to the cello-oligosaccharide 6, 13-17.

7-9

.A

10-12

. Recently, more detailed theories have been proposed that suggest the

21

It has been reported that deletions in the loop regions of the catalytic tunnel in Cel7D

22

from Phanerochaete chrysosporium (PcCel7D) created a long cellulose-binding groove,

23

which resulted in more binding sites compared to that of TrCel7A

18

. In addition,

ACS Paragon Plus Environment

3

Biomacromolecules 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 38

1

cellobiohydrolase Cel6A (E2) from the bacterium Thermobifida fusca (TfCel6A) had

2

enhanced endoglucanase activity after deletion of a surface loop

3

temperatures protein loops tend to be denatured; hence the loops are considered

4

unstable “hotspots” on proteins. As a result, the reported thermostability of

5

cellobiohydrolases, in which the catalytic tunnel is rich in loops, was relatively low

6

However, the surface of crystalline cellulose was described as so smooth that the

7

salient points could only be observed at near-atomic resolution by atomic force

8

microscopy

9

become the challenge for cellulose hydrolysis, especially at high temperatures, and

10

novel modes for crystalline cellulose degradation are currently being studied. A

11

subsequent advancement in this field of research is the discovery of processive

12

endoglucanases that can hydrolyze crystalline cellulose.

13

21

19

, and it was thus difficult for cellulases without loops to bind

. At elevated

20

.

22-23

. This has

The first processive endoglucanase TfCel9A was identified from the filamentous soil 24

14

bacterium T. fusca, and it had the ability to degrade cellulose

15

GH9 catalytic module with a family 3c CBM fused to its C-terminal end, which is

16

essential for processivity, and a family 3b CBM at the C-terminus of the intact protein,

17

which is essential for binding to crystalline cellulose. TfCel9A-68, generated from

18

TfCel9A by limited proteolysis, contained only the GH9 and CBM3c modules. Its crystal

19

structure revealed that these two modules provided an elongated flat face, to which

20

linear cellulose chains can adhere. When binding to a cellopentaose molecule, there

21

was clear electron density for six glucosyl units in two clearly separate segments in its

22

active-site architecture: a four-glucosyl segment with a reducing end that had an α-

23

configuration and a two-glucosyl segment. The six glucosyl binding sites in the cleft

. The enzyme contains a

ACS Paragon Plus Environment

4

Page 5 of 38 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

Biomacromolecules

1

were numbered from −4 to +2 from the nonreducing to the reducing end. Dozens of

2

amino acid residues that interacted with glucosyl subsites −4 to +2 of the active site

3

were mutated to analyze the mechanism of action of TfCel9A. It was found that E424

4

and D58 were the catalytic acid and base, respectively 25-26. The processivity of mutants

5

R378K and D261A, as well as the double mutant D261A/R378K, was examined and

6

found to be improved, reduced, and even more reduced, respectively. This suggests

7

that processivity requires a precise balance between the catalytic cleft binding on the

8

two sides of the cleavage site, as R378 interacted with subsites +1 and +2, while D261

9

interacted with subsites −1 to −4

26

. Interestingly, at the nonreducing end of the active-

10

site cleft of TfCel9A, there is a blockage formed by residues 245–255. This block of

11

residues leaves insufficient space for a hypothetical subsite −5 to bind in a linear cleft,

12

which leads to product yields no longer than cellotetraose

13

TfCel9A was deduced whereby the enzyme is bound to the surface of the crystalline

14

substrate via CBM2, whereas CBM3c together with GH9 adheres to the extracted single

15

chain

16

have been reported, most from cellulosome-producing Clostridium spp. (Fig. 1). For

17

example, Cel9G from C. cellulolyticum consists of GH9, CBM3c, and a dockerin domain

18

instead of a cellulose-targeting CBM3b. This and other dockeirn-bearing enzymes rely

19

on the CBM3a or CBM3b of the cellulosomal scaffoldin subunit for the cellulose-

20

targeting property. The crystal structure of C. cellulolyticum Cel9G is very similar to that

21

of TfCel9A, and further analysis revealed complementary surfaces between the flat

22

surface of crystalline cellulose and Cel9G

23

belong to noncellulosomal cellulases and cellulosomal components of this species,

24

. The productivity state of

27

. In addition to TfCel9A, several other GH9/CBM3c processive endoglucanases

28

. Cel9I and Cel9R from C. thermocellum

ACS Paragon Plus Environment

5

Biomacromolecules 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 38

29-30

. Both of them exhibit the GH9/CBM3c modular scheme, and both

1

respectively

2

release cellotetraose as their initial product.

3

Clostridium cellulosi is a thermophilic, Gram-positive obligate anaerobe, which can

4

ferment a broad range of carbon sources and generate ethanol, acetate, CO2, and

5

hydrogen

6

cell-free extracts of Geobacillus toebii DSM 14590 to the medium of the bacterium

7

utilizing lignocellulose as a sole carbon source at an optimal temperature of 60°C.

8

Zymography analysis revealed that two GHs, ORF3880 and ORF3883, play key roles in

9

the hydrolysis of lignocellulose. ORF3880 and ORF3883 are both multimodular

10

endocellulases that consist of one GH9 module and several CBMs (Fig. 1), and were

11

thus named CcCel9A and CcCel9B, respectively

12

understanding how CcCel9A with an open cleft maintains processivity by analyzing its

13

multimodular structure and active-site architecture.

14

Experimental Section

31

. C. cellulosi CS-4-4 was isolated from decayed corn stalk by addition of

32

. In this study, we focus on

15

Cloning of the gene encoding CcCel9A and its truncated mutants-Genomic DNA

16

of C. cellulosi CS-4-4 was extracted according to Andreou et al (2013) 33. The ORF3880

17

gene product (GenBank accession number LM995447) was designated CcCel9A. The

18

coding sequences for wild-type CcCel9A and its truncated mutants (Fig. 2, TM1–TM5)

19

were amplified from the genomic DNA of C. cellulosi CS-4-4 by PCR using Pyrobest

20

DNA polymerase (Takara, Dalian, Japan). Genes of mutants TM14, TM15 and TM16

21

were obtained by fusing two crude fragments amplified from the wild-type CcCel9A

22

gene. The PCR products were cloned into the pEASY-E2 vector according to the

ACS Paragon Plus Environment

6

Page 7 of 38 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

Biomacromolecules

1

protocols described by the manufacturer (TransGen Biotech, Beijing, China), and the

2

correct colonies were selected based on DNA sequencing.

3

Assay of enzymatic activity-The specific enzyme activities for carboxymethyl

4

cellulose (CMC), regenerated amorphous cellulose (RAC), and microcrystalline cellu-

5

lose (Avicel) were assayed as previously described

6

mM CaCl2 to the reaction buffer. Avicel PH-101 and CMC were purchased from Sigma

7

(Beijing, China); RAC was prepared according to Zhang et.al (2006)

8

values of CcCel9A and its mutants toward CMC, RAC, and Avicel were determined by

9

fitting the initial velocities at different substrate concentrations using the Michaelis–

10

Menten equation. The assays were performed in triplicate with at least five substrate

11

concentrations, varying from approximately 1/10 to 20-fold of the Km values, under

12

optimal conditions.

34

with a modification of adding 1

36

. The Km and kcat

13

Measurement of processivity-Processivity was determined by using a protocol

14

described by Ozdemir, et al (2012) 37 with modifications. Single discs of Whatman No. 1

15

filter paper (GE Healthcare, China) with a radius of 5 mm and total weight of 3 mg were

16

incubated with 9 µM purified proteins at 60°C in PC buffer (50 mM phosphate, 12 mM

17

citrate, 1 mM sodium azide, pH 6.5) for 16 h except for TM4 and TM5 which were

18

incubated in 50 mM sodium acetate buffer, pH 5). Then the supernatant fluids were

19

transferred to another tube, and 150 µL were taken for measuring the amount of

20

reducing sugars. The filter paper discs were washed with PC buffer or acetate buffer

21

three times. After washing, 150 µL of the same buffer were added to the discs. Two

22

hundred µL of 3,5-Dinitrosalicylic acid (DNS) reagent were added to the supernatant

23

fluids and to the filter paper tubes. Reducing sugars were determined as described in

ACS Paragon Plus Environment

7

Biomacromolecules 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 38

1

the enzyme assay methods. The ratio of soluble sugars (found in the supernatant

2

fraction) to insoluble sugars (found on the filter paper) was calculated as processivity. All

3

of the reactions were performed in triplicate.

4

Binding of CcCel9A and its truncated mutants to cellulose-The potential of wild-

5

type and mutant CcCel9A to bind to cellulose was quantitatively measured at 60°C as

6

previously described 32.

7

Thermostability assay-The thermostability of CcCel9A and its mutants was

8

determined by incubating 1 µM of each enzyme at temperatures between 50°C and

9

75°C for 1 h. Residual activity was expressed as the ratio of specific activity on Avicel to

10

the activity of proteins without incubation measured at optimal temperature.

11

Computational simulation-In order to discern the three-dimensional structure of

12

CcCel9A, homology modeling was conducted using the X-ray crystal structure of Cel9G

13

from Clostridium cellulolyticum as the template. We searched the Protein Data Bank

14

(PDB) at National Center for Biotechnology and Information using BLAST with the

15

amino acid sequence of CcCel9A. The results showed that C. cellulolyticum Cel9G had

16

the highest sequence identity (63%) with the query sequence. Its X-ray crystal structure

17

(PDB ID: 1KFG) was selected as the template for homology modeling and the three-

18

dimensional structure of CcCel9A was constructed using Discovery Studio 2.5 based on

19

the MODELER program. The structure of the modeled protein was viewed using PyMol

20

1.5. The cellohexaose structure was extracted from the crystal structure of TfCel9A

21

(PDB ID: 4TF4) and docked into the cleft of CcCel9A according to the hotspots of the

22

substrate and enzyme using Libdock. Amino acid residues within 5 Å to cellohexaose

ACS Paragon Plus Environment

8

Page 9 of 38 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

Biomacromolecules

1

were considered capable of interacting with the substrate. Then, the residues of

2

CcCel9A were aligned with TfCel9A

3

and CcCel9B as previously described 38. All highly conserved residues were selected as

4

targets in the site-directed mutagenesis test.

5

24

, CtCel9I

29

, CtCel9R

30

, CpCel9

35

, CtCel9N

29

,

Molecular dynamics (MD) simulations-MD simulations were performed using 39

40

6

Gromacs 4.5

7

force field

8

adding 12.5 Å TIP3P water in a cubic box, and sodium ions were used to neutralize the

9

system. Before free energy calculations, 5,000 steps of energy minimization followed by

10

1 ns MD simulation at constant pressure (1 atm) and temperature (300 K) were

11

performed to equilibrate the system. The pressure was regulated using the extended

12

ensemble Parrinello–Rahman approach

13

modified Berendson thermostat

14

evaluate the contribution of the long-range electrostatic interactions. A non-bonded pair

15

list cutoff of 12.0 Å was used, and the non-bonded pair list was updated every five steps.

16

All bonds involving hydrogen in proteins were constrained using the LINC algorithm

17

whereas bonds and angles of water molecules were forced by the SETTLE algorithm 49,

18

allowing a time step of 0.002 ps.

41

, with the Amber ff99SB force field

for the cello-oligomer in TIP3P water

for the protein and GLYCAM06

42

. The protein was solvated by

43-44

, and the temperature was controlled by a

45

. The Particle Mesh Ewald Method

46-47

was used to

48

,

19

Site-directed mutagenesis-Several aromatic amino acid residues located in each

20

subsite were selected (Table S2) for the deletion of their functional side chain, because

21

they had been previously identified as crucial for substrate interactions

22

directed mutagenesis, 100 ng of plasmid encoding CcCel9A was used as a template for

23

PCR amplification, and the steps were performed using KAPA HiFi HotStart polymerase

50-52

. For site-

ACS Paragon Plus Environment

9

Biomacromolecules 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 38

1

(Kapa Biosystems, MA) as follows: initial denaturation at 94°C for 3 min, followed by 18

2

cycles of 98°C for 30 s, 63°C for 30 s, and 72°C for 2.5 min (Primers are listed in Table

3

S3). The PCR product was digested using DpnI (New England Biolabs, MA) at 37°C for

4

4 h to degrade the template plasmid DNA. Further, the digestion product was

5

transformed into E. coli DH5α competent cells, and the competent cells were plated on

6

Luria-Bertani agar containing 100 µg/mL of ampicillin, using the spread-plate technique,

7

and the plates were incubated overnight at 37°C. Single colonies were propagated to

8

extract plasmids. The inserts were then sequenced (Genwiz Inc., Beijing, China) to

9

confirm the presence of the correct mutation. Gene expression and protein purification

10

were performed as previously described 32.

11

Analysis of hydrolysis products- The hydrolysis reaction mixture (150 µL)

12

consisted of 0.2 µM enzymes and 100 µL of 1% RAC in PC buffer. Aliquots of the

13

reaction mixture were removed at intervals and boiled for 10 min to inactivate the

14

enzyme. The enzymatic hydrolysis products of soluble cello-oligosaccharides were

15

analyzed by fluorescence-assisted carbohydrate electrophoresis (FACE). The gray

16

values of each band were acquired from the peaks generated using ImageJ k 1.45 and

17

extracted quantitatively53.

18

Calculation of affinity constants-The binding abilities of the wild-type and mutant 54

19

enzymes were evaluated using fluorescence spectra as previously described

20

(0.1 mg/mL) was mixed with 0–1.0% (w/v) CMC (0.1% interval) in 50 mM PBS buffer

21

(pH 6.5) for 30 min, which was then used to test the quenching of fluorescence using a

22

fluorescence spectrophotometer (Hitachi F-4600, Japan). All samples were excited at

. Protein

ACS Paragon Plus Environment

10

Page 11 of 38 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

Biomacromolecules

1

295 nm and scanned from 300–500 nm to integrate the spectra. CMC without protein

2

was used as control. Data were fit with non-linear regression: y = 1 / (a + bx)

3 (Eq. 1)

4 5

where y and x represent fluorescence intensity and CMC concentration (mM),

6

respectively. Ka (affinity constant) was calculated by b/a.

7

Results

8

Enzymatic characterization of truncated mutants of CcCel9A-Nine different

9

constructs containing the Cel9A catalytic module were cloned and expressed in

10

Escherichia coli (Fig. 2), using the primers listed in Table S1. The catalytic efficiencies of

11

mutants TM1, TM2, and TM3 that do not contain the cellulose-targeting CBM3b showed

12

higher activities towards CMC than those of CcCel9A, whereas those of mutants TM14,

13

TM15, and TM16, which possess the CBM3b, were higher on RAC and Avicel. The

14

processivities of mutants TM14, TM15 and CcCel9A with 1–3 CBMX2s were similar,

15

while the processivity of TM16 that lacks the CBMX2 was significantly lower. The Km

16

values on insoluble substrates (RAC and Avicel) of CBM3b-lacking mutants were

17

increased as the CBMX2 was deleted one by one from TM1–TM3 (Table S4).

18

Binding affinity of CcCel9A and its CBMX2-lacking mutants-To test whether the

19

CBMX2 is necessary for binding of the wild-type CcCel9A to cellulose, the binding

20

isotherms of the wild-type and CcCel9A mutants TM14, TM15, and TM16 were

21

compared. The dissociation constant (Kd) of wild-type CcCel9A was the lowest and the

22

maximum binding (qmax) was the highest, compared to the those of the mutants (Fig.

ACS Paragon Plus Environment

11

Biomacromolecules 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 38

1

3A). The partition coefficient (α), which reflects the binding ability of an enzyme, showed

2

that CcCel9A exhibited the highest binding to Avicel, whereas the binding of mutants

3

TM14 and TM15 was similar but 25% higher than that of TM16, suggesting that the

4

three CBMX2s contributed approximately 50% to the binding of insoluble substrates

5

(Fig. 3A).

6

Thermostability of CcCel9A and its CBMX2 mutants-To test whether the CBMX2 is

7

crucial for the thermostability of wild-type CcCel9A, the relative activities of each protein,

8

after incubating at different temperatures for 1 h, were determined (Fig. 3B). At

9

temperatures below 60°C, the relative activity of wild-type CcCel9A was 30%–100%

10

higher than those of the mutant TM14, TM15, and TM16 enzymes. In contrast, above

11

60°C, wild-type CcCel9A exhibited lower activity, and the thermostabilities of mutant

12

enzymes were indistinguishable from one another. These results suggest that the

13

CBMX2 does not play a key role in maintaining the thermostability of CcCel9A.

14

Key amino acid residues in the active-site architecture-A structural model of 24

and C.

15

CcCel9A TM4 was obtained by homologous modeling with TfCel9A

16

cellulolyticum Cel9G

17

catalytic cleft by molecular docking (Fig. 4A). Six glucosyl units, located in the active-

18

site architecture, were categorized as substrate-binding subsites +1/+2 and product-

19

releasing sites −1 to −4, according to their electron density. Aromatic amino acid

20

residues interacted with each catalytic subsite, and several highly conserved residues

21

were selected and ex-changed with alanine in order to remove the functional group (Fig.

22

5 & 4B). Overall, the product-releasing subsite-directed mutants exhibited similar or

23

improved activity on soluble substrates, compared to the wild type. All of the mutants

28

, with a 17-glucosyl unit single cellulose chain embedded in its

ACS Paragon Plus Environment

12

Page 13 of 38 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

Biomacromolecules

1

lost some activity on insoluble substrates, yet the substrate-binding subsite-directed

2

mutants exhibited a greater loss of activity (Table S5). However, an exception to these

3

observations was the specific activity of the −1 subsite-directed mutant on Avicel, which

4

was low, indicating that the −1 subsite was essential for degrading the crystalline

5

substrate. Importantly, though the activity and processivity of mutant Y253A remained at

6

the same level as those of the wild-type enzyme, their values decreased significantly in

7

the double mutant Y253A/D259A, indicating a sophisticated architecture of subsite −4

8

that is crucial for processivity. Analysis of the kinetic parameters revealed that the kcat

9

values of most mutants at subsites −1 to −4 were either improved or equal to that of wild

10

type, owing to their higher product-releasing velocity.

11

The processivity of double mutants at subsites −1 and +1/+2, which nearly lost activity

12

on Avicel, were particularly low, demonstrating that processivity was directly linked to

13

the ability of a cellulase to act on the crystalline substrate 6. The affinity constant for

14

Y253A was the lowest, followed by Y415A, indicating stronger interaction between the

15

substrate and subsites −4 and +1/+2 (Fig. 4D). The interaction forces between the

16

substrate and subsites (Fig.4C) suggest that subsites +1/+2 and −4 bind more tightly

17

than the others, confirming their key role in processivity.

18

Product-releasing

patterns

of

various

mutants

demons-

ACS Paragon Plus Environment

13

Biomacromolecules 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 38

1

trate that those mutated at subsites −1 to −4 were unable to secure relevant lengths of

2

cello-oligosaccharides (Fig. 6). The concentrations of released glucose (G1) and

3

cellotriose (G3) were higher than that of cellobiose (G2) for the subsite −2 mutant,

4

whereas for the subsite −3 mutant, the concentration of cellobiose (G2) was higher than

5

those of glucose (G1) and cellotriose (G3), which showed a similar pattern to that of the

6

subsite −1 mutant (Table S6). In Fig. 6E, no cellotetraose (G4) was detected after

7

hydrolysis by the subsite −4 mutant. As for the substrate-binding site mutants, G4

8

accumulated increasingly with extended reaction times, because it could not be further

9

hydrolyzed by these mutants.

10

Discussion

11

C. cellulosi is a thermophilic cellulolytic microbe that produces a series of extracellular GHs

12

and is capable of fermenting some difficult substrates that other thermophilic Clostridium spp.

13

are unable to utilize, including glycogen, inulin, mannitol, and sucrose, to produce

14

acetylmethylcarbinol

15

cellulosi CS-4-4, comprising a family 9 GH module and five CBMs 32. Recombinant CcCel9A is

16

active on CMC, RAC, Avicel, and xylan at 60°C.

31

. CcCel9A is one of the most abundant GH enzymes secreted by C.

17

Structural studies of processive glycoside hydrolases revealed that the polysaccharide-binding

18

sites form a groove or ring, suggesting that partial or complete substrate enclosure is a

19

prerequisite for processivity 5. Many cellulases contain large grooves, and, in general,

20

exoglucanases have evolved loops that fold over their active sites, resulting in an enclosed

21

substrate-binding tunnel

22

elevated temperatures. Instead, family 9 processive endoglucanases like TfCel9A do not possess

55-57

. However, these loops create structural instability for enzymes at

ACS Paragon Plus Environment

14

Page 15 of 38 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

Biomacromolecules

24

1

mobile loops, but their active sites reside in a long, open cleft

. It is assumed that a cellulose

2

chain binds to the long cleft, where the initial cleavage occurs. CcCel9A, a homologue of

3

TfCel9A, possesses higher processivity at 60°C. The absence of the loops serves to improve

4

thermostability, but in exchange, the basic element for forming a closed tunnel is lost. In order to

5

remedy the loss of affinity to the substrate in the cleft, CcCel9A has developed a unique

6

supramodular structure to achieve processivity on a crystalline substrate.

7

CcCel9A contains an unusually high number of five CBMs, each serving a different function.

8

Molecular docking results show that, in addition to the catalytic subsites +2 to −4, CBM3c

9

extends the catalytic cleft to subsite +13, in order to stabilize the single cellulose chain in the

10

open cleft (Fig. 4A). Intriguingly, several conserved residues in the CBM3c of TfCel9A were

11

also found to be well aligned to interact with the cellulose chain, leading into the catalytic cleft,

12

which would increase the potential binding sites along CBM3c to render the binding in the cleft

13

more favorable 58. Biochemical data revealed that the processivity and activity on Avicel of the

14

CBM3c-deletion mutant (TM5) were undetectable, indicating that CBM3c is crucial for

15

hydrolyzing crystalline substrates. It has been proposed that CBM3c could feed a single cellulose

16

chain into the active-site architecture of TfCel9A 59. It was also confirmed that there was a strong,

17

selective non-covalent binding interaction between GH9 and CBM3c via physical association

18

and co-crystallization of the two modules60-61.

19

The processivity of CcCel9A was improved by 40% after adding a CBMX2 to the

20

GH9/CBM3c unit and by 2.36 times for further addition of CBM3b. Adding a second CBMX2 to

21

these proteins resulted in an increase of 84% and 4.73 times in the absence and the presence of

22

the substrate-targeting CBM3b module, respectively (Fig. 2). Furthermore, the activity on

23

crystalline substrates and processivity of enzymes with the CBM3b, including the wild type and

ACS Paragon Plus Environment

15

Biomacromolecules 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 38

1

mutants TM14, TM15, and TM16, were significantly higher than related enzymes without the

2

CBM3b. The characteristics of TM1 and the wild type revealed that the CBM3b contributed

3

25%–45% toward the maintenance of processivity, demonstrating that the CBM3b is necessary

4

for processivity and hydrolysis of insoluble substrates (Fig. 2). Subfamilies of family 3 CBMs

5

had been verified to play an eminent cellulose-binding role in bacterial cellulases from early on,

6

as CBM3a and 3b were found associated with cellulosomes and/or free cellulases, respectively62-

7

66

8

lacks all three CBMX2, decreased by 70%. However, the processivity of the wild type or

9

mutants lacking one or two CBMX2s were similar, thus indicating little or no cumulative effect

10

of CBMX2 on processivity. Nevertheless, with three CBMX2s a cumulative effect is observed as

11

the binding to Avicel was improved by 64% from TM14 to the wild-type enzyme (Fig. 3A). The

12

gradual increase in Km values of CBM3b-lacking mutants (TM1–TM3) on RAC and Avicel

13

indicates the contribution of CBMX2 to substrate affinity (Table S4).

. The CBMX2 was also necessary for processivity, because the processivity of TM16, which

14

Site-directed mutants designed for each catalytic subsite were constructed and characterized.

15

The data revealed that the substrate-binding subsites play a vital role in processivity. Regardless

16

of specific activity, the processivity or catalytic efficiency of subsite +1/+2 mutants were all

17

decreased significantly, compared to those of the wild type and other subsite-directed mutants

18

(Fig. 5A and Table S5). The activities of the mutants on Avicel as substrate, including

19

Y203A/Y424A, Y253A/D259A, Y415A, H372A/Y415A, H123A, and R374A, decreased more

20

than 90%, while in each case their processivity was almost completely lost, thus revealing the

21

relationship between processivity and ability to degrade crystalline cellulose. Double mutants

22

Y203A/Y424A, Y253A/D258A, and R374A/Y415A of subsites −1, −4, and +1/+2, respectively,

23

had a clear additive effect on the characteristics of CcCel9A, because they comprise highly

ACS Paragon Plus Environment

16

Page 17 of 38 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

Biomacromolecules

1

conserved amino acids residues: Y203, D258, and R374. Most of the amino acid residues that

2

interacted with subsites −2 and −3 were not conserved, except for the aromatic residues such as

3

W257, Y202, and Y315.

4

Interestingly, the specific activity and processivity of M384A was similar or slightly higher

5

than that of the wild type, because of its weak binding to the glucan chain (Table S5). Although

6

the affinity between Y253 and the glucan chain was the strongest among the aromatic amino acid

7

residues distributed in the glucosyl subsites (Fig. 4D), especially those forming hydrogen bonds

8

with subsite −3 and stacking interactions with subsite −4 (Fig. 4B), the characteristics of Y253A

9

were barely altered (Fig.4A & Table S5). It is possible that some water molecules exist within

10

the site of Y253, as the cleft opening was relatively near to the exit, which would help retract the

11

cellulose chain, and therefore, maintain high processivity. The Ka of Y415A also decreased,

12

which highlighted the importance of substrate-binding subsites +1/+2. Taken together with the

13

strong interactive forces observed between the substrate and subsites +1/+2 (Fig. 4C), we can

14

conclude that both subsites +1 and +2 form strong hydrophobic interactions with the cellulose

15

chain and play a role similar to that of W38 and W40, which are located at the entrance of the

16

active-site tunnel of TrCel7A7, 10-11, 67.

17

The hydrolysis product patterns revealed that CcCel9A initially produced cellotetraose (G4)

18

(Fig. 6A). Compared to TrCel7A, four glucosyl units bound more tightly than cellobiose to the

19

enzyme, which prevented the cellulose chain from drawing back. Following this, G4 was

20

subsequently hydrolyzed to glucose (G1), cellobiose (G2), and cellotriose (G3) as secondary

21

reaction products. Additionally, the ability to produce smaller oligosaccharides consistent with

22

subsite −4 to −1 mutants was reduced, whereas G4 was accumulated by subsite +1/+2 mutants

ACS Paragon Plus Environment

17

Biomacromolecules 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 18 of 38

1

(Fig. 6F). Since G4 is a competitive inhibitor of CcCel9A 27, the increasing amounts of G4 could

2

not be further hydrolyzed, and therefore, repressed catalysis. In summary, subsites +1/+2 are

3

instrumental in disrupting the hydrogen bonds between the substrate and CBM, whereas subsites

4

−4 to −1 are responsible for drawing the single chain and moving it through the catalytic cleft.

5

These two processes are cooperative, allowing the enzyme to slide along the cellulose chain

6

processively.

7

Hydrolyzing recalcitrant crystalline cellulose has posed a general bottleneck for cellulase

8

application. Other than the cellobiohydrolases that possess a structural advantage for processivity,

9

the thermophilic processive endoglucanase, CcCel9A, has evolved a novel mode to accomplish

10

this process. In this context, multiple CBMs (3 × CBMX2 + CBM3b) improve the affinity to the

11

flat surface of crystalline cellulosic substrates. In addition, CBM3c extends the binding cleft of

12

GH9 module during decrystallization, and the strong interactive forces between subsites +1/+2

13

and the glucosyl rings of the substrate enable disruption of hydrogen bonds and physical

14

extension of the single cellulose chain. Then the chain is drawn to subsite −4 via interaction

15

between the glucosyl rings and the product-releasing subsites, and finally cellotetraoses are

16

released from the catalytic cleft as an initial product in order to achieve high processivity. In

17

contrast to the tunnel-threading mode of cellobiohydrolases, we propose a “wirewalking mode”

18

as a descriptive term for the GH9 processive endoglucanases that underscores their open cleft-

19

like structure at the active center (presented schematically in Fig. 7).

20

Conclusion

21

CcCel9A recruits multiple modules, whereby the GH9 and CBM3c together form an open

22

cleft and four product-releasing glucosyl sites, consisting of up to 17 glucosyl binding subsites,

ACS Paragon Plus Environment

18

Page 19 of 38 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

Biomacromolecules

1

to achieve processivity. Additional modules coordinate separate functions, whereby CBM3b

2

serves as a targeting agent that also promotes decrystallization of the cellulosic substrate, while

3

the CBMX2s assist in the stripping out of the single cellulose chain from the crystalline region.

4

In the open cleft, the substrate chain is first bound tightly by 13 subsites (i.e., subsites +1 through

5

+13) and is then drawn forwards to the location of subsite −4 by the product-releasing subsites.

6

Consequently, the enzyme “walks” 4 glucosyl units along the cellulose chain, which constitutes a

7

proposed “wirewalking mode” of action. The two sides of the cleavage site work synergistically,

8

producing cellotetraose in assembly-line fashion. Because of the open cleft formed in CcCel9A,

9

cellotetraose is further hydrolyzed to cellotriose, cellobiose, and glucose, by working in a manner

10

similar to that of a typical endoglucanase. Multiple CBMs, the long binding cleft, and larger

11

product size are all unique and advantageous features of endoglucanases to achieve processivity.

12

These findings provide a potential platform for construction of potent industrial cellulases.

ACS Paragon Plus Environment

19

Biomacromolecules 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 20 of 38

1

ASSOCIATED CONTENT

2

Supporting Information. This material is available free of charge via the Internet at

3

http://pubs.acs.org.

4

Figure S1. Amino acids sequence alignment of GH9 cellulases from different strains shown in

5

Figure 1. Table S1. PCR primers used for amplification of truncated mutants of CcCel9A.

6

Underlined portions are overlapping regions for isothermal assembly. Table S2. The amino acids

7

that interacted with the glucosyl subsites in the active-site architecture. Table S3. PCR primers

8

used for construction of site-directed mutants. Table S4. Kinetic parameters of wild-type

9

CcCel9A and its truncated mutants on CMC, RAC and Avicel. Table S5. Characteristics of wild-

10

type CcCel9A and its active-site mutants. Table S6. Gray values of the cellooligosaccharide

11

bands in Figure 6.

12

AUTHOR INFORMATION

13

Corresponding Author

14

* [email protected]

15

* [email protected]

16

ORCID

17

Kun-Di Zhang: 0000-0003-0320-5979

18

Lu-Shan Wang:0000-0002-9625-0359

19

Fu-Li Li: 0000-0003-0381-2234

20

Notes

ACS Paragon Plus Environment

20

Page 21 of 38 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

Biomacromolecules

1

The authors declare no competing financial interest.

2

ACKNOWLEDGMENT

3

We thank Dr. Ely Morag from the Weizmann Institute of Science, Israel, for helpful discussions

4

and suggestions. This work was supported by the National Science Foundation of China

5

(31600050, 31600051, 31770077 and 31770054), the Shandong Province Natural Science Funds

6

for Distinguished Young Scholar (JQ201507), and the Key Scientific and Technological Project

7

of Shandong province (2015ZDXX0403A01).

8

ABBREVIATIONS

9

GH, glycoside hydrolase; CBM, carbohydrate-binding module; TrCel7A, Cel7A from

10

Trichoderma reesei; CcCel9A, Cel9A from Clostridium cellulosi; TfCel9A, Cel9A from

11

Thermobifida fusca; CtCel9I&CtCel9R, Cel9I&Cel9R from Clostridium thermocellum; CpCel9,

12

Cel9 from Clostridium phytofermentans; CMC, carboxymethyl cellulose; RAC, regenerated

13

amorphous cellulose; DNS, 3,5-Dinitrosalicylic acid; TM, truncated mutant; MD, molecular

14

dynamics;

FACE,

fluorescence-assisted

carbohydrate

electrophoresis.

ACS Paragon Plus Environment

21

Biomacromolecules 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

1 2

Page 22 of 38

REFERENCES 1.

Varrot, A.; Hastrup, S.; Schulein, M.; Davies, G. J., Crystal structure of the catalytic core

3

domain of the family 6 cellobiohydrolase II, Cel6A, from Humicola insolens, at 1.92 A

4

resolution. Biochem J 1999, 337 ( Pt 2), 297-304.

5 6 7 8

2.

Von Hippel, P. H.; Fairfield, F. R.; Dolejsi, M. K., On the processivity of polymerases.

Ann N Y Acad Sci 1994, 726, 118-31. 3.

Hemsworth, G. R.; Dejean, G.; Davies, G. J.; Brumer, H., Learning from microbial

strategies for polysaccharide degradation. Biochem Soc Trans 2016, 44 (1), 94-108.

9

4.

Wilson, D. B.; Kostylev, M., Cellulase processivity. Methods Mol Biol 2012, 908, 93-9.

10

5.

Breyer, W. A.; Matthews, B. W., A structural basis for processivity. Protein Sci 2001, 10

11 12

(9), 1699-711. 6.

Payne, C. M.; Jiang, W.; Shirts, M. R.; Himmel, M. E.; Crowley, M. F.; Beckham, G. T.,

13

Glycoside hydrolase processivity is directly related to oligosaccharide binding free energy. J Am

14

Chem Soc 2013, 135 (50), 18831-9.

15 16 17

7.

Mulakala, C.; Reilly, P. J., Hypocrea jecorina (Trichoderma reesei) Cel7A as a molecular

machine: A docking study. Proteins 2005, 60 (4), 598-605. 8.

Cruys-Bagger, N.; Alasepp, K.; Andersen, M.; Ottesen, J.; Borch, K.; Westh, P., Rate of

18

Threading a Cellulose Chain into the Binding Tunnel of a Cellulase. J Phys Chem B 2016, 120

19

(25), 5591-600.

20

9.

Beckham, G. T.; Bomble, Y. J.; Bayer, E. A.; Himmel, M. E.; Crowley, M. F.,

21

Applications of computational science for understanding enzymatic deconstruction of cellulose.

22

Curr Opin Biotechnol 2011, 22 (2), 231-8.

23

10. Igarashi, K.; Koivula, A.; Wada, M.; Kimura, S.; Penttila, M.; Samejima, M., High speed

24

atomic

force

microscopy

visualizes

processive

movement

of

Trichoderma

25

cellobiohydrolase I on crystalline cellulose. J Biol Chem 2009, 284 (52), 36186-90.

reesei

ACS Paragon Plus Environment

22

Page 23 of 38 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

Biomacromolecules

1

11. Kont, R.; Kari, J.; Borch, K.; Westh, P.; Valjamae, P., Inter-domain Synergism Is

2

Required for Efficient Feeding of Cellulose Chain into Active Site of Cellobiohydrolase Cel7A.

3

J Biol Chem 2016, 291 (50), 26013-26023.

4

12. Donohoe, B. S.; Resch, M. G., Mechanisms employed by cellulase systems to gain access

5

through the complex architecture of lignocellulosic substrates. Curr Opin Chem Biol 2015, 29,

6

100-7.

7

13. Bu, L.; Beckham, G. T.; Shirts, M. R.; Nimlos, M. R.; Adney, W. S.; Himmel, M. E.;

8

Crowley, M. F., Probing carbohydrate product expulsion from a processive cellulase with

9

multiple absolute binding free energy methods. J Biol Chem 2011, 286 (20), 18161-9.

10 11

14. Kurasin, M.; Valjamae, P., Processivity of cellobiohydrolases is limited by the substrate. J Biol Chem 2011, 286 (1), 169-77.

12

15. Knott, B. C.; Crowley, M. F.; Himmel, M. E.; Stahlberg, J.; Beckham, G. T.,

13

Carbohydrate-protein interactions that drive processive polysaccharide translocation in enzymes

14

revealed from a computational study of cellobiohydrolase processivity. J Am Chem Soc 2014,

15

136 (24), 8810-9.

16

16. Kuusk, S.; Sorlie, M.; Valjamae, P., The predominant molecular state of bound enzyme

17

determines the strength and type of product inhibition in the hydrolysis of recalcitrant

18

polysaccharides by processive enzymes. J Biol Chem 2015, 290 (18), 11678-91.

19

17. Shang, B. Z.; Chang, R.; Chu, J. W., Systems-level modeling with molecular resolution

20

elucidates the rate-limiting mechanisms of cellulose decomposition by cellobiohydrolases. J Biol

21

Chem 2013, 288 (40), 29081-9.

22

18. Munoz, I. G.; Ubhayasekera, W.; Henriksson, H.; Szabo, I.; Pettersson, G.; Johansson, G.;

23

Mowbray, S. L.; Stahlberg, J., Family 7 cellobiohydrolases from Phanerochaete chrysosporium:

24

crystal structure of the catalytic module of Cel7D (CBH58) at 1.32 A resolution and homology

25

models of the isozymes. J Mol Biol 2001, 314 (5), 1097-111.

ACS Paragon Plus Environment

23

Biomacromolecules 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 24 of 38

1

19. Meinke, A.; Damude, H. G.; Tomme, P.; Kwan, E.; Kilburn, D. G.; Miller, R. C., Jr.;

2

Warren, R. A.; Gilkes, N. R., Enhancement of the endo-beta-1,4-glucanase activity of an

3

exocellobiohydrolase by deletion of a surface loop. J Biol Chem 1995, 270 (9), 4383-6.

4

20. Voutilainen, S. P.; Murray, P. G.; Tuohy, M. G.; Koivula, A., Expression of Talaromyces

5

emersonii cellobiohydrolase Cel7A in Saccharomyces cerevisiae and rational mutagenesis to

6

improve its thermostability and activity. Protein Eng Des Sel 2010, 23 (2), 69-79.

7

21. Baker, A. A.; Helbert, W.; Sugiyama, J.; Miles, M. J., New insight into cellulose

8

structure by atomic force microscopy shows the i(alpha) crystal phase at near-atomic resolution.

9

Biophys J 2000, 79 (2), 1139-45.

10

22. Nishiyama, Y.; Langan, P.; Chanzy, H., Crystal structure and hydrogen-bonding system

11

in cellulose Ibeta from synchrotron X-ray and neutron fiber diffraction. J Am Chem Soc 2002,

12

124 (31), 9074-82.

13

23. Nishiyama, Y.; Sugiyama, J.; Chanzy, H.; Langan, P., Crystal structure and hydrogen

14

bonding system in cellulose I(alpha) from synchrotron X-ray and neutron fiber diffraction. J Am

15

Chem Soc 2003, 125 (47), 14300-6.

16 17 18 19

24. Sakon, J.; Irwin, D.; Wilson, D. B.; Karplus, P. A., Structure and mechanism of endo/exocellulase E4 from Thermomonospora fusca. Nat Struct Biol 1997, 4 (10), 810-8. 25. Zhou, W.; Irwin, D. C.; Escovar-Kousen, J.; Wilson, D. B., Kinetic studies of Thermobifida fusca Cel9A active site mutant enzymes. Biochemistry 2004, 43 (30), 9655-63.

20

26. Li, Y.; Irwin, D. C.; Wilson, D. B., Increased crystalline cellulose activity via

21

combinations of amino acid changes in the family 9 catalytic domain and family 3c cellulose

22

binding module of Thermobifida fusca Cel9A. Appl Environ Microbiol 2010, 76 (8), 2582-8.

23

27. Kostylev, M.; Moran-Mirabal, J. M.; Walker, L. P.; Wilson, D. B., Determination of the

24

molecular states of the processive endocellulase Thermobifida fusca Cel9A during crystalline

25

cellulose depolymerization. Biotechnol Bioeng 2012, 109 (1), 295-9.

ACS Paragon Plus Environment

24

Page 25 of 38 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

Biomacromolecules

1

28. Mandelman, D.; Belaich, A.; Belaich, J. P.; Aghajari, N.; Driguez, H.; Haser, R., X-Ray

2

crystal structure of the multidomain endoglucanase Cel9G from Clostridium cellulolyticum

3

complexed with natural and synthetic cello-oligosaccharides. J Bacteriol 2003, 185 (14), 4127-

4

35.

5

29. Zverlov, V. V.; Velikodvorskaya, G. A.; Schwarz, W. H., Two new cellulosome

6

components encoded downstream of celI in the genome of Clostridium thermocellum: the non-

7

processive endoglucanase CelN and the possibly structural protein CseP. Microbiology 2003,

8

149 (Pt 2), 515-24.

9

30. Zverlov, V. V.; Schantz, N.; Schwarz, W. H., A major new component in the cellulosome

10

of Clostridium thermocellum is a processive endo-beta-1,4-glucanase producing cellotetraose.

11

FEMS Microbiol Lett 2005, 249 (2), 353-8.

12

31. He, Y. L.; Ding, Y. F.; Long, Y. Q., Two cellulolytic Clostridium species: Clostridium

13

cellulosi sp. nov. and Clostridium cellulofermentans sp. nov. Int J Syst Bacteriol 1991, 41 (2),

14

306-9.

15

32. Zhang, K.; Chen, X.; Schwarz, W. H.; Li, F., Synergism of glycoside hydrolase

16

secretomes from two thermophilic bacteria cocultivated on lignocellulose. Appl Environ

17

Microbiol 2014, 80 (8), 2592-601.

18 19

33. Andreou, L. V., Preparation of genomic DNA from bacteria. Methods Enzymol 2013, 529, 143-51.

20

34. Zhang, Y. H.; Hong, J.; Ye, X., Cellulase assays. Methods Mol Biol 2009, 581, 213-31.

21

35. Zhang, X. Z.; Sathitsuksanoh, N.; Zhang, Y. H., Glycoside hydrolase family 9 processive

22

endoglucanase from Clostridium phytofermentans: heterologous expression, characterization,

23

and synergy with family 48 cellobiohydrolase. Bioresour Technol 2010, 101 (14), 5534-8.

24

36. Zhang, Y. H.; Cui, J.; Lynd, L. R.; Kuang, L. R., A transition from cellulose swelling to

25

cellulose dissolution by o-phosphoric acid: evidence from enzymatic hydrolysis and

26

supramolecular structure. Biomacromolecules 2006, 7 (2), 644-8.

ACS Paragon Plus Environment

25

Biomacromolecules 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 26 of 38

1

37. Ozdemir, I.; Blumer-Schuette, S. E.; Kelly, R. M., S-layer homology domain proteins

2

Csac_0678 and Csac_2722 are implicated in plant polysaccharide deconstruction by the

3

extremely thermophilic bacterium Caldicellulosiruptor saccharolyticus. Appl Environ Microbiol

4

2012, 78 (3), 768-77.

5

38. Tian, L.; Liu, S.; Wang, S.; Wang, L., Ligand-binding specificity and promiscuity of the

6

main lignocellulolytic enzyme families as revealed by active-site architecture analysis. Sci Rep

7

2016, 6, 23605.

8

39. Hess, B.; Kutzner, C.; van der Spoel, D.; Lindahl, E., GROMACS 4: Algorithms for

9

Highly Efficient, Load-Balanced, and Scalable Molecular Simulation. J Chem Theory Comput

10

2008, 4 (3), 435-47.

11

40. Hornak, V.; Abel, R.; Okur, A.; Strockbine, B.; Roitberg, A.; Simmerling, C.,

12

Comparison of multiple Amber force fields and development of improved protein backbone

13

parameters. Proteins 2006, 65 (3), 712-25.

14

41. Kirschner, K. N.; Yongye, A. B.; Tschampel, S. M.; Gonzalez-Outeirino, J.; Daniels, C.

15

R.; Foley, B. L.; Woods, R. J., GLYCAM06: a generalizable biomolecular force field.

16

Carbohydrates. J Comput Chem 2008, 29 (4), 622-55.

17

42. Jorgensen, W. L.; Chandrasekhar, J.; Madura, J. D.; Impey, R. W.; Klein, M. L.,

18

Comparison of Simple Potential Functions for Simulating Liquid Water. J Chem Phys 1983, 79

19

(2), 926-935.

20 21 22 23 24 25

43. Nose, S.; Klein, M. L., Constant Pressure Molecular-Dynamics for Molecular-Systems. Mol Phys 1983, 50 (5), 1055-1076. 44. Parrinello, M.; Rahman, A., Polymorphic Transitions in Single-Crystals - a New Molecular-Dynamics Method. J Appl Phys 1981, 52 (12), 7182-7190. 45. Berendsen, H. J. C., Transport-Properties Computed by Linear Response through WeakCoupling to a Bath. Nato Adv Sci I E-App 1991, 205, 139-155.

ACS Paragon Plus Environment

26

Page 27 of 38 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

1 2 3 4

Biomacromolecules

46. Darden, T.; York, D.; Pedersen, L., Particle Mesh Ewald - an N.Log(N) Method for Ewald Sums in Large Systems. J Chem Phys 1993, 98 (12), 10089-10092. 47. Essmann, U.; Perera, L.; Berkowitz, M. L.; Darden, T.; Lee, H.; Pedersen, L. G., A Smooth Particle Mesh Ewald Method. J Chem Phys 1995, 103 (19), 8577-8593.

5

48. Hess, B.; Bekker, H.; Berendsen, H. J. C.; Fraaije, J. G. E. M., LINCS: A linear

6

constraint solver for molecular simulations. Journal of Computational Chemistry 1997, 18 (12),

7

1463-1472.

8 9

49. Miyamoto, S.; Kollman, P. A., Settle - an Analytical Version of the Shake and Rattle Algorithm for Rigid Water Models. Journal of Computational Chemistry 1992, 13 (8), 952-962.

10

50. Zakariassen, H.; Aam, B. B.; Horn, S. J.; Varum, K. M.; Sorlie, M.; Eijsink, V. G.,

11

Aromatic residues in the catalytic center of chitinase A from Serratia marcescens affect

12

processivity, enzyme activity, and biomass converting efficiency. J Biol Chem 2009, 284 (16),

13

10610-7.

14

51. Payne, C. M.; Bomble, Y. J.; Taylor, C. B.; McCabe, C.; Himmel, M. E.; Crowley, M. F.;

15

Beckham, G. T., Multiple functions of aromatic-carbohydrate interactions in a processive

16

cellulase examined with molecular simulation. J Biol Chem 2011, 286 (47), 41028-35.

17

52. Kari, J.; Kont, R.; Borch, K.; Buskov, S.; Olsen, J. P.; Cruyz-Bagger, N.; Valjamae, P.;

18

Westh, P., Anomeric Selectivity and Product Profile of a Processive Cellulase. Biochemistry

19

2017, 56 (1), 167-178.

20

53. Zhang, Q.; Zhang, X.; Wang, P.; Li, D.; Chen, G.; Gao, P.; Wang, L., Determination of

21

the action modes of cellulases from hydrolytic profiles over a time course using fluorescence-

22

assisted carbohydrate electrophoresis. Electrophoresis 2015, 36 (6), 910-7.

23

54. Wang, P. P.; Zhang, L. L.; Zhang, Q.; Wang, L. S.; Chen, G. J., Quantitative

24

Determination of Substrate Binding Affinity Influenced by a Single Amino Acid Mutation in

25

Cellulase Active-site Architecture Using Fluorescent Spectrometry. Prog Biochem Biophys 2016,

26

43 (11), 1086-1093.

ACS Paragon Plus Environment

27

Biomacromolecules 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

1 2

Page 28 of 38

55. Bodenheimer, A. M.; Meilleur, F., Crystal structures of wild-type Trichoderma reesei Cel7A catalytic domain in open and closed states. FEBS Lett 2016, 590 (23), 4429-4438.

3

56. Smith, M. A.; Rentmeister, A.; Snow, C. D.; Wu, T.; Farrow, M. F.; Mingardon, F.;

4

Arnold, F. H., A diverse set of family 48 bacterial glycoside hydrolase cellulases created by

5

structure-guided recombination. FEBS J 2012, 279 (24), 4453-65.

6

57. Silveira, R. L.; Skaf, M. S., Molecular Dynamics Simulations of Family 7

7

Cellobiohydrolase Mutants Aimed at Reducing Product Inhibition. J Phys Chem B 2015, 119

8

(29), 9295-303.

9

58. Irwin, D.; Shin, D. H.; Zhang, S.; Barr, B. K.; Sakon, J.; Karplus, P. A.; Wilson, D. B.,

10

Roles of the catalytic domain and two cellulose binding domains of Thermomonospora fusca E4

11

in cellulose hydrolysis. J Bacteriol 1998, 180 (7), 1709-14.

12 13

59. Oliveira, O. V.; Freitas, L. C.; Straatsma, T. P.; Lins, R. D., Interaction between the CBM of Cel9A from Thermobifida fusca and cellulose fibers. J Mol Recognit 2009, 22 (1), 38-45.

14

60. Burstein, T.; Shulman, M.; Jindou, S.; Petkun, S.; Frolow, F.; Shoham, Y.; Bayer, E. A.;

15

Lamed, R., Physical association of the catalytic and helper modules of a family-9 glycoside

16

hydrolase is essential for activity. FEBS Lett 2009, 583 (5), 879-84.

17

61. Petkun, S.; Rozman Grinberg, I.; Lamed, R.; Jindou, S.; Burstein, T.; Yaniv, O.; Shoham,

18

Y.; Shimon, L. J.; Bayer, E. A.; Frolow, F., Reassembly and co-crystallization of a family 9

19

processive endoglucanase from its component parts: structural and functional significance of the

20

intermodular linker. PeerJ 2015, 3, e1126.

21

62. Poole, D. M.; Morag, E.; Lamed, R.; Bayer, E. A.; Hazlewood, G. P.; Gilbert, H. J.,

22

Identification of the cellulose-binding domain of the cellulosome subunit S1 from Clostridium

23

thermocellum YS. FEMS Microbiol Lett 1992, 78 (2-3), 181-6.

24

63. Morag, E.; Lapidot, A.; Govorko, D.; Lamed, R.; Wilchek, M.; Bayer, E. A.; Shoham, Y.,

25

Expression, purification, and characterization of the cellulose-binding domain of the scaffoldin

26

subunit from the cellulosome of Clostridium thermocellum. Appl Environ Microbiol 1995, 61 (5),

27

1980-6.

ACS Paragon Plus Environment

28

Page 29 of 38 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

Biomacromolecules

1

64. Tormo, J.; Lamed, R.; Chirino, A. J.; Morag, E.; Bayer, E. A.; Shoham, Y.; Steitz, T. A.,

2

Crystal structure of a bacterial family-III cellulose-binding domain: a general mechanism for

3

attachment to cellulose. EMBO J 1996, 15 (21), 5739-51.

4

65. Ye, L.; Su, X.; Schmitz, G. E.; Moon, Y. H.; Zhang, J.; Mackie, R. I.; Cann, I. K.,

5

Molecular and biochemical analyses of the GH44 module of CbMan5B/Cel44A, a bifunctional

6

enzyme from the hyperthermophilic bacterium Caldicellulosiruptor bescii. Appl Environ

7

Microbiol 2012, 78 (19), 7048-59.

8

66. Su, X.; Mackie, R. I.; Cann, I. K., Biochemical and mutational analyses of a multidomain

9

cellulase/mannanase from Caldicellulosiruptor bescii. Appl Environ Microbiol 2012, 78 (7),

10

2230-40.

11

67. Kari, J.; Olsen, J.; Borch, K.; Cruys-Bagger, N.; Jensen, K.; Westh, P., Kinetics of

12

cellobiohydrolase (Cel7A) variants with lowered substrate affinity. J Biol Chem 2014, 289 (47),

13

32459-68.

ACS Paragon Plus Environment

29

Biomacromolecules 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 30 of 38

1 2 3

ACS Paragon Plus Environment

30

Page 31 of 38 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

Biomacromolecules

Figure 1. The schematic structures of GH9+CBM3c-bearing cellulases from different strains. Cc, Clostridium cellulosi, this study; Tf, Thermobifida fusca24; Ct, Clostridium thermocellum29-30; Cp, Clostridium phytofermentans35. SP, signal peptide; GH, glycoside hydrolase; CBM, carbohydrate-binding module; Doc, dockerin. 172x78mm (300 x 300 DPI)

ACS Paragon Plus Environment

Biomacromolecules 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

Figure 2. Schematic modular architectures and characteristics of the wild-type (WT) CcCel9A and its truncated mutants (TMs). Processivities and specific activities towards Avicel and CMC were all shown as % of those of the wild-type enzyme. 173x79mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 32 of 38

Page 33 of 38 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

Biomacromolecules

Figure 3. Characteristics of wild-type and mutants of CcCel9A. Quantitative studies of the binding of CcCel9A and truncated mutants TM14, TM15 and TM16 to Avicel at 60°C. Avicel was mixed with various concentrations of the designated proteins, and the binding activities were estimated as described in Experimental procedures. The graphs showed the binding isotherms between bound proteins (nmol/g of Avicel) and free proteins (A). Thermostability of wild-type (WT) and mutant (TM14, TM15 and TM16) proteins. Proteins were incubated for 1 h at the indicated temperatures, and the residual enzyme activity was assayed. Results represent the means and standard deviations from three independent experiments (B). 76x33mm (300 x 300 DPI)

ACS Paragon Plus Environment

Biomacromolecules 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

Figure 4. Modelled crystal structure of CcCel9A mutant TM4 with the docking of a cellulose chain in the catalytic cleft. Molecular surface of TM4 (A). The megascopic active-site architecture of the GH9 module. The cellulose chain is colored green, and amino acid residues that interact with the glucosyl ring are colored red. The 17 glucosyl binding sites in the cleft, are numbered from subsite −4 to +13, from the non-reducing to the reducing end (A), and subsites −4 to +2 constitute the active-site architecture (B). Interaction forces between substrate and subsites in the active-site architecture by molecular dynamics (C). Affinity constants (Ka) of aromatic amino acid-directed mutants compared to the wild type determined for CMC (D). 173x110mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 34 of 38

Page 35 of 38 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

Biomacromolecules

Figure 5. Processivities (% of wild type) of site-directed mutants (A), related to 6 subsites, and analysis of the conserved sites based on the amino acid sequences of enzymes shown in Fig. 1 using WebLogo (B). CS, cleavage sites. 115x77mm (300 x 300 DPI)

ACS Paragon Plus Environment

Biomacromolecules 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

Figure 6. Hydrolysis patterns of wild-type CcCel9A (A), site-directed mutants of subsites −1 (B), −2 (C), −3 (D), −4 (E) and +1/+2 (F) on RAC. Enzymes (0.3 µM) were incubated with 5 mg/ml RAC at 60°C, and samples were collected at the indicated time intervals. The subsite −1, −2, −3, −4 and +1/+2 mutants represent Y203A/Y424A, Y202A/W257A, Q308A/Y315A, Y253A/D258A and H372A/Y415A, respectively. 173x120mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 36 of 38

Page 37 of 38 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

Biomacromolecules

Figure 7. Working model of the “wirewalking” mode of CcCel9A action. The CBM3b module recognizes the gap on the surface of crystalline cellulose (A); CBM3b binds to the non-reducing end of the gap (B); CBM3b decrystallizes the substrate with the aid of the CBMX2s (C); A single cellulose chain is extracted and fed into the catalytic cleft by CBM3c (D); Subsites +1/+2 enable extension of the cellulose chain into the active-site cleft (E); The non-reducing end of the chain is drawn to the position of subsite −4 (F); Cleavage occurs between subsites +1 and −1 (G); Cellotetraose is produced and released through the exit of cleft (H). 173x116mm (300 x 300 DPI)

ACS Paragon Plus Environment

Biomacromolecules 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

Table of contents use only 70x35mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 38 of 38