Stabilization of Pseudomonas aeruginosa cytochrome c(551) by systematic amino acid substitutions based on the structure of thermophilic Hydrogenobacter thermophilus cytochrome c(552)

Heterologous expression and biochemical comparison of two homologous SoxX proteins of endosymbiontic Candidatus Vesicomyosocius okutanii and free-living Hydrogenovibrio crunogenus from deep-sea vent environments

Candidatus Vesicomyosocius okutanii is a currently uncultured endosymbiotic bacterium of the clam Pheragena okutanii, which lives in deep-sea vent environments. The genome of Ca. V. okutanii encodes a sulfur-oxidizing (Sox) enzyme complex, presumably generating biological energy for the host from inorganic sulfur compounds. Here, Ca. V. okutanii SoxX (VoSoxX), a mono-heme cytochrome c component of the Sox complex, was shown to be phylogenetically related to its homologous counterpart (HcSoxX) from a free-living deep-sea vent bacterium, Hydrogenovibrio crunogenus. Both proteins were heterologously expressed in Escherichia coli cells with co-expressing cytochrome c maturation genes. Biochemical analysis using the recombinant proteins showed that VoSoxX had a significantly lower thermal stability than HcSoxX, possibly due to structural differences. For example, the Asn-60 residue in VoSoxX may be hydrophobically disadvantageous compared with the spatially corresponding Val-73 residue in HcSoxX. This study represents the first successful case of heterologous expression of genes from Ca. V. okutanii, suggesting that the endosymbiotic VoSoxX protein does not require stabilization, unlike the free-living HcSoxX protein.

Reviewing Challenges of Predicting Protein Melting Temperature Change Upon Mutation Through the Full Analysis of a Highly Detailed Dataset with High-Resolution Structures

Predicting the effects of mutations on protein stability is a key problem in fundamental and applied biology, still unsolved even for the relatively simple case of small, soluble, globular, monomeric, two-state-folder proteins. Many articles discuss the limitations of prediction methods and of the datasets used to train them, which result in low reliability for actual applications despite globally capturing trends. Here, we review these and other issues by analyzing one of the most detailed, carefully curated datasets of melting temperature change (ΔTm) upon mutation for proteins with high-resolution structures. After examining the composition of this dataset to discuss imbalances and biases, we inspect several of its entries assisted by an online app for data navigation and structure display and aided by a neural network that predicts ΔTm with accuracy close to that of programs available to this end. We pose that the ΔTm predictions of our network, and also likely those of other programs, account only for a baseline-like general effect of each type of amino acid substitution which then requires substantial corrections to reproduce the actual stability changes. The corrections are very different for each specific case and arise from fine structural details which are not well represented in the dataset and which, despite appearing reasonable upon visual inspection of the structures, are hard to encode and parametrize. Based on these observations, additional analyses, and a review of recent literature, we propose recommendations for developers of stability prediction methods and for efforts aimed at improving the datasets used for training. We leave our interactive interface for analysis available online at http://lucianoabriata.altervista.org/papersdata/proteinstability2021/s1626navigation.html so that users can further explore the dataset and baseline predictions, possibly serving as a tool useful in the context of structural biology and protein biotechnology research and as material for education in protein biophysics.

Thermal destabilization mechanism of cytochrome c' from psychrophilic Shewanella violacea

Cytochrome c' is a nitric oxide (NO)-binding heme protein found in Gram negative bacteria. The thermal stability of psychrophilic Shewanella violacea cytochrome c' (SVCP) is lower than those of its homologues from other 2 psychrophilic Shewanella species, indicating that thermal destabilization mechanism for low-temperature adaptation accumulates in SVCP. In order to understand this mechanism at the amino acid level, here the stability and function of SVCP variants, modeled using the 2 homologues, were examined. The variants exhibited increased stability, and they bound NO similar to the wild type. The vulnerability as to the SVCP stability could be attributed to less hydrogen bond at the subunit interface, more flexible loop structure, and less salt bridge on the protein surface, which appear to be its destabilization mechanism. This study provides an example for controlling stability without spoiling function in psychrophilic proteins.

Protein surface charge effect on 3D domain swapping in cells for c-type cytochromes

Many c-type cytochromes (cyts) can form domain-swapped oligomers. The positively charged Hydrogenobacter thermophilus (HT) cytochrome (cyt) c₅₅₂ forms domain-swapped oligomers during expression in the Escherichia coli (E. coli) expression system, but the factors influencing the oligomerization remain unrevealed. Here, we found that the dimer of the negatively charged Shewanella violacea (SV) cyt c₅ exhibits a domain-swapped structure, in which the N-terminal helix is exchanged between protomers, similar to the structures of the HT cyt c₅₅₂ and Pseudomonas aeruginosa (PA) cyt c₅₅₁ domain-swapped dimers. Positively charged horse cyt c and HT cyt c₅₅₂ domain swapped during expression in E. coli, whereas negatively charged PA cyt c₅₅₁ and SV cyt c₅ did not. Oligomers were formed during expression in E. coli for HT cyt c₅₅₂ attached to either a co- or post-translational signal peptide for transportation through the cytoplasm membrane, but not for PA cyt c₅₅₁ attached to either signal peptide. HT cyt c₅₅₂ formed oligomers in E. coli in the presence and absence of rare codons. More oligomers were obtained from the in vitro folding of horse cyt c and HT cyt c₅₅₂ by the addition of negatively charged liposomes during folding, whereas the amount of oligomers for the in vitro folding of PA cyt c₅₅₁ and SV cyt c₅ did not change significantly by the addition. These results indicate that the protein surface charge affects the oligomerization of c-type cyts in cells; positively charged c-type cyts assemble on a negatively charged membrane, inducing formation of domain-swapped oligomers during folding.

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.

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.

Pseudomonas aeruginosa cytochrome c551 denaturation by five systematic urea derivatives that differ in the alkyl chain length

Reversible denaturation of Pseudomonas aeruginosa cytochrome c551 (PAc551) could be followed using five systematic urea derivatives that differ in the alkyl chain length, i.e. urea, N-methylurea (MU), N-ethylurea (EU), N-propylurea (PU), and N-butylurea (BU). The BU concentration was the lowest required for the PAc551 denaturation, those of PU, EU, MU, and urea being gradually higher. Furthermore, the accessible surface area difference upon PAc551 denaturation caused by BU was found to be the highest, those by PU, EU, MU, and urea being gradually lower. These findings indicate that urea derivatives with longer alkyl chains are stronger denaturants. In this study, as many as five systematic urea derivatives could be applied for the reversible denaturation of a single protein, PAc551, for the first time, and the effects of the alkyl chain length on protein denaturation were systematically verified by means of thermodynamic parameters.

Structural and functional insights into thermally stable cytochrome c' from a thermophile

Thermophilic Hydrogenophilus thermoluteolus cytochrome c′ (PHCP) exhibits higher thermal stability than a mesophilic counterpart, Allochromatium vinosum cytochrome c′ (AVCP), which has a homo‐dimeric structure and ligand‐binding ability. To understand the thermal stability mechanism and ligand‐binding ability of the thermally stable PHCP protein, the crystal structure of PHCP was first determined. It formed a homo‐dimeric structure, the main chain root mean square deviation (rmsd) value between PHCP and AVCP being 0.65 Å. In the PHCP structure, six specific residues appeared to strengthen the heme‐related and subunit–subunit interactions, which were not conserved in the AVCP structure. PHCP variants having altered subunit–subunit interactions were more severely destabilized than ones having altered heme‐related interactions. The PHCP structure further revealed a ligand‐binding channel and a penta‐coordinated heme, as observed in the AVCP protein. A spectroscopic study clearly showed that some ligands were bound to the PHCP protein. It is concluded that the dimeric PHCP from the thermophile is effectively stabilized through heme‐related and subunit–subunit interactions with conservation of the ligand‐binding ability.

Brief Summary

We report the X‐ray crystal structure of cytochrome c′ (PHCP) from thermophilic Hydrogenophilus thermoluteolus. The high thermal stability of PHCP was attributed to heme‐related and subunit–subunit interactions, which were confirmed by a mutagenesis study. The ligand‐binding ability of PHCP was examined by spectrophotometry. PHCP acquired the thermal stability with conservation of the ligand‐binding ability. This study furthers the understanding of the stability and function of cytochromes c.

PDB Code(s): 5B3I

Domain swapping oligomerization of thermostable c-type cytochrome in E. coli cells

Knowledge on domain swapping in vitro is increasing, but domain swapping may not occur regularly in vivo, and its information in cells is limited. Herein, we show that domain-swapped oligomers of a thermostable c-type cytochrome, Hydrogenobacter thermophilus cyt c₅₅₂, are formed in E. coli which expresses cyt c₅₅₂. The region containing the N-terminal α-helix and heme was domain-swapped between protomers in the dimer formed in E. coli. The amount of cyt c₅₅₂ oligomers increased in E. coli as the cyt c₅₅₂ concentration was increased, whereas that of high-order oligomers decreased in the order of decrease in protein stability, indicating that domain swapping decreases in cells when the protein stability decreases. Apo cyt c₅₅₂ was detected in the cyt c₅₅₂ oligomer formed in E. coli, but not in that of the A5F/M11V/Y32F/Y41E/I76V mutant. The cyt c₅₅₂ oligomer containing its apo protein may form at the periplasm, since the apo protein detected by mass measurements did not contain the signal peptide. These results show that domain-swapped cyt c₅₅₂ oligomers were formed in E. coli, owing to the stability of the transient oligomer containing the apo protein before heme attachment. This is an indication that exceedingly stable proteins may have disadvantages forming domain-swapped oligomers in cells.

Comparison of the backbone dynamics of wild-type Hydrogenobacter thermophilus cytochrome c(552) and its b-type variant

Cytochrome c552 from the thermophilic bacterium Hydrogenobacter thermophilus is a typical c-type cytochrome which binds heme covalently via two thioether bonds between the two heme vinyl groups and two cysteine thiol groups in a CXXCH sequence motif. This protein was converted to a b-type cytochrome by substitution of the two cysteine residues by alanines (Tomlinson and Ferguson in Proc Natl Acad Sci USA 97:5156-5160, 2000a). To probe the significance of the covalent attachment of the heme in the c-type protein, (15)N relaxation and hydrogen exchange studies have been performed for the wild-type and b-type proteins. The two variants share very similar backbone dynamic properties, both proteins showing high (15)N order parameters in the four main helices, with reduced values in an exposed loop region (residues 18-21), and at the C-terminal residue Lys80. Some subtle changes in chemical shift and hydrogen exchange protection are seen between the wild-type and b-type variant proteins, not only for residues at and neighbouring the mutation sites, but also for some residues in the heme binding pocket. Overall, the results suggest that the main role of the covalent linkages between the heme group and the protein chain must be to increase the stability of the protein.

Domain-swapped dimer of Pseudomonas aeruginosa cytochrome c551: structural insights into domain swapping of cytochrome c family proteins

Cytochrome c (cyt c) family proteins, such as horse cyt c, Pseudomonas aeruginosa cytochrome c₅₅₁ (PA cyt c₅₅₁), and Hydrogenobacter thermophilus cytochrome c₅₅₂ (HT cyt c₅₅₂), have been used as model proteins to study the relationship between the protein structure and folding process. We have shown in the past that horse cyt c forms oligomers by domain swapping its C-terminal helix, perturbing the Met–heme coordination significantly compared to the monomer. HT cyt c₅₅₂ forms dimers by domain swapping the region containing the N-terminal α-helix and heme, where the heme axial His and Met ligands belong to different protomers. Herein, we show that PA cyt c₅₅₁ also forms domain-swapped dimers by swapping the region containing the N-terminal α-helix and heme. The secondary structures of the M61A mutant of PA cyt c₅₅₁ were perturbed slightly and its oligomer formation ability decreased compared to that of the wild-type protein, showing that the stability of the protein secondary structures is important for domain swapping. The hinge loop of domain swapping for cyt c family proteins corresponded to the unstable region specified by hydrogen exchange NMR measurements for the monomer, although the swapping region differed among proteins. These results show that the unstable loop region has a tendency to become a hinge loop in domain-swapped proteins.

Inversion of the stereochemistry around the sulfur atom of the axial methionine side chain through alteration of amino acid side chain packing in Hydrogenobacter thermophilus cytochrome C552 and its functional consequences

In cytochrome c, the coordination of the axial Met Sδ atom to the heme Fe atom occurs in one of two distinctly different stereochemical manners, i.e., R and S configurations, depending upon which of the two lone pairs of the Sδ atom is involved in the bond; hence, the Fe-coordinated Sδ atom becomes a chiral center. In this study, we demonstrated that an alteration of amino acid side chain packing induced by the mutation of a single amino acid residue, i.e., the A73V mutation, in Hydrogenobacter thermophilus cytochrome c552 (HT) forces the inversion of the stereochemistry around the Sδ atom from the R configuration [Travaglini-Allocatelli, C., et al. (2005) J. Biol. Chem. 280, 25729-25734] to the S configuration. Functional comparison between the wild-type HT and the A73V mutant possessing the R and S configurations as to the stereochemistry around the Sδ atom, respectively, demonstrated that the redox potential (Em) of the mutant at pH 6.00 and 25 °C exhibited a positive shift of ∼20 mV relative to that of the wild-type HT, i.e., 245 mV, in an entropic manner. Because these two proteins have similar enthalpically stabilizing interactions, the difference in the entropic contribution to the Em value between them is likely to be due to the effect of the conformational alteration of the axial Met side chain associated with the inversion of the stereochemistry around the Sδ atom due to the effect of mutation on the internal mobility of the loop bearing the axial Met. Thus, the present study demonstrated that the internal mobility of the loop bearing the axial Met, relevant to entropic control of the redox function of the protein, is affected quite sensitively by the contextual stereochemical packing of amino acid side chains in the proximity of the axial Met.

From protein engineering to immobilization: promising strategies for the upgrade of industrial enzymes

Enzymes found in nature have been exploited in industry due to their inherent catalytic properties in complex chemical processes under mild experimental and environmental conditions. The desired industrial goal is often difficult to achieve using the native form of the enzyme. Recent developments in protein engineering have revolutionized the development of commercially available enzymes into better industrial catalysts. Protein engineering aims at modifying the sequence of a protein, and hence its structure, to create enzymes with improved functional properties such as stability, specific activity, inhibition by reaction products, and selectivity towards non-natural substrates. Soluble enzymes are often immobilized onto solid insoluble supports to be reused in continuous processes and to facilitate the economical recovery of the enzyme after the reaction without any significant loss to its biochemical properties. Immobilization confers considerable stability towards temperature variations and organic solvents. Multipoint and multisubunit covalent attachments of enzymes on appropriately functionalized supports via linkers provide rigidity to the immobilized enzyme structure, ultimately resulting in improved enzyme stability. Protein engineering and immobilization techniques are sequential and compatible approaches for the improvement of enzyme properties. The present review highlights and summarizes various studies that have aimed to improve the biochemical properties of industrially significant enzymes.

An atomistic view on human hemoglobin carbon monoxide migration processes

A significant amount of work has been devoted to obtaining a detailed atomistic knowledge of the human hemoglobin mechanism. Despite this impressive research, to date, the ligand diffusion processes remain unclear and controversial. Using recently developed computational techniques, PELE, we are capable of addressing the ligand migration processes. First, the methodology was tested on myoglobin's CO migration, and the results were compared with the wealth of theoretical and experimental studies. Then, we explored both hemoglobin tense and relaxed states and identified the differences between the α-and β-subunits. Our results indicate that the proximal site, equivalent to the Xe1 cavity in myoglobin, is never visited. Furthermore, strategically positioned residues alter the diffusion processes within hemoglobin's subunits and suggest that multiple pathways exist, especially diversified in the α-globins. A significant dependency of the ligand dynamics on the tertiary structure is also observed.

Enhancement of the thermostability of Hydrogenobacter thermophilus cytochrome c(552) through introduction of an extra methylene group into its hydrophobic protein interior

Careful scrutiny of the protein interior of Hydrogenobacter thermophilus cytochrome c(552) (HT) on the basis of its X-ray structure [Travaglini-Allocatelli, C., Gianni, S., Dubey, V. K., Borgia, A., Di Matteo, A., Bonivento, D., Cutruzzola, F., Bren, K. L., and Brunori, M. (2005) J. Biol. Chem. 280, 25729-25734] indicated that a void space, which is large enough to accommodate a methyl group, exists in the hydrophobic protein interior near the heme. We tried to reduce the void space through the replacement of a Val by Ile or Leu (Val/Ile or Val/Leu mutation), and then the structural and functional consequences of these two mutations were characterized in order to elucidate the relationship between the nature of the packing of hydrophobic residues and the functional properties of the protein. The study demonstrated striking differences in the structural and functional consequences between the two mutations. The Val/Ile mutation was found to cause further enhancement of the thermostability of the oxidized HT, as reflected in the increase of the denaturation temperature (T(m)) value by ∼ 3 deg, whereas the thermostability of the reduced form was essentially unaffected. As a result, the redox potential (E(m)) of the Val/Ile mutant exhibited a negative shift of ∼ 50 mV relative to that of the wild-type protein in an enthalpic manner, this being consistent with our previous finding that a protein with higher stability in its oxidized form exhibits a lower E(m) value [Terui, N., Tachiiri, N., Matsuo, H., Hasegawa, J., Uchiyama, S., Kobayashi, Y., Igarashi, Y., Sambongi, Y., and Yamamoto, Y. (2003) J. Am. Chem. Soc. 125, 13650-13651]. In contrast, the Val/Leu mutation led to a decrease in thermostability of both the redox forms of the protein, as reflected in the decreases of the T(m) values of the oxidized and reduced proteins by ∼ 3 and ∼ 5 deg, respectively, and the E(m) value of the Val/Leu mutant happened to be similar to that of the Val/Ile one. The E(m) value of the Val/Leu mutant could be reasonably interpreted in terms of the different effects of the mutation on the stabilities of the two different redox forms of the protein. Thus, the present study demonstrated that the stability of the protein is affected quite sensitively by the contextual stereochemical packing of hydrophobic residues in the protein interior and that the structural properties of the hydrophobic core in the protein interior are crucial for control of the redox function of the protein. These findings provide novel insights as to functional control of a protein, which could be utilized for tuning of the T(m) and E(m) values of the protein by means of protein engineering.

Stability of the heme Fe-N-terminal amino group coordination bond in denatured cytochrome c

In the denatured states of Hydrogenobacter thermophilus cytochrome c(552) (HT) and Pseudomonas aeruginosa cytochrome c(551) (PA), and their mutants, the N-terminal amino group of the polypeptide chain is coordinated to heme Fe in place of the axial Met, the His-N(term) form being formed. The coordination of the N-terminal amino group to heme Fe leads to loop formation by the N-terminal stretch preceding the first Cys residue bound to the heme, and the N-terminal stretches of HT and PA are different from each other in terms of both the sequence and the number of constituent amino acid residues. The His-N(term) form was shown to be rather stable, and hence it can influence the stability of the denatured state. We have investigated the heme Fe coordination structures and stabilities of the His-N(term) forms emerging upon guanidine hydrochloric acid-induced unfolding of the oxidized forms of the proteins. The Fe-N(term) coordination bond in the His-N(term) form with a 9-residue N-terminal stretch of HT proteins was found to be tilted to some extent away from the heme normal, as reflected by the great heme methyl proton shift spread. On the other hand, the small heme methyl proton shift spread of the His-N(term) form with an 11-residue stretch of PA proteins indicated that its Fe-N(term) bond is nearly parallel with the heme normal. The stability of the His-N(term) form was found to be affected by the structural properties of the N-terminal stretch, such as its length and the N-terminal residue. With a given N-terminal residue, the stability of the His-N(term) form is higher for a 9-residue N-terminal stretch than an 11-residue one. In addition, with a given length of the N-terminal stretch, the His-N(term) form with an N-terminal Glu is stabilized by a few kJ mol(-1) relative to that with an N-terminal Asn. These results provide a novel insight into the stabilizing interactions in the denatured cyts c that will facilitate elucidation of the folding/unfolding mechanisms of the proteins.

Isotopic labeling of c-type multiheme cytochromes overexpressed in E. coli

Progresses made in bacterial genome sequencing show a remarkable profusion of multiheme c-type cytochromes in many bacteria, highlighting the importance of these proteins in different cellular events. However, the characterization of multiheme cytochromes has been significantly retarded by the numerous experimental challenges encountered by researchers who attempt to overexpress these proteins, especially if isotopic labeling is required. Here we describe a methodology for isotopic labeling of multiheme cytochromes c overexpressed in Escherichia coli, using the triheme cytochrome PpcA from Geobacter sulfurreducens as a model protein. By combining different strategies previously described and using E. coli cells containing the gene coding for PpcA and the cytochrome c maturation gene cluster, an experimental labeling methodology was developed that is based on two major aspects: (i) use of a two-step culture growth procedure, where cell growth in rich media was followed by transfer to minimal media containing (15)N-labeled ammonium chloride, and (ii) incorporation of the heme precursor delta-aminolevulinic acid in minimal culture media. The yields of labeled protein obtained were comparable to those obtained for expression of PpcA in rich media. Proper protein folding and labeling were confirmed by UV-visible and NMR spectroscopy. To our knowledge, this is the first report of a recombinant multiheme cytochrome labeling and it represents a major breakthrough for functional and structural studies of multiheme cytochromes.

Roles of a short connecting disulfide bond in the stability and function of psychrophilic Shewanella violacea cytochrome c (5)

Cys-59 and Cys-62, forming a disulfide bond in the four-residue loop of Shewanella violacea cytochrome c (5) (SV cytc (5)), contribute to protein stability but not to redox function. These Cys residues were substituted with Ala in SV cytc (5), and the structural and functional properties of the resulting C59A/C62A variant were determined and compared with those of the wild-type. The variant had similar features to those of the wild-type in absorption, circular dichroic, and paramagnetic (1)H NMR spectra. In addition, the redox potentials of the wild-type and variant were essentially the same, indicating that removal of the disulfide bond from SV cytc (5) does not affect the redox function generated in the vicinity of heme. However, calorimetric analysis of the wild-type and variant showed that the mutations caused a drastic decrease in the protein stability through enthalpy, but not entropy. Four residues are encompassed by the SV cytc (5) disulfide bond, which is the shortest one that has been proved to affect protein stability. The protein stability of SV cytc (5) can be controlled without changing the redox function, providing a new strategy for regulating the stability and function of cytochrome c.

Residue cluster additivity of thermodynamic stability in the hydrophobic core of mesophile vs. hyperthermophile rubredoxins

The branched sidechain residues 24 and 33 in the hydrophobic core of rubredoxin differ between the Clostridium pasteurianum (Cp) and Pyrococcus furiosus (Pf) sequences. Their X-ray structures indicate that these two sidechains are in van der Waals contact with each other, while neither appears to significantly interact with the other nonconserved residues. The simultaneous interchange of residues 24 and 33 between the Cp and Pf rubredoxin sequences yield a complementary pair of hybrid proteins for which the sum of their thermodynamic stabilities equals that of the parental rubredoxins. The 1.2 kcal/mol change arising from this two residues interchange accounts for 21% of the differential thermodynamic stability between the mesophile and hyperthermophile proteins. The additional interchange of the sole nonconserved aromatic residue in the hydrophobic core yields a 0.78 kcal/mol deviation from thermodynamic additivity.

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.

Relationship between local structural entropy and protein thermostability

We developed a technique to compute structural entropy directly from protein sequences. We explored the possibility of using structural entropy to identify residues involved in thermal stabilization of various protein families. Examples include methanococcal adenylate kinase, Ribonuclease HI and holocytochrome c(551). Our results show that the positions of the largest structural entropy differences between wild type and mutant usually coincide with the residues relevant to thermostability. We also observed a good linear relationship between the average structural entropy and the melting temperatures for adenylate kinase and its chimeric constructs. To validate this linear relationship, we compiled a large dataset comprised of 1153 sequences and found that most protein families still display similar linear relationships. Our results suggest that the multitude of interactions involved in thermal stabilization may be generalized into the tendency of proteins to maintain local structural conservation. The linear relationship between structural entropy and protein thermostability should be useful in the study of protein thermal stabilization.

Unveiling a hidden folding intermediate in c-type cytochromes by protein engineering

Several investigators have highlighted a correlation between the basic features of the folding process of a protein and its topology, which dictates the folding pathway. Within this conceptual framework we proposed that different members of the cytochrome c (cyt c) family share the same folding mechanism, involving a consensus partially structured state. Pseudomonas aeruginosa cyt c(551) (Pa cyt c(551)) folds via an apparent two-state mechanism through a high energy intermediate. Here we present kinetic evidence demonstrating that it is possible to switch its folding mechanism from two to three state, stabilizing the high energy intermediate by rational mutagenesis. Characterization of the folding kinetics of one single-site mutant of the Pa cyt c(551) (Phe(7) to Ala) indeed reveals an additional refolding phase and a fast unfolding process which are explained by the accumulation of a partially folded species. Further kinetic analysis highlights the presence of two parallel processes both leading to the native state, suggesting that the above mentioned species is a non obligatory on-pathway intermediate. Determination of the crystallographic structure of F7A shows the presence of an extended internal cavity, which hosts three "bound" water molecules and a H-bond in the N-terminal helix, which is shorter than in the wild type protein. These two features allow us to propose a detailed structural interpretation for the stabilization of the native and especially the intermediate states induced by a single crucial mutation. These results show how protein engineering, x-ray crystallography and state-of-the-art kinetics concur to unveil a folding intermediate and the structural determinants of its stability.

Unexpected elevated production of Aquifex aeolicus cytochrome c555 in Escherichia coli cells lacking disulfide oxidoreductases

Mutant strains of Escherichia coli lacking DsbA, DsbB, or DsbD (proteins required for disulfide bond formation in the periplasm) did not produce mitochondrial or chloroplast cytochromes c, as previously observed for bacterial ones. Unexpectedly, however, cytochrome c(555) (AA c(555)) from a hyperthermophile, Aquifex aeolicus, was produced in the E. coli periplasm without Dsb proteins, three times more than with them. These results indicate that the Dsb proteins are not necessarily required for AA c(555) production in E. coli, possibly because of hyperthermophilic origin compared with the others.

Additivity in Both Thermodynamic Stability and Thermal Transition Temperature for Rubredoxin Chimeras via Hybrid Native Partitioning

Given any operational definition of pairwise interaction, the set of residues that differ between two structurally homologous proteins can be uniquely partitioned into subsets of clusters for which no such interactions occur between clusters. Although hybrid protein sequences that preserve such clustering are consistent with tertiary structures composed of only parental native-like interactions, the stability of such predicted structures will depend upon the physical robustness of the assumed interaction potential. A simple distance cutoff criterion was applied to the most thermostable protein known to predict such a seven-residue cluster in the metal binding site region of Pyrococcus furiosus rubredoxin and a mesophile homolog. Both conformational stability and thermal transition temperature measurements demonstrate that 39% of the differential stability arises from these seven residues.

Cloning, expression, crystallization and preliminary X-ray characterization of cytochrome c552 from a moderate thermophilic bacterium, Hydrogenophilus thermoluteolus

The amino-acid sequence of cytochrome c552 (PH c552) from a moderately thermophilic bacterium, Hydrogenophilus thermoluteolus, was more than 50% identical to that of cytochrome c from an extreme thermophile, Hydrogenobacter thermophilus (HT c552), and from a mesophile, Pseudomonas aeruginosa (PA c551). The PH c552 gene was overexpressed as a correctly processed holoprotein in the Escherichia coli periplasm. The overexpressed PH c552 has been crystallized by vapour diffusion from polyethylene glycol 4000 pH 6.5. The crystals belong to space group C222(1), with unit-cell parameters a = 48.98, b = 57.99, c = 56.20 A. The crystals diffract X-rays to around 2.1 A resolution.

Amino acid coupling patterns in thermophilic proteins

Structural analysis is useful in elucidating structural features responsible for enhanced thermal stability of proteins. However, due to the rapid increase of sequenced genomic data, there are far more protein sequences than the corresponding three-dimensional (3D) structures. The usual sequence-based amino acid composition analysis provides useful but simplified clues about the amino acid types related to thermal stability of proteins. In this work, we developed a statistical approach to identify the significant amino acid coupling sequence patterns in thermophilic proteins. The amino acid coupling sequence pattern is defined as any 2 types of amino acids separated by 1 or more amino acids. Using this approach, we construct the rho profiles for the coupling patterns. The rho value gives a measure of the relative occurrence of a coupling pattern in thermophiles compared with mesophiles. We found that thermophiles and mesophiles exhibit significant bias in their amino acid coupling patterns. We showed that such bias is mainly due to temperature adaptation instead of species or GC content variations. Though no single outstanding coupling pattern can adequately account for protein thermostability, we can use a group of amino acid coupling patterns having strong statistical significance (p values < 10(-7)) to distinguish between thermophilic and mesophilic proteins. We found a good correlation between the optimal growth temperatures of the genomes and the occurrences of the coupling patterns (the correlation coefficient is 0.89). Furthermore, we can separate the thermophilic proteins from their mesophilic orthologs using the amino acid coupling patterns. These results may be useful in the study of the enhanced stability of proteins from thermophiles-especially when structural information is scarce. Proteins 2005. (c) 2005 Wiley-Liss, Inc.

Five amino acid residues responsible for the high stability of Hydrogenobacter thermophilus cytochrome c552: reciprocal mutation analysis

Five amino acid residues responsible for extreme stability have been identified in cytochrome c(552) (HT c(552)) from a thermophilic bacterium, Hydrogenobacter thermophilus. The five residues, which are spatially distributed in three regions of HT c(552), were replaced with the corresponding residues in the homologous but less stable cytochrome c(551) (PA c(551)) from Pseudomonas aeruginosa. The quintuple HT c(552) variant (A7F/M13V/Y34F/Y43E/I78V) showed the same stability against guanidine hydrochloride denaturation as that of PA c(551), suggesting that the five residues in HT c(552) necessarily and sufficiently contribute to the overall stability. In the three HT c(552) variants carrying mutations in each of the three regions, the Y34F/Y43E mutations resulted in the greatest destabilization, by -13.3 kJ mol(-1), followed by A7F/M13V (-3.3 kJ mol(-1)) and then I78V (-1.5 kJ mol(-1)). The order of destabilization in HT c(552) was the same as that of stabilization in PA c(551) with reverse mutations such as F34Y/E43Y, F7A/V13M, and V78I (13.4, 10.3, and 0.3 kJ mol(-1), respectively). The results of guanidine hydrochloride denaturation were consistent with those of thermal denaturation for the same variants. The present study established a method for reciprocal mutation analysis. The effects of side-chain contacts were experimentally evaluated by swapping the residues between the two homologous proteins that differ in stability. A comparative study of the two proteins was a useful tool for assessing the amino acid contribution to the overall stability.

Rational engineering of enzyme stability

During the past 15 years there has been a continuous flow of reports describing proteins stabilized by the introduction of mutations. These reports span a period from pioneering rational design work on small enzymes such as T4 lysozyme and barnase to protein design, and directed evolution. Concomitantly, the purification and characterization of naturally occurring hyperstable proteins has added to our understanding of protein stability. Along the way, many strategies for rational protein stabilization have been proposed, some of which (e.g. entropic stabilization by introduction of prolines or disulfide bridges) have reasonable success rates. On the other hand, comparative studies and efforts in directed evolution have revealed that there are many mutational strategies that lead to high stability, some of which are not easy to define and rationalize. Recent developments in the field include increasing awareness of the importance of the protein surface for stability, as well as the notion that normally a very limited number of mutations can yield a large increase in stability. Another development concerns the notion that there is a fundamental difference between the "laboratory stability" of small pure proteins that unfold reversibly and completely at high temperatures and "industrial stability", which is usually governed by partial unfolding processes followed by some kind of irreversible inactivation process (e.g. aggregation). Provided that one has sufficient knowledge of the mechanism of thermal inactivation, successful and efficient rational stabilization of enzymes can be achieved.

Do Proteins Always Benefit from a Stability Increase? Relevant and Residual Stabilisation in a Three-state Protein by Charge Optimisation

The vast majority of our knowledge on protein stability arises from the study of simple two-state models. However, proteins displaying equilibrium intermediates under certain conditions abound and it is unclear whether the energetics of native/intermediate equilibria is well represented in current knowledge. We consider here that the overall conformational stability of three-state proteins is made of a "relevant" term and a "residual" one, corresponding to the free energy differences of the native to intermediate (N-to-I) and intermediate to denatured (I-to-D) equilibria, respectively. The N-to-I free energy difference is considered to be the relevant stability because protein-unfolding intermediates are likely devoid of biological activity. We use surface charge optimisation to first increase the overall (N-to-D) stability of a model three-state protein (apoflavodoxin) and then investigate whether the stabilisation obtained is realised into relevant or into residual stability. Most of the mutations designed from electrostatic calculations or from simple sequence conservation analysis produce large increases in the overall stability of the protein. However, in most cases, this simply leads to similarly large increases of the residual stability. Two mutations, nevertheless, show a different trend and increase the relevant stability of the protein substantially. When all the mutations are mapped onto the structure of the apoflavodoxin thermal-unfolding intermediate (obtained independently by equilibrium phi-analysis and NMR) they cluster perfectly so that the mutations increasing the relevant stability appear in the small unstructured region of the intermediate and the others in the native-like region. This illustrates the need for specific investigation of N-to-I equilibria and the structure of protein intermediates, and indicates that it is possible to rationally stabilise a protein against partial unfolding once the structure of the intermediate conformation is known, even if at low resolution.

Large improvement in the thermal stability of a tetrameric malate dehydrogenase by single point mutations at the dimer-dimer interface

The stability of tetrameric malate dehydrogenase from the green phototrophic bacterium Chloroflexus aurantiacus (CaMDH) is at least in part determined by electrostatic interactions at the dimer-dimer interface. Since previous studies had indicated that the thermal stability of CaMDH becomes lower with increasing pH, attempts were made to increase the stability by removal of (excess) negative charge at the dimer-dimer interface. Mutation of Glu165 to Gln or Lys yielded a dramatic increase in thermal stability at pH 7.5 (+23.6 -- + 23.9 degrees C increase in apparent t(m)) and a more moderate increase at pH 4.4 (+4.6 -- + 5.4 degrees C). The drastically increased stability at neutral pH was achieved without forfeiture of catalytic performance at low temperatures. The crystal structures of the two mutants showed only minor structural changes close to the mutated residues, and indicated that the observed stability effects are solely due to subtle changes in the complex network of electrostatic interactions in the dimer-dimer interface. Both mutations reduced the concentration dependency of thermal stability, suggesting that the oligomeric structure had been reinforced. Interestingly, the two mutations had similar effects on stability, despite the charge difference between the introduced side-chains. Together with the loss of concentration dependency, this may indicate that both E165Q and E165K stabilize CaMDH to such an extent that disruption of the inter-dimer electrostatic network around residue 165 no longer limits kinetic thermal stability.

Effects of axial methionine coordination on the in-plane asymmetry of the heme electronic structure of cytochrome c

The paramagnetic susceptibility ( chi) tensors of the oxidized forms of thermophile Hydrogenobacter thermophilus cytochrome c(552) (Ht cyt c(552)) and a quintuple mutant (F7A/V13 M/F34Y/E43Y/V78I; qm) of mesophile Pseudomonas aeruginosa cytochrome c(551) (Pa cyt c(551)) have been determined on the basis of the redox-dependent (1)H NMR shift changes of the main-chain NH and C(alpha)H proton resonances of non-coordinated amino acid residues and the NMR structures of the reduced forms of the corresponding proteins (J. Hasegawa, T. Yoshida, T. Yamazaki, Y. Sambongi, Y. Yu, Y. Igarashi, T. Kodama, K. Yamazaki, Y. Kyogoku, Y. Kobayashi (1998) Biochemistry 37:9641-9649; J. Hasegawa, S. Uchiyama, Y. Tanimoto, M. Mizutani, Y. Kobayashi, Y. Sambongi,Y. Igarashi (2000) J Biol Chem 275:37824-37828). From the chi tensors determined, we obtained the contact shifts for heme methyl proton resonances, which provided the heme electronic structures of the oxidized forms of Ht cyt c(552) and qm. We also characterized the heme electronic structure of the cyanide adducts of the proteins, where the axial Met was replaced by an exogenous cyanide ion, through the analysis of (1)H NMR spectra. The results indicated that the heme electronic structures of both the proteins in their oxidized forms with axial His and Met coordination are largely different to each other, while those in their cyanide adducts are similar to each other. These results demonstrated that the orientation of the axial Met sulfur lone pair, with respect to heme, predominantly contributes to the spin delocalization into the porphyrin-pi system of heme in the oxidized proteins with axial His and Met coordination.

Modulation of mammalian cell growth and death by prokaryotic and eukaryotic cytochrome c

Cytochrome c(551), an 8,685-Da haem-containing protein, is known to be involved in electron transfer during dissimilative denitrification by Pseudomonas aeruginosa. Both cytochrome c(551) and copper-containing redox protein azurin have been used in vitro as partners in electron transfer. Azurin has been reported to induce apoptosis in macrophages and cancer cells. We now report that, unlike azurin, cytochrome c(551), termed Cyt c(551), has very little ability to induce apoptosis in J774 cell line-derived macrophages but demonstrates significant inhibition of cell cycle progression in such cells. A mutant form of Cyt c(551), V23DI59E, has significantly reduced ability to inhibit cell cycle progression but demonstrates a higher level of apoptosis-inducing activity in macrophages, compared with WT Cyt c(551). Interestingly, the WT Cyt c(551), but not the mutant form, significantly enhances the level of tumor suppressor protein p16(Ink4a), a known inhibitor of cell cycle progression whereas the mutant form seems to form a complex with tumor suppressor protein p53, thereby enhancing its intracellular level to some extent. Eukaryotic cytochromes such as horse and bovine cytochrome c have also been shown to induce apoptosis but not inhibition of cell cycle progression in J774 cells, thus signifying a role of this redox protein in entry to, and in the induction of, cell death in mammalian cells.

Induction of apoptosis in macrophages by Pseudomonas aeruginosa azurin: tumour-suppressor protein p53 and reactive oxygen species, but not redox activity, as critical elements in cytotoxicity

Azurin is a copper-containing protein involved in electron transfer during denitrification. We reported recently that purified azurin demonstrates cytotoxicity to macrophages by forming a complex with the tumour-suppressor protein p53, thereby stabilizing it and enhancing its function as an inducer of proapoptotic activity (Yamada, T., Goto, M., Punj, V., Zaborina, O., Kimbara, K., Das Gupta, T. K., and Chakrabarty, A. M. 2002, Infect Immun70: 7054-7062). It is, however, not known whether the oxidoreductase (redox) activity of azurin or the involvement of copper is important for its cytotoxicity. We have isolated apo-azurin devoid of copper and site-directed mutants that are redox negative because of either replacement of a cysteine residue (Cys-112) involved in co-ordination with copper or mutational replacement of two methionine residues (Met-44 and Met-64) that are present in the hydrophobic patch of azurin and allow interaction of azurin with its redox partner cytochrome c551. We demonstrate that, although the wild type (wt) and the Cys-112 Asp mutant azurin can form complexes with the tumour-suppressor protein p53 and generate high levels of reactive oxygen species (ROS), the redox-negative Met-44LysMet-64Glu mutant azurin is defective in complex formation with p53, generates low levels of ROS and lacks appreciable cytotoxicity towards macrophages. Thus, complex formation with p53 and ROS generation, rather than azurin redox activity, are important in the cytotoxic action of azurin towards macrophages.

The bacterial redox protein azurin induces apoptosis in J774 macrophages through complex formation and stabilization of the tumor suppressor protein p53

Two redox proteins, azurin and cytochrome c(551) elaborated by Pseudomonas aeruginosa, demonstrate significant cytotoxic activity towards macrophages. Azurin can enter macrophages, localize in the cytosol and nuclear fractions, and induce apoptosis. Two redox-negative mutants of azurin have less cytotoxicity than does wild-type (wt) azurin. Azurin has been shown to form a complex with the tumor suppressor protein p53, a known inducer of apoptosis, thereby stabilizing it and enhancing its intracellular level. A higher level of reactive oxygen species (ROS), generated during treatment of macrophages with wt azurin, correlates with its cytotoxicity. Treatment with some ROS-removing antioxidants greatly reduces azurin-mediated cytotoxicity, thus demonstrating a novel virulence property of this bacterial redox protein.

Cytochrome c from a thermophilic bacterium has provided insights into the mechanisms of protein maturation, folding, and stability

Cytochrome c is widely distributed in bacterial species, from mesophiles to thermophiles, and is one of the best-characterized redox proteins in terms of biogenesis, folding, structure, function, and evolution. Experimental molecular biology techniques (gene cloning and expression) have become applicable to cytochrome c, enabling its engineering and manipulation. Heterologous expression systems for cytochromes c in bacteria, for use in mutagenesis studies, have been established by extensive investigation of the biological process by which the functional structure is formed. Mutagenesis and structure analyses based on comparative studies using a thermophile Hydrogenobacter thermophilus cytochrome c-552 and its mesophilic counterpart have provided substantial clues to the mechanism underlying protein stability at the amino-acid level. The molecular mechanisms underlying protein maturation, folding, and stability in bacterial cytochromes c are beginning to be understood.