The rotamer of the second-sphere histidine in AA9 lytic polysaccharide monooxygenase is pH dependent.

Lytic polysaccharide monooxygenases (LPMOs) catalyze a reaction that is crucial for the biological decomposition of various biopolymers and for the industrial conversion of plant biomass. Despite the importance of LPMOs, the exact molecular-level nature of the reaction mechanism is still debated today. Here, we investigated the pH-dependent conformation of a second-sphere histidine (His) that we call the stacking histidine, which is conserved in fungal AA9 LPMOs and is speculated to assist catalysis in several of the LPMO reaction pathways. Using constant-pH and accelerated molecular dynamics simulations, we monitored the dynamics of the stacking His in different protonation states for both the resting Cu(II) and active Cu(I) forms of two fungal LPMOs. Consistent with experimental crystallographic and neutron diffraction data, our calculations suggest that the side chain of the protonated and positively charged form is rotated out of the active site toward the solvent. Importantly, only one of the possible neutral states of histidine (HIE state) is observed in the stacking orientation at neutral pH or when bound to cellulose. Our data predict that, in solution, the stacking His may act as a stabilizer (via hydrogen bonding) of the Cu(II)-superoxo complex after the LPMO-Cu(I) has reacted with O2 in solution, which, in fine, leads to H2O2 formation. Also, our data indicate that the HIE-stacking His is a poor acid/base catalyst when bound to the substrate and, in agreement with the literature, may play an important stabilizing role (via hydrogen bonding) during the peroxygenase catalysis. Our study reveals the pH titration midpoint values of the pH-dependent orientation of the stacking His should be considered when modeling and interpreting LPMO reactions, whether it be for classical LPMO kinetics or in industry-oriented enzymatic cocktails, and for understanding LPMO behavior in slightly acidic natural processes such as fungal wood decay.

How Do Branched Detergents Stabilize GPCRs in Micelles?

The structural and functional properties of G protein-coupled receptors (GPCRs) are often studied in a detergent micellar environment, but many GPCRs tend to denature or aggregate in short alkyl chain detergents. In our previous work [Lee, S., et al. (2016) J. Am. Chem. Soc. 138, 15425-15433], we showed that GPCRs in alkyl glucosides were highly dynamic, resulting in the penetration of detergent molecules between transmembrane α-helices, which is the initial step in receptor denaturation. Although this was not observed for GPCRs in dodecyl maltoside (DDM, also known as lauryl maltoside), even this detergent is not mild enough to preserve the integrity of many GPCRs during purification. Lauryl maltose neopentylglycol (LMNG) detergents have been found to have significant advantages for purifying GPCRs in a native state as they impart more stability to the receptor than DDM. To gain insights into how they stabilize GPCRs, we used atomistic molecular dynamics simulations of wild type adenosine A2A receptor (WT-A2AR), thermostabilized A2AR (tA2AR), and wild type ß2-adrenoceptor (ß2AR) in a variety of detergents (LMNG, DMNG, OGNG, and DDM). Analysis of molecular dynamics simulations of tA2AR in LMNG, DMNG, and OGNG showed that this series of detergents exhibited behavior very similar to that of an analogous series of detergents DDM, DM, and OG in our previous study. However, there was a striking difference upon comparison of the behavior of LMNG to that of DDM. LMNG showed considerably less motion than DDM, which resulted in the enhanced density of the aliphatic chains around the hydrophobic regions of the receptor and considerably more hydrogen bond formation between the head groups. This contributed to enhanced interaction energies between both detergent molecules and between the receptor and detergent, explaining the enhanced stability of GPCRs purified in this detergent. Branched detergents occlude between transmembrane helices and reduce their flexibility. Our results provide a rational foundation to develop detergent variants for stabilizing membrane proteins.

Prediction of Conformation Specific Thermostabilizing Mutations for Class A G Protein-Coupled Receptors.

G protein-coupled receptors (GPCRs) are highly flexible and prone to denaturation during protein extraction in detergents and purification. This poses a huge challenge to purify a conformationally homogeneous solution of GPCRs. Thermostabilizing mutations have been used widely to purify and obtain crystal structures of several GPCRs. However, identifying thermostabilizing mutations for GPCRs remains a tedious and expensive task as they are not transferable even among closely related GPCRs. Additionally, the mutations stabilizing one conformational state of a GPCR do not always stabilize other conformational state(s) of the same GPCR. Previously we developed a computational method, LiticonDesign, for rapid prediction of thermostabilizing mutations for a specific GPCR conformation. In this study, we have used LiticonDesign to predict thermostabilizing mutations for the agonist bound active-intermediate state of the human adenosine receptor (A2AR) using the structure of the inactive state of the same GPCR and vice versa. Our study shows that the thermostable mutation predictions using LiticonDesign, for an active-intermediate state of a GPCR (A2AR in our case), requires a homology model that is derived from an active/active-intermediate state GPCR structure as a template. Similarly, the homology models derived from inactive state GPCR conformations are better in predicting the thermostable mutations for the inactive state of A2AR. Overall, LiticonDesign method is not only efficient in predicting thermostabilizing mutations for a given GPCR sequence but also can recover conformation specific mutations for a state of interest, if a suitable starting structure of desired conformation is chosen.

Engineering Salt Bridge Networks between Transmembrane Helices Confers Thermostability in G-Protein-Coupled Receptors.

Introduction of specific point mutations has been an effective strategy in enhancing the thermostability of G-protein-coupled receptors (GPCRs). Our previous work showed that a specific residue position on transmembrane helix 6 (TM6) in class A GPCRs consistently yields thermostable mutants. The crystal structure of human chemokine receptor CCR5 also showed increased thermostability upon mutation of two positions, A233D6.33 and K303E7.59. With the goal of testing the transferability of these two thermostabilizing mutations in other chemokine receptors, we tested the mutations A237D6.33 and R307E7.59 in human CCR3 for thermostability and aggregation properties in detergent solution. Interestingly, the double mutant exhibited a 6-10-fold decrease in the aggregation propensity of the wild-type protein. This is in stark contrast to the two single mutants whose aggregation properties resemble the wild type (WT). Moreover, unlike in CCR5, the two single mutants separately showed no increase in thermostability compared to the wild-type CCR3, while the double-mutant A237D6.33/R307E7.59 confers an increase of 2.6 °C in the melting temperature compared to the WT. Extensive all-atom molecular dynamics (MD) simulations in detergent micelles show that a salt bridge network between transmembrane helices TM3, TM6, and TM7 that is absent in the two single mutants confers stability in the double mutant. The free energy surface of the double mutant shows conformational homogeneity compared to the single mutants. An annular n-dodecyl maltoside detergent layer packs tighter to the hydrophobic surface of the double-mutant CCR3 compared to the single mutants providing additional stability. The purification of other C-C chemokine receptors lacking such stabilizing residues may benefit from the incorporation of these two point mutations.

Correlation of structure, function and protein dynamics in GH7 cellobiohydrolases from Trichoderma atroviride, T. reesei and T. harzianum.

BACKGROUND: The ascomycete fungus Trichoderma reesei is the predominant source of enzymes for industrial conversion of lignocellulose. Its glycoside hydrolase family 7 cellobiohydrolase (GH7 CBH) TreCel7A constitutes nearly half of the enzyme cocktail by weight and is the major workhorse in the cellulose hydrolysis process. The orthologs from Trichoderma atroviride (TatCel7A) and Trichoderma harzianum (ThaCel7A) show high sequence identity with TreCel7A, ~ 80%, and represent naturally evolved combinations of cellulose-binding tunnel-enclosing loop motifs, which have been suggested to influence intrinsic cellobiohydrolase properties, such as endo-initiation, processivity, and off-rate. RESULTS: The TatCel7A, ThaCel7A, and TreCel7A enzymes were characterized for comparison of function. The catalytic domain of TatCel7A was crystallized, and two structures were determined: without ligand and with thio-cellotriose in the active site. Initial hydrolysis of bacterial cellulose was faster with TatCel7A than either ThaCel7A or TreCel7A. In synergistic saccharification of pretreated corn stover, both TatCel7A and ThaCel7A were more efficient than TreCel7A, although TatCel7A was more sensitive to thermal inactivation. Structural analyses and molecular dynamics (MD) simulations were performed to elucidate important structure/function correlations. Moreover, reverse conservation analysis (RCA) of sequence diversity revealed divergent regions of interest located outside the cellulose-binding tunnel of Trichoderma spp. GH7 CBHs. CONCLUSIONS: We hypothesize that the combination of loop motifs is the main determinant for the observed differences in Cel7A activity on cellulosic substrates. Fine-tuning of the loop flexibility appears to be an important evolutionary target in Trichoderma spp., a conclusion supported by the RCA data. Our results indicate that, for industrial use, it would be beneficial to combine loop motifs from TatCel7A with the thermostability features of TreCel7A. Furthermore, one region implicated in thermal unfolding is suggested as a primary target for protein engineering.

Improving the thermal stability of cellobiohydrolase Cel7A from Hypocrea jecorina by directed evolution.

Secreted mixtures of Hypocrea jecorina cellulases are able to efficiently degrade cellulosic biomass to fermentable sugars at large, commercially relevant scales. H. jecorina Cel7A, cellobiohydrolase I, from glycoside hydrolase family 7, is the workhorse enzyme of the process. However, the thermal stability of Cel7A limits its use to processes where temperatures are no higher than 50 °C. Enhanced thermal stability is desirable to enable the use of higher processing temperatures and to improve the economic feasibility of industrial biomass conversion. Here, we enhanced the thermal stability of Cel7A through directed evolution. Sites with increased thermal stability properties were combined, and a Cel7A variant (FCA398) was obtained, which exhibited a 10.4 °C increase in Tm and a 44-fold greater half-life compared with the wild-type enzyme. This Cel7A variant contains 18 mutated sites and is active under application conditions up to at least 75 °C. The X-ray crystal structure of the catalytic domain was determined at 2.1 Å resolution and showed that the effects of the mutations are local and do not introduce major backbone conformational changes. Molecular dynamics simulations revealed that the catalytic domain of wild-type Cel7A and the FCA398 variant exhibit similar behavior at 300 K, whereas at elevated temperature (475 and 525 K), the FCA398 variant fluctuates less and maintains more native contacts over time. Combining the structural and dynamic investigations, rationales were developed for the stabilizing effect at many of the mutated sites.

Aromatic-Mediated Carbohydrate Recognition in Processive Serratia marcescens Chitinases.

Microorganisms use a host of enzymes, including processive glycoside hydrolases, to deconstruct recalcitrant polysaccharides to sugars. Processive glycoside hydrolases closely associate with polymer chains and repeatedly cleave glycosidic linkages without dissociating from the crystalline surface after each hydrolytic step; they are typically the most abundant enzymes in both natural secretomes and industrial cocktails by virtue of their significant hydrolytic potential. The ubiquity of aromatic residues lining the enzyme catalytic tunnels and clefts is a notable feature of processive glycoside hydrolases. We hypothesized that these aromatic residues have uniquely defined roles, such as substrate chain acquisition and binding in the catalytic tunnel, that are defined by their local environment and position relative to the substrate and the catalytic center. Here, we investigated this hypothesis with variants of Serratia marcescens family 18 processive chitinases ChiA and ChiB. We applied molecular simulation and free energy calculations to assess active site dynamics and ligand binding free energies. Isothermal titration calorimetry provided further insight into enthalpic and entropic contributions to ligand binding free energy. Thus, the roles of six aromatic residues, Trp-167, Trp-275, and Phe-396 in ChiA, and Trp-97, Trp-220, and Phe-190 in ChiB, have been examined. We observed that point mutation of the tryptophan residues to alanine results in unfavorable changes in the free energy of binding relative to wild-type. The most drastic effects were observed for residues positioned at the "entrances" of the deep substrate-binding clefts and known to be important for processivity. Interestingly, phenylalanine mutations in ChiA and ChiB had little to no effect on chito-oligomer binding, in accordance with the limited effects of their removal on chitinase functionality.

Processivity, Substrate Positioning, and Binding: The Role of Polar Residues in a Family 18 Glycoside Hydrolase.

The enzymatic degradation of recalcitrant polysaccharides such as cellulose (ß-1,4-linked glucose) and chitin (ß-1,4-linked N-acetylglucosamine) by glycoside hydrolases (GHs) is of significant biological and economical importance. In nature, depolymerization is primarily accomplished by processive GHs, which remain attached to the substrate between subsequent hydrolytic reactions. Recent computational efforts have suggested that the processive ability of a GH is directly linked to the ligand binding free energy. The contribution of individual aromatic residues in the active site of these enzymes has been extensively studied. In this study, we offer the first experimental evidence confirming correlation of binding free energy and degree of processivity and evidence that polar residues are essential for maintaining processive ability. Exchanging Thr(276) with Ala in substrate binding subsite -2 in the processive ChiA of Serratia marcescens results in a decrease in both the enthalpy (2.6 and 3.8 kcal/mol) and free energy (0.5 and 2.2 kcal/mol) for the binding to the substrate (GlcNAc)6 and the inhibitor allosamidin, respectively, compared to that of the wild type. Moreover, the initial apparent processivity as measured by [(GlcNAc)2]/[GlcNAc] ratios (17.1 ± 0.4) and chitin degradation efficiency (20%) are greatly reduced for ChiA-T276A versus those of the wild type (30.1 ± 1.5 and 75%, respectively). Mutation of Arg(172) to Ala reduces the level of recognition and positioning of the substrate into the active site. Molecular dynamics simulations indicate ChiA-R172A behaves like the wild type, but the dynamics of ChiA-T276A are greatly influenced by mutation, which is reflective of their influence on processivity.

Sequencing, biochemical characterization, crystal structure and molecular dynamics of cellobiohydrolase Cel7A from Geotrichum candidum 3C.

The ascomycete Geotrichum candidum is a versatile and efficient decay fungus that is involved, for example, in biodeterioration of compact discs; notably, the 3C strain was previously shown to degrade filter paper and cotton more efficiently than several industrial enzyme preparations. Glycoside hydrolase (GH) family 7 cellobiohydrolases (CBHs) are the primary constituents of industrial cellulase cocktails employed in biomass conversion, and feature tunnel-enclosed active sites that enable processive hydrolytic cleavage of cellulose chains. Understanding the structure-function relationships defining the activity and stability of GH7 CBHs is thus of keen interest. Accordingly, we report the comprehensive characterization of the GH7 CBH secreted by G. candidum (GcaCel7A). The bimodular cellulase consists of a family 1 cellulose-binding module (CBM) and linker connected to a GH7 catalytic domain that shares 64% sequence identity with the archetypal industrial GH7 CBH of Hypocrea jecorina (HjeCel7A). GcaCel7A shows activity on Avicel cellulose similar to HjeCel7A, with less product inhibition, but has a lower temperature optimum (50 °C versus 60-65 °C, respectively). Five crystal structures, with and without bound thio-oligosaccharides, show conformational diversity of tunnel-enclosing loops, including a form with partial tunnel collapse at subsite -4 not reported previously in GH7. Also, the first O-glycosylation site in a GH7 crystal structure is reported--on a loop where the glycan probably influences loop contacts across the active site and interactions with the cellulose surface. The GcaCel7A structures indicate higher loop flexibility than HjeCel7A, in accordance with sequence modifications. However, GcaCel7A retains small fluctuations in molecular simulations, suggesting high processivity and low endo-initiation probability, similar to HjeCel7A. DATABASE: Structural data are available in the Protein Data Bank under the accession numbers 5AMP, 4ZZV, 4ZZW, 4ZZT, and 4ZZU. The Geotrichum candidum GH family 7 cellobiohydrolase nucleotide sequence is available in GenBank under accession number KJ958925. ENZYMES: Glycoside hydrolase family 7 reducing end acting cellobiohydrolase.

Thermodynamic Relationships with Processivity in Serratia marcescens Family 18 Chitinases.

The enzymatic degradation of recalcitrant polysaccharides is accomplished by synergistic enzyme cocktails of glycoside hydrolases (GHs) and accessory enzymes. Many GHs are processive which means that they remain attached to the substrate in between subsequent hydrolytic reactions. Chitinases are GHs that catalyze the hydrolysis of chitin (ß-1,4-linked N-acetylglucosamine). Previously, a relationship between active site topology and processivity has been suggested while recent computational efforts have suggested a link between the degree of processivity and ligand binding free energy. We have investigated these relationships by employing computational (molecular dynamics (MD)) and experimental (isothermal titration calorimetry (ITC)) approaches to gain insight into the thermodynamics of substrate binding to Serratia marcescens chitinases ChiA, ChiB, and ChiC. We show that increased processive ability indeed corresponds to more favorable binding free energy and that this likely is a general feature of GHs. Moreover, ligand binding in ChiB is entropically driven; in ChiC it is enthalpically driven, and the enthalpic and entropic contributions to ligand binding in ChiA are equal. Furthermore, water is shown to be especially important in ChiA-binding. This work provides new insight into oligosaccharide binding, getting us one step closer to understand how GHs efficiently degrade recalcitrant polysaccharides.