Flexibility of cold- and heat-adapted subtilisin-like serine proteinases evaluated with fluorescence quenching and molecular dynamics

The subtilisin-like serine proteinases, VPR, from a psychrotrophic Vibrio species and aqualysin I (AQUI) from the thermophile Thermus aquaticus, are structural homologues, but differ significantly with respect to stability and catalytic properties. It has been postulated that the higher catalytic activity of cold adapted enzymes when compared to homologues from thermophiles, reflects their higher molecular flexibility. To assess a potential difference in molecular flexibility between the two homologous proteinases, we have measured their Trp fluorescence quenching by acrylamide at different temperatures. We also investigated protein dynamics of VPR and AQUI at an atomic level by molecular dynamics simulations. VPR contains four Trp residues, three of which are at corresponding sites in the structure of AQUI. To aid in the comparison, a Tyr at the fourth corresponding site in AQUI was mutated to Trp (Y191W). A lower quenching effect of acrylamide on the intrinsic fluorescence of the thermophilic AQUI_Y191W was observed at all temperatures measured (10–55 °C), suggesting that it possesses a more rigid structure than VPR. The MD analysis (Cα rmsf profiles) showed that even though VPR and AQUI have similar flexibility profiles, the cold adapted VPR displays higher flexibility in most regions of the protein structure. Some of these regions contain or are in proximity to some of the Trp residues (Trp6, Trp114 and Trp208) in the proteins. Thus, we observe an overall agreement between the fluorescence quenching data and the flexibility profiles obtained from the MD simulations to different flexibilities of specific regions in the proteins.     

Sequence,Structure, and Evolution of Cellulases in Glycoside Hydrolase Family 48

Currently, the cost of cellulase enzymes remains a key economic impediment to commercialization of biofuels (1). Enzymes from glycoside hydrolase family 48 (GH48) are a critical component of numerous natural lignocellulose-degrading systems. Although computational mining of large genomic data sets is a promising new approach for identifying novel cellulolytic activities, current computational methods are unable to distinguish between cellulases and enzymes with different substrate specificities that belong to the same protein family. We show that by using a robust computational approach supported by experimental studies, cellulases and non-cellulases can be effectively identified within a given protein family. Phylogenetic analysis of GH48 showed non-monophyletic distribution, an indication of horizontal gene transfer. Enzymatic function of GH48 proteins coded by horizontally transferred genes was verified experimentally, which confirmed that these proteins are cellulases. Computational and structural studies of GH48 enzymes identified structural elements that define cellulases and can be used to computationally distinguish them from non-cellulases. We propose that the structural element that can be used for in silico discrimination between cellulases and non-cellulases belonging to GH48 is an ω-loop located on the surface of the molecule and characterized by highly conserved rare amino acids. These markers were used to screen metagenomics data for “true” cellulases.     

Comparative thermal unfolding study of psychrophilic and mesophilic subtilisin-like serine proteases by molecular dynamics simulations

Molecular dynamics (MD) simulations of a subtilisin-like serine protease VPR from the psychrophilic marine bacterium Vibrio sp. PA-44 and its mesophilic homologue, proteinase K (PRK), have been performed for 20 ns at four different temperatures (300, 373, 473, and 573 K). The comparative analyses of MD trajectories reveal that at almost all temperatures, VPR exhibits greater structural fluctuations/deviations, more unstable regular secondary structural elements, and higher global flexibility than PRK. Although these two proteases follow similar unfolding pathways at high temperatures, VPR initiates unfolding at a lower temperature and unfolds faster at the same high temperatures than PRK. These observations collectively indicate that VPR is less stable and more heat-labile than PRK. Analyses of the structural/geometrical properties reveal that, when compared to PRK, VPR has larger radius of gyration (Rg), less intramolecular contacts and hydrogen bonds (HBs), more protein-solvent HBs, and smaller burial of nonpolar area and larger exposure of polar area. These suggest that the increased flexibility of VPR would be most likely caused by its reduced intramolecular interactions and more favourable protein-solvent interactions arising from the larger exposure of the polar area, whereas the enhanced stability of PRK could be ascribed to its increased intramolecular interactions arising from the better optimized hydrophobicity. The factors responsible for the significant differences in local flexibility between these two proteases were also analyzed and ascertained. This study provides insights into molecular basis of thermostability of homologous serine proteases adapted to different temperatures.     

Comparison of the synergistic action of two thermostable xylanases from GH families 10 and 11 with thermostable cellulases in lignocellulose hydrolysis

Recombinant xylanase preparations from Nonomuraea flexuosa (Nf Xyn, GH11) and Thermoascus aurantiacus (Ta Xyn, GH10) were evaluated for their abilities to hydrolyze hydrothermally pretreated wheat straw. The GH family 10 enzyme Ta Xyn was clearly more efficient in solubilizing xylan from pretreated wheat straw. Improvement of the hydrolysis of hydrothermally pretreated wheat straw by addition of the thermostable xylanase preparations to thermostable cellulases was evaluated. Clear synergistic enhancement of hydrolysis of cellulose was observed when cellulases were supplemented even with a low amount of pure xylanases. Xylobiose was the main hydrolysis product from xylan. It was found that the hydrolysis of cellulose increased nearly linearly with xylan removal during the enzymatic hydrolysis. The results also showed that the xylanase preparation from T. aurantiacus, belonging to GH family 10 always showed better hydrolytic capacity of solubilizing xylan and acting synergistically with thermostable cellulases in the hydrolysis of hydrothermally pretreated wheat straw.     

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.     

Large-Scale Analyses of Glycosylation in Cellulases

Cellulases are important glycosyl hydrolases (GHs) that hydrolyze cellulose polymers into smaller oligosaccharides by breaking the cellulose β (1→4) bonds,and they are widely used to produce cellulosic ethanol from the plant biomass.N-linked and O-linked glycosylations were proposed to impact the catalytic efficiency,cellulose binding affinity and the stability of cellulases based on observations of individual cellulases.As far as we know,there has not been any systematic analysis of the distributions of N-...     

Computer simulations of the flexibility of a series of synthetic cyclic peptide analogues

The flexibility of a series of cyclic peptides derived from the epitope of a snake toxin is investigated using computer simulation techniques. Molecular dynamics (MD) simulations and vibrational analyses are performed on chemically constrained peptides modeled on the parent protein loop. In the 50 ps MD simulations, large variations in the atomic fluctuations are seen between the peptides, and can be related to the nature of the chemical constraints present in the molecules. Normal mode analyses are performed on energy-minimized configurations derived from the dynamics trajectories. The atomic fluctuations calculated from the normal modes are about 30% of those of the molecular dynamics for the more flexible peptides and 70% for the more constrained molecules. The calculated differences in flexibility between the molecules are much less significant in the harmonic approximation. © 1993 John Wiley & Sons, Inc.     

Conformational flexibility of two RNA trimers explored by computational tools and database search

Two RNA sequences, AAA and AUG, were studied by the conformational search program CICADA and by molecular dynamics (MD) in the framework of the AMBER force field, and also via thorough PDB database search. CICADA was used to provide detailed information about conformers and conformational interconversions on the energy surfaces of the above molecules. Several conformational families were found for both sequences. Analysis of the results shows differences, especially between the energy of the single families, and also in flexibility and concerted conformational movement. Therefore, several MD trajectories (altogether 16 ns) were run to obtain more details about both the stability of conformers belonging to different conformational families and about the dynamics of the two systems. Results show that the trajectories strongly depend on the starting structure. When the MD start from the global minimum found by CICADA, they provide a stable run, while MD starting from another conformational family generates a trajectory where several different conformational families are visited. The results obtained by theoretical methods are compared with the thorough database search data. It is concluded that all except for the highest energy conformational families found in theoretical result also appear in experimental data. Registry numbers: adenylyl-(3' --> 5')-adenylyl-(3' --> 5')-adenosine [917-44-2] adenylyl-(3' --> 5')-uridylyl-(3' --> 5')-guanosine [3494-35-7].     

Altering the linker in processive GH5 endoglucanase 1 modulates lignin binding and catalytic properties

Background

The non-productive adsorption of cellulases onto lignin in biomass is a key issue for the biofuel process economy. It would be helpful to reduce the inhibitory effect of lignin on enzymatic hydrolysis by engineering weak lignin-binding cellulases. Cellulase linkers are highly divergent in their lengths, compositions, and glycosylations. Numerous studies have revealed that linkers can facilitate optimal interactions between structured domains. Recently, efforts have focused on the contributions and mechanisms of carbohydrate-binding modules and catalytic domains that affect lignin affinity and processivity of cellulases, but our understanding of the effects of the linker regions on lignin adsorption and processivity of GH5 processive endoglucanases is still limited.

Results

Eight GH5 endoglucanase 1 variants of varying length, flexibility, and sequence in the linker region were constructed. Their characteristics were then compared to the wild-type enzyme (EG1). Remarkably, significant differences in the lignin adsorption profiles and processivities were observed for EG1 and other variants. Our studies suggest that either the length or the specific amino acid composition of the linker has a prominent influence on the lignin-binding affinity of the enzymes. Comparatively, the processivity may depend primarily on the length of the linker and less so on the specific amino acid composition. EG1-ApCel5A, a variant with better performance in enzymatic hydrolysis in the presence of lignin, was obtained by replacing a longer, flexible linker. In total, up to between 28.2 and 30.1% more reducing sugars were generated from filter paper by EG1-ApCel5A in the presence of lignin compared to EG1.

Conclusions

Our results highlight the relevance of the linker region in the lignin adsorption and processivity of a processive endoglucanase. Our findings suggest that the linker region may be used as a target for the design of more active and weaker lignin-binding cellulases.

     

Cel9D, an atypical 1,4-beta-D-glucan glucohydrolase from Fibrobacter succinogenes: characteristics, catalytic residues, and synergistic interactions with other cellulases

The increasing demands of renewable energy have led to the critical emphasis on novel enzymes to enhance cellulose biodegradation for biomass conversion. To identify new cellulases in the ruminal bacterium Fibrobacter succinogenes, a cell extract of cellulose-grown cells was separated by ion-exchange chromatography and cellulases were located by zymogram analysis and identified by peptide mass fingerprinting. An atypical family 9 glycoside hydrolase (GH9), Cel9D, with less than 20% identity to typical GH9 cellulases, was identified. Purified recombinant Cel9D enhanced the production of reducing sugar from acid swollen cellulose (ASC) and Avicel by 1.5- to 4-fold when mixed separately with each of four other glucanases, although it had low activity on these substrates. Cel9D degraded ASC and cellodextrins with a degree of polymerization higher than 2 to glucose with no apparent endoglucanase activity, and its activity was restricted to beta-1-->4-linked glucose residues. It catalyzed the hydrolysis of cellulose by an inverting mode of reaction, releasing glucose from the nonreducing end. Unlike many GH9 cellulases, calcium ions were not required for its function. Cel9D had increased kcat/Km values for cello-oligosaccharides with higher degrees of polymerization. The kcat/Km value for cellohexaose was 2,300 times higher than that on cellobiose. This result indicates that Cel9D is a 1,4-beta-D-glucan glucohydrolase (EC 3.2.1.74) in the GH9 family. Site-directed mutagenesis of Cel9D identified Asp166 and Glu612 as the candidate catalytic residues, while Ser168, which is not present in typical GH9 cellulases, has a crucial structural role. This enzyme has an important role in crystalline cellulose digestion by releasing glucose from accessible cello-oligosaccharides.     

Characterization of All Family-9 Glycoside Hydrolases Synthesized by the Cellulosome-producing Bacterium Clostridium cellulolyticum

The genome of Clostridium cellulolyticum encodes 13 GH9 enzymes that display seven distinct domain organizations. All but one contain a dockerin module and were formerly detected in the cellulosomes, but only three of them were previously studied (Cel9E, Cel9G, and Cel9M). In this study, the 10 uncharacterized GH9 enzymes were overproduced in Escherichia coli and purified, and their activity pattern was investigated in the free state or in cellulosome chimeras with key cellulosomal cellulases. The newly purified GH9 enzymes, including those that share similar organization, all exhibited distinct activity patterns, various binding capacities on cellulosic substrates, and different synergies with pivotal cellulases in mini-cellulosomes. Furthermore, one enzyme (Cel9X) was characterized as the first genuine endoxyloglucanase belonging to this family, with no activity on soluble and insoluble celluloses. Another GH9 enzyme (Cel9V), whose sequence is 78% identical to the cellulosomal cellulase Cel9E, was found inactive in the free and complexed states on all tested substrates. The sole noncellulosomal GH9 (Cel9W) is a cellulase displaying a broad substrate specificity, whose engineered form bearing a dockerin can act synergistically in minicomplexes. Finally, incorporation of all GH9 cellulases in trivalent cellulosome chimera containing Cel48F and Cel9G generated a mixture of heterogeneous mini-cellulosomes that exhibit more activity on crystalline cellulose than the best homogeneous tri-functional complex. Altogether, our data emphasize the importance of GH9 diversity in bacterial cellulosomes, confirm that Cel9G is the most synergistic GH9 with the major endoprocessive cellulase Cel48F, but also identify Cel9U as an important cellulosomal component during cellulose depolymerization.     

Molecular Dynamics Study of the Structure,Flexibility and Dynamics of Thermostable L1 Lipase at High Temperatures

Molecular Dynamics (MD) simulations have been used to understand how protein structure, dynamics, and flexibility are affected by adaptation to high temperature for several years. We report here the results of the high temperature MD simulations of Bacillus stearothermophilus L1 (L1 lipase). We found that the N-terminal moiety of the enzyme showed a high flexibility and dynamics during high temperature simulations which preceded and followed by clear structural changes in two specific regions; the small domain and the main catalytic domain or core domain of the enzyme. These two domains interact with each other through a Zn²⁺-binding coordination with Asp-61 and Asp-238 from the core domain and His-81 and His-87 from the small domain. Interestingly, the His-81 and His-87 were among the highly fluctuated and mobile residues at high temperatures. The results appear to suggest that tight interactions of Zn²⁺-binding coordination with specified residues became weak at high temperature which suggests the contribution of this region to the thermostability of the enzyme. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Conformational Flexibility of Two RNA Trimers Explored by Computational Tools and Database Search

Abstract

Two RNA sequences, AAA and AUG, were studied by the conformational search program CICADA and by molecular dynamics (MD) in the framework of the AMBER force field, and also via thorough PDB database search. CICADA was used to provide detailed information about conformers and conformational interconversions on the energy surfaces of the above molecules. Several conformational families were found for both sequences. Analysis of the results shows differences, especially between the energy of the single families, and also in flexibility and concerted conformational movement. Therefore, several MD trajectories (altogether 16 ns) were run to obtain more details about both the stability of conformers belonging to different conformational families and about the dynamics of the two systems. Results show that the trajectories strongly depend on the starting structure. When the MD start from the global minimum found by CICADA, they provide a stable run, while MD starting from another conformational family generates a trajectory where several different conformational families are visited. The results obtained by theoretical methods are compared with the thorough database search data. It is concluded that all except for the highest energy conformational families found in theoretical result also appear in experimental data.

Registry numbers:

adenylyl-(3′ →5′)-adenylyl-(3′ →5′)-adenosine [917-44-2]

adenylyl-(3′ →5′)-uridylyl-(3′ →5′)-guanosine [3494-35-7]     

Glycosylation effects on the structure and dynamics of a full-length Cel7A cellulase

Fungi cellulases are used to degrade cellulose-containing biomass for bioethanol production. Industrial cellulases such as Cel7A from Trichoderma reesei (TrCel7A) are critical in this process. Thus, the understanding of structure and dynamics is crucial for engineering variants with improved cellulolytic activity. This cellulase consists of two domains connected by a flexible and highly glycosylated linker. However, the linker flexibility has hindered the determination of Cel7A complete structure. Herein, based on atomic and sparse data, we applied integrative modelling to build a model of the complete enzyme structure. Next, through simulations, we studied the glycosylation effects on the structure and dynamics of a solubilized TrCel7A. Essential dynamics analysis showed that O-glycosylation in the linker led to the stabilization of protein overall dynamics. O-linked glycans seem to restrict protein dihedral angles distribution in this region, selecting more elongated conformations. Besides the reduced flexibility, functional interdomain motions occurred in a more concerted way in the glycosylated system. In contrast, in the absence of glycosylation, we observed vast conformational plasticity with the functional domains frequently collapsing. We report here evidence that targeting Cel7A linker flexibility by point mutations including modification of glycosylation sites could be a promising design strategy to improve cellulase activity.     

Atomistic and coarse-grained analysis of double spectrin repeat units: the molecular origins of flexibility

Spectrin is an ubiquitous protein in metazoan cells, and its flexibility is one of the keys to maintaining cellular structure and organization. Both alpha-spectrin and beta-spectrin polypeptides consist primarily of triple coiled-coil modular repeat units, and two important factors that determine spectrin flexibility are the bending flexibility between two consecutive repeat units and the conformational flexibility of individual repeat units. Atomistic molecular dynamics (MD) simulations are used here to study double spectrin repeat units (DSRUs) from the human erythrocyte beta-spectrin (HEbeta89) and the chicken brain alpha-spectrin (CBalpha1617). From the results of MD simulations, a highly conserved Trp residue in the A-helix of most repeat units that has been suggested to be important in conferring stability to the coiled-coil structures is found not to have a significant effect on the conformational flexibility of individual repeat units. Characterization of the bending flexibility for two consecutive repeats of spectrin via atomistic simulations and coarse-grained (CG) modeling indicate that the bending flexibility is governed by the interactions between the AB-loop of the first repeat unit, the BC-loop of the second repeat unit and the linker region. Specifically, interactions between residues in these regions can lead to a strong directionality in the bending behavior of two repeat units. The biological implications of these finding are discussed.