Comparative study on dihydrofolate reductases from Shewanella species living in deep-sea and ambient atmospheric-pressure environments

Contribution of a surface salt bridge to the protein stability of deep-sea Shewanella benthica cytochrome c'

Two homologous cytochromes c', SBCP and SVCP, from deep-sea Shewanella benthica and Shewanella violacea respectively exhibit only nine surface amino acid substitutions, along with one at the N-terminus. Despite the small sequence difference, SBCP is thermally more stable than SVCP. Here, we examined the thermal stability of SBCP variants, each containing one of the nine substituted residues in SVCP, and found that the SBCP K87V variant was the most destabilized. We then determined the X-ray crystal structure of the SBCP K87V variant at a resolution of 2.1 Å. The variant retains a four-helix bundle structure similar to the wild-type, but notable differences are observed in the hydration structure around the mutation site. Instead of forming of the intrahelical salt bridge between Lys-87 and Asp-91 in the wild-type, a clathrate-like hydration around Val-87 through a hydrogen bond network with the nearby amino acid residues is observed. This network potentially enhances the ordering of surrounding water molecules, leading to an entropic destabilization of the protein. These results suggest that the unfavorable hydrophobic hydration environment around Val-87 and the inability to form the Asp-91-mediated salt bridge contribute to the observed difference in stability between SBCP and SVCP. These findings will be useful in future protein engineering for controlling protein stability through the manipulation of surface intrahelical salt bridges.

Study on spoilage potential and its molecular basis of Shewanella putrefaciens in response to cold conditions by Label-free quantitative proteomic analysis

To elucidate how Shewanella putrefaciens survives and produces spoilage products in response to cold conditions, the metabolic and protease activity of S. putrefaciens DSM6067 cultured at three different temperatures (30 °C, 10 °C, and 4 °C) was studied by determining the bacterial growth, total volatile basic nitrogen (TVB-N), biogenic amines, extracellular protease activity, as well as the differential expressed proteins via Label-free quantitative proteomics analysis. The lag phase of the strain cultured at 10 °C and 4 °C was about 20 h and 120 h longer than at 30 °C, respectively. The TVB-N increased to 89.23 mg N/100 g within 28 h at 30 °C, and it needed at least 72 h and 224 h at 10 °C and 4 °C, respectively. Cold temperatures (10 °C and 4 °C) also inhibited the yield factors and the extracellular protease activity per cell at the lag phase. However, the protease activity per cell and the yield factors of the sample cultivated at 10 °C and 4 °C well recovered, especially at the mid and latter stages of the log phase. The further quantitative proteomic analysis displayed a complex biological network to tackle cold stress: cold stress responses, nutrient uptake, and energy conservation strategy. It was observed that the protease and peptidase were upregulated, so as to the degradation pathways of serine, arginine, and aspartate, which might lead to the accumulation of spoilage products. This study highlighted the spoilage potential of S. putrefaciens still should be concerned even at low temperatures.

Pressure Adaptations in Deep-Sea Moritella Dihydrofolate Reductases: Compressibility versus Stability

Simple Summary

Deep-sea organisms must have proteins that function under high hydrostatic pressure to survive. Adaptations used in proteins from “pressure-loving” piezophiles may include greater compressibility or greater stability against pressure-induced destabilization. However, while greater compressibility can be accomplished by greater void volume, larger cavities in a protein have been associated with greater destabilization and even unfolding as pressure is increased. Here, computer simulations of dihydrofolate reductase from a moderate piezophile and a hyperpiezophile were performed to understand the balance between adaptations for greater compressibility and those against pressure destabilization and unfolding. The results indicate that while compressibility appears to be important for deep-sea microbes, adaptation for the greatest depths may be to prevent water penetration into the interior.

Abstract

Proteins from “pressure-loving” piezophiles appear to adapt by greater compressibility via larger total cavity volume. However, larger cavities in proteins have been associated with lower unfolding pressures. Here, dihydrofolate reductase (DHFR) from a moderate piezophile Moritella profunda (Mp) isolated at ~2.9 km in depth and from a hyperpiezophile Moritella yayanosii (My) isolated at ~11 km in depth were compared using molecular dynamics simulations. Although previous simulations indicate that MpDHFR is more compressible than a mesophile DHFR, here the average properties and a quasiharmonic analysis indicate that MpDHFR and MyDHFR have similar compressibilities. A cavity analysis also indicates that the three unique mutations in MyDHFR are near cavities, although the cavities are generally similar in size in both. However, while a cleft overlaps an internal cavity, thus forming a pathway from the surface to the interior in MpDHFR, the unique residue Tyr103 found in MyDHFR forms a hydrogen bond with Leu78, and the sidechain separates the cleft from the cavity. Thus, while Moritella DHFR may generally be well suited to high-pressure environments because of their greater compressibility, adaptation for greater depths may be to prevent water entry into the interior cavities.

The Molecular Basis for Life in Extreme Environments

Sampling and genomic efforts over the past decade have revealed an enormous quantity and diversity of life in Earth's extreme environments. This new knowledge of life on Earth poses the challenge of understandingits molecular basis in such inhospitable conditions, given that such conditions lead to loss of structure and of function in biomolecules from mesophiles. In this review, we discuss the physicochemical properties of extreme environments. We present the state of recent progress in extreme environmental genomics. We then present an overview of our current understanding of the biomolecular adaptation to extreme conditions. As our current and future understanding of biomolecular structure-function relationships in extremophiles requires methodologies adapted to extremes of pressure, temperature, and chemical composition, advances in instrumentation for probing biophysical properties under extreme conditions are presented. Finally, we briefly discuss possible future directions in extreme biophysics.

Adaptations for Pressure and Temperature Effects on Loop Motion in Escherichia coli and Moritella profunda Dihydrofolate Reductase

Determining how enzymes in piezophilic microbes function at high pressure can give insights into how life adapts to living at high pressure. Here, the effects of pressure and temperature on loop motions are compared Escherichia coli (Ec) and Moritella profunda (Mp) dihydrofolate reductase (DHFR) via molecular dynamics simulations at combinations of the growth temperature and pressure of the two organisms. Analysis indicates that a flexible CD loop in MpDHFR is an adaptation for cold because it makes the adenosine binding subdomain more flexible. Also, analysis indicates that the Thr113-Glu27 hydrogen bond in MpDHFR is an adaptation for high pressure because it provides flexibility within the loop subdomain compared to the very strong Thr113-Asp27 hydrogen bond in EcDHFR, and affects the correlation of the Met20 and GH loops. In addition, the results suggest that temperature might affect external loops more strongly while pressure might affect motion between elements within the protein.

Differences in biochemical properties of two 5'-nucleotidases from deep- and shallow-sea Shewanella species under various harsh conditions

Deep-sea Shewanella violacea 5'-nucleotidase (SVNTase) activity exhibited higher NaCl tolerance than that of a shallow-sea Shewanella amazonensis homologue (SANTase), the sequence identity between them being 70.4%. Here, SVNTase exhibited higher activity than SANTase with various inorganic salts, similar to the difference in their NaCl tolerance. In contrast, SVNTase activity decreased with various organic solvents, while SANTase activity was retained with the same concentrations of the solvents. Therefore, SVNTase is more robust than SANTase with inorganic salts, but more vulnerable with organic solvents. As to protein stability, SANTase was more stable against organic solvents and heat than SVNTase, which correlated with the differences in their enzymatic activities. We also found that SANTase retained higher activity for three weeks than SVNTase did in the presence of glycerol. These findings will facilitate further application of these enzymes as appropriate biological catalysts under various harsh conditions. Abbreviations: NTase: 5'-nucleotidase; SANTase: Shewanella amazonensis 5'-nucleotidase; SVNTase: Shewanella violacea 5'-nucleotidase; CD: circular dichroism.

Enzymes from piezophiles

The discovery of microbial communities in extreme conditions that would seem hostile to life leads to the question of how the molecules making up these microbes can maintain their structure and function. While microbes that live under extremes of temperature have been heavily studied, those that live under extremes of pressure, or "piezophiles", are now increasingly being studied because of advances in sample collection and high-pressure cells for biochemical and biophysical measurements. Here, adaptations of enzymes in piezophiles against the effects of pressure are discussed in light of recent experimental and computational studies. However, while concepts from studies of enzymes from temperature extremophiles can provide frameworks for understanding adaptations by piezophile enzymes, the effects of temperature and pressure on proteins differ in significant ways. Thus, the state of the knowledge of adaptation in piezophile enzymes is still in its infancy and many more experiments and computational studies on different enzymes from a variety of piezophiles are needed.

Similar structural stabilities of 3-isopropylmalate dehydrogenases from the obligatory piezophilic bacterium Shewanella benthica strain DB21MT-2 and its atmospheric congener S. oneidensis strain MR-1

We previously found that the enzymatic activity of 3-isopropylmalate dehydrogenase from the obligatory piezophilic bacterium Shewanella benthica strain DB21MT-2 (SbIPMDH) was pressure-tolerant up to 100 MPa, but that from its atmospheric congener S. oneidensis strain MR-1 (SoIPMDH) was pressure-sensitive. Such characteristics were determined by only one amino acid residue at position 266, serine (SoIPMDH) or alanine (SbIPMDH) [Y. Hamajima et al. Extremophiles 20: 177, 2016]. In this study, we investigated the structural stability of these enzymes. At pH 7.6, SoIPMDH was slightly more stable against hydrostatic pressure than SbIPMDH, contrary to the physiological pressures of their normal environments. Pressure unfolding of these IPMDHs followed a two-state unfolding model between a native dimer and two unfolded monomers, and the dimer structure was pressure-tolerant up to 200 MPa, employing a midpoint pressure of 245.3 ± 0.1 MPa and a volume change of -225 ± 24 mL mol^-1 for the most unstable mutant, SbIPMDH A266S. Thus, their pressure-dependent activity did not originate from structural perturbations such as unfolding or dimer dissociation. Conversely, urea-induced unfolding of these IPMDHs followed a three-state unfolding model, including a dimer intermediate. Interestingly, the first transition was strongly pH-dependent but pressure-independent; however, the second transition showed the opposite pattern. Obtained volume changes due to urea-induced unfolding were almost equal for both IPMDHs, approximately +10 and -30 mL mol^-1 for intermediate formation and dimer dissociation, respectively. These results indicated that both IPMDHs have similar structural stability, and a pressure-adaptation mechanism was provided for only the enzymatic activity of SbIPMDH.

Commonly stabilized cytochromes c from deep-sea Shewanella and Pseudomonas

Two cytochromes c5 (SBcytc and SVcytc) have been derived from Shewanella living in the deep-sea, which is a high pressure environment, so it could be that these proteins are more stable at high pressure than at atmospheric pressure, 0.1 MPa. This study, however, revealed that SBcytc and SVcytc were more stable at 0.1 MPa than at higher pressure. In addition, at 0.1-150 MPa, the stability of SBcytc and SVcytc was higher than that of homologues from atmospheric-pressure Shewanella, which was due to hydrogen bond formation with the heme in the former two proteins. This study further revealed that cytochrome c551 (PMcytc) of deep-sea Pseudomonas was more stable than a homologue of atmospheric-pressure Pseudomonas aeruginosa, and that specific hydrogen bond formation with the heme also occurred in the former. Although SBcytc and SVcytc, and PMcytc are phylogenetically very distant, these deep-sea cytochromes c are commonly stabilized through hydrogen bond formation.

Structure-based analysis of Bacilli and plasmid dihydrofolate reductase evolution

Dihydrofolate reductase (DHFR), a key enzyme in tetrahydrofolate-mediated biosynthetic pathways, has a structural motif known to be highly conserved over a wide range of organisms. Given its critical role in purine and amino acid synthesis, DHFR is a well established therapeutic target for treating a wide range of prokaryotic and eukaryotic infections as well as certain types of cancer. Here we present a structural-based computer analysis of bacterial (Bacilli) and plasmid DHFR evolution. We generated a structure-based sequence alignment using 7 wild-type DHFR x-ray crystal structures obtained from the RCSB Protein Data Bank and 350 chromosomal and plasmid homology models we generated from sequences obtained from the NCBI Protein Database. We used these alignments to compare active site and non-active site conservation in terms of amino acid residues, secondary structure and amino acid residue class. With respect to amino acid sequences and residue classes, active-site positions in both plasmid and chromosomal DHFR are significantly more conserved than non-active site positions. Secondary structure conservation was similar for active site and non-active site positions. Plasmid-encoded DHFR proteins have greater degree of sequence and residue class conservation, particularly in sequence positions associated with a network of concerted protein motions, than chromosomal-encoded DHFR proteins. These structure-based were used to build DHFR specific phylogenetic trees from which evidence for horizontal gene transfer was identified.

Comparative study on stabilization mechanism of monomeric cytochrome c5 from deep-sea piezophilic Shewanella violacea

Monomeric cytochrome c5 from deep-sea piezophilic Shewanella violacea (SVcytc5) was stable against heat and denaturant compared with the homologous protein from shallow-sea piezo-sensitive Shewanella livingstonensis (SLcytc5). Here, the SVcytc5 crystal structure revealed that the Lys-50 side chain on the flexible loop formed a hydrogen bond with heme whereas that of corresponding hydrophobic Leu-50 could not form such a bond in SLcytc5, which appeared to be one of possible factors responsible for the difference in stability between the two proteins. This structural insight was confirmed by a reciprocal mutagenesis study on the thermal stability of these two proteins. As SVcytc5 was isolated from a deep-sea piezophilic bacterium, the present comparative study indicates that adaptation of monomeric SVcytc5 to high pressure environments results in stabilization against heat.

Pressure adaptation of 3-isopropylmalate dehydrogenase from an extremely piezophilic bacterium is attributed to a single amino acid substitution

3-Isopropylmalate dehydrogenase (IPMDH) from the extreme piezophile Shewanella benthica (SbIPMDH) is more pressure-tolerant than that from the atmospheric pressure-adapted Shewanella oneidensis (SoIPMDH). To understand the molecular mechanisms of this pressure tolerance, we analyzed mutated enzymes. The results indicate that only a single mutation at position 266, corresponding to Ala (SbIPMDH) and Ser (SoIPMDH), essentially affects activity under higher-pressure conditions. Structural analyses of SoIPMDH suggests that penetration of three water molecules into the cleft around Ser266 under high-pressure conditions could reduce the activity of the wild-type enzyme; however, no water molecule is observed in the Ala266 mutant.

Effects of salt on the structure, stability, and function of a halophilic dihydrofolate reductase from a hyperhalophilic archaeon, Haloarcula japonica strain TR-1

The effects of salt on the structure, stability, and enzymatic function of a novel dihydrofolate reductase (HjDHFR P1) from a hyperhalophilic archaeon, Haloarcula japonica strain TR-1 living in a Japanese saltern, were studied using ultraviolet absorption, circular dichroism (CD), and fluorescence spectroscopy. HjDHFR P1 had a partial structure at pH 8.0 in the absence of NaCl, and the addition of NaCl (0-500 mM concentration) induced significant structural formation to HjDHFR P1. The addition of NADPH, which is a coenzyme for its catalytic reaction, and lowering the pH from 8 to 6 also induced the same CD change, indicating the formation of the NADPH-binding site in HjDHFR P1. The NaCl dependence of thermal and urea-induced unfolding measurements suggested that protein stability increased depending on NaCl concentration regardless of structural formation, and HjDHFR P1 achieved the same stability as Escherichia coli DHFR at 750 mM NaCl. Halophilic characteristics were also observed for enzymatic function, although its structure had already formed under the conditions that enzymatic activity was measured at due to the presence of NADPH. These results suggest that the halophilic mechanism on structural stability and function was caused by factors other than structural formation, which are suggested to be the contributions of preferential interactions between the protein and salt ions and the specific binding of salt ions.

Environmental Adaptation of Dihydrofolate Reductase from Deep-Sea Bacteria

In order to elucidate the molecular adaptation mechanisms of enzymes to the high hydrostatic pressure of the deep sea, we cloned, purified, and characterized more than ten dihydrofolate reductases (DHFRs) from bacteria living in deep-sea and ambient atmospheric pressure environments. The nucleotide and amino acid sequences of these DHFRs indicate the deep-sea bacteria are adapted to their environments after the differentiation of their genus from ancestors inhabiting atmospheric pressure environments. In particular, the backbone structure of the deep-sea DHFR from Moritella profunda (mpDHFR) almost overlapped with the normal homolog from Escherichia coli (ecDHFR). Thus, those of other DHFRs would also overlap on the basis of their sequence similarities. However, the structural stability of both DHFRs was quite different: compared to ecDHFR, mpDHFR was more thermally stable but less stable against urea and pressure unfolding. The smaller volume changes due to unfolding suggest that the native structure of mpDHFR has a smaller cavity and/or enhanced hydration compared to ecDHFR. High hydrostatic pressure reduced the enzymatic activity of many DHFRs, but three deep-sea DHFRs and the D27E mutant of ecDHFR exhibited pressure-dependent activation. The inverted activation volumes from positive to negative values indicate the modification of their structural dynamics, conversion of the rate-determining step of the enzymatic reaction, and different contributions of the cavity and hydration to the transition-state structure. Since the cavity and hydration depend on amino acid side chains, DHFRs would adapt to the deep-sea environment by regulating the cavity and hydration by substituting their amino acid side chains without altering their backbone structure. The results of this study clearly indicate that the cavity and hydration play important roles in the adaptation of enzymes to the deep-sea environment.

Thermodynamic and functional characteristics of deep-sea enzymes revealed by pressure effects

Hydrostatic pressure analysis is an ideal approach for studying protein dynamics and hydration. The development of full ocean depth submersibles and high pressure biological techniques allows us to investigate enzymes from deep-sea organisms at the molecular level. The aim of this review was to overview the thermodynamic and functional characteristics of deep-sea enzymes as revealed by pressure axis analysis after giving a brief introduction to the thermodynamic principles underlying the effects of pressure on the structural stability and function of enzymes.

Solvent environments significantly affect the enzymatic function of Escherichia coli dihydrofolate reductase: comparison of wild-type protein and active-site mutant D27E

To investigate the contribution of solvent environments to the enzymatic function of Escherichia coli dihydrofolate reductase (DHFR), the salt-, pH-, and pressure-dependence of the enzymatic function of the wild-type protein were compared with those of the active-site mutant D27E in relation to their structure and stability. The salt concentration-dependence of enzymatic activity indicated that inorganic cations bound to and inhibited the activity of wild-type DHFR at neutral pH. The BaCl2 concentration-dependence of the (1)H-(15)N HSQC spectra of the wild-type DHFR-folate binary complex showed that the cation-binding site was located adjacent to the Met20 loop. The insensitivity of the D27E mutant to univalent cations, the decreased optimal pH for its enzymatic activity, and the increased Km and Kd values for its substrate dihydrofolate suggested that the substrate-binding cleft of the mutant was slightly opened to expose the active-site side chain to the solvent. The marginally increased fluorescence intensity and decreased volume change due to unfolding of the mutant also supported this structural change or the modified cavity and hydration. Surprisingly, the enzymatic activity of the mutant increased with pressurization up to 250MPa together with negative activation volumes of -4.0 or -4.8mL/mol, depending on the solvent system, while that of the wild-type was decreased and had positive activation volumes of 6.1 or 7.7mL/mol. These results clearly indicate that the insertion of a single methylene at the active site could substantially change the enzymatic reaction mechanism of DHFR, and solvent environments play important roles in the function of this enzyme.

The respiratory system of the piezophile Photobacterium profundum SS9 grown under various pressures

It is known that the facultative piezophile Shewanella violacea DSS12 alters its respiratory components under the influence of hydrostatic pressure during growth. This can be considered one of the mechanisms of bacterial adaptation to high pressure. In this study, we investigated the respiratory system of another well-studied piezophile, Photobacterium profundum SS9. We analyzed cytochrome contents, the expression of genes encoding respiratory components in P. profundum SS9 grown under various conditions, and the pressure dependency of the terminal oxidase activities. Activity was more tolerant of relatively high pressures, such as 125 MPa when the cells were grown under high pressure as compared with cells grown under atmospheric pressure. Such properties observed are similar to the case of S. violacea. However, the contents of the cytochromes and expression of the respiratory genes were not influenced by growth pressure in P. profundum SS9, inconsistent with the case of S. violacea. We suggest that the mechanism of the piezoadaptation of the respiratory system of P. profundum SS9 differs from that of S. violacea, as described above, and that each strain chooses its own strategy.

Pressure dependence of activity and stability of dihydrofolate reductases of the deep-sea bacterium Moritella profunda and Escherichia coli

To understand the pressure-adaptation mechanism of deep-sea enzymes, we studied the effects of pressure on the enzyme activity and structural stability of dihydrofolate reductase (DHFR) of the deep-sea bacterium Moritella profunda (mpDHFR) in comparison with those of Escherichia coli (ecDHFR). mpDHFR exhibited optimal enzyme activity at 50MPa whereas ecDHFR was monotonically inactivated by pressure, suggesting inherent pressure-adaptation mechanisms in mpDHFR. The secondary structure of apo-mpDHFR was stable up to 80°C, as revealed by circular dichroism spectra. The free energy changes due to pressure and urea unfolding of apo-mpDHFR, determined by fluorescence spectroscopy, were smaller than those of ecDHFR, indicating the unstable structure of mpDHFR against pressure and urea despite the three-dimensional crystal structures of both DHFRs being almost the same. The respective volume changes due to pressure and urea unfolding were -45 and -53ml/mol at 25°C for mpDHFR, which were smaller (less negative) than the corresponding values of -77 and -85ml/mol for ecDHFR. These volume changes can be ascribed to the difference in internal cavity and surface hydration of each DHFR. From these results, we assume that the native structure of mpDHFR is loosely packed and highly hydrated compared with that of ecDHFR in solution.

The role of large-scale motions in catalysis by dihydrofolate reductase

Dihydrofolate reductase has long been used as a model system to study the coupling of protein motions to enzymatic hydride transfer. By studying environmental effects on hydride transfer in dihydrofolate reductase (DHFR) from the cold-adapted bacterium Moritella profunda (MpDHFR) and comparing the flexibility of this enzyme to that of DHFR from Escherichia coli (EcDHFR), we demonstrate that factors that affect large-scale (i.e., long-range, but not necessarily large amplitude) protein motions have no effect on the kinetic isotope effect on hydride transfer or its temperature dependence, although the rates of the catalyzed reaction are affected. Hydrogen/deuterium exchange studies by NMR-spectroscopy show that MpDHFR is a more flexible enzyme than EcDHFR. NMR experiments with EcDHFR in the presence of cosolvents suggest differences in the conformational ensemble of the enzyme. The fact that enzymes from different environmental niches and with different flexibilities display the same behavior of the kinetic isotope effect on hydride transfer strongly suggests that, while protein motions are important to generate the reaction ready conformation, an optimal conformation with the correct electrostatics and geometry for the reaction to occur, they do not influence the nature of the chemical step itself; large-scale motions do not couple directly to hydride transfer proper in DHFR.

Eicosapentaenoic acid plays a role in stabilizing dynamic membrane structure in the deep-sea piezophile Shewanella violacea: a study employing high-pressure time-resolved fluorescence anisotropy measurement

Shewanella violacea DSS12 is a psychrophilic piezophile that optimally grows at 30MPa. It contains a substantial amount of eicosapentaenoic acid (EPA) in the membrane. Despite evidence linking increased fatty acid unsaturation and bacterial growth under high pressure, little is known of how the physicochemical properties of the membrane are modulated by unsaturated fatty acids in vivo. By means of the newly developed system performing time-resolved fluorescence anisotropy measurement under high pressure (HP-TRFAM), we demonstrate that the membrane of S. violacea is highly ordered at 0.1MPa and 10°C with the order parameter S of 0.9, and the rotational diffusion coefficient D(w) of 5.4μs(-1) for 1-[4-(trimethylamino)pheny]-6-phenyl-1,3,5-hexatriene in the membrane. Deletion of pfaA encoding the omega-3 polyunsaturated fatty acid synthase caused disorder of the membrane and enhanced the rotational motion of acyl chains, in concert with a 2-fold increase in the palmitoleic acid level. While the wild-type membrane was unperturbed over a wide range of pressures with respect to relatively small effects of pressure on S and D(w), the ΔpfaA membrane was disturbed judging from the degree of increased S and decreased D(w). These results suggest that EPA prevents the membrane from becoming hyperfluid and maintains membrane stability against significant changes in pressure. Our results counter the generally accepted concept that greater fluidity is a membrane characteristic of microorganisms that inhabit cold, high-pressure environments. We suggest that retaining a certain level of membrane physical properties under high pressure is more important than conferring membrane fluidity alone.

Thermal stability of cytochrome c₅ of pressure-sensitive Shewanella livingstonensis

Cytochrome c₅ of pressure-sensitive Shewanella livingstonensis (SL cytc₅) exhibits lower thermal stability than a highly homologous counterpart of pressure-tolerant Shewanella violacea. This stability difference is due to an enthalpic effect that can be attributed to the amino acid residue at position 50 (Leu or Lys). These cytc₅ proteins are appropriate materials for understanding the protein stability mechanism.