Adopting Selected Hydrogen Bonding and Ionic Interactions from Aspergillus fumigatus Phytase Structure Improves the Thermostability of Aspergillus niger PhyA Phytase

Although it has been widely used as a feed supplement to reduce manure phosphorus pollution of swine and poultry, Aspergillus niger PhyA phytase is unable to withstand heat inactivation during feed pelleting. Crystal structure comparisons with its close homolog, the thermostable Aspergillus fumigatus phytase (Afp), suggest associations of thermostability with several key residues (E35, S42, R168, and R248) that form a hydrogen bond network in the E35-to-S42 region and ionic interactions between R168 and D161 and between R248 and D244. In this study, loss-of-function mutations (E35A, R168A, and R248A) were introduced singularly or in combination into seven mutants of Afp. All seven mutants displayed decreases in thermostability, with the highest loss (25% [P < 0.05]) in the triple mutant (E35A R168A R248A). Subsequently, a set of corresponding substitutions were introduced into nine mutants of PhyA to strengthen the hydrogen bonding and ionic interactions. While four mutants showed improved thermostability, the best response came from the quadruple mutant (A58E P65S Q191R T271R), which retained 20% greater (P < 0.05) activity after being heated at 80°C for 10 min and had a 7°C higher melting temperature than that of wild-type PhyA. This study demonstrates the functional importance of the hydrogen bond network and ionic interaction in supporting the high thermostability of Afp and the feasibility of adopting these structural units to improve the thermostability of a homologous PhyA phytase.     

Improvement of the Thermostability and Activity of a Pectate Lyase by Single Amino Acid Substitutions, Using a Strategy Based on Melting-Temperature-Guided Sequence Alignment

In the vast number of random mutagenesis experiments that have targeted protein thermostability, single amino acid substitutions that increase the apparent melting temperature (T_m) of the enzyme more than 1 to 2°C are rare and often require the creation of a large library of mutated genes. Here we present a case where a single beneficial mutation (R236F) of a hemp fiber-processing pectate lyase of Xanthomonas campestris origin (PL_Xc) produced a 6°C increase in T_m and a 23-fold increase in the half-life at 45°C without compromising the enzyme's catalytic efficiency. This success was based on a variation of sequence alignment strategy where a mesophilic amino acid sequence is matched with the sequences of its thermophilic counterparts that have established T_m values. Altogether, two-thirds of the nine targeted single amino acid substitutions were found to have effects either on the thermostability or on the catalytic activity of the enzyme, evidence of a high success rate of mutation without the creation of a large gene library and subsequent screening of clones. Combination of R236F with another beneficial mutation (A31G) resulted in at least a twofold increase in specific activity while preserving the improved T_m value. To understand the structural basis for the increased thermal stability or activity, the variant R236F and A31G R236F proteins and wild-type PL_Xc were purified and crystallized. By structure analysis and computational methods, hydrophobic desolvation was found to be the driving force for the increased stability with R236F.     

Modeling Reduction of Uranium U(VI) under Variable Sulfate Concentrations by Sulfate-Reducing Bacteria

The kinetics for the reduction of sulfate alone and for concurrent uranium [U(VI)] and sulfate reduction, by mixed and pure cultures of sulfate-reducing bacteria (SRB) at 21 ± 3°C were studied. The mixed culture contained the SRB Desulfovibrio vulgaris along with a Clostridium sp. determined via 16S ribosomal DNA analysis. The pure culture was Desulfovibrio desulfuricans (ATCC 7757). A zero-order model best fit the data for the reduction of sulfate from 0.1 to 10 mM. A lag time occurred below cell concentrations of 0.1 mg (dry weight) of cells/ml. For the mixed culture, average values for the maximum specific reaction rate, V_max, ranged from 2.4 ± 0.2 μmol of sulfate/mg (dry weight) of SRB · h⁻¹) at 0.25 mM sulfate to 5.0 ± 1.1 μmol of sulfate/mg (dry weight) of SRB · h⁻¹ at 10 mM sulfate (average cell concentration, 0.52 mg [dry weight]/ml). For the pure culture, V_max was 1.6 ± 0.2 μmol of sulfate/mg (dry weight) of SRB · h⁻¹ at 1 mM sulfate (0.29 mg [dry weight] of cells/ml). When both electron acceptors were present, sulfate reduction remained zero order for both cultures, while uranium reduction was first order, with rate constants of 0.071 ± 0.003 mg (dry weight) of cells/ml · min⁻¹ for the mixed culture and 0.137 ± 0.016 mg (dry weight) of cells/ml · min⁻¹ (U₀ = 1 mM) for the D. desulfuricans culture. Both cultures exhibited a faster rate of uranium reduction in the presence of sulfate and no lag time until the onset of U reduction in contrast to U alone. This kinetics information can be used to design an SRB-dominated biotreatment scheme for the removal of U(VI) from an aqueous source.     

In Silico Rational Design and Systems Engineering of Disulfide Bridges in the Catalytic Domain of an Alkaline α-Amylase from Alkalimonas amylolytica To Improve Thermostability

High thermostability is required for alkaline α-amylases to maintain high catalytic activity under the harsh conditions used in textile production. In this study, we attempted to improve the thermostability of an alkaline α-amylase from Alkalimonas amylolytica through in silico rational design and systems engineering of disulfide bridges in the catalytic domain. Specifically, 7 residue pairs (P35-G426, Q107-G167, G116-Q120, A147-W160, G233-V265, A332-G370, and R436-M480) were chosen as engineering targets for disulfide bridge formation, and the respective residues were replaced with cysteines. Three single disulfide bridge mutants—P35C-G426C, G116C-Q120C, and R436C-M480C—of the 7 showed significantly enhanced thermostability. Combinational mutations were subsequently assessed, and the triple mutant P35C-G426C/G116C-Q120C/R436C-M480C showed a 6-fold increase in half-life at 60°C and a 5.2°C increase in melting temperature compared with the wild-type enzyme. Interestingly, other biochemical properties of this mutant also improved: the optimum temperature increased from 50°C to 55°C, the optimum pH shifted from 9.5 to 10.0, the stable pH range extended from 7.0 to 11.0 to 6.0 to 12.0, and the catalytic efficiency (k_cat/K_m) increased from 1.8 × 10⁴ to 2.4 × 10⁴ liters/g · min. The possible mechanism responsible for these improvements was explored through comparative analysis of the model structures of wild-type and mutant enzymes. The disulfide bridge engineering strategy used in this work may be applied to improve the thermostability of other industrial enzymes.     

Engineering a Disulfide Bond in the Lid Hinge Region of Rhizopus chinensis Lipase: Increased Thermostability and Altered Acyl Chain Length Specificity

The key to enzyme function is the maintenance of an appropriate balance between molecular stability and structural flexibility. The lid domain which is very important for “interfacial activation” is the most flexible part in the lipase structure. In this work, rational design was applied to explore the relationship between lid rigidity and lipase activity by introducing a disulfide bond in the hinge region of the lid, in the hope of improving the thermostability of R. chinensis lipase through stabilization of the lid domain without interfering with its catalytic performance. A disulfide bridge between F95C and F214C was introduced into the lipase from R. chinensis in the hinge region of the lid according to the prediction of the “Disulfide by Design” algorithm. The disulfide variant showed substantially improved thermostability with an eleven-fold increase in the t _1/2 value at 60°C and a 7°C increase of T _m compared with the parent enzyme, probably contributed by the stabilization of the geometric structure of the lid region. The additional disulfide bond did not interfere with the catalytic rate (k _cat) and the catalytic efficiency towards the short-chain fatty acid substrate, however, the catalytic efficiency of the disulfide variant towards pNPP decreased by 1.5-fold probably due to the block of the hydrophobic substrate channel by the disulfide bond. Furthermore, in the synthesis of fatty acid methyl esters, the maximum conversion rate by RCLCYS reached 95% which was 9% higher than that by RCL. This is the first report on improving the thermostability of the lipase from R. chinensis by introduction of a disulfide bond in the lid hinge region without compromising the catalytic rate.     

Affinity for phosphatidylinositol 4,5-bisphosphate determines muscarinic agonist sensitivity of Kv7 K+ channels

Kv7 K⁺-channel subunits differ in their apparent affinity for PIP₂ and are differentially expressed in nerve, muscle, and epithelia in accord with their physiological roles in those tissues. To investigate how PIP₂ affinity affects the response to physiological stimuli such as receptor stimulation, we exposed homomeric and heteromeric Kv7.2, 7.3, and 7.4 channels to a range of concentrations of the muscarinic receptor agonist oxotremorine-M (oxo-M) in a heterologous expression system. Activation of M₁ receptors by oxo-M leads to PIP₂ depletion through G_q and phospholipase C (PLC). Chinese hamster ovary cells were transiently transfected with Kv7 subunits and M₁ receptors and studied under perforated-patch voltage clamp. For Kv7.2/7.3 heteromers, the EC₅₀ for current suppression was 0.44 ± 0.08 µM, and the maximal inhibition (Inhib_max) was 74 ± 3% (n = 5–7). When tonic PIP₂ abundance was increased by overexpression of PIP 5-kinase, the EC₅₀ was shifted threefold to the right (1.2 ± 0.1 µM), but without a significant change in Inhib_max (73 ± 4%, n = 5). To investigate the muscarinic sensitivity of Kv7.3 homomers, we used the A315T pore mutant (Kv7.3^T) that increases whole-cell currents by 30-fold without any change in apparent PIP₂ affinity. Kv7.3^T currents had a slightly right-shifted EC₅₀ as compared with Kv7.2/7.3 heteromers (1.0 ± 0.8 µM) and a strongly reduced Inhib_max (39 ± 3%). In contrast, the dose–response curve of homomeric Kv7.4 channels was shifted considerably to the left (66 ± 8 nM), and Inhib_max was slightly increased (81 ± 6%, n = 3–4). We then studied several Kv7.2 mutants with altered apparent affinities for PIP₂ by coexpressing them with Kv7.3^T subunits to boost current amplitudes. For the lower affinity (Kv7.2 (R463Q)/Kv7.3^T) or higher affinity (Kv7.2 (R463E)/Kv7.3^T) channels, the EC₅₀ and Inhib_max were similar to Kv7.4 or Kv7.3^T homomers (0.12 ± 0.08 µM and 79 ± 6% [n = 3–4] and 0.58 ± 0.07 µM and 27 ± 3% [n = 3–4], respectively). The very low-affinity Kv7.2 (R452E, R459E, and R461E) triple mutant was also coexpressed with Kv7.3^T. The resulting heteromer displayed a very low EC₅₀ for inhibition (32 ± 8 nM) and a slightly increased Inhib_max (83 ± 3%, n = 3–4). We then constructed a cellular model that incorporates PLC activation by oxo-M, PIP₂ hydrolysis, PIP₂ binding to Kv7-channel subunits, and K⁺ current through Kv7 tetramers. We were able to fully reproduce our data and extract a consistent set of PIP₂ affinities.     

Mutagenesis of Trichoderma reesei endoglucanase I: impact of expression host on activity and stability at elevated temperatures

Background

Trichoderma reesei is a key cellulase source for economically saccharifying cellulosic biomass for the production of biofuels. Lignocellulose hydrolysis at temperatures above the optimum temperature of T. reesei cellulases (~50°C) could provide many significant advantages, including reduced viscosity at high-solids loadings, lower risk of microbial contamination during saccharification, greater compatibility with high-temperature biomass pretreatment, and faster rates of hydrolysis. These potential advantages motivate efforts to engineer T. reesei cellulases that can hydrolyze lignocellulose at temperatures ranging from 60–70°C.

Results

A B-factor guided approach for improving thermostability was used to engineer variants of endoglucanase I (Cel7B) from T. reesei (TrEGI) that are able to hydrolyze cellulosic substrates more rapidly than the recombinant wild-type TrEGI at temperatures ranging from 50–70°C. When expressed in T. reesei, TrEGI variant G230A/D113S/D115T (G230A/D113S/D115T Tr_TrEGI) had a higher apparent melting temperature (3°C increase in T_m) and improved half-life at 60°C (t_1/2 = 161 hr) than the recombinant (T. reesei host) wild-type TrEGI (t_1/2 = 74 hr at 60°C, Tr_TrEGI). Furthermore, G230A/D113S/D115T Tr_TrEGI showed 2-fold improved activity compared to Tr_TrEGI at 65°C on solid cellulosic substrates, and was as efficient in hydrolyzing cellulose at 60°C as Tr_TrEGI was at 50°C. The activities and stabilities of the recombinant TrEGI enzymes followed similar trends but differed significantly in magnitude depending on the expression host (Escherichia coli cell-free, Saccharomyces cerevisiae, Neurospora crassa, or T. reesei). Compared to N.crassa-expressed TrEGI, S. cerevisiae-expressed TrEGI showed inferior activity and stability, which was attributed to the lack of cyclization of the N-terminal glutamine in Sc_TrEGI and not to differences in glycosylation. N-terminal pyroglutamate formation in TrEGI expressed in S. cerevisiae was found to be essential in elevating its activity and stability to levels similar to the T. reesei or N. crassa-expressed enzyme, highlighting the importance of this ubiquitous modification in GH7 enzymes.

Conclusion

Structure-guided evolution of T. reesei EGI was used to engineer enzymes with increased thermal stability and activity on solid cellulosic substrates. Production of TrEGI enzymes in four hosts highlighted the impact of the expression host and the role of N-terminal pyroglutamate formation on the activity and stability of TrEGI enzymes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12896-015-0118-z) contains supplementary material, which is available to authorized users.     

Improving the Thermostability and Activity of a Thermophilic Subtilase by Incorporating Structural Elements of Its Psychrophilic Counterpart

The incorporation of the structural elements of thermostable enzymes into their less stable counterparts is generally used to improve enzyme thermostability. However, the process of engineering enzymes with both high thermostability and high activity remains an important challenge. Here, we report that the thermostability and activity of a thermophilic subtilase were simultaneously improved by incorporating structural elements of a psychrophilic subtilase. There were 64 variable regions/residues (VRs) in the alignment of the thermophilic WF146 protease, mesophilic sphericase, and psychrophilic S41. The WF146 protease was subjected to systematic mutagenesis, in which each of its VRs was replaced with those from S41 and sphericase. After successive rounds of combination and screening, we constructed the variant PBL5X with eight amino acid residues from S41. The half-life of PBL5X at 85°C (57.1 min) was approximately 9-fold longer than that of the wild-type (WT) WF146 protease (6.3 min). The substitutions also led to an increase in the apparent thermal denaturation midpoint temperature (T_m) of the enzyme by 5.5°C, as determined by differential scanning calorimetry. Compared to the WT, PBL5X exhibited high caseinolytic activity (25 to 95°C) and high values of K_m and k_cat (25 to 80°C). Our study may provide a rational basis for developing highly stable and active enzymes, which are highly desired in industrial applications.     

A Model for the Action of Vinblastine in Vivo

A model for the action of vinblastine (VLB) on cells multiplying exponentially in vivo with a generation time, T_G, has been derived. It is based on the assumption that cells attempting to pass through mitosis in the presence of VLB lose their proliferative capacity and that this lethal effect occurs only when the cells are exposed to a concentration of VLB which is above a critical value, C_k. The model leads to two predictions. First, that the percentage of cells surviving at any time after exposure to a dose, D, of VLB is 100% if D < D_k and decreases to 0% after a time, T_G, following a dose D ≥ D_k·2^{T G/T1/2}, where D_k represents the dose of VLB required to produce the concentration C_k, and T_1/2 is the half-life of the VLB in vivo. Second, that the time, T_G, at which the percentage of cells surviving an exposure to VLB, at doses greater than D_k·2^{U G/T1/2}, decreases to zero should be equal to the generation time of the cells. Both of these predictions were confirmed experimentally which indicates that the model adequately explains the action of VLB in vivo.     

Enhancing the Thermostability of Serratia plymuthica Sucrose Isomerase Using B-Factor-Directed Mutagenesis

The sucrose isomerase of Serratia plymuthica AS9 (AS9 PalI) was expressed in Escherichia coli BL21(DE3) and characterized. The half-life of AS9 PalI was 20 min at 45°C, indicating that it was unstable. In order to improve its thermostability, six amino acid residues with higher B-factors were selected as targets for site-directed mutagenesis, and six mutants (E175N, K576D, K174D, G176D, S575D and N577K) were designed using the RosettaDesign server. The E175N and K576D mutants exhibited improved thermostability in preliminary experiments, so the double mutant E175N/K576D was constructed. These three mutants (E175N, K576D, E175N/K576D) were characterized in detail. The results indicate that the three mutants exhibit a slightly increased optimal temperature (35°C), compared with that of the wild-type enzyme (30°C). The mutants also share an identical pH optimum of 6.0, which is similar to that of the wild-type enzyme. The half-lives of the E175N, K576D and E175N/K576D mutants were 2.30, 1.78 and 7.65 times greater than that of the wild-type enzyme at 45°C, respectively. Kinetic studies showed that the K_m values for the E175N, K576D and E175N/K576D mutants decreased by 6.6%, 2.0% and 11.0%, respectively, and their k_cat/K_m values increased by 38.2%, 4.2% and 19.4%, respectively, compared with those of the wild-type enzyme. After optimizing the conditions for isomaltulose production at 45°C, we found that the E175N, K576D and E175N/K576D mutants displayed slightly improved isomaltulose yields, compared with the wild-type enzyme. Therefore, the mutants produced in this study would be more suitable for industrial biosynthesis of isomaltulose.     

DCM-Related Tropomyosin Mutants E40K/E54K Over-Inhibit the Actomyosin Interaction and Lead to a Decrease in the Number of Cycling Cross-Bridges

Two DCM mutants (E40K and E54K) of tropomyosin (Tm) were examined using the thin-filament extraction/reconstitu­tion technique. The effects of the Ca²⁺, ATP, phos­phate (Pi), and ADP concentrations on isometric tension and its transients were studied at 25°C, and the results were com­pared to those for the WT protein. Our results indicate that both E40K and E54K have a significantly lower T _HC (high Ca²⁺ ten­sion at pCa 4.66) (E40K: 1.21±0.06 T _a, ±SEM, N = 34; E54K: 1.24±0.07 T _a, N = 28), a significantly lower T _LC (low- Ca²⁺ tension at pCa 7.0) (E40K: 0.07±0.02 T _a, N = 34; E54K: 0.06±0.02 T _a, N = 28), and a significantly lower T _act (Ca²⁺ activatable tension) (T _act = T _HC–T_LC, E40K: 1.15±0.08 T _a, N = 34; E54K: 1.18±0.06 T _a, N = 28) than WT (T _HC = 1.53±0.07 T _a, T _LC = 0.12±0.01 T _a, T _act = 1.40±0.07 T _a, N = 25). All tensions were normalized to T _a ( = 13.9±0.8 kPa, N = 57), the ten­sion of actin-filament reconstituted cardiac fibers (myocardium) under the standard activating conditions. The Ca²⁺ sensitivity (pCa₅₀) of E40K (5.23±0.02, N = 34) and E54K (5.24±0.03, N = 28) was similar to that of the WT protein (5.26±0.03, N = 25). The cooper­a­tivity increased significantly in E54K (3.73±0.25, N = 28) compared to WT (2.80±0.17, N = 25). Seven kinetic constants were deduced using sinusoidal analysis at pCa 4.66. These results enabled us to calculate the cross-bridge distribution in the strongly attached states, and thereby deduce the force/cross-bridge. The results indicate that the force/cross-bridge is ∼15% less in E54K than WT, but remains similar to that of the WT protein in the case of E40K. We conclude that over-inhibition of the actomyosin interaction by E40K and E54K Tm mutants leads to a decreased force-generating ability at systole, which is the main mechanism underlying the early pathogenesis of DCM.     

Measurement of Pulmonary Flow Reserve and Pulmonary Index of Microcirculatory Resistance for Detection of Pulmonary Microvascular Obstruction

Background

The pulmonary microcirculation is the chief regulatory site for resistance in the pulmonary circuit. Despite pulmonary microvascular dysfunction being implicated in the pathogenesis of several pulmonary vascular conditions, there are currently no techniques for the specific assessment of pulmonary microvascular integrity in humans. Peak hyperemic flow assessment using thermodilution-derived mean transit-time (T_mn) facilitate accurate coronary microcirculatory evaluation, but remain unvalidated in the lung circulation. Using a high primate model, we aimed to explore the use of T_mn as a surrogate of pulmonary blood flow for the purpose of measuring the novel indices Pulmonary Flow Reserve [PFR = (maximum hyperemic)/(basal flow)] and Pulmonary Index of Microcirculatory Resistance [PIMR = (maximum hyperemic distal pulmonary artery pressure)×(maximum hyperemic T_mn)]. Ultimately, we aimed to investigate the effect of progressive pulmonary microvascular obstruction on PFR and PIMR.

Methods and Results

Temperature- and pressure-sensor guidewires (TPSG) were placed in segmental pulmonary arteries (SPA) of 13 baboons and intravascular temperature measured. T_mn and hemodynamics were recorded at rest and following intra-SPA administration of the vasodilator agents adenosine (10–400 µg/kg/min) and papaverine (3–24 mg). Temperature did not vary with intra-SPA sensor position (0.010±0.009 v 0.010±0.009°C; distal v proximal; p = 0.1), supporting T_mn use in lung for the purpose of hemodynamic indices derivation. Adenosine (to 200 µg/kg/min) & papaverine (to 24 mg) induced dose-dependent flow augmentations (40±7% & 35±13% T_mn reductions v baseline, respectively; p<0.0001). PFR and PIMR were then calculated before and after progressive administration of ceramic microspheres into the SPA. Cumulative microsphere doses progressively reduced PFR (1.41±0.06, 1.26±0.19, 1.17±0.07 & 1.01±0.03; for 0, 10⁴, 10⁵ & 10⁶ microspheres; p = 0.009) and increased PIMR (5.7±0.6, 6.3±1.0, 6.8±0.6 & 7.6±0.6 mmHg.sec; p = 0.0048).

Conclusions

Thermodilution-derived mean transit time can be accurately and reproducibly measured in the pulmonary circulation using TPSG. Mean transit time-derived PFR and PIMR can be assessed using a TPSG and adenosine or papaverine as hyperemic agents. These novel indices detect progressive pulmonary microvascular obstruction and thus have with a potential role for pulmonary microcirculatory assessment in humans.     

Molecular Structural Basis for the Cold Adaptedness of the Psychrophilic β-Glucosidase BglU in Micrococcus antarcticus

Psychrophilic enzymes play crucial roles in cold adaptation of microbes and provide useful models for studies of protein evolution, folding, and dynamic properties. We examined the crystal structure (2.2-Å resolution) of the psychrophilic β-glucosidase BglU, a member of the glycosyl hydrolase 1 (GH1) enzyme family found in the cold-adapted bacterium Micrococcus antarcticus. Structural comparison and sequence alignment between BglU and its mesophilic and thermophilic counterpart enzymes (BglB and GlyTn, respectively) revealed two notable features distinct to BglU: (i) a unique long-loop L3 (35 versus 7 amino acids in others) involved in substrate binding and (ii) a unique amino acid, His299 (Tyr in others), involved in the stabilization of an ordered water molecule chain. Shortening of loop L3 to 25 amino acids reduced low-temperature catalytic activity, substrate-binding ability, the optimal temperature, and the melting temperature (T_m). Mutation of His299 to Tyr increased the optimal temperature, the T_m, and the catalytic activity. Conversely, mutation of Tyr301 to His in BglB caused a reduction in catalytic activity, thermostability, and the optimal temperature (45 to 35°C). Loop L3 shortening and H299Y substitution jointly restored enzyme activity to the level of BglU, but at moderate temperatures. Our findings indicate that loop L3 controls the level of catalytic activity at low temperatures, residue His299 is responsible for thermolability (particularly heat lability of the active center), and long-loop L3 and His299 are jointly responsible for the psychrophilic properties. The described structural basis for the cold adaptedness of BglU will be helpful for structure-based engineering of new cold-adapted enzymes and for the production of mutants useful in a variety of industrial processes at different temperatures.     

Role of Key Salt Bridges in Thermostability of G. thermodenitrificans EstGtA2: Distinctive Patterns within the New Bacterial Lipolytic Enzyme Family XV

Bacterial lipolytic enzymes were originally classified into eight different families defined by Arpigny and Jaeger (families I-VIII). Recently, the discovery of new lipolytic enzymes allowed for extending the original classification to fourteen families (I-XIV). We previously reported that G. thermodenitrificans EstGtA2 (access no. AEN92268) belonged to a novel group of bacterial lipolytic enzymes. Here we propose a 15^th family (family XV) and suggest criteria for the assignation of protein sequences to the N’ subfamily. Five selected salt bridges, hallmarks of the N’ subfamily (E3/R54, E12/R37, E66/R140, D124/K178 and D205/R220) were disrupted in EstGtA2 using a combinatorial alanine-scanning approach. A set of 14 (R/K→A) mutants was produced, including five single, three double, three triple and three quadruple mutants. Despite a high tolerance to non-conservative mutations for folding, all the alanine substitutions were destabilizing (decreasing T _m by 5 to 14°C). A particular combination of four substitutions exceeded this tolerance and prevents the correct folding of EstGtA2, leading to enzyme inactivation. Although other mutants remain active at low temperatures, the accumulation of more than two mutations had a dramatic impact on EstGtA2 activity at high temperatures suggesting an important role of these conserved salt bridge-forming residues in thermostability of lipolytic enzymes from the N’ subfamily. We also identified a particular interloop salt bridge in EstGtA2 (D194/H222), located at position i -2 and i -4 residues from the catalytic Asp and His respectively which is conserved in other related bacterial lipolytic enzymes (families IV and XIII) with high tolerance to mutations and charge reversal. We investigated the role of residue identity at position 222 in controlling stability-pH dependence in EstGtA2. The introduction of a His to Arg mutation led to increase thermostability under alkaline pH. Our results suggest primary targets for optimization of EstGtA2 for specific biotechnological purposes.     

Subtropical mouse-tailed bats use geothermally heated caves for winter hibernation

We report that two species of mouse-tailed bats (Rhinopoma microphyllum and R. cystops) hibernate for five months during winter in geothermally heated caves with stable high temperature (20°C). While hibernating, these bats do not feed or drink, even on warm nights when other bat species are active. We used thermo-sensitive transmitters to measure the bats’ skin temperature in the natural hibernacula and open flow respirometry to measure torpid metabolic rate at different ambient temperatures (T_a, 16–35°C) and evaporative water loss (EWL) in the laboratory. Bats average skin temperature at the natural hibernacula was 21.7 ± 0.8°C, and no arousals were recorded. Both species reached the lowest metabolic rates around natural hibernacula temperatures (20°C, average of 0.14 ± 0.01 and 0.16 ± 0.04 ml O₂ g⁻¹ h⁻¹ for R. microphyllum and R. cystops, respectively) and aroused from torpor when T_a fell below 16°C. During torpor the bats performed long apnoeas (14 ± 1.6 and 16 ± 1.5 min, respectively) and had a very low EWL. We hypothesize that the particular diet of these bats is an adaptation to hibernation at high temperatures and that caves featuring high temperature and humidity during winter enable these species to survive this season on the northern edge of their world distribution.     

Isomaltulose Synthase from Klebsiella sp. Strain LX3: Gene Cloning and Characterization and Engineering of Thermostability

The gene (palI) encoding isomaltulose synthase (PalI) from a soil bacterial isolate, Klebsiella sp. strain LX3, was cloned and characterized. PalI converts sucrose into isomaltulose, trehalulose, and trace amounts of glucose and fructose. Sequence domain analysis showed that PalI contains an α-amylase domain and (β/α)₈-barrel structures, suggesting that it belongs to the α-amylase family. Sequence alignment indicated that the five amino acid residues of catalytic importance in α-amylases and glucosyltransferases (Asp²⁴¹, Glu²⁹⁵, Asp³⁶⁹, His¹⁴⁵, and His³⁶⁸) are conserved in PalI. Purified recombinant PalI displayed high catalytic efficiency, with a K_m of 54.6 ± 1.7 mM for sucrose, and maximum activity (approximately 328.0 ± 2.5 U/mg) at pH 6.0 and 35°C. PalI activity was strongly inhibited by Fe³⁺ and Hg²⁺ and was enhanced by Mn²⁺ and Mg²⁺. The half-life of PalI was 1.8 min at 50°C. Replacement of selected amino acid residues by proline significantly increased the thermostability of PalI. Simultaneous replacement of Glu⁴⁹⁸ and Arg³¹⁰ with proline resulted in an 11-fold increase in the half-life of PalI at 50°C.     

Structural basis for uracil DNA glycosylase interaction with uracil: NMR study

Two dimensional (2D) NMR and molecular dynamics simulations have been used to determine the three dimensional (3D) structure of a hairpin DNA, d-CTA-GAGGATCC-TUTT-GGATCCT (22mer; abbreviated as U2-hairpin), which has uracil at the second position from the 5′ end of the tetraloop. The ¹H resonances of this hairpin have been assigned almost completely. NMR restrained molecular dynamics and energy minimization procedures have been used to describe the 3D structure of U2-hairpin. This study establishes that the stem of the hairpin adopts a right-handed B-DNA conformation, while the T₁₂ and T₁₅ nucleotides stack upon 3′ and 5′ ends of the stem, respectively. Further, T₁₄ stacks upon both T₁₂ and T₁₅. Though U₁₃ partially stacks upon T₁₄, no stacking interaction is observed between U₁₃ and T₁₂. All the individual nucleotide bases belonging to the stem and T₁₂ and T₁₅ of the loop adopt ‘anti’ conformation with respect to their sugar moiety, while the U₁₃ and T₁₄ of the loop are in ‘syn’ conformation. The turning phosphate in the loop is located between T₁₃ and T₁₄. This study and a concurrent NMR structural study on yet another hairpin DNA d-CTAGAGGAATAA-TTTU-GGATCCT (22mer; abbreviated as U4-hairpin), with uracil at the fourth position from the 5′ end of the tetraloop throw light upon various interactions which have been reported between Escherichia coli uracil DNA glycosylase (UDG) and uracil containing DNA. The of T₁₂ and α, β, γ, and ζ of U₁₃ and γ of T₁₄, which partially influence the local conformation of U₁₃ in U2-hairpin are all locked in ‘trans’ conformation. Such stretched out backbone conformation in the vicinity of U₁₃ could be the reason as to why the U2-hairpin is found to be the poor substrate for its interaction with UDG compared to the other substrates in which the uracil is at first, third and fourth positions of the tetraloop from its 5′ end, as reported earlier by Vinay and Varshney. This study shows that UDG actively promotes the flipping of uracil from a stacked conformation and rules out the possibility of UDG recognizing the flipped out uracil bases.     

Cyclic AMP- and (Rp)-cAMPS-induced Conformational Changes in a Complex of the Catalytic and Regulatory (RIα) Subunits of Cyclic AMP-dependent Protein Kinase

We took a discovery approach to explore the actions of cAMP and two of its analogs, one a cAMP mimic ((S_p)-adenosine cyclic 3′:5′-monophosphorothioate ((S_p)-cAMPS)) and the other a diastereoisomeric antagonist ((R_p)-cAMPS), on a model system of the type Iα cyclic AMP-dependent protein kinase holoenzyme, RIα(91–244)·C-subunit, by using fluorescence spectroscopy and amide H/²H exchange mass spectrometry. Specifically, for the fluorescence experiments, fluorescein maleimide was conjugated to three cysteine single residue substitution mutants, R92C, T104C, and R239C, of RIα(91–244), and the effects of cAMP, (S_p)-cAMPS, and (R_p)-cAMPS on the kinetics of R-C binding and the time-resolved anisotropy of the reporter group at each conjugation site were measured. For the amide exchange experiments, ESI-TOF mass spectrometry with pepsin proteolytic fragmentation was used to assess the effects of (R_p)-cAMPS on amide exchange of the RIα(91–244)·C-subunit complex. We found that cAMP and its mimic perturbed at least parts of the C-subunit interaction Sites 2 and 3 but probably not Site 1 via reduced interactions of the linker region and αC of RIα(91–244). Surprisingly, (R_p)-cAMPS not only increased the affinity of RIα(91–244) toward the C-subunit by 5-fold but also produced long range effects that propagated through both the C- and R-subunits to produce limited unfolding and/or enhanced conformational flexibility. This combination of effects is consistent with (R_p)-cAMPS acting by enhancing the internal entropy of the R·C complex. Finally, the (R_p)-cAMPS-induced increase in affinity of RIα(91–244) toward the C-subunit indicates that (R_p)-cAMPS is better described as an inverse agonist because it decreases the fractional dissociation of the cyclic AMP-dependent protein kinase holoenzyme and in turn its basal activity.Cyclic AMP-dependent protein kinase (PKA)¹ plays a crucial role in a plethora of cellular functions. All isoforms of PKA are composed of two catalytic (C) subunits and homodimeric regulatory (R) subunits (1–3). As the name implies, cAMP is a major PKA regulator (4). Much progress has been made in the last decade in delineating the molecular basis of action of cAMP. An important tactic in this endeavor has been through the comparison of the effects of cAMP with those of two phosphorothioate cAMP analogs: (S_p)-cAMPS (a cAMP mimic) and (R_p)-cAMPS (an antagonist and a diastereoisomer of (S_p)-cAMPS). Although the importance of geometry of the sulfur substitution is critical in determining the pharmacological properties of the two phosphorothioate cAMP analogs, the molecular basis for this behavior is not fully understood. To date, these comparisons have only been made using either wild-type or truncated mutants of the type Iα regulatory subunit (RIα) that are free in solution, not complexed to the C-subunit. X-ray spectroscopic examination of ligand-bound RIα(92–379) complexes reveals few differences between ligand-bound complexes, but the (R_p)-cAMPS complex is structurally “looser” with higher thermal factors than complexes formed with either cAMP or (S_p)-cAMPS (5). This is consistent with the observation that both cAMP and (S_p)-cAMPS, but not (R_p)-cAMPS, raise the urea concentration required for wild-type RIα unfolding (6). Further insight into the structural basis of cAMP action stems from NMR spectroscopic comparison of the effects of (R_p)-cAMPS, cAMP, and (S_p)-cAMPS on chemical shifts and ¹⁵N relaxation of the RIα(119–244) mutant (7). In addition to producing fewer significant chemical shift changes than either cAMP or (S_p)-cAMPS, (R_p)-cAMPS binding is associated with enhanced millisecond to microsecond time scale backbone motions of a β-turn (β2,3 loop) and around the phosphate-binding cassette (PBC) (7).Further insight into the molecular basis of actions of cAMP and its analogs should come from the analysis of ligand-bound R·C complexes. Unfortunately, the large size of even the heterodimeric R·C complex (∼95 kDa) and the difficulty of preparing (R_p)-cAMPS·R·C-subunit crystals currently preclude the use of both NMR spectroscopy and x-ray crystallography. Consequently, we took two alternative lower resolution approaches to this issue. One approach involves the use of site-directed labeling combined with fluorescence spectroscopy to examine both the effects of cAMP and its analogs on R-C subunit binding kinetics and on the conformational dynamics of RIα(91–244). RIα(91–244) includes the “A” cyclic nucleotide binding (CNB) domain, the pseudosubstrate, and linker domains and represents the minimal segments necessary for high affinity C-subunit binding (Fig. 1) (8). The other approach involves an examination of the effects of cAMP and its analogs on solvent exposure/conformational flexibility of RIα(91–244)·C-subunit complex using H/²H amide exchange measured with a combination of mass spectrometry (ESI-Q-TOF) and proteolytic fragmentation. In the first approach, fluorescein maleimide (FM) was conjugated to three cysteine substitution mutants with the substitution sites located near or within the pseudosubstrate sequence, the linker domain, or αC (R92C, T104C, and R239C, respectively) of RIα(91–244) (Fig. 1). The time-resolved fluorescence anisotropy results suggest that cAMP and (S_p)-cAMPS reduce the interaction of the RIα linker domain and αC with the two peripheral R-C interaction sites on the C-subunit (so-called Sites 2 and 3) without affecting the interaction of the pseudosubstrate sequence with the active site cleft (so-called Site 1). Because of limitations of the amide H/²H exchange experiments, only the effects of (R_p)-cAMPS on H/²H amide exchange in RIα(91–244)·C-subunit complex could be investigated. The results showed that (R_p)-cAMPS induces a relatively widespread increase in amide exchange, indicating limited unfolding and/or enhanced conformational flexibility that is propagated almost globally through the C-subunit and, at least, part of RIα. These conformational changes were accompanied by a 5-fold increase in the affinity of RIα(91–244) toward C-subunit, suggesting that, at least, some of the (R_p)-cAMPS effects are mediated by an increase in internal entropy. Finally, the (R_p)-cAMPS-induced increase in R-C affinity indicates that (R_p)-cAMPS is better described as an inverse agonist because the basal activity of the PKA holoenzyme should be decreased by (R_p)-cAMPS.Open in a separate windowFig. 1.Overview of PKA structure and cAMP analogs. A, domain organization of RIα showing the domain boundaries of RIα(91–244) where the pseudosubstrate in green is connected to CNB-A domain in blue by a linker segment. B, structure of RIα(91–244) in the C-subunit-bound conformation (Protein Data Bank code 1U7E (23)) showing the pseudosubstrate in green, linker in yellow, and helical subdomain comprising helices αN, αA, αB, and αC in blue and β-subdomain in tan. The PBC is in red. C, structure of the C·RIα(91–244) holoenzyme showing the C-subunit in tan and RIα(91–244) in blue. Sites for introduction of cysteines by site-directed mutagenesis are represented by red circles. The cAMP binding site (PBC) is in red. D, structure of cAMP showing the 2′-OH group and 3′–5′ phosphodiester bond. The exocyclic oxygens upon replacement with sulfur atoms to generate the (S_p)-cAMPS and (R_p)-cAMPS diastereomers are highlighted.     

Evidence from Liposome Encapsulation for Transport-Limited Microbial Metabolism of Solid Alkanes

The recalcitrance of xenobiotics may be caused by an absence of transforming enzymes or by their inability to enter microbial cells. A nondestructive method for differentiating between these two possibilities is described. The solid n-alkanes octadecane (C₁₈) and hexatriacontane (C₃₆) were encapsulated into phosphatidylcholine bilayers (liposomes). The uptake and metabolism rates of encapsulated and unencapsulated substrates were then compared. During 1 h at 25°C, a Pseudomonas isolate took up 1.3% of radiolabeled and unencapsulated C₁₈ (solid state) versus 23.5% of labeled and encapsulated C₁₈. Growth at 25°C occurred with an apparent k_s of 2453 ± 148 mg/liter. Liposome encapsulation decreased this K_s to 60 ± 12 mg/liter. At 34°C, growth on C₁₈ (liquid state) occurred with an apparent K_s of 819 ± 83 mg/liter and on the readily available carbon source succinate, K_s values were 80 ± 10 and 13 ± 7 mg/liter at 25 and 34°C, respectively. At 25°C, the isolate grew on C₃₆ with an apparent K_s of 2,698 ± 831 mg/liter. Liposome encapsulation decreased the K_s more than 60-fold to 41 ± 7 mg/liter, resulting in the complete utilization of 400 mg of C₃₆ per liter in 16 h. Since controls excluded the metabolic utilization of phosphatidylcholine, the results clearly identify transport limitation as the cause for C₃₆ recalcitrance.     

Thermal niche for in situ seed germination by Mediterranean mountain streams: model prediction and validation for Rhamnus persicifolia seeds

Background and Aims

Mediterranean mountain species face exacting ecological conditions of rainy, cold winters and arid, hot summers, which affect seed germination phenology. In this study, a soil heat sum model was used to predict field emergence of Rhamnus persicifolia, an endemic tree species living at the edge of mountain streams of central eastern Sardinia.

Methods

Seeds were incubated in the light at a range of temperatures (10–25 and 25/10 °C) after different periods (up to 3 months) of cold stratification at 5 °C. Base temperatures (T_b), and thermal times for 50 % germination (θ₅₀) were calculated. Seeds were also buried in the soil in two natural populations (Rio Correboi and Rio Olai), both underneath and outside the tree canopy, and exhumed at regular intervals. Soil temperatures were recorded using data loggers and soil heat sum (°Cd) was calculated on the basis of the estimated T_b and soil temperatures.

Key Results

Cold stratification released physiological dormancy (PD), increasing final germination and widening the range of germination temperatures, indicative of a Type 2 non-deep PD. T_b was reduced from 10·5 °C for non-stratified seeds to 2·7 °C for seeds cold stratified for 3 months. The best thermal time model was obtained by fitting probit germination against log °Cd. θ₅₀ was 2·6 log °Cd for untreated seeds and 2·17–2·19 log °Cd for stratified seeds. When θ₅₀ values were integrated with soil heat sum estimates, field emergence was predicted from March to April and confirmed through field observations.

Conclusions

T_b and θ₅₀ values facilitated model development of the thermal niche for in situ germination of R. persicifolia. These experimental approaches may be applied to model the natural regeneration patterns of other species growing on Mediterranean mountain waterways and of physiologically dormant species, with overwintering cold stratification requirement and spring germination.