An obligatory intermediate in the folding pathway of cytochrome c552 from Hydrogenobacter thermophilus

Enantioselective Synthesis of α-Trifluoromethyl Amines via Biocatalytic N-H Bond Insertion with Acceptor-Acceptor Carbene Donors

The biocatalytic toolbox has recently been expanded to include enzyme-catalyzed carbene transfer reactions not occurring in Nature. Herein, we report the development of a biocatalytic strategy for the synthesis of enantioenriched α-trifluoromethyl amines through an asymmetric N–H carbene insertion reaction catalyzed by engineered variants of cytochrome c₅₅₂ from Hydrogenobacter thermophilus. Using a combination of protein and substrate engineering, this metalloprotein scaffold was redesigned to enable the synthesis of chiral α-trifluoromethyl amino esters with up to >99% yield and 95:5 er using benzyl 2-diazotrifluoropropanoate as the carbene donor. When the diazo reagent was varied, the enantioselectivity of the enzyme could be inverted to produce the opposite enantiomers of these products with up to 99.5:0.5 er. This methodology is applicable to a broad range of aryl amine substrates, and it can be leveraged to obtain chemoenzymatic access to enantioenriched β-trifluoromethyl-β-amino alcohols and halides. Computational analyses provide insights into the interplay of protein- and reagent-mediated control on the enantioselectivity of this reaction. This work introduces the first example of a biocatalytic N–H carbenoid insertion with an acceptor–acceptor carbene donor, and it offers a biocatalytic solution for the enantioselective synthesis of α-trifluoromethylated amines as valuable synthons for medicinal chemistry and the synthesis of bioactive molecules.

Use of 3D domain swapping in constructing supramolecular metalloproteins

Supramolecules, which are formed by assembling multiple molecules by noncovalent intermolecular interactions instead of covalent bonds, often show additional properties that cannot be exhibited by a single molecule. Supramolecules have evolved into molecular machines in the field of chemistry, and various supramolecular proteins are responsible for life activities in the field of biology. The design and creation of supramolecular proteins will lead to development of new enzymes, functional biomaterials, drug delivery systems, etc.; thus, the number of studies on the regulation of supramolecular proteins is increasing year by year. Several methods, including disulfide bond, metal coordination, and surface-surface interaction, have been utilized to construct supramolecular proteins. In nature, proteins have been shown to form oligomers by 3D domain swapping (3D-DS), a phenomenon in which a structural region is exchanged between molecules of the same protein. We have been studying the mechanism of 3D-DS and utilizing 3D-DS to construct supramolecular metalloproteins. Cytochrome c forms cyclic oligomers and polymers by 3D-DS, whereas other metalloproteins, such as various c-type cytochromes and azurin form small oligomers and myoglobin forms a compact dimer. We have also utilized 3D-DS to construct heterodimers with different active sites, a protein nanocage encapsulating a Zn-SO₄ cluster in the internal cavity, and a tetrahedron with a designed building block protein. Protein oligomer formation was controlled for the 3D-DS dimer of a dimer-monomer transition protein. This article reviews our research on supramolecular metalloproteins.

Insight into the Folding and Dimerization Mechanisms of the N-Terminal Domain from Human TDP-43

TAR DNA-binding protein 43 (TDP-43) is a 414-residue long nuclear protein whose deposition into intraneuronal insoluble inclusions has been associated with the onset of amyotrophic lateral sclerosis (ALS) and other diseases. This protein is physiologically a homodimer, and dimerization occurs through the N-terminal domain (NTD), with a mechanism on which a full consensus has not yet been reached. Furthermore, it has been proposed that this domain is able to affect the formation of higher molecular weight assemblies. Here, we purified this domain and carried out an unprecedented characterization of its folding/dimerization processes in solution. Exploiting a battery of biophysical approaches, ranging from FRET to folding kinetics, we identified a head-to-tail arrangement of the monomers within the dimer. We found that folding of NTD proceeds through the formation of a number of conformational states and two parallel pathways, while a subset of molecules refold slower, due to proline isomerism. The folded state appears to be inherently prone to form high molecular weight assemblies. Taken together, our results indicate that NTD is inherently plastic and prone to populate different conformations and dimeric/multimeric states, a structural feature that may enable this domain to control the assembly state of TDP-43.

Structural characterization of an on-pathway intermediate and transition state in the folding of the N-terminal SH2 domain from SHP2

Src Homology 2 (SH2) domains are a class of protein domains that present a conserved three-dimensional structure and possess a crucial role in mediating protein-protein interactions. Despite their importance and abundance in the proteome, knowledge about the folding properties of SH2 domain is limited. Here we present an extensive mutational analysis (Φ value analysis) of the folding pathway of the N-SH2 domain of the Src homology region 2 domain-containing phosphatase-2 (SHP2) protein, a 104 residues domain that presents the classical SH2 domain fold (two α-helices flanking a central β-sheet composed of 3-5 antiparallel β-strands), with a fundamental role in mediating the interaction of SHP2 with its substrates and triggering key metabolic pathways in the cell. By analysing folding kinetic data we demonstrated that the folding pathway of N-SH2 presents an obligatory on-pathway intermediate that accumulates during the folding reaction. The production of 24 conservative site-directed variants allowed us to perform a Φ value analysis, by which we could fully characterize the intermediate and the transition state native-like interactions, providing a detailed quantitative analysis of the folding pathway of N-SH2. Results highlight the presence of a hydrophobic nucleus that stabilizes the intermediate, leading to a higher degree of native-like interactions in the transition state. Data are discussed and compared with previous works on SH2 domains.

Influence of heme c attachment on heme conformation and potential

Heme c is characterized by its covalent attachment to a polypeptide. The attachment is typically to a CXXCH motif in which the two Cys form thioether bonds with the heme, "X" can be any amino acid other than Cys, and the His serves as a heme axial ligand. Some cytochromes c, however, contain heme attachment motifs with three or four intervening residues in a CX₃CH or CX₄CH motif. Here, the impacts of these variations in the heme attachment motif on heme ruffling and electronic structure are investigated by spectroscopically characterizing CX₃CH and CX₄CH variants of Hydrogenobacter thermophilus cytochrome c₅₅₂. In addition, a novel CXCH variant is studied. ¹H and ¹³C NMR, EPR, and resonance Raman spectra of the protein variants are analyzed to deduce the extent of ruffling using previously reported relationships between these spectral data and heme ruffling. In addition, the reduction potentials of these protein variants are measured using protein film voltammetry. The CXCH and CX₄CH variants are found to have enhanced heme ruffling and lower reduction potentials. Implications of these results for the use of these noncanonical motifs in nature, and for the engineering of novel heme peptide structures, are discussed.

From protein engineering to artificial enzymes - biological and biomimetic approaches towards sustainable hydrogen production

Hydrogen gas is used extensively in industry today and is often put forward as a suitable energy carrier due its high energy density. Currently, the main source of molecular hydrogen is fossil fuels via steam reforming. Consequently, novel production methods are required to improve the sustainability of hydrogen gas for industrial processes, as well as paving the way for its implementation as a future solar fuel. Nature has already developed an elaborate hydrogen economy, where the production and consumption of hydrogen gas is catalysed by hydrogenase enzymes. In this review we summarize efforts on engineering and optimizing these enzymes for biological hydrogen gas production, with an emphasis on their inorganic cofactors. Moreover, we will describe how our understanding of these enzymes has been applied for the preparation of bio-inspired/-mimetic systems for efficient and sustainable hydrogen production.

Unveiling the folding mechanism of the Bromodomains

Bromodomains (BRDs) are small protein domains often present in large multidomain proteins involved in transcriptional regulation in eukaryotic cells. They currently represent valuable targets for the development of inhibitors of aberrant transcriptional processes in a variety of human diseases. Here we report urea-induced equilibrium unfolding experiments monitored by circular dichroism (CD) and fluorescence on two structurally similar BRDs: BRD2(2) and BRD4(1), showing that BRD4(1) is more stable than BRD2(2). Moreover, we report a description of their kinetic folding mechanism, as obtained by careful analysis of stopped-flow and temperature-jump data. The presence of a high energy intermediate for both proteins, suggested by the non-linear dependence of the folding rate on denaturant concentration in the millisec time regime, has been experimentally observed by temperature-jump experiments. Quantitative global analysis of all the rate constants obtained over a wide range of urea concentrations, allowed us to propose a common, three-state, folding mechanism for these two BRDs. Interestingly, the intermediate of BRD4(1) appears to be more stable and structurally native-like than that populated by BRD2(2). Our results underscore the role played by structural topology and sequence in determining and tuning the folding mechanism.

•

A three-state mechanism for the folding of two representative Bromodomains is proposed.

•

Global analyses of BRD2(2) and BRD4(1) folding kinetics highlights the presence of an on-pathway, folding intermediate.

•

The folding intermediate of BRD4(1) is proposed to be more native-like than that of BRD2(2).

Identification and Structural Characterization of an Intermediate in the Folding of the Measles Virus X Domain

Although most proteins fold by populating intermediates, the transient nature of such states makes it difficult to characterize their structures. In this work we identified and characterized the structure of an intermediate of the X domain of phosphoprotein (P) of measles virus. We obtained this result by a combination of equilibrium and kinetic measurements and NMR chemical shifts used as structural restraints in replica-averaged metadynamics simulations. The structure of the intermediate was then validated by rationally designing four mutational variants predicted to affect the stability of this state. These results provide a detailed view of an intermediate state and illustrate the opportunities offered by a synergistic use of experimental and computational methods to describe non-native states at atomic resolution.

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.

Effects of protein structure on iron-polypeptide vibrational dynamic coupling in cytochrome c

Cytochrome c (Cyt c) has a heme covalently bound to the polypeptide via a Cys-X-X-Cys-His (CXXCH) linker that is located in the interface region for protein-protein interactions. To determine whether the polypeptide matrix influences iron vibrational dynamics, nuclear resonance vibrational spectroscopy (NRVS) measurements were performed on (57)Fe-labeled ferric Hydrogenobacter thermophilus cytochrome c-552, and variants M13V, M13V/K22M, and A7F, which have structural modifications that alter the composition or environment of the CXXCH pentapeptide loop. Simulations of the NRVS data indicate that the 150-325 cm(-1) region is dominated by NHis-Fe-SMet axial ligand and polypeptide motions, while the 325-400 cm(-1) region shows dominant contributions from ν(Fe-NPyr) (Pyr = pyrrole) and other heme-based modes. Diagnostic spectral signatures that directly relate to structural features of the heme active site are identified using a quantum chemistry-centered normal coordinate analysis (QCC-NCA). In particular, spectral features that directly correlate with CXXCH loop stiffness, the strength of the Fe-His interaction, and the degree of heme distortion are identified. Cumulative results from our investigation suggest that compared to the wild type (wt), variants M13V and M13V/K22M have a more rigid CXXCH pentapeptide segment, a stronger Fe-NHis interaction, and a more ruffled heme. Conversely, the A7F variant has a more planar heme and a weaker Fe-NHis bond. These results are correlated to the observed changes in reduction potential between wt protein and the variants studied here. Implications of these results for Cyt c biogenesis and electron transfer are also discussed.

Change in structure and ligand binding properties of hyperstable cytochrome c555 from Aquifex aeolicus by domain swapping

Cytochrome c555 from hyperthermophilic bacteria Aquifex aeolicus (AA cyt c555 ) is a hyperstable protein belonging to the cyt c protein family, which possesses a unique long 310 -α-310 helix containing the heme-ligating Met61. Herein, we show that AA cyt c555 forms dimers by swapping the region containing the extra 310 -α-310 helix and C-terminal α-helix. The asymmetric unit of the crystal of dimeric AA cyt c555 contained two dimer structures, where the structure of the hinge region (Val53-Lys57) was different among all four protomers. Dimeric AA cyt c555 dissociated to monomers at 92 ± 1°C according to DSC measurements, showing that the dimer was thermostable. According to CD measurements, the secondary structures of dimeric AA cyt c555 were maintained at pH 2.2-11.0. CN(-) and CO bound to dimeric AA cyt c555 in the ferric and ferrous states, respectively, owing to the flexibility of the hinge region close to Met61 in the dimer, whereas these ligands did not bind to the monomer under the same conditions. In addition, CN(-) and CO bound to the oxidized and reduced dimer at neutral pH and a wide range of pH (pH 2.2-11.0), respectively, in a wide range of temperature (25-85°C), owing to the thermostability and pH tolerance of the dimer. These results show that the ligand binding character of hyperstable AA cyt c555 changes upon dimerization by domain swapping.

Oligomerization enhancement and two domain swapping mode detection for thermostable cytochrome c552via the elongation of the major hinge loop

High-order oligomers increased whereas N-terminal domain swapping and C-terminal domain swapping were elucidated by the insertion of Gly residues at the major hinge loop of cytochrome c₅₅₂.

The influence of heme ruffling on spin densities in ferricytochromes c probed by heme core 13C NMR

The heme in cytochromes c undergoes a conserved out-of-plane distortion known as ruffling. For cytochromes c from the bacteria Hydrogenobacter thermophilus and Pseudomonas aeruginosa , NMR and EPR spectra have been shown to be sensitive to the extent of heme ruffling and to provide insights into the effect of ruffling on the electronic structure. Through the use of mutants of each of these cytochromes that differ in the amount of heme ruffling, NMR characterization of the low-spin (S = ½) ferric proteins has confirmed and refined the developing understanding of how ruffling influences the spin distribution on heme. The chemical shifts of the core heme carbons were obtained through site-specific labeling of the heme via biosynthetic incorporation of (13)C-labeled 5-aminolevulinic acid derivatives. Analysis of the contact shifts of these core heme carbons allowed Fermi contact spin densities to be estimated and changes upon ruffling to be evaluated. The results allow a deconvolution of the contributions to heme hyperfine shifts and a test of the influence of heme ruffling on the electronic structure and hyperfine shifts. The data indicate that as heme ruffling increases, the spin densities on the β-pyrrole carbons decrease while the spin densities on the α-pyrrole carbons and meso carbons increase. Furthermore, increased ruffling is associated with stronger bonding to the heme axial His ligand.

A new cytoplasmic monoheme cytochrome c from Acidithiobacillus ferrooxidans involved in sulfur oxidation

Acidithiobacillus ferrooxidans can obtain energy from the oxidation of various reduced inorganic sulfur compounds (RISCs, e.g., sulfur) and ferrous iron in bioleaching so has multiple branched respiratory pathways with a diverse range of electron transporters, especially cytochrome c proteins. A cytochrome c family gene, afe1130, which has never been reported before, was found by screening the whole genome of A. ferrooxidans. Here we report the differential gene transcription, bioinformatics analysis, and molecular modeling of the protein encoded by the afe1130 gene (AFE1130). The differential transcription of the target afe1130 gene versus the reference rrs gene in the A. ferrooxidans, respectively, on the culture conditions of sulfur and ferrous energy sources was performed through quantitative reverse transcription polymerase chain reaction (qRT-PCR) with a SYBR green-based assay according to the standard curves method. The qRT-PCR results showed that the afe1130 gene in sulfur culture condition was obviously more transcribed than that in ferrous culture condition. Bioinformatics analysis indicated that the AFE1130 was affiliated to the subclass ID of class I of cytochrome c and located in cytoplasm. Molecular modeling results exhibited that the AFE1130 protein consisted of 5 alpha-helices harboring one heme c group covalently bonded by Cys13 and Cys16 and ligated by His17 and Met62 and owned a big raised hydrophobic surface responsible for attaching to inner cytomembrane. So the AFE1130 in A. ferrooxidans plays a role in the RISCs oxidation in bioleaching in cytoplasm bound to inner membrane.

Structural characterization of nitrosomonas europaea cytochrome c-552 variants with marked differences in electronic structure

Nitrosomonas europaea cytochrome c-552 (Ne c-552) variants with the same His/Met axial ligand set but with different EPR spectra have been characterized structurally, to aid understanding of how molecular structure determines heme electronic structure. Visible light absorption, Raman, and resonance Raman spectroscopy of the protein crystals was performed along with structure determination. The structures solved are those of Ne c-552, which displays a "HALS" (or highly anisotropic low-spin) EPR spectrum, and of the deletion mutant Ne N64Δ, which has a rhombic EPR spectrum. Two X-ray crystal structures of wild-type Ne c-552 are reported; one is of the protein isolated from N. europaea cells (Ne c-552n, 2.35 Å resolution), and the other is of recombinant protein expressed in Escherichia coli (Ne c-552r, 1.63 Å resolution). Ne N64Δ crystallized in two different space groups, and two structures are reported [monoclinic (2.1 Å resolution) and hexagonal (2.3 Å resolution)]. Comparison of the structures of the wild-type and mutant proteins reveals that heme ruffling is increased in the mutant; increased ruffling is predicted to yield a more rhombic EPR spectrum. The 2.35 Å Ne c-552n structure shows 18 molecules in the asymmetric unit; analysis of the structure is consistent with population of more than one axial Met configuration, as seen previously by NMR. Finally, the mutation was shown to yield a more hydrophobic heme pocket and to expel water molecules from near the axial Met. These structures reveal that heme pocket residue 64 plays multiple roles in regulating the axial ligand orientation and the interaction of water with the heme. These results support the hypothesis that more ruffled hemes lead to more rhombic EPR signals in cytochromes c with His/Met axial ligation.

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.

Heme-protein vibrational couplings in cytochrome c provide a dynamic link that connects the heme-iron and the protein surface

The active site of cytochrome c (Cyt c) consists of a heme covalently linked to a pentapeptide segment (Cys-X-X-Cys-His), which provides a link between the heme and the protein surface, where the redox partners of Cyt c bind. To elucidate the vibrational properties of heme c, nuclear resonance vibrational spectroscopy (NRVS) measurements were performed on (57)Fe-labeled ferric Hydrogenobacter thermophilus cytochrome c(552), including (13)C(8)-heme-, (13)C(5)(15)N-Met-, and (13)C(15)N-polypeptide (pp)-labeled samples, revealing heme-based vibrational modes in the 200- to 450-cm(-1) spectral region. Simulations of the NRVS spectra of H. thermophilus cytochrome c(552) allowed for a complete assignment of the Fe vibrational spectrum of the protein-bound heme, as well as the quantitative determination of the amount of mixing between local heme vibrations and pp modes from the Cys-X-X-Cys-His motif. These results provide the basis to propose that heme-pp vibrational dynamic couplings play a role in electron transfer (ET) by coupling vibrations of the heme directly to vibrations of the pp at the protein-protein interface. This could allow for the direct transduction of the thermal (vibrational) energy from the protein surface to the heme that is released on protein/protein complex formation, or it could modulate the heme vibrations in the protein/protein complex to minimize reorganization energy. Both mechanisms lower energy barriers for ET. Notably, the conformation of the distal Met side chain is fine-tuned in the protein to localize heme-pp mixed vibrations within the 250- to 400-cm(-1) spectral region. These findings point to a particular orientation of the distal Met that maximizes ET.

Modulation of ligand-field parameters by heme ruffling in cytochromes c revealed by EPR spectroscopy

Electron paramagnetic resonance (EPR) spectra of variants of Hydrogenobacter thermophilus cytochrome c(552) (Ht c-552) and Pseudomonas aeruginosa cytochrome c(551) (Pa c-551) are analyzed to determine the effect of heme ruffling on ligand-field parameters. Mutations introduced at positions 13 and 22 in Ht c-552 were previously demonstrated to influence hydrogen bonding in the proximal heme pocket and to tune reduction potential (E(m)) over a range of 80 mV [Michel, L. V.; Ye, T.; Bowman, S. E. J.; Levin, B. D.; Hahn, M. A.; Russell, B. S.; Elliott, S. J.; Bren, K. L. Biochemistry 2007, 46, 11753-11760]. These mutations are shown here to also increase heme ruffling as E(m) decreases. The primary effect on electronic structure of increasing heme ruffling is found to be a decrease in the axial ligand-field term Δ/λ, which is proposed to arise from an increase in the energy of the d(xy) orbital. Mutations at position 7, previously demonstrated to influence heme ruffling in Pa c-551 and Ht c-552, are utilized to test this correlation between molecular and electronic structure. In conclusion, the structure of the proximal heme pocket of cytochromes c is shown to play a role in determining heme conformation and electronic structure.

Methionine ligand lability in bacterial monoheme cytochromes c: an electrochemical study

The direct electrochemical analysis of adsorbed redox active proteins has proven to be a powerful technique in biophysical chemistry, frequently making use of the electrode material pyrolytic "edge-plane" graphite. However, many heme-bearing proteins such as cytochromes c have been also examined systematically at alkanethiol-modified gold surfaces, and previously we reported the characterization of the redox properties of a series of bacterial cytochromes c in a side-by-side comparison of carbon and gold electrode materials. In our prior findings, we reported an unanticipated, low potential (E(m) ∼ -100 mV vs SHE) redox couple that could be analogously observed when a variety of monoheme cytochromes c are adsorbed onto carbon-based electrodes. Here we demonstrate that our prior phenomological data can be understood quantitatively in the loss of the methionine ligand of the heme iron, using the cytochrome c from Hydrogenbacter thermophilum as a model system. Through the comparison of wild-type protein with M61H and M61A mutants, in direct electrochemical analyses conducted as a function of temperature and exogenous ligand concentration, we are able to show that Met-ligated cytochromes c have a propensity to lose their Met ligand at graphite surfaces, and that energetics of this process (6.3 ± 0.2 kJ/mol) is similar to the energies associated with "foldons" of known protein folding pathways.

Temperature dependent equilibrium native to unfolded protein dynamics and properties observed with IR absorption and 2D IR vibrational echo experiments

Dynamic and structural properties of carbonmonoxy (CO)-coordinated cytochrome c(552) from Hydrogenobacter thermophilus (Ht-M61A) at different temperatures under thermal equilibrium conditions were studied with infrared absorption spectroscopy and ultrafast two-dimensional infrared (2D IR) vibrational echo experiments using the heme-bound CO as the vibrational probe. Depending on the temperature, the stretching mode of CO shows two distinct bands corresponding to the native and unfolded proteins. As the temperature is increased from low temperature, a new absorption band for the unfolded protein grows in and the native band decreases in amplitude. Both the temperature-dependent circular dichroism and the IR absorption area ratio R(A)(T), defined as the ratio of the area under the unfolded band to the sum of the areas of the native and unfolded bands, suggest a two-state transition from the native to the unfolded protein. However, it is found that the absorption spectrum of the unfolded protein increases its inhomogeneous line width and the center frequency shifts as the temperature is increased. The changes in line width and center frequency demonstrate that the unfolding does not follow simple two-state behavior. The temperature-dependent 2D IR vibrational echo experiments show that the fast dynamics of the native protein are virtually temperature independent. In contrast, the fast dynamics of the unfolded protein are slower than those of the native protein, and the unfolded protein fast dynamics and at least a portion of the slower dynamics of the unfolded protein change significantly, becoming faster as the temperature is raised. The temperature dependence of the absorption spectrum and the changes in dynamics measured with the 2D IR experiments confirm that the unfolded ensemble of conformers continuously changes its nature as unfolding proceeds, in contrast to the native state, which displays a temperature-independent distribution of structures.

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.

The denatured state dictates the topology of two proteins with almost identical sequence but different native structure and function

The protein folding problem is often studied by comparing the mechanisms of proteins sharing the same structure but different sequence. The recent design of the two proteins G(A)88 and G(B)88, displaying different structures and functions while sharing 88% sequence identity (49 out of 56 amino acids), allows the unique opportunity for a complementary approach. At which stage of its folding pathway does a protein commit to a given topology? Which residues are crucial in directing folding mechanisms to a given structure? By using a combination of biophysical and computational techniques, we have characterized the folding of both G(A)88 and G(B)88. We show that, contrary to expectation, G(B)88, characterized by a native α+β fold, displays in the denatured state a content of native-like helical structure greater than G(A)88, which is all-α in its native state. Both experiments and simulations indicate that such residual structure may be tuned by changing pH. Thus, despite the high sequence identity, the folding pathways for these two proteins appear to diverge as early as in the denatured state. Our results suggest a mechanism whereby protein topology is committed very early along the folding pathway, being imprinted in the residual structure of the denatured state.

NMR and DFT investigation of heme ruffling: functional implications for cytochrome c

Out-of-plane (OOP) deformations of the heme cofactor are found in numerous heme-containing proteins and the type of deformation tends to be conserved within functionally related classes of heme proteins. We demonstrate correlations between the heme ruffling OOP deformation and the (13)C and (1)H nuclear magnetic resonance (NMR) hyperfine shifts of heme aided by density functional theory (DFT) calculations. The degree of ruffling in the heme cofactor of Hydrogenobacter thermophilus cytochrome c(552) has been modified by a single amino acid mutation in the second coordination sphere of the cofactor. The (13)C and (1)H resonances of the cofactor have been assigned using one- and two-dimensional NMR spectroscopy aided by selective (13)C-enrichment of the heme. DFT has been used to predict the NMR hyperfine shifts and electron paramagnetic resonance (EPR) g-tensor at several points along the ruffling deformation coordinate. The DFT-predicted NMR and EPR parameters agree with the experimental observations, confirming that an accurate theoretical model of the electronic structure and its response to ruffling has been established. As the degree of ruffling increases, the heme methyl (1)H resonances move upfield while the heme methyl and meso (13)C resonances move downfield. These changes are a consequence of altered overlap of the Fe 3d and porphyrin pi orbitals, which destabilizes all three occupied Fe 3d-based molecular orbitals and decreases the positive and negative spin density on the beta-pyrrole and meso carbons, respectively. Consequently, the heme ruffling deformation decreases the electronic coupling of the cofactor with external redox partners and lowers the reduction potential of heme.

Variation and analysis of second-sphere interactions and axial histidinate character in c-type cytochromes

The electron-donating properties of the axial His ligand to heme iron in cytochromes c (cyts c) are found to be correlated with the midpoint reduction potential (E(m)) in variants of Hydrogenobacter thermophilus cytochrome c(552) (Ht cyt c(552)) in which mutations have been made in and near the Cys-X-X-Cys-His (CXXCH) heme-binding motif. To probe the strength of the His-Fe(III) interaction, we have measured (13)C nuclear magnetic resonance (NMR) chemical shifts for (13)CN(-) bound to heme iron trans to the axial His in Ht Fe(III) cyt c(552) variants. We observe a linear relationship between these (13)C chemical shifts and E(m), indicating that the His-Fe(III) bond strength correlates with E(m). To probe a conserved hydrogen bonding interaction between the axial His Hdelta1 and the backbone carbonyl of a Pro residue, we measured the pK(a) of the axial His Hdelta1 proton (pK(a(2))), which we propose to relate to the His-Fe(III) interaction, reduction potential, and local electrostatic effects. The observed linear relationship between the axial His (13)Cbeta chemical shift and E(m) is proposed to reflect histidinate (anionic) character of the ligand. A linear relationship also is seen between the average heme methyl (1)H chemical shift and E(m) which may reflect variation in axial His electron-donating properties or in the ruffling distortion of the heme plane. In summary, chemical shifts of the axial His and exogenous CN(-) bound trans to His are shown to be sensitive probes of the His-Fe(III) interaction in variants of Ht cyt c(552) and display trends that correlate with E(m).

What lessons can be learned from studying the folding of homologous proteins?

The studies of the folding of structurally related proteins have proved to be a very important tool for investigating protein folding. Here we review some of the insights that have been gained from such studies. Our highlighted studies show just how such an investigation should be designed and emphasise the importance of the synergy between experiment and theory. We also stress the importance of choosing the right system carefully, exploiting the excellent structural and sequence databases at our disposal.

Review: studies of ferric heme proteins with highly anisotropic/highly axial low spin (S = 1/2) electron paramagnetic resonance signals with bis-histidine and histidine-methionine axial iron coordination

Six-coordinated heme groups are involved in a large variety of electron transfer reactions because of their ability to exist in both the ferrous (Fe(2+)) and ferric (Fe(3+)) state without any large differences in structure. Our studies on hemes coordinated by two histidines (bis-His) and hemes coordinated by histidine and methionine (His-Met) will be reviewed. In both of these coordination environments, the heme core can exhibit ferric low spin (electron paramagnetic resonance EPR) signals with large g(max) values (also called Type I, highly anisotropic low spin, or highly axial low spin, HALS species) as well as rhombic EPR (Type II) signals. In bis-His coordinated hemes rhombic and HALS envelopes are related to the orientation of the His groups with respect to each other such that (i) parallel His planes results in a rhombic signal and (ii) perpendicular His planes results in a HALS signal. Correlation between the structure of the heme and its ligands for heme with His-Met axial ligation and ligand-field parameters, as derived from a large series of cytochrome c variants, show, however, that for such a combination of axial ligands there is no clear-cut difference between the large g(max) and the "small g-anisotropy" cases as a result of the relative Met-His arrangements. Nonetheless, a new linear correlation links the average shift delta of the heme methyl groups with the g(max) values.

Early events, kinetic intermediates and the mechanism of protein folding in cytochrome C

Kinetic studies of the early events in cytochrome c folding are reviewed with a focus on the evidence for folding intermediates on the submillisecond timescale. Evidence from time-resolved absorption, circular dichroism, magnetic circular dichroism, fluorescence energy and electron transfer, small-angle X-ray scattering and amide hydrogen exchange studies on the t ≤ 1 ms timescale reveals a picture of cytochrome c folding that starts with the ~ 1-μs conformational diffusion dynamics of the unfolded chains. A fractional population of the unfolded chains collapses on the 1 – 100 μs timescale to a compact intermediate I_C containing some native-like secondary structure. Although the existence and nature of I_C as a discrete folding intermediate remains controversial, there is extensive high time-resolution kinetic evidence for the rapid formation of I_C as a true intermediate, i.e., a metastable state separated from the unfolded state by a discrete free energy barrier. Final folding to the native state takes place on millisecond and longer timescales, depending on the presence of kinetic traps such as heme misligation and proline mis-isomerization. The high folding rates observed in equilibrium molten globule models suggest that I_C may be a productive folding intermediate. Whether it is an obligatory step on the pathway to the high free energy barrier associated with millisecond timescale folding to the native state, however, remains to be determined.

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.

Cloning, expression and purification of cytochrome c(6) from the brown alga Hizikia fusiformis and complete X-ray diffraction analysis of the structure

The primary sequence of cytochrome c(6) from the brown alga Hizikia fusiformis has been determined by cDNA cloning and the crystal structure has been solved at 1.6 A resolution. The crystal belonged to the tetragonal space group P4(1)2(1)2, with unit-cell parameters a = b = 84.58, c = 232.91 A and six molecules per asymmetric unit. The genome code, amino-acid sequence and crystal structure of H. fusiformis cytochrome c(6) were most similar to those of red algal cytochrome c(6). These results support the hypothesis that brown algae acquired their chloroplasts via secondary endosymbiosis involving a red algal endosymbiont and a eukaryote host.

Heme attachment motif mobility tunes cytochrome c redox potential

Hydrogen exchange (HX) rates and midpoint potentials (Em) of variants of cytochrome c from Pseudomonas aeruginosa (Pa cyt c551) and Hydrogenobacter thermophilus (Ht cyt c552) have been characterized in an effort to develop an understanding of the impact of properties of the Cys-X-X-Cys-His pentapeptide c-heme attachment (CXXCH) motif on heme redox potential. Despite structural conservation of the CXXCH motif, Ht cyt c552 exhibits a low level of protection from HX for amide protons within this motif relative to Pa cyt c551. Site-directed mutants have been prepared to determine the structural basis for and functional implications of these variations on HX behavior. The double mutant Ht-M13V/K22M displays suppressed HX within the CXXCH motif as well as a decreased Em (by 81 mV), whereas the corresponding double mutant of Pa cyt c551 (V13M/M22K) exhibits enhanced HX within the CXXCH pentapeptide and a modest increase in Em (by 30 mV). The changes in Em correlate with changes in axial His chemical shifts in the ferric proteins reflecting the extent of histidinate character. Thus, the mobility of the CXXCH pentapeptide is found to impact the His-Fe(III) interaction and therefore the heme redox potential.

Characterization of N-terminal amino group–heme ligation emerging upon guanidine hydrochloric acid induced unfolding of Hydrogenobacter thermophilus ferricytochrome c 552

Nonnative heme coordination structures emerging upon guanidine hydrochloric acid (GdnHCl) induced unfolding of Hydrogenobacter thermophilus ferricytochrome c552 were characterized by means of paramagnetic NMR. The heme coordination structure possessing the N-terminal amino group of the peptide chain in place of axial Met (His-Nterm form) was determined in the presence of GdnHCl concentrations in excess of 1.5 M at neutral pH. The stability of the His-Nterm form at pH 7.0 was found to be comparable with that of the bis-His form which has been recognized as a major nonnative heme coordination structure in cytochrome c folding/unfolding. Consequently, in addition to the bis-His form, the His-Nterm form is a substantial intermediate which affects the pathway and kinetics of the folding/unfolding of cytochromes c, of which the N-terminal amino groups are not acetylated.

Local Conformational Transition of Hydrogenobacter thermophilus Cytochrome c552 Relevant to Its Redox Potential,

In order to elucidate the molecular mechanisms responsible for the apparent nonlinear behavior of the temperature dependence of the redox potential of Hydrogenobacter thermophilus cytochrome c552 [Takahashi, Y., Sasaki, H., Takayama, S. J., Mikami, S., Kawano, S., Mita, H., Sambongi, Y., and Yamamoto, Y. (2006) Biochemistry 45, 11005-11011], its heme active site structure has been characterized using variable-temperature and -pressure NMR techniques. The study revealed a temperature-dependent conformational transition between protein structures, which slightly differ in the conformation of the loop bearing the Fe-bound axial Met residue. The heme environment in the protein structure which arises at lower temperature was found to be more polar, as a result of the altered orientation of the loop with respect to the heme due to its conformational change, than that arising at higher temperature. The present study demonstrated the importance of the structural and dynamic properties of the polypeptide chain in close proximity to the heme for redox regulation of the protein.

Plasticity of the protein folding landscape: switching between on- and off-pathway intermediates

Proteins may fold via parallel routes partitioned by the relative effect of solvent conditions on the relevant transition states. Thus, intermediates may or may not necessarily be obligatory species accumulated during the folding process, but rather kinetic traps due to the ruggedness of the folding landscape. Implicit in this view is the notion of plasticity of the folding pathway: proteins can be rerouted through the energy landscape by mutational, topological or solvent perturbations. Our work was specifically aimed to the experimental identification of a switch in the folding mechanism of a c-type cytochrome from the thermophilic bacterium Hydrogenobacter thermophilus (HT cyt c(552)) induced by acidic conditions. We present evidence that, by destabilizing the relevant transition state, the native state of HT cyt c(552) can be reached along alternative folding routes, which may involve an off-pathway intermediate.