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.     

Machine learning reveals sequence-function relationships in family 7 glycoside hydrolases

Family 7 glycoside hydrolases (GH7) are among the principal enzymes for cellulose degradation in nature and industrially. These enzymes are often bimodular, including a catalytic domain and carbohydrate-binding module (CBM) attached via a flexible linker, and exhibit an active site that binds cello-oligomers of up to ten glucosyl moieties. GH7 cellulases consist of two major subtypes: cellobiohydrolases (CBH) and endoglucanases (EG). Despite the critical importance of GH7 enzymes, there remain gaps in our understanding of how GH7 sequence and structure relate to function. Here, we employed machine learning to gain data-driven insights into relationships between sequence, structure, and function across the GH7 family. Machine-learning models, trained only on the number of residues in the active-site loops as features, were able to discriminate GH7 CBHs and EGs with up to 99% accuracy, demonstrating that the lengths of loops A4, B2, B3, and B4 strongly correlate with functional subtype across the GH7 family. Classification rules were derived such that specific residues at 42 different sequence positions each predicted the functional subtype with accuracies surpassing 87%. A random forest model trained on residues at 19 positions in the catalytic domain predicted the presence of a CBM with 89.5% accuracy. Our machine learning results recapitulate, as top-performing features, a substantial number of the sequence positions determined by previous experimental studies to play vital roles in GH7 activity. We surmise that the yet-to-be-explored sequence positions among the top-performing features also contribute to GH7 functional variation and may be exploited to understand and manipulate function.     

Biofeedback as a means to alter electromyographic activity in a total knee replacement patient

This paper presents a single case controlled study of a 75-year-old male having bilateral total knee replacement. Baseline EMG recordings demonstrated differential levels of vastus medialis and vastus lateralis muscle activity in both knees during exercise, with increased vastus lateralis activity compared to vastus medialis activity. The purpose of the study was to use electromyographic (EMG) biofeedback training to train the patient to equalize vastus medialis and vastus lateralis EMG activity during exercise. After 11 and 13 training sessions for the left and right knees, respectively, differences between vastus medialis and vastus lateralis activity had markedly decreased. Following the termination of biofeedback training, EMG activity during exercise showed a return toward baseline levels. Several concomitant changes in psychological and physical function were noted. These results suggested that EMG biofeedback can be used to train vastus medialis and vastus lateralis activity in total knee replacement patients, and that biofeedback training may produce positive benefits in other functional areas.The authors wish to express thanks to Dr. Karen Gil for her helpful comments on a draft of this article, and Christianne Herman and Allison Roodman for their help in data collection and entry.     

Status of RNAs, localized in Xenopus laevis oocytes, in the frogs Rana pipiens and Eleutherodactylus coqui

Early development in the frog model, Xenopus laevis, is governed by RNAs, localized to the vegetal cortex of the oocyte. These RNAs include Xdazl RNA, which is involved in primordial germ cell formation, and VegT RNA, which specifies the mesoderm and endoderm. In order to determine whether orthologues of these RNAs are localized and have similar functions in other frogs, we cloned RpDazl and RpVegT from Rana pipiens, a frog that is phylogenetically distant from X. laevis. RNAs from both genes are localized to the vegetal cortex of the R. pipiens oocyte, indicating that the vegetal localization is likely the basal state. The animal location of EcVegT RNA in Eleutherodactylus coqui that we found previously (Beckham et al., 2003) is then a derived state, probably due to the great increase in egg size required for direct development of this species. To answer the question of function, we injected RpVegT or EcVegT RNAs into X. laevis embryos, and assayed animal caps for gene expression. Both of these RNAs induced the expression of endodermal, mesodermal, and organizer genes, showing that the function of RpVegT and EcVegT as meso-endodermal determinants is conserved in frogs. The RNA localizations and the function of VegT orthologues in germ layer specification may be synapomorphies for anuran amphibians.     

The O-glycosylated linker from the Trichoderma reesei Family 7 cellulase is a flexible, disordered protein

Fungi and bacteria secrete glycoprotein cocktails to deconstruct cellulose. Cellulose-degrading enzymes (cellulases) are often modular, with catalytic domains for cellulose hydrolysis and carbohydrate-binding modules connected by linkers rich in serine and threonine with O-glycosylation. Few studies have probed the role that the linker and O-glycans play in catalysis. Since different expression and growth conditions produce different glycosylation patterns that affect enzyme activity, the structure-function relationships that glycosylation imparts to linkers are relevant for understanding cellulase mechanisms. Here, the linker of the Trichoderma reesei Family 7 cellobiohydrolase (Cel7A) is examined by simulation. Our results suggest that the Cel7A linker is an intrinsically disordered protein with and without glycosylation. Contrary to the predominant view, the O-glycosylation does not change the stiffness of the linker, as measured by the relative fluctuations in the end-to-end distance; rather, it provides a 16 Å extension, thus expanding the operating range of Cel7A. We explain observations from previous biochemical experiments in the light of results obtained here, and compare the Cel7A linker with linkers from other cellulases with sequence-based tools to predict disorder. This preliminary screen indicates that linkers from Family 7 enzymes from other genera and other cellulases within T. reesei may not be as disordered, warranting further study.     

Virus-like particles of the Ty3 retrotransposon assemble in association with P-body components

Retroviruses and retrotransposons assemble intracellular immature core particles around a RNA genome, and nascent particles collect in association with membranes or as intracellular clusters. How and where genomic RNA are identified for retrovirus and retrotransposon assembly, and how translation and assembly processes are coordinated is poorly understood. To understand this process, the subcellular localization of Ty3 RNA and capsid proteins and virus-like particles was investigated. We demonstrate that mRNAs, proteins, and virus-like particles of the yeast Ty3 retrotransposon accumulate in association with cytoplasmic P-bodies, which are sites of mRNA translation repression, storage, and degradation. Deletions of genes encoding P-body proteins decreased Ty3 transposition and caused changes in the pattern of Ty3 foci, underscoring the biological significance of the association of Ty3 virus-like protein components and P-bodies. These results suggest the hypothesis that P-bodies may serve to segregate translation and assembly functions of the Ty3 genomic RNA to promote assembly of virus-like particles. Because Ty3 has features of a simple retrovirus and P-body functions are conserved between yeast and metazoan organisms, these findings may provide insights into host factors that facilitate retrovirus assembly.     

Binding preferences, surface attachment, diffusivity, and orientation of a family 1 carbohydrate-binding module on cellulose

Cellulase enzymes often contain carbohydrate-binding modules (CBMs) for binding to cellulose. The mechanisms by which CBMs recognize specific surfaces of cellulose and aid in deconstruction are essential to understand cellulase action. The Family 1 CBM from the Trichoderma reesei Family 7 cellobiohydrolase, Cel7A, is known to selectively bind to hydrophobic surfaces of native cellulose. It is most commonly suggested that three aromatic residues identify the planar binding face of this CBM, but several recent studies have challenged this hypothesis. Here, we use molecular simulation to study the CBM binding orientation and affinity on hydrophilic and hydrophobic cellulose surfaces. Roughly 43 μs of molecular dynamics simulations were conducted, which enables statistically significant observations. We quantify the fractions of the CBMs that detach from crystal surfaces or diffuse to other surfaces, the diffusivity along the hydrophobic surface, and the overall orientation of the CBM on both hydrophobic and hydrophilic faces. The simulations demonstrate that there is a thermodynamic driving force for the Cel7A CBM to bind preferentially to the hydrophobic surface of cellulose relative to hydrophilic surfaces. In addition, the simulations demonstrate that the CBM can diffuse from hydrophilic surfaces to the hydrophobic surface, whereas the reverse transition is not observed. Lastly, our simulations suggest that the flat faces of Family 1 CBMs are the preferred binding surfaces. These results enhance our understanding of how Family 1 CBMs interact with and recognize specific cellulose surfaces and provide insights into the initial events of cellulase adsorption and diffusion on cellulose.     

Computational investigation of glycosylation effects on a family 1 carbohydrate-binding module

Carbohydrate-binding modules (CBMs) are ubiquitous components of glycoside hydrolases, which degrade polysaccharides in nature. CBMs target specific polysaccharides, and CBM binding affinity to cellulose is known to be proportional to cellulase activity, such that increasing binding affinity is an important component of performance improvement. To ascertain the impact of protein and glycan engineering on CBM binding, we use molecular simulation to quantify cellulose binding of a natively glycosylated Family 1 CBM. To validate our approach, we first examine aromatic-carbohydrate interactions on binding, and our predictions are consistent with previous experiments, showing that a tyrosine to tryptophan mutation yields a 2-fold improvement in binding affinity. We then demonstrate that enhanced binding of 3-6-fold over a nonglycosylated CBM is achieved by the addition of a single, native mannose or a mannose dimer, respectively, which has not been considered previously. Furthermore, we show that the addition of a single, artificial glycan on the anterior of the CBM, with the native, posterior glycans also present, can have a dramatic impact on binding affinity in our model, increasing it up to 140-fold relative to the nonglycosylated CBM. These results suggest new directions in protein engineering, in that modifying glycosylation patterns via heterologous expression, manipulation of culture conditions, or introduction of artificial glycosylation sites, can alter CBM binding affinity to carbohydrates and may thus be a general strategy to enhance cellulase performance. Our results also suggest that CBM binding studies should consider the effects of glycosylation on binding and function.     

Structural characterisation of Tpx from Yersinia pseudotuberculosis reveals insights into the binding of salicylidene acylhydrazide compounds

Thiol peroxidase, Tpx, has been shown to be a target protein of the salicylidene acylhydrazide class of antivirulence compounds. In this study we present the crystal structures of Tpx from Y. pseudotuberculosis (ypTpx) in the oxidised and reduced states, together with the structure of the C61S mutant. The structures solved are consistent with previously solved atypical 2-Cys thiol peroxidases, including that for "forced" reduced states using the C61S mutant. In addition, by investigating the solution structure of ypTpx using small angle X-ray scattering (SAXS), we have confirmed that reduced state ypTpx in solution is a homodimer. The solution structure also reveals flexibility around the dimer interface. Notably, the conformational changes observed between the redox states at the catalytic triad and at the dimer interface have implications for substrate and inhibitor binding. The structural data were used to model the binding of two salicylidene acylhydrazide compounds to the oxidised structure of ypTpx. Overall, the study provides insights into the binding of the salicylidene acylhydrazides to ypTpx, aiding our long-term strategy to understand the mode of action of this class of compounds.