Role of loop dynamics in thermal stability of mesophilic and thermophilic adenylosuccinate synthetase: a molecular dynamics and normal mode analysis study

Enzymes from thermophiles are poorly active at temperatures at which their mesophilic homologs exhibit high activity and attain corresponding active states at high temperatures. In this study, comparative molecular dynamics (MD) simulations, supplemented by normal mode analysis, have been performed on an enzyme Adenylosuccinate synthetase (AdSS) from E. coli (mesophilic) and P. horikoshii (thermophilic) systems to understand the effects of loop dynamics on thermal stability of AdSS. In mesophilic AdSS, both ligand binding and catalysis are facilitated through the coordinated movement of five loops on the protein. The simulation results suggest that thermophilic P. horikoshii preserves structure and catalytic function at high temperatures by using the movement of only a subset of loops (two out of five) for ligand binding and catalysis unlike its mesophilic counterpart in E. coli. The pre-arrangement of the catalytic residues in P. horikoshii is well-preserved and salt bridges remain stable at high temperature (363K). The simulations suggest a general mechanism (including pre-arrangement of catalytic residues, increased polar residue content, stable salt bridges, increased rigidity, and fewer loop movements) used by thermophilic enzymes to preserve structure and be catalytically active at elevated temperatures.     

An electrostatic basis for the stability of thermophilic proteins

Two factors provide key contributions to the stability of thermophilic proteins relative to their mesophilic homologues: electrostatic interactions of charged residues in the folded state and the dielectric response of the folded protein. The dielectric response for proteins in a "thermophilic series" globally modulates the thermal stability of its members, with the calculated dielectric constant for the protein increasing from mesophiles to hyperthermophiles. This variability results from differences in the distribution of charged residues on the surface of the protein, in agreement with structural and genetic observations. Furthermore, the contribution of electrostatic interactions to the stability of the folded state is more favorable for thermophilic proteins than for their mesophilic homologues. This leads to the conclusion that electrostatic interactions play an important role in determining the stability of proteins at high temperatures. The interplay between electrostatic interactions and dielectric response also provides further rationalization for the enhanced stability of thermophilic proteins with respect to cold-denaturation. Taken together, the distribution of charged residues and their fluctuations have been shown to be factors in modulating protein stability over the entire range of biologically relevant temperatures.     

Targeted random mutagenesis to identify functionally important residues in the D2 protein of photosystem II in Synechocystis sp. strain PCC 6803

To identify important residues in the D2 protein of photosystem II (PSII) in the cyanobacterium Synechocystis sp. strain PCC 6803, we randomly mutagenized a region of psbDI (coding for a 96-residue-long C-terminal part of D2) with sodium bisulfite. Mutagenized plasmids were introduced into a Synechocystis sp. strain PCC 6803 mutant that lacks both psbD genes, and mutants with impaired PSII function were selected. Nine D2 residues were identified that are important for PSII stability and/or function, as their mutation led to impairment of photoautotrophic growth. Five of these residues are likely to be involved in the formation of the Q(A)-binding niche; these are Ala249, Ser254, Gly258, Ala260, and His268. Three others (Gly278, Ser283, and Gly288) are in transmembrane alpha-helix E, and their alteration leads to destabilization of PSII but not to major functional alterations of the remaining centers, indicating that they are unlikely to interact directly with cofactors. In the C-terminal lumenal tail of D2, only one residue (Arg294) was identified as functionally important for PSII. However, from the number of mutants generated it is likely that most or all of the 70 residues that are susceptible to bisulfite mutagenesis have been altered at least once. The fact that mutations in most of these residues have not been picked up by our screening method suggests that these mutations led to a normal photoautotrophic phenotype. A novel method of intragenic complementation in Synechocystis sp. strain PCC 6803 was developed to facilitate genetic analysis of psbDI mutants containing several amino acid changes in the targeted domain. Recombination between genome copies in the same cell appears to be much more prevalent in Synechocystis sp. strain PCC 6803 than was generally assumed.     

Engineering a selectable marker for hyperthermophiles

Limited thermostability of antibiotic resistance markers has restricted genetic research in the field of extremely thermophilic Archaea and bacteria. In this study, we used directed evolution and selection in the thermophilic bacterium Thermus thermophilus HB27 to find thermostable variants of a bleomycin-binding protein from the mesophilic bacterium Streptoalloteichus hindustanus. In a single selection round, we identified eight clones bearing five types of double mutated genes that provided T. thermophilus transformants with bleomycin resistance at 77 degrees C, while the wild-type gene could only do so up to 65 degrees C. Only six different amino acid positions were altered, three of which were glycine residues. All variant proteins were produced in Escherichia coli and analyzed biochemically for thermal stability and functionality at high temperature. A synthetic mutant resistance gene with low GC content was designed that combined four substitutions. The encoded protein showed up to 17 degrees C increased thermostability and unfolded at 85 degrees C in the absence of bleomycin, whereas in its presence the protein unfolded at 100 degrees C. Despite these highly thermophilic properties, this mutant was still able to function normally at mesophilic temperatures in vivo. The mutant protein was co-crystallized with bleomycin, and the structure of the binary complex was determined to a resolution of 1.5 A. Detailed structural analysis revealed possible molecular mechanisms of thermostabilization and enhanced antibiotic binding, which included the introduction of an intersubunit hydrogen bond network, improved hydrophobic packing of surface indentations, reduction of loop flexibility, and alpha-helix stabilization. The potential applicability of the thermostable selection marker is discussed.     

The thermostability of DNA-binding protein HU from mesophilic, thermophilic, and extreme thermophilic bacteria

Based on primary structure comparison between four highly homologous DNA-binding proteins (HUs) displaying differential thermostability, we have employed in vitro site-directed mutagenesis to decipher their thermostability mechanism at the molecular level. The contribution of the 11 amino acids that differ between the thermophilic HUBst from Bacillus stearothermophilus (Tm = 61.6 degrees C) and the mesophilic HUBsu from Bacillus subtilis (Tm = 39.7 degrees C) was evaluated by replacing these amino acids in HUBst with their mesophilic counterparts. Among 11 amino acids, three residues, Gly-15, Glu-34, and Val-42, which are highly conserved in the thermophilic HUs, have been found to be responsible for the thermostability of HUBst. These amino acids in combination (HUBst-G15E/E34D/V42I) reduce the thermostability of the protein (Tm = 45.1 degrees C) at the level of its mesophilic homologue HUBsu. By replacing these amino acids in HUBsu with their thermophilic counterparts, the HUBsu-E15G/D34E/142V mutant was generated with thermostability (Tm = 57.8 degrees C) at the level of thermophilic HUBst. Employing the same strategy, we generated several mutants in the extremely thermophilic HUTmar from Thermotoga maritima (Tm = 80.5 degrees C), and obtained data consistent with the previous results. The triplet mutant HUTmar-G15E/E34D/V421 (Tm = 35.9 degrees C) converted the extremely thermophilic protein HUTmar to mesophilic. The various forms of HU proteins were overproduced in Escherichia coli, highly purified, and the thermostability of the mutants confirmed by circular dichroism spectroscopy. The results presented here were elucidated on the basis of the X-ray structure of HUBst and HUTmar (our unpublished results), and their mechanism was proposed at the molecular level. The results clearly show that three individual local interactions located at the helix-turn-helix part of the protein are responsible for the stability of HU proteins by acting cooperatively in a common mechanism for thermostability.     

Comparing Residue Clusters from Thermophilic and Mesophilic Enzymes Reveals Adaptive Mechanisms

Understanding how proteins adapt to function at high temperatures is important for deciphering the energetics that dictate protein stability and folding. While multiple principles important for thermostability have been identified, we lack a unified understanding of how internal protein structural and chemical environment determine qualitative or quantitative impact of evolutionary mutations. In this work we compare equivalent clusters of spatially neighboring residues between paired thermophilic and mesophilic homologues to evaluate adaptations under the selective pressure of high temperature. We find the residue clusters in thermophilic enzymes generally display improved atomic packing compared to mesophilic enzymes, in agreement with previous research. Unlike residue clusters from mesophilic enzymes, however, thermophilic residue clusters do not have significant cavities. In addition, anchor residues found in many clusters are highly conserved with respect to atomic packing between both thermophilic and mesophilic enzymes. Thus the improvements in atomic packing observed in thermophilic homologues are not derived from these anchor residues but from neighboring positions, which may serve to expand optimized protein core regions.     

Solution structure and backbone dynamics of the K18G/R82E Alicyclobacillus acidocaldarius thioredoxin mutant: a molecular analysis of its reduced thermal stability

No general strategy for thermostability has been yet established, because the extra stability of thermophiles appears to be the sum of different cumulative stabilizing interactions. In addition, the increase of conformational rigidity observed in many thermophilic proteins, which in some cases disappears when mesophilic and thermophilic proteins are compared at their respective physiological temperatures, suggests that evolutionary adaptation tends to maintain corresponding states with respect to conformational flexibility. In this study, we accomplished a structural analysis of the K18G/R82E Alicyclobacillus acidocaldarius thioredoxin (BacTrx) mutant, which has reduced heat resistance with respect to the thermostable wild-type. Furthermore, we have also achieved a detailed study, carried out at 25, 45, and 65 degrees C, of the backbone dynamics of both the BacTrx and its K18G/R82E mutant. Our findings clearly indicate that the insertion of the two mutations causes a loss of energetically favorable long-range interactions and renders the secondary structure elements of the double mutants more similar to those of the mesophilic Escherichia coli thioredoxin. Moreover, protein dynamics analysis shows that at room temperature the BacTrx, as well as the double mutant, are globally as rigid as the mesophilic thioredoxins; differently, at 65 degrees C, which is in the optimal growth temperature range of A. acidocaldarius, the wild-type retains its rigidity while the double mutant is characterized by a large increase of the amplitude of the internal motions. Finally, our research interestingly shows that fast motions on the pico- to nanosecond time scale are not detrimental to protein stability and provide an entropic stabilization of the native state. This study further confirms that protein thermostability is reached through diverse stabilizing interactions, which have the key role to maintain the structural folding stable and functional at the working temperature.     

Mesophilic and thermophilic biotreatment of BTEX-polluted air in reactors

This study compares the removal of a mixture of benzene, toluene, ethylbenzene, and all three xylene isomers (BTEX) in mesophilic and thermophilic (50 degrees C) bioreactors. In the mesophilic reactor fungi became dominant after long-term operation, while bacteria dominated in the thermophilic unit. Microbial acclimation was achieved by exposing the biofilters to initial BTEX loads of 2-15 g m(-3) h(-1), at an empty bed residence time of 96 s. After adaptation, the elimination capacities ranged from 3 to 188 g m(-3) h(-1), depending on the inlet load, for the mesophilic biofilter with removal efficiencies reaching 96%. On the other hand, in the thermophilic reactor the average removal efficiency was 83% with a maximum elimination capacity of 218 g m(-3) h(-1). There was a clear positive relationship between temperature gradients as well as CO(2) production and elimination capacities across the biofilters. The gas phase was sampled at different depths along the reactors observing that the percentage pollutant removal in each section was strongly dependant on the load applied. The fate of individual alkylbenzene compounds was checked, showing the unusually high biodegradation rate of benzene at high loads under thermophilic conditions (100%) compared to its very low removal in the mesophilic reactor at such load (<10%). Such difference was less pronounced for the other pollutants. After 210 days of operation, the dry biomass content for the mesophilic and thermophilic reactors were 0.300 and 0.114 g g(-1) (support), respectively, reaching higher removals under thermophilic conditions with a lower biomass accumulation, that is, lower pressure drop.     

Two types of functionally distinct NAD(P)H dehydrogenases in Synechocystis sp. strain PCC6803

The ndhD gene encodes a membrane protein component of NAD(P)H dehydrogenase. The genome of Synechocystis sp. PCC6803 contains 6 ndhD genes. Three mutants were constructed by disrupting highly homologous ndhD genes in pairs. Only the DeltandhD1/DeltandhD2 (DeltandhD1/D2) mutant was unable to grow under photoheterotrophic conditions and exhibited low respiration rate, although the mutant grew normally under photoautotrophic conditions in air. The DeltandhD3/DeltandhD4 (DeltandhD3/D4) mutant grew very slowly in air and did not take up CO(2). The results demonstrated the presence of two types of functionally distinct NAD(P)H dehydrogenases in Synechocystis PCC6803 cells. TheDeltandhD5/DeltandhD6 (DeltandhD5/D6) mutant grew like the wild-type strain. Under far-red light (>710 nm), the level of P700(+) was high in DeltandhD1/D2 and M55 (ndhB-less mutant) at low intensities. The capacity of Q(A) (tightly bound plastoquinone) reduction by plastoquinone pool, as measured by the fluorescence increase in darkness upon addition of KCN, was much less in DeltandhD1/D2 and M55 than in DeltandhD3/D4 and DeltandhD5/D6. We conclude that electrons from NADPH are transferred to the plastoquinone pool mainly by the NdhD1.NdhD2 type of NAD(P)H dehydrogenases.     

Conferring thermostability to mesophilic proteins through optimized electrostatic surfaces

Recently, there have been several experimental reports of proteins displaying appreciable stability gains through mutation of one or two amino acid residues. Here, we employ a simple theoretical model to quickly screen mutant structures for increased thermostability through optimization of the protein's electrostatic surface. Our results are able to reproduce the experimental observation that elimination of like-charge repulsions and creation of opposite-charge attractions on the protein surface is an efficient method to confer thermostability to a mesophilic protein. Using Poisson-Boltzmann electrostatics, we calculate relative protein stabilities for the exhaustive surface mutagenesis of the cold shock, RNase T1, and CheY proteins. Comparison with 25 experimentally characterized cold shock protein mutants reveals an average correlation of 0.86. The model is also quantitatively accurate when reproducing the experimental D49A and D49H mutant stabilities of RNase T1. This work represents the first comprehensive in silico screening of mutant candidates likely to confer thermostability to mesophilic proteins through optimization of surface electrostatics. Systematic single mutant, followed by double mutant, screening yields a limited number of mutant structures displaying significant stability gains suitable for subsequent experimental characterization.     

Site-specific mutations of the N-terminal threonines in the D1 protein affects photoautotrophic growth but not D1 protein stability in Synechocystis 6803

The cyanobacterium Synechocystis 6803 was engineered to produce a D1 protein where one or more of the N-terminal threonines at positions 2, 3 and 4 were replaced by other amino acid residues. No phenotypic effects were found for the T2S or T2L mutations, whereas the T2V, T2L;T4V and T2V;T3V;T4V mutations resulted in reduced photoautotrophic growth rate and oxygen evolving activity. The mutant strain T2V;T3V;T4V exhibited an oxygen evolution activity that was only half of that for the wild-type strain. Despite of that, both accumulation and stability of the D1 protein in the thylakoid membrane appeared unaffected in the mutant.     

Selected mutations in a mesophilic cytochrome c confer the stability of a thermophilic counterpart

Mesophilic cytochrome c(551) of Pseudomonas aeruginosa (PA c(551)) became as stable as its thermophilic counterpart, Hydrogenobacter thermophilus cytochrome c(552) (HT c(552)), through only five amino acid substitutions. The five residues, distributed in three spatially separated regions, were selected and mutated with reference to the corresponding residues in HT c(552) through careful structure comparison. Thermodynamic analysis indicated that the stability of the quintuple mutant of PA c(551) could be partly attained through an enthalpic factor. The solution structure of the mutant showed that, as in HT c(552), there were tighter side chain packings in the mutated regions. Furthermore, the mutant had an increased total accessible surface area, resulting in great negative hydration free energy. Our results provide a novel example of protein stabilization in that limited amino acid substitutions can confer the overall stability of a natural highly thermophilic protein upon a mesophilic molecule.     

Two exposed amino acid residues confer thermostability on a cold shock protein

Thermophilic organisms produce proteins of exceptional stability. To understand protein thermostability at the molecular level we studied a pair of cold shock proteins, one of mesophilic and one of thermophilic origin, by systematic mutagenesis. Although the two proteins differ in sequence at 12 positions, two surface-exposed residues are responsible for the increase in stability of the thermophilic protein (by 15.8 kJ mol-1 at 70 degrees C). 11.5 kJ mol-1 originate from a predominantly electrostatic contribution of Arg 3 and 5.2 kJ mol-1 from hydrophobic interactions of Leu 66 at the carboxy terminus. The mesophilic protein could be converted to a highly thermostable form by changing the Glu residues at positions 3 and 66 to Arg and Leu, respectively. The variation of surface residues may thus provide a simple and powerful approach for increasing the thermostability of a protein.     

Site-directed mutagenesis of glutamate residues in the large extrinsic loop of the photosystem II protein CP 43 affects oxygen-evolving activity and PS II assembly

The psbC gene encodes the intrinsic chlorophyll protein CP 43, a component of photosystem II in higher plants, green algae, and cyanobacteria. Oligonucleotide-directed mutagenesis was used to introduce mutations into the portion of psbC that encodes the large extrinsic loop E of CP 43 in the cyanobacterium Synechocystis 6803. Three mutations, E293Q, E339Q, and E352Q, each produced a strain with impaired photosystem II activity. The E293Q mutant strain grew photoautotrophically at rates comparable to the control strain. Immunological analyses of several PS II components indicated that this mutant accumulated normal quantities of PS II proteins. However, this mutant evolved oxygen to only 56% of control rates at saturating light intensities. Measurements of total variable fluorescence yield indicated that this mutant assembled approximately 60% of the fully functional PS II centers found in the control strain. The E339Q mutant grew photoautotrophically at a severely reduced rate. Both immunological analysis and variable fluorescence yield experiments indicated that E339Q assembled a normal complement of PS II centers. However, this mutant was capable of evolving oxygen to only 20% of control rates. Variable fluorescence yield experiments demonstrated that this mutant was inefficient at using water as an electron donor. Both E293Q and E339Q strains exhibited an increased (approximately 2-fold) sensitivity to photoinactivation. The E352Q mutant was the most severely affected. This mutant failed to grow photoautotrophically and exhibited essentially no capacity for oxygen evolution. Measurements of total variable fluorescence yield indicated that this mutant assembled no functional PS II centers. Immunological analysis of isolated thylakoid membranes from E352Q revealed a complete absence of CP 43 and reduced levels of both the D1 and manganese-stabilizing proteins. These results suggest that the mutations E293Q and E339Q each produce a defect associated with the oxygen-evolving complex of photosystem II. The E352Q mutation appears to affect the stability of the PS II complex. This is the first report showing that alteration of negatively charged residues in the CP 43 large extrinsic loop results in mutations affecting PS II assembly/function.     

A single proline substitution is critical for the thermostabilization of Clostridium beijerinckii alcohol dehydrogenase

Analysis of the three-dimensional structures of three closely related mesophilic, thermophilic, and hyperthermophilic alcohol dehydrogenases (ADHs) from the respective microorganisms Clostridium beijerinckii (CbADH), Entamoeba histolytica (EhADH1), and Thermoanaerobacter brockii (TbADH) suggested that a unique, strategically located proline residue (Pro100) might be crucial for maintaining the thermal stability of EhADH1. To determine whether proline substitution at this position in TbADH and CbADH would affect thermal stability, we used site-directed mutagenesis to replace the complementary residues in both enzymes with proline. The results showed that replacing Gln100 with proline significantly enhanced the thermal stability of the mesophilic ADH: DeltaT(1/2) (60 min) = + 8 degrees C (temperature of 50% inactivation after incubation for 60 min), DeltaT(1/2) (CD) = +11.5 degrees C (temperature at which 50% of the original CD signal at 218 nm is lost upon heating between 30 degrees and 98 degrees C). A His100 --> Pro substitution in the thermophilic TbADH had no effect on its thermostability. An analysis of the three-dimensional structure of the crystallized thermostable mutant Q100P-CbADH suggested that the proline residue at position 100 stabilized the enzyme by reinforcing hydrophobic interactions and by reducing the flexibility of a loop at this strategic region.     

Stabilization of mesophilic Allochromatium vinosum cytochrome c′ through specific mutations modeled by a thermophilic homologue

AVCP cytochrome c′ from mesophilic Allochromatium vinosum exhibits lower stability than a thermophilic counterpart, Hydrogenophilus thermoluteolus cytochrome c′ (PHCP), in which the six specific amino acid residues that are not conserved in AVCP are responsible for its stability. Here we measured the stability of AVCP variants carrying these specific residues instead of the original AVCP ones. Among the six single AVCP variants, all of which formed a dimeric structure similar to that of the wild-type, three were successfully stabilized compared with the wild-type, while one showed lower stability than the wild-type. In addition, the most stabilized and destabilized AVCP variants could bind CO, similar to the wild-type. These results indicated that mesophilic AVCP could be stabilized through specific three mutations modeled by the thermophilic counterpart, PHCP, without changing the CO binding ability.     

Different packing of external residues can explain differences in the thermostability of proteins from thermophilic and mesophilic organisms

MOTIVATION: Understanding the basis of protein stability in thermophilic organisms raises a general question: what structural properties of proteins are responsible for the higher thermostability of proteins from thermophilic organisms compared to proteins from mesophilic organisms? RESULTS: A unique database of 373 structurally well-aligned protein pairs from thermophilic and mesophilic organisms is constructed. Comparison of proteins from thermophilic and mesophilic organisms has shown that the external, water-accessible residues of the first group are more closely packed than those of the second. Packing of interior parts of proteins (residues inaccessible to water molecules) is the same in both cases. The analysis of amino acid composition of external residues of proteins from thermophilic organisms revealed an increased fraction of such amino acids as Lys, Arg and Glu, and a decreased fraction of Ala, Asp, Asn, Gln, Thr, Ser and His. Our theoretical investigation of folding/unfolding behavior confirms the experimental observations that the interactions that differ in thermophilic and mesophilic proteins form only after the passing of the transition state during folding. Thus, different packing of external residues can explain differences in thermostability of proteins from thermophilic and mesophilic organisms. AVAILABILITY: The database of 373 structurally well-aligned protein pairs is available at http://phys.protres.ru/resources/termo_meso_base.html. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.     

Thermodynamic differences among homologous thermophilic and mesophilic proteins.

Here, we analyze the thermodynamic parameters and their correlations in families containing homologous thermophilic and mesophilic proteins which show reversible two-state folding <--> unfolding transitions between the native and the denatured states. For the proteins in these families, the melting temperatures correlate with the maximal protein stability change (between the native and the denatured states) as well as with the enthalpic and entropic changes at the melting temperature. In contrast, the heat capacity change is uncorrelated with the melting temperature. These and additional results illustrate that higher melting temperatures are largely obtained via an upshift and broadening of the protein stability curves. Both thermophilic and mesophilic proteins are maximally stable around room temperature. However, the maximal stabilities of thermophilic proteins are considerably greater than those of their mesophilic homologues. At the living temperatures of their respective source organisms, homologous thermophilic and mesophilic proteins have similar stabilities. The protein stability at the living temperature of the source organism does not correlate with the living temperature of the protein. We tie thermodynamic observations to microscopics via the hydrophobic effect and a two-state model of the water structure. We conclude that, to achieve higher stability and greater resistance to high and low temperatures, specific interactions, particularly electrostatic, should be engineered into the protein. The effect of these specific interactions is largely reflected in an increased enthalpy change at the melting temperature.     

Lys13 plays a crucial role in the functional adaptation of the thermophilic triose-phosphate isomerase from Bacillus stearothermophilus to high temperatures.

The thermophilic triose-phosphate isomerases (TIMs) of Bacillus stearothermophilus (bTIM) and Thermotoga maritima (tTIM) have been found to possess a His12-Lys13 pair instead of the Asn12-Gly13 pair normally present in mesophilic TIMs. His12 in bTIM was proposed to prevent deamidation at high temperature, while the precise role of Lys13 is unknown. To investigate the role of the His12 and Lys13 pair in the enzyme's thermoadaptation, we reintroduced the "mesophilic residues" Asn and Gly into both thermophilic TIMs. Neither double mutant displayed diminished structural stability, but the bTIM double mutant showed drastically reduced catalytic activity. No similar behavior was observed with the tTIM double mutant, suggesting that the presence of the His12 and Lys13 cannot be systematically correlated to thermoadaptation in TIMs. We determined the crystal structure of the bTIM double mutant complexed with 2-phosphoglycolate to 2.4-A resolution. A molecular dynamics simulation showed that upon substitution of Lys13 to Gly an increase of the flexibility of loop 1 is observed, causing an incorrect orientation of the catalytic Lys10. This suggests that Lys13 in bTIM plays a crucial role in the functional adaptation of this enzyme to high temperature. Analysis of bTIM single mutants supports this assumption.