The esterase from the thermophilic eubacterium Bacillus acidocaldarius: structural-functional relationship and comparison with the esterase from the hyperthermophilic archaeon Archaeoglobus fulgidus

Computer-aided rational design strategy based on protein surface charge to improve the thermal stability of a novel esterase from Geobacillus jurassicus

OBJECTIVES

Although Geobacillus are significant thermophilic bacteria source, there are no reports of thermostable esterase gene in Geobacillus jurassicus or rational design strategies to increase the thermal stability of esterases.

RESULTS

Gene gju768 showed a highest similarity of 15.20% to esterases from Geobacillus sp. with detail enzymatic properties. Using a combination of Gibbs Unfolding Free Energy (∆∆G) calculator and the distance from the mutation site to the catalytic site (DsCα-Cα) to screen suitable mutation sites with elimination of negative surface charge, the mutants (D24N, E221Q, and E253Q) displayed stable mutants with higher thermal stability than the wild-type (WT). Mutant E253Q exhibited the best thermal stability, with a half-life (T1/2) at 65 °C of 32.4 min, which was 1.8-fold of the WT (17.9 min).

CONCLUSION

Cloning of gene gju768 and rational design based on surface charge engineering contributed to the identification of thermostable esterase from Geobacillus sp. and the exploration of evolutionary strategies for thermal stability.

Identification, Characterization, and Computer-Aided Rational Design of a Novel Thermophilic Esterase from Geobacillus subterraneus, and Application in the Synthesis of Cinnamyl Acetate

Investigation of a novel thermophilic esterase gene from Geobacillus subterraneus DSMZ 13552 indicated a high amino acid sequence similarity of 25.9% to a reported esterase from Geobacillus sp. A strategy that integrated computer-aided rational design tools was developed to select mutation sites. Six mutants were selected from four criteria based on the simulated saturation mutation (including 19 amino acid residues) results. Of these, the mutants Q78Y and G119A were found to retain 87% and 27% activity after incubation at 70 °C for 20 min, compared with the 19% activity for the wild type. Subsequently, a double-point mutant (Q78Y/G119A) was obtained and identified with optimal temperature increase from 65 to 70 °C and a 41.51% decrease in Km. The obtained T1/2 values of 42.2 min (70 °C) and 16.9 min (75 °C) for Q78Y/G119A showed increases of 340% and 412% compared with that in the wild type. Q78Y/G119A was then employed as a biocatalyst to synthesize cinnamyl acetate, for which the conversion rate reached 99.40% with 0.3 M cinnamyl alcohol at 60 °C. The results validated the enhanced enzymatic properties of the mutant and indicated better prospects for industrial application as compared to that in the wild type. This study reported a method by which an enzyme could evolve to achieve enhanced thermostability, thereby increasing its potential for industrial applications, which could also be expanded to other esterases.

A comparative molecular dynamics study of thermophilic and mesophilic β-fructosidase enzymes

Arabidopsis thaliana cell wall invertase 1 (AtcwINV1) and Thermotoga maritima β-fructosidase (BfrA) are among the best structurally studied members of the glycoside hydrolase family 32. Both enzymes hydrolyze sucrose as the main substrate but differ strongly in their thermal stability. Mesophilic AtcwINV1 and thermophilic BfrA have divergent sequence similarities in the N-terminal five bladed β-propeller catalytic domain (31 %) and the C-terminal β-sandwich domain (15 %) of unknown function. The two enzymes were subjected to 200 ns molecular dynamics simulations at 300 K (27 °C) and 353 K (80 °C). Regular secondary structure regions, but not loops, in AtcwINV1 and BfrA showed no significant fluctuation differences at both temperatures. BfrA was more rigid than AtcwINV1 at 300 K. The simulation at 353 K did not alter the structural stability of BfrA, but did increase the overall flexibility of AtcwINV1 exhibiting the most fluctuating regions in the β-propeller domain. The simulated heat treatment also increased the gyration radius and hydrophobic solvent accessible surface area of the plant enzyme, consistent with the initial steps of an unfolding process. The preservation of the conformational rigidity of BfrA at 353 K is linked to the shorter size of the protein loops. Shortening of BfrA loops appears to be a key mechanism for thermostability.

Thermostabilization of the β1-adrenergic receptor correlates with increased entropy of the inactive state

The dynamic nature of GPCRs is a major hurdle in their purification and crystallization. Thermostabilization can facilitate GPCR structure determination, as has been shown by the structure of the thermostabilized β1-adrenergic receptor (β1AR) mutant, m23-β1AR, which has been thermostabilized in the inactive state. However, it is unclear from the structure how the six thermostabilizing mutations in m23-β1AR affect receptor dynamics. We have used molecular dynamics simulations in explicit solvent to compare the conformational ensembles for both wild type β1AR (wt-β1AR) and m23-β1AR. Thermostabilization results in an increase in the number of accessible microscopic conformational states within the inactive state ensemble, effectively increasing the side chain entropy of the inactive state at room temperature, while suppressing large-scale main chain conformational changes that lead to activation. We identified several diverse mechanisms of thermostabilization upon mutation. These include decrease of long-range correlated movement between residues in the G-protein coupling site to the extracellular region (Y227A(5.58), F338M(7.48)), formation of new hydrogen bonds (R68S), and reduction of local stress (Y227(5.58), F327(7.37), and F338(7.48)). This study provides insights into microscopic mechanisms underlying thermostability that leads to an understanding of the effect of these mutations on the structure of the receptor.

Using empirical phase diagrams to understand the role of intramolecular dynamics in immunoglobulin G stability

Understanding the relationship between protein dynamics and stability is of paramount importance to the fields of biology and pharmaceutics. Clarifying this relationship is complicated by the large amount of experimental data that must be generated and analyzed if motions that exist over the wide range of timescales are to be included. To address this issue, we propose an approach that utilizes a multidimensional vector-based empirical phase diagram (EPD) to analyze a set of dynamic results acquired across a temperature-pH perturbation plane. This approach is applied to a humanized immunoglobulin G1 (IgG1), a protein of major biological and pharmaceutical importance whose dynamic nature is linked to its multiple biological roles. Static and dynamic measurements are used to characterize the IgG and to construct both static and dynamic EPDs. Between pH 5 and 8, a single, pH-dependent transition is observed that corresponds to thermal unfolding of the IgG. Under more acidic conditions, evidence exists for the formation of a more compact, aggregation resistant state of the immunoglobulin, known as A-form. The dynamics-based EPD presents a considerably more detailed pattern of apparent phase transitions over the temperature-pH plane. The utility and potential applications of this approach are discussed.

Rapid, specific and sensitive electrochemical detection of foodborne bacteria

Electrochemical biochips are an emerging tool for point-of-care diagnostic systems in medicine, food and environmental monitoring. In the current study, a thermostable reporter enzyme, esterase 2 (EST2) from Alicyclobacillus acidocaldarius, is used for specific and sensitive detection of bacteria by one-step rRNA/DNA hybridization between a bacterium-specific capture oligodeoxynucleotide (ODN), bacterial 16S rRNA and an uniform EST2-ODN reporter conjugate. The detection limit corresponds to approximately 500 colony forming units (cfu) Escherichia coli. Beside high sensitivity, the application of electrochemical biochips allows discrimination of two gram-negative and two gram-positive bacteria demonstrating the specificity and the potential for parallel detection of microorganisms. The feasibility of identification of foodborne bacteria was studied with meat juice contaminated with E. coli. This detection system has the capability to be applied for monitoring of bacterial food contamination.

A comparative molecular dynamics study of thermophilic and mesophilic ribonuclease HI enzymes

We studied a pair of homologous thermophilic and mesophilic ribonuclease HI enzymes by molecular dynamics simulations. Each protein was subjected to three 5 ns simulations in explicit water at both 310 K and 340 K. The thermophilic enzyme showed larger overall positional fluctuations at both temperatures, while only the mesophilic enzyme at the higher temperature showed significant instability. When the temperature is changed, the relative flexibility of different local segments on the two proteins changed differently. Principal component analysis showed that the simulations of the two proteins explored largely overlapping regions in the conformational space. However, at 340 K, the collective structure variations of the thermophilic protein are different from those of the mesophilic protein. Our results, although not in accordance with the view that hyperthermostability of proteins may originate from their conformational rigidity, are consistent with several recent experimental and simulation studies which showed that thermophilic proteins may be conformationally more flexible than their mesophilic counterparts. The decorrelation between conformational rigidity and hyperthermostability may be attributed to the temperature dependence and long range nature of electrostatic interactions that play more important roles in the structural stability of thermophilic proteins.

D-Trehalose/D-maltose-binding protein from the hyperthermophilic archaeon Thermococcus litoralis: The binding of trehalose and maltose results in different protein conformational states

In this work, we used fluorescence spectroscopy, molecular dynamics simulation, and Fourier transform infrared spectroscopy for investigating the effect of trehalose binding and maltose binding on the structural properties and the physical parameters of the recombinant D-trehalose/D-maltose binding protein (TMBP) from the hyperthermophilic archaeon Thermococcus litoralis. The binding of the two sugars to TMBP was studied in the temperature range 20 degrees-100 degrees C. The results show that TMBP possesses remarkable temperature stability and its secondary structure does not melt up to 90 degrees C. Although both the secondary structure itself and the sequence of melting events were not significantly affected by the sugar binding, the protein assumes different conformations with different physical properties depending whether maltose or trehalose is bound to the protein. At low and moderate temperatures, TMBP possesses a structure that is highly compact both in the absence and in the presence of two sugars. At about 90 degrees C, the structure of the unliganded TMBP partially relaxes whereas both the TMBP/maltose and the TMBP/trehalose complexes remain in the compact state. In addition, Fourier transform infrared results show that the population of alpha-helices exposed to the solvent was smaller in the absence than in the presence of the two sugars. The spectroscopic results are supported by molecular dynamics simulations. Our data on dynamics and stability of TMBP can contribute to a better understanding of transport-related functions of TMBP and constitute ground for targeted modifications of this protein for potential biotechnological applications.

Pressure effect on the stability and the conformational dynamics of the D-Galactose/D-Glucose-binding protein from Escherichia coli

The effect of the pressure on the structure and stability of the D-Galactose/D-Glucose binding protein (GGBP) from Escherichia coli was studied by steady-state and time-resolved fluorescence spectroscopy, and the ability of glucose ligand to stabilize the GGBP structure was also investigated. Steady-state fluorescence experiments showed a marked quenching of fluorescence emission of GGBP in the absence of glucose. Instead, the presence of glucose seems to stabilize the structure of GGBP at low and moderate pressure values. Time-resolved fluorescence measurements showed that the GGBP taumean in the absence of glucose varies significantly up to 600 bar, while in the presence of the ligand it is almost unaffected by pressure increase up to 600 bar. The effect of the pressure on GGBP was also studied by molecular dynamics simulations. The simulation data support the spectroscopic results and confirm that the presence of glucose is able to contrast the negative effects of pressure on the protein structure. Taken together, the spectroscopic and computer simulation studies suggest that at pressure values up to 2000 bar the structure of GGBP in the absence of glucose remains folded, but a significant perturbation of the protein secondary structures can be detected. The binding of glucose reduces the negative effect of pressure on protein structure and confers protection from perturbation especially at moderate pressure values.

The role of calcium in the conformational dynamics and thermal stability of the D-galactose/D-glucose-binding protein from Escherichia coli

We have characterized stability and conformational dynamics of the calcium depleted D-galactose/D-glucose-binding protein (GGBP) from Escherichia coli. The structural stability of the protein was investigated by steady state and time resolved fluorescence, and far-UV circular dichroism in the temperature range from 20 degrees C to 70 degrees C. We have found that the absence of the Ca(2+) ion results in a significant destabilization of the C-terminal domain of the protein. In particular, the melting temperature decreases by about 10 degrees C with the simultaneous loss of the melting cooperativity. Time resolved fluorescence quenching revealed significant loosening of the protein when highly shielded Trp residue(s) became accessible to acrylamide at higher temperatures. We have documented a significant stabilizing effect of glucose that mostly reverts the effect of calcium, that is, the thermal stability of the protein increases by about 10 degrees C and the melting cooperativity is restored. Moreover, the protein structure remains compact with low amplitude of the segmental mobility up to high temperatures. We have used molecular dynamics to identify the structural feature responsible for changes in the temperature stability. Disintegration of the Ca(2+)-binding loop seems to be responsible for the loss of the stability in the absence of calcium. The new insights on the structural properties and temperature stability of the calcium depleted GGBP contribute to better understanding of the protein function and constitute important information for the development of new biotechnological applications of this class of proteins.

Secondary structure, conformational stability and glycosylation of a recombinant Candida rugosa lipase studied by Fourier-transform infrared spectroscopy

The secondary structure of lipase 1 from Candida rugosa, a model system for large monomeric enzymes, has been studied by FTIR (Fourier-transform infrared) spectroscopy in water and 2H2O. The secondary structure content, determined by the analysis of the amide I band absorption through second derivative and curve fitting procedures, is in agreement with that estimated by X-ray data and predicts, in addition, the existence of two classes of alpha-helices. We have also investigated the enzyme stability and aggregation at high temperature by following the protein unfolding. The thermal stability determined by FTIR is in excellent agreement with the temperature dependence of the lipase activity. Furthermore, new insights on the glycosylation of the recombinant protein produced in Pichia pastoris and on its heterogeneity related to different fermentation batches were obtained by the analysis of the IR absorption in the 1200-900 cm(-1) carbohydrate region. A drastic reduction of the intensity of this band was found after enzymic deglycosylation of the protein. To confirm that the FTIR absorption in the 1200-900 cm(-1) region depends on the carbohydrate content and glycoform distribution, we performed an MS analysis of the protein sugar moieties. Glycosidic structures of the high mannose type were found, with mannoses ranging from 8 to 25 residues.

Structural and thermal stability characterization of Escherichia coli D-galactose/D-glucose-binding protein

The effect of temperature and glucose binding on the structure of the galactose/glucose-binding protein from Escherichia coli was investigated by circular dichroism, Fourier transform infrared spectroscopy, and steady-state and time-resolved fluorescence. The data showed that the glucose binding induces a moderate change of the secondary structure content of the protein and increases the protein thermal stability. The infrared spectroscopy data showed that some protein stretches, involved in alpha-helices and beta strand conformations, are particularly sensitive to temperature. The fluorescence studies showed that the intrinsic tryptophanyl fluorescence of the protein is well represented by a three-exponential model and that in the presence of glucose the protein adopts a structure less accessible to the solvent. The new insights on the structural properties of the galactose/glucose-binding protein can contribute to a better understanding of the protein functions and represent fundamental information for the development of biotechnological applications of the protein.

Analysis of Thermal Adaptation in the HSL Enzyme Family

The recently solved three-dimensional (3D) structures of two thermostable members of the carboxylesterase/lipase HSL family, namely the Alicyclobacillus (formerly Bacillus) acidocaldarius and Archaeoglobus fulgidus carboxylesterases (EST2 and AFEST, respectively) were compared with that of the mesophilic homologous counterpart Brefeldine A esterase from Bacillus subtilis. Since the 3D homology models of other members of the HSL family were also available, we performed a structural alignment with all these sequences. The resulting alignment was used to assess the amino acid "traffic rule" in the HSL family. Quite surprisingly, the data were in very good agreement with those recently reported from two independent groups and based on the comparison of a huge number of homologous sequences from the genus Bacillus, Methanococcus and Deinococcus/Thermus. Taken as a whole, the data point to the statistical meaning of defined amino acid conversions going from psychrophilic to hyperthermophilic sequences. We identified and mapped several such changes onto the EST2 structure and observed that such mutations were localized mostly in loops regions or alpha-helices and were mostly excluded from the active site. A site-directed mutagenesis of two of the identified residues confirmed they were involved in thermal stability.

A substrate-induced switch in the reaction mechanism of a thermophilic esterase: kinetic evidences and structural basis

The reaction mechanism of the esterase 2 (EST2) from Alicyclobacillus acidocaldarius was studied at the kinetic and structural level to shed light on the mechanism of activity and substrate specificity increase previously observed in its double mutant M211S/R215L. In particular, the values of kinetic constants (k1, k(-1), k2, and k3) along with activation energies (E1, E(-1), E2, and E3) were measured for wild type and mutant enzyme. The previously suggested substrate-induced switch in the reaction mechanism from kcat=k3 with a short acyl chain substrate (p-nitrophenyl hexanoate) to kcat=k2 with a long acyl chain substrate (p-nitrophenyl dodecanoate) was validated. The inhibition afforded by an irreversible inhibitor (1-hexadecanesulfonyl chloride), structurally related to p-nitrophenyl dodecanoate, was studied by kinetic analysis. Moreover the three-dimensional structure of the double mutant bound to this inhibitor was determined, providing essential information on the enzyme mechanism. In fact, structural analysis explained the observed substrate-induced switch because of an inversion in the binding mode of the long acyl chain derivatives with respect to the acyl- and alcohol-binding sites.

Thermotoga neapolitana adenylate kinase is highly active at 30 degrees C

The adenylate kinase (AK) gene from Thermotoga neapolitana, a hyperthermophilic bacterium, was cloned and overexpressed in Escherichia coli, and the recombinant enzyme was biochemically characterized. The T. neapolitana AK (TNAK) sequence indicates that this enzyme belongs to the long bacterial AKs. TNAK contains the four cysteine residues that bind Zn(2+) in all Gram-positive AKs and in a few other Zn(2+)-containing bacterial AKs. Atomic emission spectroscopy and titration data indicate a content of 1 mol of Zn(2+)/mol of recombinant TNAK. The EDTA-treated enzyme has a melting temperature (T (m)=93.5 degrees C) 6.2 degrees C below that of the holoenzyme (99.7 degrees C), identifying Zn(2+) as a stabilizing feature in TNAK. TNAK is a monomeric enzyme with a molecular mass of approx. 25 kDa. TNAK displays V (max) and K (m) values at 30 degrees C identical with those of the E. coli AK at 30 degrees C, and displays very high activity at 80 degrees C, with a specific activity above 8000 units/mg. The unusually high activity of TNAK at 30 degrees C makes it an interesting model to test the role of enzyme flexibility in activity.

Enantiomeric resolution of 2-aryl propionic esters with hyperthermophilic and mesophilic esterases: contrasting thermodynamic mechanisms

The enantiomeric resolution of 2-aryl propionic esters by hyperthermophilic and mesophilic esterases was found to be governed by contrasting thermodynamic mechanisms. Entropic contributions predominated for mesophilic esterases from Candida rugosa and Rhizomucor miehei, while enthalpic forces controlled this resolution by the esterase from the extremely thermoacidophilic archaeon, Sulfolobus solfataricus P1. This disparity in thermodynamic mechanism can be attributed to the differences in conformational flexibility of mesophilic and thermophilic enzymes as they relate to the temperature range (4-70 degrees C) examined.

Molecular dynamics study of a hyperthermophilic and a mesophilic rubredoxin

In recent years, increased interest in the origin of protein thermal stability has gained attention both for its possible role in understanding the forces governing the folding of a protein and for the design of new highly stable engineered biocatalysts. To study the origin of thermostability, we have performed molecular dynamics simulations of two rubredoxins, from the mesophile Clostridium pasteurianum and from the hyperthermophile Pyrococcus furiosus. The simulations were carried out at two temperatures, 300 and 373 K, for each molecule. The length of the simulations was within the range of 6-7.2 ns. The rubredoxin from the hyperthermophilic organism was more flexible than its mesophilic counterpart at both temperatures; however, the overall flexibility of both molecules at their optimal growth temperature was the same, despite 59% sequence homology. The conformational space sampled by both molecules was larger at 300 K than at 373 K. The essential dynamics analysis showed that the principal overall motions of the two molecules are significantly different. On the contrary, each molecule showed similar directions of motion at both temperatures.

Residues at the active site of the esterase 2 from Alicyclobacillus acidocaldarius involved in substrate specificity and catalytic activity at high temperature

The recently solved three-dimensional structure of the thermophilic esterase 2 from Alicyclobacillus acidocaldarius allowed us to have a snapshot of an enzyme-sulfonate complex, which mimics the second stage of the catalytic reaction, namely the covalent acyl-enzyme intermediate. The aim of this work was to design, by structure-aided analysis and to generate by site-directed and saturation mutagenesis, EST2 variants with changed substrate specificity in the direction of preference for monoacylesters whose acyl-chain length is greater than eight carbon atoms. Positions 211 and 215 of the polypeptide chain were chosen to introduce mutations. Among five variants with single and double amino acid substitutions, three were obtained, M211S, R215L, and M211S/R215L, that changed the catalytic efficiency profile in the desired direction. Kinetic characterization of mutants and wild type showed that this change was achieved by an increase in k(cat) and a decrease in K(m) values with respect to the parental enzyme. The M211S/R215L specificity constant for p-nitrophenyl decanoate substrate was 6-fold higher than the wild type. However, variants M211T, M211S, and M211V showed strikingly increased activity as well as maximal activity with monoacylesters with four carbon atoms in the acyl chain, compared with the wild type. In the case of mutant M211T, the k(cat) for p-nitrophenyl butanoate was 2.4-fold higher. Overall, depending on the variant and on the substrate, we observed improved catalytic activity at 70 degrees C with respect to the wild type, which was a somewhat unexpected result for an enzyme with already high k(cat) values at high temperature. In addition, variants with altered specificity toward the acyl-chain length were obtained. The results were interpreted in the context of the EST2 three-dimensional structure and a proposed catalytic mechanism in which k(cat), e.g. the limiting step of the reaction, was dependent on the acyl chain length of the ester substrate.