Stability and apoptotic activity of recombinant human cytochrome c

An efficient system for producing human cytochrome c variants is important to help us understand the roles of this protein in biological processes relevant to human diseases including apoptosis and oxidative stress. Here, we describe an Escherichia coli expression system for producing recombinant human cytochrome c. We also characterize the structure, stability, and function of the protein and show its utility for studying apoptosis. Yields of greater than 8 mg of pure protein per liter culture were attained. Circular dichroism spectropolarimetry studies show that the secondary and tertiary structures of the human protein are nearly identical to those of the horse protein, but the human protein is more stable than other eukaryotic cytochromes c. Furthermore, recombinant human cytochrome c is capable of inducing caspase-3 activity in a cell-free caspase activation assay. We use data from this assay along with data from the literature to define the apaf-1 binding site on human cytochrome c.     

Folding barrier in horse cytochrome c: support for a classical folding pathway

Native-state structures and conformations of ferrocytochrome c, nitrosylcytochrome c, and carbonmonoxycytochrome c are very similar. They are, however, immensely different from each other in terms of thermodynamic stability. The dramatic destabilization of ferrocytochrome c to the extent of 12 kcal mol(-1) produces no effect on the folding rate, and this is so in spite of the fact that all three test-tube variants fold in an apparent two-state manner. For all three proteins the folding barrier is early in time, sizable in energy, and is of the same magnitude (approximately 6.5 kcal mol(-1)). These results raise some challenges to the "new view" of protein folding. An early transition state, the search for which consumes most of the observed folding time, is suggested.     

Exploring the energy landscape of protein folding using replica-exchange and conventional molecular dynamics simulations

Two independent replica-exchange molecular dynamics (REMD) simulations with an explicit water model were performed of the Trp-cage mini-protein. In the first REMD simulation, the replicas started from the native conformation, while in the second they started from a nonnative conformation. Initially, the first simulation yielded results qualitatively similar to those of two previously published REMD simulations: the protein appeared to be over-stabilized, with the predicted melting temperature 50-150K higher than the experimental value of 315K. However, as the first REMD simulation progressed, the protein unfolded at all temperatures. In our second REMD simulation, which starts from a nonnative conformation, there was no evidence of significant folding. Transitions from the unfolded to the folded state did not occur on the timescale of these simulations, despite the expected improvement in sampling of REMD over conventional molecular dynamics (MD) simulations. The combined 1.42 micros of simulation time was insufficient for REMD simulations with different starting structures to converge. Conventional MD simulations at a range of temperatures were also performed. In contrast to REMD, the conventional MD simulations provide an estimate of Tm in good agreement with experiment. Furthermore, the conventional MD is a fraction of the cost of REMD and continuous, realistic pathways of the unfolding process at atomic resolution are obtained.     

Direct observation of microscopic reversibility in single-molecule protein folding

Both folded and unfolded conformations should be observed for a protein at its melting temperature (T(m)), where DeltaG between these states is zero. In an all-atom molecular dynamics simulation of chymotrypsin inhibitor 2 (CI2) at its experimental T(m), the protein rapidly loses its low-temperature native structure; it then unfolds before refolding to a stable, native-like conformation. The initial unfolding follows the unfolding pathway described previously for higher-temperature simulations: the hydrophobic core is disrupted, the beta-sheet pulls apart and the alpha-helix unravels. The unfolded state reached under these conditions maintains a kernel of structure in the form of a non-native hydrophobic cluster. Refolding simply reverses this path, the side-chain interactions shift, the helix refolds, and the native packing and hydrogen bonds are recovered. The end result of this refolding is not the initial crystal structure; it contains the proper topology and the majority of the native contacts, but the structure is expanded and the contacts are long. We believe this to be the native state at elevated temperature, and the change in volume and contact lengths is consistent with experimental studies of other native proteins at elevated temperature and the chemical denaturant equivalent of T(m).     

Fast folding kinetics and stabilization of apo-cytochrome c

It is generally accepted that in the c-type cytochromes the covalently bound heme plays a primary role in the acquisition of the folded state. Here, we show that a stabilized site-directed variant of apo-cyt c551 from Pseudomonas aeruginosa (Pa-apocyt F7A/W77F) retains native-like features in the presence of sodium sulfate even in the absence of heme. By time-resolved intrinsic fluorescence, we have evidence that Pa-apocyt F7A/W77F may acquire a compact, native-like conformation within microseconds. These results challenge current thinking about the role of the heme group in the folding of c-type cytochromes.     

Yeast cytochrome c is a sequence-specific DNA-binding protein

A protein binding to the alcohol oxidase 2 upstream activation sequence (AOX2UAS) of the methylotropic yeast, Pichia pastoris, has been purified and identified as cytochrome c (cyt c). Cyt c purified from P. pastoris or Saccharomyces cerevisiae binds to AOX2UAS. Specific point mutations in AOX2UAS abolish cyt c binding. We conclude that yeast cyt c is a sequence-specific DNA-binding protein and may have a regulatory role in the nucleus.     

Parallel pathways in cytochrome c(551) folding

The folding of cytochrome c(551) from Pseudomonas aeruginosa was previously thought to follow a simple sequential mechanism, consistent with the lack of histidine residues, other than the native His16 heme ligand, that can give rise to mis-coordinated species. However, further kinetic analysis reveals complexities indicative of a folding mechanism involving parallel pathways. Double-jump interrupted refolding experiments at low pH indicate that approximately 50% of the unfolded cytochrome c(551) population can reach the native state via a fast (10 ms) folding track, while the rest follows a slower folding path with populated intermediates. Stopped-flow experiments using absorbance at 695 nm to monitor refolding confirm the presence of a rapidly folding species containing the native methionine-iron bond while measurements on carboxymethylated cytochrome c(551) (which lacks the Met-Fe coordination bond) indicate that methionine ligation occurs late during folding along the fast folding track, which appears to be dominant at physiological pH. Continuous-flow measurements of tryptophan-heme energy transfer, using a capillary mixer with a dead time of about 60 micros, show evidence for a rapid chain collapse within 100 micros preceding the rate-limiting folding phase on the milliseconds time scale. A third process with a time constant in the 10-50 ms time range is consistent with a minor population of molecules folding along a parallel channel, as confirmed by quantitative kinetic modeling. These findings indicate the presence of two or more slowly inter-converting ensembles of denatured states that give rise to pH-dependent partitioning among fast and slow-folding pathways.     

Cooperative folding in a multi-domain protein

Most protein domains are found in multi-domain proteins, yet most studies of protein folding have concentrated on small, single-domain proteins or on isolated domains from larger proteins. Spectrin domains are small (106 amino acid residues), independently folding domains consisting of three long alpha-helices. They are found in multi-domain proteins with a number of spectrin domains in tandem array. Structural studies have shown that in these arrays the last helix of one domain forms a continuous helix with the first helix of the following domain. It has been demonstrated that a number of spectrin domains are stabilised by their neighbours. Here we investigate the molecular basis for cooperativity between adjacent spectrin domains 16 and 17 from chicken brain alpha-spectrin (R16 and R17). We show that whereas the proteins unfold as a single cooperative unit at 25 degrees C, cooperativity is lost at higher temperatures and in the presence of stabilising salts. Mutations in the linker region also cause the cooperativity to be lost. However, the cooperativity does not rely on specific interactions in the linker region alone. Most mutations in the R17 domain cause a decrease in cooperativity, whereas proteins with mutations in the R16 domain still fold cooperatively. We propose a mechanism for this behaviour.     

5-Fluorotryptophan as dual probe for ground-state heterogeneity and excited-state dynamics in apoflavodoxin

The apoflavodoxin protein from Azotobacter vinelandii harboring three tryptophan (Trp) residues, was biosynthetically labeled with 5-fluorotryptophan (5-FTrp). 5-FTrp has the advantage that chemical differences in its microenvironment can be sensitively visualized via ¹⁹F NMR. Moreover, it shows simpler fluorescence decay kinetics. The occurrence of FRET was earlier observed via the fluorescence anisotropy decay of WT apoflavodoxin and the anisotropy decay parameters are in excellent agreement with distances between and relative orientations of all Trp residues. The anisotropy decay in 5-FTrp apoflavodoxin demonstrates that the distances and orientations are identical for this protein. This work demonstrates the added value of replacing Trp by 5-FTrp to study structural features of proteins via ¹⁹F NMR and fluorescence spectroscopy.     

Intimate view of a kinetic protein folding intermediate: residue-resolved structure, interactions, stability, folding and unfolding rates, homogeneity

A cytochrome c kinetic folding intermediate was studied by hydrogen exchange (HX) pulse labeling. Advances in the technique and analysis made it possible to define the structured and unstructured regions, equilibrium stability, and kinetic opening and closing rates, all at an amino acid-resolved level. The entire N-terminal and C-terminal helices are formed and docked together at their normal native positions. They fray in both directions from the interaction region, due to a progression in both unfolding and refolding rates, leading to the surprising suggestion that helix propagation may proceed very slowly in the condensed milieu. Several native-like beta turns are formed. Some residues in the segment that will form the native 60s helix are protected but others are not, suggesting energy minimization to some locally non-native conformation in the transient intermediate. All other regions are unprotected, presumably dynamically disordered. The intermediate resembles a partially constructed native state. It is early, on-pathway, and all of the refolding molecules pass through it. These and related results consistently point to distinct, homogeneous, native-like intermediates in a stepwise sequential pathway, guided by the same factors that determine the native structure. Previous pulse labeling efforts have always assumed EX2 exchange during the labeling pulse, often leading to the suggestion of heterogeneous intermediates in alternative parallel pathways. The present work reveals a dominant role for EX1 exchange in the high pH labeling pulse, which will mimic heterogeneous behavior when EX2 exchange is assumed. The general problem of homogeneous versus heterogeneous intermediates and pathways is discussed.     

Production and preliminary characterization of a recombinant triheme cytochrome c7 from Geobacter sulfurreducens in Escherichia coli

Multiheme cytochromes c have been found in a number of sulfate- and metal ion-reducing bacteria. Geobacter sulfurreducens is one of a family of microorganisms that oxidize organic compounds, with Fe(III) oxide as the terminal electron acceptor. A triheme 9.6 kDa cytochrome c₇ from G. sulfurreducens is a part of the metal ion reduction pathway. We cloned the gene for cytochrome c₇ and expressed it in Escherichiacoli together with the cytochrome c maturation gene cluster, ccmABCDEFGH, on a separate plasmid. We designed two constructs, with and without an N-terminal His-tag. The untagged version provided a good yield (up to 6 mg/l of aerobic culture) of the fully matured protein, with all three hemes attached, while the N-terminal His-tag appeared to be detrimental for proper heme incorporation. The recombinant protein (untagged) is properly folded, it has the same molecular weight and displays the same absorption spectra, both in reduced and in oxidized forms, as the protein isolated from G. sulfurreducens and it is capable of reducing metal ions in vitro. The shape parameters for the recombinant cytochrome c₇ determined by small angle X-ray scattering are in good agreement with the ones calculated from a homologous cytochrome c₇ of known structure.     

Nucleotide binding to nucleoside diphosphate kinases: X-ray structure of human NDPK-A in complex with ADP and comparison to protein kinases

NDPK-A, product of the nm23-H1 gene, is one of the two major isoforms of human nucleoside diphosphate kinase. We analyzed the binding of its nucleotide substrates by fluorometric methods. The binding of nucleoside triphosphate (NTP) substrates was detected by following changes of the intrinsic fluorescence of the H118G/F60W variant, a mutant protein engineered for that purpose. Nucleoside diphosphate (NDP) substrate binding was measured by competition with a fluorescent derivative of ADP, following the fluorescence anisotropy of the derivative. We also determined an X-ray structure at 2.0A resolution of the variant NDPK-A in complex with ADP, Ca(2+) and inorganic phosphate, products of ATP hydrolysis. We compared the conformation of the bound nucleotide seen in this complex and the interactions it makes with the protein, with those of the nucleotide substrates, substrate analogues or inhibitors present in other NDP kinase structures. We also compared NDP kinase-bound nucleotides to ATP bound to protein kinases, and showed that the nucleoside monophosphate moieties have nearly identical conformations in spite of the very different protein environments. However, the beta and gamma-phosphate groups are differently positioned and oriented in the two types of kinases, and they bind metal ions with opposite chiralities. Thus, it should be possible to design nucleotide analogues that are good substrates of one type of kinase, and poor substrates or inhibitors of the other kind.     

FRET-based in vivo screening for protein folding and increased protein stability

Fluorescence resonance energy transfer (FRET) was used to establish a novel in vivo screening system that allows rapid detection of protein folding and protein variants with increased thermodynamic stability in the cytoplasm of Escherichia coli. The system is based on the simultaneous fusion of the green fluorescent protein (GFP) to the C terminus of a protein X of interest, and of blue-fluorescent protein (BFP) to the N terminus of protein X. Efficient FRET from BFP to GFP in the ternary fusion protein is observed in vivo only when protein X is folded and brings BFP and GFP into close proximity, while FRET is lost when BFP and GFP are far apart due to unfolding or intracellular degradation of protein X. The screening system was validated by identification of antibody V(L) intradomains with increased thermodynamic stabilities from expression libraries after random mutagenesis, bacterial cell sorting, and colony screening.     

Cooperative omega loops in cytochrome c: role in folding and function

Hydrogen exchange experiments under slow exchange conditions show that an omega loop in cytochrome c (residues 40-57) acts as a cooperative unfolding/refolding unit under native conditions. This unit behavior accounts for an initial step on the unfolding pathway, a final step in refolding, and a number of other structural, functional and evolutionary properties.     

Structure analysis of the protein translocating channel TatA in membranes using a multi-construct approach

The twin-arginine-translocase (Tat) can transport proteins in their folded state across bacterial or thylakoid membranes. In Bacillus subtilis the Tat-machinery consists of only two integral (inner) membrane proteins, TatA and TatC. Multiple copies of TatA are supposed to form the transmembrane channel, but little structural data is available on this 70-residue component. We used a multi-construct approach for expressing several characteristic fragments of TatA_d, to determine their individual structures and to cross-validate them comprehensively within the architecture of the full-length protein. Here, we report the design, high-yield expression, detergent-aided purification and lipid-reconstitution of five constructs of TatA_d, overcoming difficulties associated with the very different hydrophobicities and sizes of these membrane protein fragments. Circular dichroism (CD) and oriented CD (OCD) were used to determine their respective conformations and alignments in suitable, negatively charged phospholipid bilayers. CD spectroscopy showed an N-terminal α-helix, a central helical stretch, and an unstructured C-terminus, thus proving the existence of these secondary structures in TatA_d for the first time. The OCD spectra demonstrated a transmembrane orientation of the N-terminal α-helix and a surface alignment of the central amphiphilic helix in lipid bilayers, thus supporting the postulated topology model and function of TatA as a transmembrane channel.     

First characterization of co-chaperonin protein 10 from hyper-thermophilic Aquifex aeolicus

All known co-chaperonin protein 10 (cpn10) molecules are heptamers of seven identical subunits that are linked together by beta-strand interactions. Here, we report the first characterization of a cpn10 protein from a thermophilic organism: Aquifex aeolicus. Primary-structure alignment of A. aeolicus cpn10 (Aaecpn10) shows high homology with mesophilic cpn10 sequences, except for a unique 25-residue C-terminal extension not found in any other cpn10. Recombinant Aaecpn10 adopts a heptameric structure in solution at pH values above 4 (20 degrees C). Both monomers and heptamers are folded at 20 degrees C, although the thermal stability of the monomers (pH 3; Tm approximately 58 degrees C) is lower than that of the heptamers (pH 7; Tm approximately 115 degrees C). Aaecpn10 functions in a GroEL-dependent in vitro activity assay. Taken together, Aaecpn10 appears similar in secondary, tertiary, and quaternary structure, as well as in many biophysical features, to its mesophilic counterparts despite a functional temperature of 90 degrees C.     

Analysis of the stability of cytochrome c(6) with an improved stopped-flow protocol

This work presents an improved stopped-flow protocol for the simultaneous measurement of thermodynamic and kinetic protein stability data from a single experiment, along with a formalism for the global analysis of the data. The method was applied to the comparison of the stabilities of cytochrome c(6) from Anabaena sp. PCC 7119 and one of its mutants (D72K). Compared to the wild type the mutant was found to have a significantly reduced thermodynamic (deltadeltaG(U0)=2.7 kJ mol(-1)) and kinetic stability, as well as a more pronounced shift in transition state structure upon destabilization.     

The effects of ATP and sodium chloride on the cytochrome c-cardiolipin interaction: the contrasting behavior of the horse heart and yeast proteins

In cells a portion of cytochrome c (cyt c) (15–20%) is tightly bound to cardiolipin (CL), one of the phospholipids constituting the mitochondrial membrane. The CL-bound protein, which has nonnative tertiary structure, altered heme pocket, and disrupted Fe(III)-M80 axial bond, is thought to play a role in the apoptotic process. This has attracted considerable interest in order to clarify the mechanisms governing the cyt c–CL interaction. Herein we have investigated the binding reaction of CL with the c-type cytochromes from horse heart and yeast. Although the two proteins possess a similar tertiary architecture, yeast cyt c displays lower stability and, contrary to the equine protein, it does not bind ATP and lacks pro-apoptotic activity. The study has been performed in the absence and in the presence of ATP and NaCl, two compounds that influence the (horse cyt c)-CL binding process and, thus, the pro-apoptotic activity of the protein. The two proteins behave differently: while CL interaction with horse cyt c is strongly influenced by the two effectors, no effect is observed for yeast cyt c. It is noteworthy that NaCl induces dissociation of the (horse cyt c)–CL complex but has no influence on that of yeast cyt c. The differences found for the two proteins highlight that specific structural factors, such as the different local structure conformation of the regions involved in the interactions with either CL or ATP, can significantly affect the behavior of cyt c in its reaction with liposomes and the subsequent pro-apoptotic action of the protein.     

Protein disulfide isomerase assisted protein folding in malaria parasites

In eukaryotes, the formation of protein disulfide bonds among cysteine residues is mediated by protein disulfide isomerases and occurs in the highly oxidised environment of the endoplasmic reticulum. This process is poorly understood in malaria parasites. In this paper, we report the gene isolation, sequence and phylogenetic comparisons, protein structure and thioredoxin-domain analyses of nine protein disulfide isomerases-like molecules from five species of malaria parasites including Plasmodium falciparum and Plasmodium vivax (human), Plasmodium knowlesi (simian) and Plasmodium berghei and Plasmodium yoelii (murine). Four of the studied protein disulfide isomerases belong to P. falciparum malaria and have been named PfPDI-8, PfPDI-9, PfPDI-11 and PfPDI-14, based on their chromosomal location. Among these, PfPDI-8 bears the closest similarity to a prototype PDI molecule with two thioredoxin domains (containing CGHC active sites) and a C-terminal Endoplasmic reticulum retrieval signal, SEEL. PfPDI-8 is expressed during all stages of parasite life cycle and is highly conserved (82-96% identity at amino acid level) in the other four Plasmodium species studied. Detailed biochemical analysis of PfPDI-8 revealed that this molecule is a potent oxido-reductase enzyme that facilitated the disulfide-dependent conformational folding of EBA-175, a leading malaria vaccine candidate. These studies open the avenues to understand the process of protein folding and secretory pathway in malaria parasites that in turn might aid in the production of superior recombinant vaccines and provide novel drug targets.