A single amino acid mutation enhances the thermal stability of Escherichia coli malate dehydrogenase

A conserved arginine residue in a staphylococcal anti-sigma factor is required to preserve its kinase activity, structure, and stability

RsbW, σ^B, and RsbV, encoded by Staphylococcus aureus and related bacteria, act as an anti-sigma factor, an sigma factor, and an anti-anti-sigma factor, respectively. The interaction between RsbW and σ^B blocks the transcription initiation activity of the latter protein. RsbW also functions as a serine kinase and phosphorylates RsbV in the presence of ATP. Our modeling study indicates that the RsbW-RsbV complex is stabilized by twenty-four intermolecular non-covalent bonds. Of the bond-forming RsbW residues, Arg 23, and Glu 49 are conserved residues. To understand the roles of Arg 23 in RsbW, rRsbW[R23A], a recombinant S. aureus RsbW (rRsbW) harboring Arg to Ala change at position 23, was investigated using various probes. The results reveal that rRsbW[R23A], like rRsbW, exists as the dimers in the aqueous solution. However, rRsbW[R23A], unlike rRsbW, neither interacted with a chimeric RsbV (rRsbV) nor formed the phosphorylated rRsbV in the presence of ATP. Furthermore, the tertiary structure and hydrophobic surface area of rRsbW[R23A] matched little with those of rRsbW. Conversely, both rRsbW[R23A] and rRsbW showed interaction with a recombinant σ^B (rσ^B). rRsbW and rRsbW[R23A] were also unfolded via the formation of at least one intermediate in the presence of urea. However, the thermodynamic stability of rRsbW significantly differed from that of rRsbW[R23A]. Our molecular dynamics (MD) simulation study also reveals the substantial change of structure, dimension, and stability of RsbW due to the above mutation. The ways side chain of critical Arg 23 contributes to maintaining the tertiary structure, and stability of RsbW was elaborately discussed.Communicated by Ramaswamy H. Sarma.

Investigation of the proton relay system operative in human cystosolic aminopeptidase P

Aminopeptidase P, a metalloprotease, targets Xaa-Proline peptides for cleavage [1-4]. There are two forms of human AMPP, a membrane-bound form (hmAMPP) and a soluble cytosolic form (hcAMPP)[5]. Similar to the angiotensin-I-converting enzyme, AMPP plays an important role in the catabolism of inflammatory and vasoactive peptides, known as kinins. The plasma kinin, bradykinin, was used as the substrate to conduct enzymatic activity analyses and to determine the Michaelis constant (Km) of 174 μM and the catalytic rate constant (kcat) of 10.8 s-1 for hcAMPP. Significant differences were observed in the activities of Y527F and R535A hcAMPP mutants, which displayed a 6-fold and 13.5-fold for decrease in turnover rate, respectively. Guanidine hydrochloride restored the activity of R535A hcAMPP, increasing the kcat/Km 20-fold, yet it had no impact on the activities of the wild-type or Y527F mutant hcAMPPs. Activity restoration by guanidine derivatives followed the order guanidine hydrochloride >> methyl-guanidine > amino-guanidine > N-ethyl-guanidine. Overall, the results indicate the participation of R535 in the hydrogen bond network that forms a proton relay system. The quaternary structure of hcAMPP was determined by using analytical ultracentrifugation (AUC). The results show that alanine replacement of Arg535 destabilizes the hcAMPP dimer and that guanidine hydrochloride restores the native monomer-dimer equilibrium. It is proposed that Arg535 plays an important role in hcAMMP catalysis and in stabilization of the catalytically active dimeric state.

Function, kinetic properties, crystallization, and regulation of microbial malate dehydrogenase*

Malate dehydrogenase (MDH) is an enzyme widely distributed among living organisms and is a key protein in the central oxidative pathway. It catalyzes the interconversion between malate and oxaloacetate using NAD+ or NADP+ as a cofactor. Surprisingly, this enzyme has been extensively studied in eukaryotes but there are few reports about this enzyme in prokaryotes. It is necessary to review the relevant information to gain a better understanding of the function of this enzyme. Our review of the data generated from studies in bacteria shows much diversity in their molecular properties, including weight, oligomeric states, cofactor and substrate binding affinities, as well as differences in the direction of the enzymatic reaction. Furthermore, due to the importance of its function, the transcription and activity of this enzyme are rigorously regulated. Crystal structures of MDH from different bacterial sources led to the identification of the regions involved in substrate and cofactor binding and the residues important for the dimer-dimer interface. This structural information allows one to make direct modifications to improve the enzyme catalysis by increasing its activity, cofactor binding capacity, substrate specificity, and thermostability. A comparative analysis of the phylogenetic reconstruction of MDH reveals interesting facts about its evolutionary history, dividing this superfamily of proteins into two principle clades and establishing relationships between MDHs from different cellular compartments from archaea, bacteria, and eukaryotes.

Engineering of a metagenome derived lipase toward thermal tolerance: effect of asparagine to lysine mutation on the protein surface

A highly thermostable mutant lipase was generated and characterized. Mutant enzyme demonstrated 144 fold enhanced thermostability over the wild type enzyme at 60°C. Interestingly, the overall catalytic efficiency (k(cat/)K(m)) of mutant was also enhanced (~20 folds). Circular dichroism spectroscopy, studied as function of temperature, demonstrated that the mutant lipase retained its secondary structure up to 70-80°C, whereas wild type protein structure was completely distorted above 35°C. Additionally, the intrinsic tryptophan fluorescence (a probe for the tertiary structure) also displayed difference in the conformation of two enzymes during temperature dependent unfolding. Furthermore, mutation N355K resulted in extensive H-bonding (Lys355 HZ1OE2 Glu284) with a distance 2.44 Å. In contrast to this, Wt enzyme has not shown such H-bonding interaction.

The role of rhodopsin glycosylation in protein folding, trafficking, and light-sensitive retinal degeneration

Several mutations in the N terminus of the G-protein-coupled receptor rhodopsin disrupt NXS/T consensus sequences for N-linked glycosylation (located at N2 and N15) and cause sector retinitis pigmentosa in which the inferior retina preferentially degenerates. Here we examined the role of rhodopsin glycosylation in biosynthesis, trafficking, and retinal degeneration (RD) using transgenic Xenopus laevis expressing glycosylation-defective human rhodopsin mutants. Although mutations T4K and T4N caused RD, N2S and T4V did not, demonstrating that glycosylation at N2 was not required for photoreceptor viability. In contrast, similar mutations eliminating glycosylation at N15 (N15S and T17M) caused rod death. Expression of T17M was more toxic than T4K to transgenic photoreceptors, further suggesting that glycosylation at N15 plays a more important physiological role than glycosylation at N2. Together, these results indicate that the structure of the rhodopsin N terminus must be maintained by an appropriate amino acid sequence surrounding N2 and may require a carbohydrate moiety at N15. The mutant rhodopsins were rendered less toxic in their dark inactive states, because RD was abolished or significantly reduced when transgenic tadpoles expressing T4K, T17M, and N2S/N15S were protected from light exposure. Regardless of their effect on rod viability, all of the mutants primarily localized to the outer segment and Golgi and showed little or no endoplasmic reticulum accumulation. Thus, glycosylation was not crucial for rhodopsin biosynthesis or trafficking. Interestingly, expression of similar bovine rhodopsin mutants did not cause rod cell death, possibly attributable to greater stability of bovine rhodopsin.

Thermotolerance of L3 Ostertagia (Teladorsagia) circumcincta and some of its enzymes

Parasitic nematodes of ruminants can be expected to experience temperatures in excess of 40 degrees C in faeces on pasture and, perhaps, in the host. L3 Ostertagia (Teladorsagia) circumcincta survived for at least 90 min at 45 degrees C in vitro in water, but the larvae were inactivated rapidly by only slightly higher temperatures. The glycolytic enzymes hexokinase and pyruvate kinase were inactivated in a similar temperature range, whereas malate dehydrogenase maintained its activity at temperatures in excess of 50 degrees C. These data imply that the loss of glycolytic activity might explain the loss of larval motility at temperatures between 45 degrees C and 50 degrees C.

Comparative study of the stability of poplar plastocyanin isoforms

The stability of the two isoforms of poplar plastocyanin (PCa and PCb) was studied with differential scanning calorimetry (DSC) technique. It was shown that the thermal unfolding of both isoforms is an irreversible process with two endothermic and one exothermic peaks. The melting temperature of PCb was found to be 1.3+/-0.2 K degrees higher than of PCa, which indicates that PCb is more stable. The enthalpy of unfolding was estimated from the heat capacity curves and was found to be significantly higher for PCb at salt concentration I=0.1 M. In addition, PCb unfolding enthalpy and melting temperature are much more sensitive to the changes in the salt concentration as found in the experiments done at different ionic strength. The experiments were complemented with numerical calculations. The salt effect on the stability was modeled using the X-ray structure of PCa and a homology modeled structure of PCb. It was found, in agreement with the experimental data, that the stability of PCb changes by 4.7 kJ more than PCa, as the salt concentration increases from zero to 0.1 M. Thus, the differences in only 12 amino acid positions between "a" and "b" isoforms result in a measurable difference in the folding enthalpy and a significant difference in the salt dependence. The optimization of the electrostatic energies of PCa and PCb were studied and it was shown that PCb is better electrostatically optimized.

Contribution of engineered electrostatic interactions to the stability of cytosolic malate dehydrogenase

Protein engineering is a promising tool to obtain stable proteins. Comparison between homologous thermophilic and mesophilic enzymes from a given structural family can reveal structural features responsible for the enhanced stability of thermophilic proteins. Structures from pig heart cytosolic and Thermus flavus malate dehydrogenases (cMDH, Tf MDH), two proteins showing a 55% sequence homology, were compared with the aim of increasing cMDH stability using features from the Thermus flavus enzyme. Three potential salt bridges from Tf MDH were selected on the basis of their location in the protein (surface R176-D200, inter-subunit E57-K168 and intrasubunit R149-E275) and implemented on cMDH using site-directed mutagenesis. Mutants containing E275 were not produced in any detectable amount, which shows that the energy penalty of introducing a charge imbalance in a region that was not exposed to solvent was too unfavourable to allow proper folding of the protein. The salt bridge R149-E275, if formed, would not enhance stability enough to overcome this effect. The remaining mutants were expressed and active and no differences from wild-type other than stability were found. Of the mutants assayed, Q57E/L168K led to a stability increase of 0.4 kcal/mol, as determined by either guanidinium chloride denaturalization or thermal inactivation experiments. This results in a 15 degrees C shift in the optimal temperature, thus confirming that the inter-subunit salt bridge initially present in the T.flavus enzyme was formed in the cMDH structure and that the extra energy obtained is transformed into an increase in protein stability. These results indicate that the use of structural features of thermophilic enzymes, revealed by a detailed comparison of three-dimensional structures, is a valid strategy to improve the stability of mesophilic malate dehydrogenases.

Sulfolobus acidocaldarius inorganic pyrophosphatase: structure, thermostability, and effect of metal ion in an archael pyrophosphatase

The first crystal structure of an inorganic pyrophosphatase (S-PPase) from an archaebacterium, the thermophile Sulfolobus acidocaldarius, has been solved by molecular replacement and refined to an R-factor of 19.7% at 2.7 A. S-PPase is a D3 homohexameric protein with one Mg2+ per active site in a position similar to, but not identical with, the first activating metal in mesophilic pyrophosphatases (PPase). In mesophilic PPases, Asp65, Asp70, and Asp102 coordinate the Mg2+, while only Asp65 and Asp102 do in S-PPase, and the Mg2+ moves by 0.7 A. S-PPase may therefore be deactivated at low temperature by mispositioning a key metal ion. The monomer S-PPase structure is very similar to that of Thermus thermophilus (T-PPase) and Escherichia coli (E-PPase), root-mean-square deviations around 1 A/Calpha. But the hexamer structures of S- and T-PPase are more tightly packed and more similar to each other than they are to that of E-PPase, as shown by the increase in surface area buried upon oligomerization. In T-PPase, Arg116 creates an interlocking ionic network to both twofold and threefold related monomers; S-PPase has hydrophilic interactions to threefold related monomers absent in both E- and T-PPase. In addition, the thermostable PPases have about 7% more hydrogen bonds per monomer than E-PPase, and, especially in S-PPase, additional ionic interactions anchor the C-terminus to the rest of the protein. Thermostability in PPases is thus due to subtle improvements in both monomer and oligomer interactions.

Analysis of the dynamics of rhizomucor miehei lipase at different temperatures

The dynamics of Rhizomucor miehei lipase has been studied by molecular dynamics simulations at temperatures ranging from 200-500K. Simulations carried out in periodic boundary conditions and using explicit water molecules were performed for 400 ps at each temperature. Our results indicate that conformational changes and internal motions in the protein are significantly influenced by the temperature increase. With increasing temperature, the number of internal hydrogen bonds decreases, while surface accessibility, radius of gyration and the number of residues in random coil conformation increase. In the temperature range studied, the motions can be described in a low dimensional subspace, whose dimensionality decreases with increasing temperature. Approximately 80% of the total motion is described by the first (i) 80 eigenvectors at T=200K, (ii) 30 eigenvectors at T=300K and (iii) 10 eigenvectors at T=400K. At high temperature, the alpha-helix covering the active site in the native Rhizomucor miehei lipase, the helix at which end the active site is located, and in particular, the loop (Gly35-Lys50) show extensive flexibility.

Mutational analysis of the major loop of Bacillus 1,3-1,4-beta-D-glucan 4-glucanohydrolases. Effects on protein stability and substrate binding

The carbohydrate-binding cleft of Bacillus licheniformis 1,3-1, 4-beta-D-glucan 4-glucanohydrolase is partially covered by the surface loop between residues 51 and 67, which is linked to beta-strand-(87-95) of the minor beta-sheet III of the protein core by a single disulfide bond at Cys61-Cys90. An alanine scanning mutagenesis approach has been applied to analyze the role of loop residues from Asp51 to Arg64 in substrate binding and stability by means of equilibrium urea denaturation, enzyme thermotolerance, and kinetics. The DeltaDeltaGU between oxidized and reduced forms is approximately constant for all mutants, with a contribution of 5.3 +/- 0.2 kcal.mol-1 for the disulfide bridge to protein stability. A good correlation is observed between DeltaGU values by reversible unfolding and enzyme thermotolerance. The N57A mutant, however, is more thermotolerant than the wild-type enzyme, whereas it is slightly less stable to reversible urea denaturation. Mutants with a <2-fold increase in Km correspond to mutations at residues not involved in substrate binding, for which the reduction in catalytic efficiency (kcat/Km) is proportional to the loss of stability relative to the wild-type enzyme. Y53A, N55A, F59A, and W63A, on the other hand, show a pronounced effect on catalytic efficiency, with Km > 2-fold and kcat < 5% of the wild-type values. These mutated residues are directly involved in substrate binding or in hydrophobic packing of the loop. Interestingly, the mutation M58A yields an enzyme that is more active than the wild-type enzyme (7-fold increase in kcat), but it is slightly less stable.

Malate dehydrogenase from the mesophile Chlorobium vibrioforme and from the mild thermophile Chlorobium tepidum: molecular cloning, construction of a hybrid, and expression in Escherichia coli

The genes (mdh) encoding malate dehydrogenase (MDH) from the mesophile Chlorobium vibrioforme and the moderate thermophile C. tepidum were cloned and sequenced, and the complete amino acid sequences were deduced. When the region upstream of mdh was analyzed, a sequence with high homology to an operon encoding ribosomal proteins from Escherichia coli was found. Each mdh gene consists of a 930-bp open reading frame and encodes 310 amino acid residues, corresponding to a subunit weight of 33,200 Da for the dimeric enzyme. The amino acid sequence identity of the two MDHs is 86%. Homology searches using the primary structures of the two MDHs revealed significant sequence similarity to lactate dehydrogenases. A hybrid mdh was constructed from the 3' part of mdh from C. tepidum and the 5' part of mdh from C. vibrioforme. The thermostabilities of the hybrid enzyme and of MDH from C. vibrioforme and C. tepidum were compared.

Tetrameric malate dehydrogenase from a thermophilic Bacillus: cloning, sequence and overexpression of the gene encoding the enzyme and isolation and characterization of the recombinant enzyme

The gene encoding the tetrameric malate dehydrogenase (MDH) in a thermophilic Bacillus species (BI) has been cloned in an Escherichia coli plasmid. The nucleotide sequence of the gene, the first to be elucidated for a tetrameric MDH, shows the MDH subunit to contain 312 amino acids and have a molecular mass of 33648 Da, which confirms the experimentally determined value of about 35 kDa. Like the genomic DNA of BI, the MDH gene is relatively AT-rich; this contrasts with the generally GC-rich nature of the DNA of thermophilic Bacillus species. Comparison of amino acid sequences reveals that BI MDH bears greater structural similarity to lactate dehydrogenases (LDHs) than to other (dimeric) MDHs. MDHs and LDHs resemble each other in catalytic mechanism and several other respects. However, whereas MDHs in the majority of organisms are dimers, the tetrameric structure is favoured among LDHs. The stronger structural resemblance that BI MDH has to LDHs than to the dimeric MDHs provides some explanation as to why Bacillus MDH, unlike most other MDHs, is tetrameric. A 1 kb fragment containing the BI MDH gene, produced in a PCR, has been cloned into a high-expression E. coli plasmid vector. BI MDH synthesized from this clone constitutes about 47% of the total protein in cell extracts of the E. coli strain carrying the clone. MDH purified from BI and that purified from the E. coli strain carrying the MDH gene clone appear to be identical proteins by several criteria. A number of characteristics of the MDH have been elucidated, including the molecular masses of the native enzyme and the subunit, N-terminal amino acid sequence, isoelectric point, pH optimum for activity, thermostability, stability to pH, urea and guanidinium chloride and several kinetic parameters. Whereas the MDH is a stable tetramer in the pH range 5-7, it appears to be converted into a stable dimer at pH 3.5. This suggests that the dimer is a stable intermediate in the dissociation of the tetramer to monomers at low pH.

A Bacillus subtilis malate dehydrogenase gene

A Bacillus subtilis gene for malate dehydrogenase (citH) was found downstream of genes for citrate synthase and isocitrate dehydrogenase. Disruption of citH caused partial auxotrophy for aspartate and a requirement for aspartate during sporulation. In the absence of aspartate, citH mutant cells were blocked at a late stage of spore formation.

The optimization of protein-solvent interactions: thermostability and the role of hydrophobic and electrostatic interactions

Protein-solvent interactions were analyzed using an optimization parameter based on the ratio of the solvent-accessible area in the native and the unfolded protein structure. The calculations were performed for a set of 183 nonhomologous proteins with known three-dimensional structure available in the Protein Data Bank. The dependence of the total solvent-accessible surface area on the protein molecular mass was analyzed. It was shown that there is no difference between the monomeric and oligomeric proteins with respect to the solvent-accessible area. The results also suggested that for proteins with molecular mass above some critical mass, which is about 28 kDa, a formation of domain structure or subunit aggregation into oligomers is preferred rather than a further enlargement of a single domain structure. An analysis of the optimization of both protein-solvent and charge-charge interactions was performed for 14 proteins from thermophilic organisms. The comparison of the optimization parameters calculated for proteins from thermophiles and mesophiles showed that the former are generally characterized by a high degree of optimization of the hydrophobic interactions or, in cases where the optimization of the hydrophobic interactions is not sufficiently high, by highly optimized charge-charge interactions.