Ugur Sezerman

Improved Thermostability of Clostridium thermocellum Endoglucanase Cel8A by Using Consensus-Guided Mutagenesis

ABSTRACT The use of thermostable cellulases is advantageous for the breakdown of lignocellulosic biomass toward the commercial production of biofuels. Previously, we have demonstrated the engineering of an enhanced thermostable family 8 cellulosomal endoglucanase (EC 3.2.1.4), Cel8A, from Clostridium thermocellum, using random error-prone PCR and a combination of three beneficial mutations, dominated by an intriguing serine-to-glycine substitution (M. Anbar, R. Lamed, E. A. Bayer, ChemCatChem 2:997–1003, 2010). In the present study, we used a bioinformatics-based approach involving sequence alignment of homologous family 8 glycoside hydrolases to create a library of consensus mutations in which residues of the catalytic module are replaced at specific positions with the most prevalent amino acids in the family. One of the mutants (G283P) displayed a higher thermal stability than the wild-type enzyme. Introducing this mutation into the previously engineered Cel8A triple mutant resulted in an optimized enzyme, increasing the half-life of activity by 14-fold at 85°C. Remarkably, no loss of catalytic activity was observed compared to that of the wild-type endoglucanase. The structural changes were simulated by molecular dynamics analysis, and specific regions were identified that contributed to the observed thermostability. Intriguingly, most of the proteins used for sequence alignment in determining the consensus residues were derived from mesophilic bacteria, with optimal temperatures well below that of C. thermocellum Cel8A.

Triticum durum Metallothionein ISOLATION OF THE GENE AND STRUCTURAL CHARACTERIZATION OF THE PROTEIN USING SOLUTION SCATTERING AND MOLECULAR MODELING

A novel gene sequence, with two exons and one intron, encoding a metallothionein (MT) has been identified in durum wheat Triticum durum cv. Balcali85 genomic DNA. Multiple alignment analyses on the cDNA and the translated protein sequences showed that T. durum MT (dMT) can be classified as a type 1 MT. dMT has three Cys-X-Cys motifs in each of the N- and C-terminal domains and a 42-residue-long hinge region devoid of cysteines. dMT was overexpressed in Escherichia coli as a fusion protein (GSTdMT), and bacteria expressing the fusion protein showed increased tolerance to cadmium in the growth medium compared with controls. Purified GSTdMT was characterized by SDS- and native-PAGE, size exclusion chromatography, and matrix-assisted laser desorption ionization time-of-flight mass spectrometry. It was shown that the recombinant protein binds 4 ± 1 mol of cadmium/mol of protein and has a high tendency to form stable oligomeric structures. The structure of GSTdMT and dMT was investigated by synchrotron x-ray solution scattering and computational methods. X-ray scattering measurements indicated a strong tendency for GSTdMT to form dimers and trimers in solution and yielded structural models that were compatible with a stable dimeric form in which dMT had an extended conformation. Results of homology modeling and ab initio solution scattering approaches produced an elongated dMT structure with a long central hinge region. The predicted model and those obtained from x-ray scattering are in agreement and suggest that dMT may be involved in functions other than metal detoxification.

Novel Hydrophobins from Trichoderma Define a New Hydrophobin Subclass: Protein Properties, Evolution, Regulation and Processing

Hydrophobins are small proteins, characterised by the presence of eight positionally conserved cysteine residues, and are present in all filamentous asco- and basidiomycetes. They are found on the outer surfaces of cell walls of hyphae and conidia, where they mediate interactions between the fungus and the environment. Hydrophobins are conventionally grouped into two classes (class I and II) according to their solubility in solvents, hydropathy profiles and spacing between the conserved cysteines. Here we describe a novel set of hydrophobins from Trichoderma spp. that deviate from this classification in their hydropathy, cysteine spacing and protein surface pattern. Phylogenetic analysis shows that they form separate clades within ascomycete class I hydrophobins. Using T. atroviride as a model, the novel hydrophobins were found to be expressed under conditions of glucose limitation and to be regulated by differential splicing.

A homologous production system for Trichoderma reesei secreted proteins in a cellulase-free background

Recent demands for the production of biofuels from lignocellulose led to an increased interest in engineered cellulases from Trichoderma reesei or other fungal sources. While the methods to generate such mutant cellulases on DNA level are straightforward, there is often a bottleneck in their production since a correct posttranslational processing of these enzymes is needed to obtain highly active enzymes. Their production and subsequent enzymatic analysis in the homologous host T. reesei is, however, often disturbed by the concomitant production of other endogenous cellulases. As a useful alternative, we tested the production of cellulases in T. reesei in a genetic background where cellulase formation has been impaired by deletion of the major cellulase transcriptional activator gene xyr1. Three cellulase genes (cel7a, cel7b, and cel12a) were expressed under the promoter regions of the two highly expressed genes tef1 (encoding translation elongation factor 1-alpha) or cdna1 (encoding the hypothetical protein Trire2:110879). When cultivated on d-glucose as carbon source, the Δxyr1 strain secreted all three cellulases into the medium. Related to the introduced gene copy number, the cdna1 promoter appeared to be superior to the tef1 promoter. No signs of proteolysis were detected, and the individual cellulases could be assayed over a background essentially free of other cellulases. Hence this system can be used as a vehicle for rapid and high-throughput testing of cellulase muteins in a homologous background.

Effect of codon optimization on the expression of Trichoderma reesei endoglucanase 1 in Pichia pastoris

Trichoderma reesei cellulases are important biocatalysts for a wide range of industrial applications that include the paper, feed, and textile industries. T. reesei endoglucanase 1 (egl1) was successfully expressed as an active and stable catalyst in Pichia pastoris for the first time. Codon optimization was applied to egl1 of T. reesei to enhance its expression levels in P. pastoris. When compared with the originally cloned egl1 gene of T. reesei, the synthetic codon optimized egl1 gene (egl1s) was expressed at a higher level in P. pastoris. Batch fermentations of both clones with the same copy number under controlled conditions indicated that codon optimized EGI enzyme activity increased to 1.24 fold after 72 h of methanol induction. Our research indicated that P. pastoris is a suitable host for cellulase production.

Classification of GPCRs Using Family Specific Motifs

The classification of G-Protein Coupled Receptor (GPCR) sequences is an important problem that arises from the need to close the gap between the large number of orphan receptors and the relatively small number of annotated receptors. Equally important is the characterization of GPCR Class A subfamilies and gaining insight into the ligand interaction since GPCR Class A encompasses a very large number of drug-targeted receptors. In this work, we propose a method for Class A subfamily classification using sequence-derived motifs which characterizes the subfamilies by discovering receptor-ligand interaction sites. The motifs that best characterize a subfamily are selected by the Distinguishing Power Evaluation (DPE) technique we propose. The experiments performed on GPCR sequence databases show that our method outperforms state-of-the-art classification techniques for GPCR Class A subfamily prediction. An important contribution of our work is to discover key receptor-ligand interaction sites which is very important for drug design.

Molecular characterization of cDNA encoding resistance gene-like sequences in Buchloe dactyloides

Current knowledge of resistance (R) genes and their use for genetic improvement in buffalograss (Buchloe dactyloides [Nutt.] Engelm.) lag behind most crop plants. This study was conducted to clone and characterize cDNA encoding R gene-like (RGL) sequences in buffalograss. This report is the first to clone and characterize of buffalograss RGLs. Degenerate primers designed from the conserved motifs of known R genes were used to amplify RGLs and fragments of expected size were isolated and cloned. Sequence analysis of cDNA clones and analysis of putative translation products revealed that most encoded amino acid sequences shared the similar conserved motifs found in the cloned plant disease resistance genes PRS2, MLA6, L6, RPMI, and Xa1. These results indicated diversity of the R gene candidate sequences in buffalograss. Analysis of 5′ rapid amplification of cDNA ends (RACE), applied to investigate upstream of RGLs, indicated that regulatory sequences such as TATA box were conserved among the RGLs identified. The cloned RGL in this study will further enhance our knowledge on organization, function, and evolution of R gene family in buffalo grass. With the sequences of the primers and sizes of the markers provided, these RGL markers are readily available for use in a genomics-assisted selection in buffalograss.

Modulating the thermostability of Endoglucanase I from Trichoderma reesei using computational approaches

In the last decades, effective cellulose degradation became a major point of interest due to the properties of cellulose as a renewable energy source and the widespread application of cellulases (the cellulose degrading enzymes) in many industrial processes. Effective bioconversion of lignocellulosic biomass into soluble sugars for ethanol production requires use of thermostable and highly active cellulases. The library of current cellulases includes enzymes that can work at acidic and neutral pH in a wide temperature range. However, only few cellulases are reported to be thermostable. In order to alleviate this, we have performed a hybrid approach for the thermostabilization of a key cellulase, Endoglucanase I (EGI) from Trichoderma reesei. We combined in silico and in vitro experiments to modulate the thermostability of EGI. Four different predictive algorithms were used to set up a library of mutations. Three thermostabilizer mutations (Q126F, K272F, Q274V) were selected and molecular dynamics simulations at room temperature and high temperatures were performed to analyze the effect of the mutations on enzyme structure and stability. The mutations were then introduced into the endoglucanase 1 gene, using site-directed mutagenesis, and the effect of the mutations on enzyme structure and stability were determined. MD simulations supported the fact that Q126F, K272F and Q274V mutations have a thermostabilizing effect on the protein structure. Experimental studies validated that all of the mutants exhibited higher thermostability compared with native EGI albeit with a decrease in specific activity.

Modifying the catalytic preference of tributyrin in Bacillus thermocatenulatus lipase through in-silico modeling of enzyme-substrate complex

In this study, rational design for Bacillus thermocatenulatus lipase (BTL2) was carried out to lower the activation barrier for hydrolysis of short-chain substrates. In this design, we used computational models for the enzyme-substrate (ES) complexes of tributyrin (C4) and tricaprylin (C8), which were generated through docking and molecular dynamics (MD) simulations. These ES complexes were employed in steered MD (SMD) simulations with Jarzynskis equality to estimate their relative binding free energies. Potential mutation sites for modifying the chain-length selectivity of BTL2 were found by inspecting the SMD trajectories and fine-tuning the volume and hydrophobicity of the cleft. Seven mutations (F17A, L57F, V175A, V175F, I320A, I320F and L360F) were performed to cover three binding pockets for sn-1, sn-2 and sn-3 acyl chains. The relative binding free energies of the mutant ES complexes formed by C4 and C8 ligands were calculated similarly. The experimental routines of protein engineering including site-directed mutagenesis, heterologous protein expression and purification were performed for all lipases. Steady-state specific activities towards C4 and C8 were determined for wild-type and mutant lipases, which gave an estimate of the relative change in the binding free energy of transition state complex (ES(‡)). The chain-length selectivity of mutants was determined from the relative changes in the activation barrier of hydrolysis of C4 and C8 triacylglycerol with respect to wild-type using computational and experimental findings. The most promising mutant for C4 over C8 preference was found to be L360F. We suggest that L360F may be at a critical position to lower the activation barrier for C4 and elevate it for C8 hydrolysis.

Clustering of protein families into functional subtypes using Relative Complexity Measure with reduced amino acid alphabets

BackgroundPhylogenetic analysis can be used to divide a protein family into subfamilies in the absence of experimental information. Most phylogenetic analysis methods utilize multiple alignment of sequences and are based on an evolutionary model. However, multiple alignment is not an automated procedure and requires human intervention to maintain alignment integrity and to produce phylogenies consistent with the functional splits in underlying sequences. To address this problem, we propose to use the alignment-free Relative Complexity Measure (RCM) combined with reduced amino acid alphabets to cluster protein families into functional subtypes purely on sequence criteria. Comparison with an alignment-based approach was also carried out to test the quality of the clustering.ResultsWe demonstrate the robustness of RCM with reduced alphabets in clustering of protein sequences into families in a simulated dataset and seven well-characterized protein datasets. On protein datasets, crotonases, mandelate racemases, nucleotidyl cyclases and glycoside hydrolase family 2 were clustered into subfamilies with 100% accuracy whereas acyl transferase domains, haloacid dehalogenases, and vicinal oxygen chelates could be assigned to subfamilies with 97.2%, 96.9% and 92.2% accuracies, respectively.ConclusionsThe overall combination of methods in this paper is useful for clustering protein families into subtypes based on solely protein sequence information. The method is also flexible and computationally fast because it does not require multiple alignment of sequences.