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.     

Inosine-Mediated Modulation of RNA Sensing by Toll-Like Receptor 7 (TLR7) and TLR8

RNA-specific adenosine deaminase (ADAR)-mediated adenosine-to-inosine (A-to-I) editing is a critical arm of the antiviral response. However, mechanistic insights into how A-to-I RNA editing affects viral infection are lacking. We posited that inosine incorporation into RNA facilitates sensing of nonself RNA by innate immune sensors and accordingly investigated the impact of inosine-modified RNA on Toll-like receptor 7 and 8 (TLR7/8) sensing. Inosine incorporation into synthetic single-stranded RNA (ssRNA) potentiated tumor necrosis factor alpha (TNF-α) or alpha interferon (IFN-α) production in human peripheral blood mononuclear cells (PBMCs) in a sequence-dependent manner, indicative of TLR7/8 recruitment. The effect of inosine incorporation on TLR7/8 sensing was restricted to immunostimulatory ssRNAs and was not seen with inosine-containing short double-stranded RNAs or with a deoxy-inosine-modified ssRNA. Inosine-mediated increase of self-secondary structure of an ssRNA resulted in potentiated IFN-α production in human PBMCs through TLR7 recruitment, as established through the use of a TLR7 antagonist and Tlr7-deficient cells. There was a correlation between hyperediting of influenza A viral ssRNA and its ability to stimulate TNF-α, independent of 5′-triphosphate residues, and involving Adar-1. Furthermore, A-to-I editing of viral ssRNA directly enhanced mouse Tlr7 sensing, when present in proportions reproducing biologically relevant levels of RNA editing. Thus, we demonstrate for the first time that inosine incorporation into immunostimulatory ssRNA can potentiate TLR7/8 activation. Our results suggest a novel function of A-to-I RNA editing, which is to facilitate TLR7/8 sensing of phagocytosed viral RNA.     

Metabolic engineering of Pseudomonas putida for increased polyhydroxyalkanoate production from lignin

Microbial conversion offers a promising strategy for overcoming the intrinsic heterogeneity of the plant biopolymer, lignin. Soil microbes that natively harbour aromatic-catabolic pathways are natural choices for chassis strains, and Pseudomonas putida KT2440 has emerged as a viable whole-cell biocatalyst for funnelling lignin-derived compounds to value-added products, including its native carbon storage product, medium-chain-length polyhydroxyalkanoates (mcl-PHA). In this work, a series of metabolic engineering targets to improve mcl-PHA production are combined in the P. putida chromosome and evaluated in strains growing in a model aromatic compound, p-coumaric acid, and in lignin streams. Specifically, the PHA depolymerase gene phaZ was knocked out, and the genes involved in β-oxidation (fadBA1 and fadBA2) were deleted. Additionally, to increase carbon flux into mcl-PHA biosynthesis, phaG, alkK, phaC1 and phaC2 were overexpressed. The best performing strain – which contains all the genetic modifications detailed above – demonstrated a 53% and 200% increase in mcl-PHA titre (g l⁻¹) and a 20% and 100% increase in yield (g mcl-PHA per g cell dry weight) from p-coumaric acid and lignin, respectively, compared with the wild type strain. Overall, these results present a promising strain to be employed in further process development for enhancing mcl-PHA production from aromatic compounds and lignin.     

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.     

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.