An iterative computational design approach to increase the thermal endurance of a mesophilic enzyme.

BACKGROUND: Strategies for maximizing the microbial production of bio-based chemicals and fuels include eliminating branched points to streamline metabolic pathways. While this is often achieved by removing key enzymes, the introduction of nonnative enzymes can provide metabolic shortcuts, bypassing branched points to decrease the production of undesired side-products. Pyruvate decarboxylase (PDC) can provide such a shortcut in industrially promising thermophilic organisms; yet to date, this enzyme has not been found in any thermophilic organism. Incorporating nonnative enzymes into host organisms can be challenging in cases such as this, where the enzyme has evolved in a very different environment from that of the host. RESULTS: In this study, we use computational protein design to engineer the Zymomonas mobilis PDC to resist thermal denaturation at the growth temperature of a thermophilic host. We generate thirteen PDC variants using the Rosetta protein design software. We measure thermal stability of the wild-type PDC and PDC variants using circular dichroism. We then measure and compare enzyme endurance for wild-type PDC with the PDC variants at an elevated temperature of 60 °C (thermal endurance) using differential interference contrast imaging. CONCLUSIONS: We find that increases in melting temperature (Tm) do not directly correlate with increases in thermal endurance at 60 °C. We also do not find evidence that any individual mutation or design approach is the major contributor to the most thermostable PDC variant. Rather, remarkable cooperativity among sixteen thermostabilizing mutations is key to rationally designing a PDC with significantly enhanced thermal endurance. These results suggest a generalizable iterative computational protein design approach to improve thermal stability and endurance of target enzymes.

Natural diversity of glycoside hydrolase family 48 exoglucanases: insights from structure.

Glycoside hydrolase (GH) family 48 is an understudied and increasingly important exoglucanase family found in the majority of bacterial cellulase systems. Moreover, many thermophilic enzyme systems contain GH48 enzymes. Deletion of GH48 enzymes in these microorganisms results in drastic reduction in biomass deconstruction. Surprisingly, given their importance for these microorganisms, GH48s have intrinsically low cellulolytic activity but even in low ratios synergize greatly with GH9 endoglucanases. In this study, we explore the structural and enzymatic diversity of these enzymes across a wide range of temperature optima. We have crystallized one new GH48 module from Bacillus pumilus in a complex with cellobiose and cellohexaose (BpumGH48). We compare this structure to other known GH48 enzymes in an attempt to understand GH48 structure/function relationships and draw general rules correlating amino acid sequences and secondary structures to thermostability in this GH family.

Development of an Optical Zn2+ Probe Based on a Single Fluorescent Protein.

Various fluorescent probes have been developed to reveal the biological functions of intracellular labile Zn2+. Here, we present Green Zinc Probe (GZnP), a novel genetically encoded Zn2+ sensor design based on a single fluorescent protein (single-FP). The GZnP sensor is generated by attaching two zinc fingers (ZF) of the transcription factor Zap1 (ZF1 and ZF2) to the two ends of a circularly permuted green fluorescent protein (cpGFP). Formation of ZF folds induces interaction between the two ZFs, which induces a change in the cpGFP conformation, leading to an increase in fluorescence. A small sensor library is created to include mutations in the ZFs, cpGFP and linkers between ZF and cpGFP to improve signal stability, sensor brightness and dynamic range based on rational protein engineering, and computational design by Rosetta. Using a cell-based library screen, we identify sensor GZnP1, which demonstrates a stable maximum signal, decent brightness (QY = 0.42 at apo state), as well as specific and sensitive response to Zn2+ in HeLa cells (Fmax/Fmin = 2.6, Kd = 58 pM, pH 7.4). The subcellular localizing sensors mito-GZnP1 (in mitochondria matrix) and Lck-GZnP1 (on plasma membrane) display sensitivity to Zn2+ (Fmax/Fmin = 2.2). This sensor design provides freedom to be used in combination with other optical indicators and optogenetic tools for simultaneous imaging and advancing our understanding of cellular Zn2+ function.

Comparing Residue Clusters from Thermophilic and Mesophilic Enzymes Reveals Adaptive Mechanisms.

Understanding how proteins adapt to function at high temperatures is important for deciphering the energetics that dictate protein stability and folding. While multiple principles important for thermostability have been identified, we lack a unified understanding of how internal protein structural and chemical environment determine qualitative or quantitative impact of evolutionary mutations. In this work we compare equivalent clusters of spatially neighboring residues between paired thermophilic and mesophilic homologues to evaluate adaptations under the selective pressure of high temperature. We find the residue clusters in thermophilic enzymes generally display improved atomic packing compared to mesophilic enzymes, in agreement with previous research. Unlike residue clusters from mesophilic enzymes, however, thermophilic residue clusters do not have significant cavities. In addition, anchor residues found in many clusters are highly conserved with respect to atomic packing between both thermophilic and mesophilic enzymes. Thus the improvements in atomic packing observed in thermophilic homologues are not derived from these anchor residues but from neighboring positions, which may serve to expand optimized protein core regions.

New perspective on glycoside hydrolase binding to lignin from pretreated corn stover.

BACKGROUND: Non-specific binding of cellulases to lignin has been implicated as a major factor in the loss of cellulase activity during biomass conversion to sugars. It is believed that this binding may strongly impact process economics through loss of enzyme activities during hydrolysis and enzyme recycling scenarios. The current model suggests glycoside hydrolase activities are lost though non-specific/non-productive binding of carbohydrate-binding domains to lignin, limiting catalytic site access to the carbohydrate components of the cell wall. RESULTS: In this study, we have compared component enzyme affinities of a commercial Trichoderma reesei cellulase formulation, Cellic CTec2, towards extracted corn stover lignin using sodium dodecyl sulfate-polyacrylamide gel electrophoresis and p-nitrophenyl substrate activities to monitor component binding, activity loss, and total protein binding. Protein binding was strongly affected by pH and ionic strength. ß-d-glucosidases and xylanases, which do not have carbohydrate-binding modules (CBMs) and are basic proteins, demonstrated the strongest binding at low ionic strength, suggesting that CBMs are not the dominant factor in enzyme adsorption to lignin. Despite strong adsorption to insoluble lignin, ß-d-glucosidase and xylanase activities remained high, with process yields decreasing only 4-15 % depending on lignin concentration. CONCLUSION: We propose that specific enzyme adsorption to lignin from a mixture of biomass-hydrolyzing enzymes is a competitive affinity where ß-d-glucosidases and xylanases can displace CBM interactions with lignin. Process parameters, such as temperature, pH, and salt concentration influence the individual enzymes' affinity for lignin, and both hydrophobic and electrostatic interactions are responsible for this binding phenomenon. Moreover, our results suggest that concern regarding loss of critical cell wall degrading enzymes to lignin adsorption may be unwarranted when complex enzyme mixtures are used to digest biomass.

The catalytic mechanism and unique low pH optimum of Caldicellulosiruptor bescii family 3 pectate lyase.

The unique active site of the Caldicellulosiruptor bescii family 3 pectate lyase (PL3) enzyme has been thoroughly characterized using a series of point mutations, X-ray crystallography, pK(a) calculations and biochemical assays. The X-ray structures of seven PL3 active-site mutants, five of them in complex with intact trigalacturonic acid, were solved and characterized structurally, biochemically and computationally. The results confirmed that Lys108 is the catalytic base, but there is no clear candidate for the catalytic acid. However, the reaction mechanism can also be explained by an antiperiplanar trans-elimination reaction, in which Lys108 abstracts a proton from the C5 atom without the help of simultaneous proton donation by an acidic residue. An acidified water molecule completes the anti ß-elimination reaction by protonating the O4 atom of the substrate. Both the C5 hydrogen and C4 hydroxyl groups of the substrate must be orientated in axial configurations, as for galacturonic acid, for this to be possible. The wild-type C. bescii PL3 displays a pH optimum that is lower than that of Bacillus subtilis PL1 according to activity measurements, indicating that C. bescii PL3 has acquired a lower pH optimum by utilizing lysine instead of arginine as the catalytic base, as well as by lowering the pK(a) of the catalytic base in a unique active-site environment.

Predicting enzyme adsorption to lignin films by calculating enzyme surface hydrophobicity.

The inhibitory action of lignin on cellulase cocktails is a major challenge to the biological saccharification of plant cell wall polysaccharides. Although the mechanism remains unclear, hydrophobic interactions between enzymes and lignin are hypothesized to drive adsorption. Here we evaluate the role of hydrophobic interactions in enzyme-lignin binding. The hydrophobicity of the enzyme surface was quantified using an estimation of the clustering of nonpolar atoms, identifying potential interaction sites. The adsorption of enzymes to lignin surfaces, measured using the quartz crystal microbalance, correlates to the hydrophobic cluster scores. Further, these results suggest a minimum hydrophobic cluster size for a protein to preferentially adsorb to lignin. The impact of electrostatic contribution was ruled out by comparing the isoelectric point (pI) values to the adsorption of proteins to lignin surfaces. These results demonstrate the ability to predict enzyme-lignin adsorption and could potentially be used to design improved cellulase cocktails, thus lowering the overall cost of biofuel production.

Computationally designed peptide inhibitors of the ubiquitin E3 ligase SCF(Fbx4).

A structure-based computational approach was used to rationally design peptide inhibitors that can target an E3 ligase (SCF(Fbx4) )-substrate (TRF1) interface and subsequent ubiquitylation. Characterization of the inhibitors demonstrates that our sequence-optimization protocol results in an increase in peptide-TRF1 affinity without compromising peptide-protein specificity.

Cellulase linkers are optimized based on domain type and function: insights from sequence analysis, biophysical measurements, and molecular simulation.

Cellulase enzymes deconstruct cellulose to glucose, and are often comprised of glycosylated linkers connecting glycoside hydrolases (GHs) to carbohydrate-binding modules (CBMs). Although linker modifications can alter cellulase activity, the functional role of linkers beyond domain connectivity remains unknown. Here we investigate cellulase linkers connecting GH Family 6 or 7 catalytic domains to Family 1 or 2 CBMs, from both bacterial and eukaryotic cellulases to identify conserved characteristics potentially related to function. Sequence analysis suggests that the linker lengths between structured domains are optimized based on the GH domain and CBM type, such that linker length may be important for activity. Longer linkers are observed in eukaryotic GH Family 6 cellulases compared to GH Family 7 cellulases. Bacterial GH Family 6 cellulases are found with structured domains in either N to C terminal order, and similar linker lengths suggest there is no effect of domain order on length. O-glycosylation is uniformly distributed across linkers, suggesting that glycans are required along entire linker lengths for proteolysis protection and, as suggested by simulation, for extension. Sequence comparisons show that proline content for bacterial linkers is more than double that observed in eukaryotic linkers, but with fewer putative O-glycan sites, suggesting alternative methods for extension. Conversely, near linker termini where linkers connect to structured domains, O-glycosylation sites are observed less frequently, whereas glycines are more prevalent, suggesting the need for flexibility to achieve proper domain orientations. Putative N-glycosylation sites are quite rare in cellulase linkers, while an N-P motif, which strongly disfavors the attachment of N-glycans, is commonly observed. These results suggest that linkers exhibit features that are likely tailored for optimal function, despite possessing low sequence identity. This study suggests that cellulase linkers may exhibit function in enzyme action, and highlights the need for additional studies to elucidate cellulase linker functions.

Transmembrane peptides used to investigate the homo-oligomeric interface and binding hotspot of latent membrane protein 1.

Epstein-Barr virus (EBV), a human γ-herpesvirus, establishes lifelong infection by targeting the adaptive immune system of the host through memory B cells. Although normally benign, EBV contributes to lymphoid malignancies and lymphoproliferative syndromes in immunocompromised individuals. The viral oncoprotein latent membrane protein 1 (LMP-1) is essential for B lymphocyte immortalization by EBV. The constitutive signaling activity of LMP-1 is dependent on homo-oligomerization of its six-spanning hydrophobic transmembrane domain (TMD). However, the mechanism driving LMP-1 intermolecular interaction is poorly understood. Here, we show that the fifth transmembrane helix (TM5) of LMP-1 strongly self-associates, forming a homotrimeric complex mediated by a polar residue embedded in the membrane, D150. Replacement of this aspartic acid residue with alanine disrupts TM5 self-association in detergent micelles and bacterial cell membranes. A full-length LMP-1 variant harboring the D150A substitution is deficient in NFκB activation, supporting the key role of the fifth transmembrane helix in constitutive activation of signaling by this oncoprotein.

Computational design of the sequence and structure of a protein-binding peptide.

The de novo design of protein-binding peptides is challenging because it requires the identification of both a sequence and a backbone conformation favorable for binding. We used a computational strategy that iterates between structure and sequence optimization to redesign the C-terminal portion of the RGS14 GoLoco motif peptide so that it adopts a new conformation when bound to Gα(i1). An X-ray crystal structure of the redesigned complex closely matches the computational model, with a backbone root-mean-square deviation of 1.1 Å.

Structural determinants of affinity enhancement between GoLoco motifs and G-protein alpha subunit mutants.

GoLoco motif proteins bind to the inhibitory G(i) subclass of G-protein α subunits and slow the release of bound GDP; this interaction is considered critical to asymmetric cell division and neuro-epithelium and epithelial progenitor differentiation. To provide protein tools for interrogating the precise cellular role(s) of GoLoco motif/Gα(i) complexes, we have employed structure-based protein design strategies to predict gain-of-function mutations that increase GoLoco motif binding affinity. Here, we describe fluorescence polarization and isothermal titration calorimetry measurements showing three predicted Gα(i1) point mutations, E116L, Q147L, and E245L; each increases affinity for multiple GoLoco motifs. A component of this affinity enhancement results from a decreased rate of dissociation between the Gα mutants and GoLoco motifs. For Gα(i1)(Q147L), affinity enhancement was seen to be driven by favorable changes in binding enthalpy, despite reduced contributions from binding entropy. The crystal structure of Gα(i1)(Q147L) bound to the RGS14 GoLoco motif revealed disorder among three peptide residues surrounding a well defined Leu-147 side chain. Monte Carlo simulations of the peptide in this region showed a sampling of multiple backbone conformations in contrast to the wild-type complex. We conclude that mutation of Glu-147 to leucine creates a hydrophobic surface favorably buried upon GoLoco peptide binding, yet the hydrophobic Leu-147 also promotes flexibility among residues 511-513 of the RGS14 GoLoco peptide.

Computational design of second-site suppressor mutations at protein-protein interfaces.

The importance of a protein-protein interaction to a signaling pathway can be established by showing that amino acid mutations that weaken the interaction disrupt signaling, and that additional mutations that rescue the interaction recover signaling. Identifying rescue mutations, often referred to as second-site suppressor mutations, controls against scenarios in which the initial deleterious mutation inactivates the protein or disrupts alternative protein-protein interactions. Here, we test a structure-based protocol for identifying second-site suppressor mutations that is based on a strategy previously described by Kortemme and Baker. The molecular modeling software Rosetta is used to scan an interface for point mutations that are predicted to weaken binding but can be rescued by mutations on the partner protein. The protocol typically identifies three types of specificity switches: knob-in-to-hole redesigns, switching hydrophobic interactions to hydrogen bond interactions, and replacing polar interactions with nonpolar interactions. Computational predictions were tested with two separate protein complexes; the G-protein Galpha(i1) bound to the RGS14 GoLoco motif, and UbcH7 bound to the ubiquitin ligase E6AP. Eight designs were experimentally tested. Swapping a buried hydrophobic residue with a polar residue dramatically weakened binding affinities. In none of these cases were we able to identify compensating mutations that returned binding to wild-type affinity, highlighting the challenges inherent in designing buried hydrogen bond networks. The strongest specificity switches were a knob-in-to-hole design (20-fold) and the replacement of a charge-charge interaction with nonpolar interactions (55-fold). In two cases, specificity was further tuned by including mutations distant from the initial design. Proteins 2010. (c) 2009 Wiley-Liss, Inc.

A peptide antagonist of the TLR4-MD2 interaction.

Toll-like receptors are an integral part of innate immunity in the central nervous system (CNS); they orchestrate a robust defense in response to both exogenous and endogenous danger signals. Recently, toll-like receptor 4 (TLR4) has emerged as a therapeutic target for the treatment of CNS-related diseases such as sepsis and chronic pain. We herein report a chemical biology approach by using a rationally designed peptide inhibitor to disrupt the TLR4-MD2 association, thereby blocking TLR4 signaling.

Structure-based protocol for identifying mutations that enhance protein-protein binding affinities.

The ability to manipulate protein binding affinities is important for the development of proteins as biosensors, industrial reagents, and therapeutics. We have developed a structure-based method to rationally predict single mutations at protein-protein interfaces that enhance binding affinities. The protocol is based on the premise that increasing buried hydrophobic surface area and/or reducing buried hydrophilic surface area will generally lead to enhanced affinity if large steric clashes are not introduced and buried polar groups are not left without a hydrogen bond partner. The procedure selects affinity enhancing point mutations at the protein-protein interface using three criteria: (1) the mutation must be from a polar amino acid to a non-polar amino acid or from a non-polar amino acid to a larger non-polar amino acid, (2) the free energy of binding as calculated with the Rosetta protein modeling program should be more favorable than the free energy of binding calculated for the wild-type complex and (3) the mutation should not be predicted to significantly destabilize the monomers. The performance of the computational protocol was experimentally tested on two separate protein complexes; Galpha(i1) from the heterotrimeric G-protein system bound to the RGS14 GoLoco motif, and the E2, UbcH7, bound to the E3, E6AP from the ubiquitin pathway. Twelve single-site mutations that were predicted to be stabilizing were synthesized and characterized in the laboratory. Nine of the 12 mutations successfully increased binding affinity with five of these increasing binding by over 1.0 kcal/mol. To further assess our approach we searched the literature for point mutations that pass our criteria and have experimentally determined binding affinities. Of the eight mutations identified, five were accurately predicted to increase binding affinity, further validating the method as a useful tool to increase protein-protein binding affinities.

Normal metabolism but different physical properties of myelin from mice deficient in proteolipid protein.

Proteolipid protein (PLP) is the primary protein component of CNS myelin, yet myelin from the PLP(null) mouse has only minor ultrastructural abnormalities. Might compensation for a potentially unstable structure involve increased myelin synthesis and turnover? This was not the case; neither accumulation nor in vivo synthesis rates for the myelin-specific lipid cerebroside was altered in PLP(null) mice relative to wild-type (wt) animals. However, the yield of myelin from PLP(null) mice, assayed as levels of cerebroside, was only about 55% of wt control levels. Loss of myelin occurred during initial centrifugation of brain homogenate at 20,000g for 20 min, which is sufficient to sediment almost all myelin from wt mice. Cerebroside-containing fragments from PLP(null) mice remaining in the supernatant could be sedimented by more stringent centrifugation, 100,000g for 60 min. Both the rapidly and the more slowly sedimenting cerebroside-containing membranes banded at the 0.85/0.32 M sucrose interface of a density gradient, as did myelin from wt mice. These results suggest at least some myelin from PLP(null) mice differs from wt myelin with respect to physical stability (fragmented into smaller particles during dispersion) and/or density. Alternatively, slowly sedimenting cerebroside-containing particles could be myelin precursor membranes that, lacking PLP, were retarded in their processing toward mature myelin and thus differ from mature myelin in physical properties. If this is so, recently synthesized cerebroside should be preferentially found in these "slower-sedimenting" myelin precursor fragments. Metabolic tracer experiments showed this was not the case. We conclude that PLP(null) myelin is physically less stable and/or less dense than wt myelin.

Alterations in metabolism and gene expression in brain regions during cuprizone-induced demyelination and remyelination.

Exposure of mice to the copper chelator, cuprizone, results in CNS demyelination. There is remyelination after removal of the metabolic insult. We present brain regional studies identifying corpus callosum as particularly severely affected; 65% of cerebroside is lost after 6 weeks of exposure. We examined recovery of cerebroside and ability to synthesize cerebroside and cholesterol following removal of the toxicant. The temporal pattern for concentration of myelin basic protein resembled that of cerebroside. We applied Affymetrix GeneChip technology to corpus callosum to identify temporal changes in levels of mRNAs during demyelination and remyelination. Genes coding for myelin structural components were greatly down-regulated during demyelination and up-regulated during remyelination. Genes related to microglia/macrophages appeared in a time-course (peaking at 6 weeks) correlating with phagocytosis of myelin and repair of lesions. mRNAs coding for many cytokines had peak expression at 4 weeks, compatible with intercellular signaling roles. Of interest were other genes with temporal patterns correlating with one of the three above patterns, but of function not obviously related to demyelination/remyelination. The ability to correlate gene expression with known pathophysiological events should help in elucidating further function of such genes as related to demyelination/remyelination.