Studies of Pseudomonas aeruginosa azurin mutants: cavities in beta-barrel do not affect refolding speed

Significance of the Disulfide Bridge in the Structure and Stability of Metalloprotein Azurin

Metalloproteins make up a class of proteins that incorporate metal ions into their structures, enabling them to perform essential functions in biological systems, such as catalysis and electron transport. Azurin is one such metalloprotein with copper cofactor, having a β-barrel structure with exceptional thermal stability. The copper metal ion is coordinated at one end of the β-barrel structure, and there is a disulfide bond at the opposite end. In this study, we explore the effect of this disulfide bond in the high thermal stability of azurin by analyzing both the native S-S bonded and S-S nonbonded (S-S open) forms using temperature replica exchange molecular dynamics (REMD). Similar to experimental observations, we find a 35 K decrease in denaturation temperature for S-S open azurin compared to that of the native holo form (420 K). As observed in the case of native holo azurin, the unfolding process of the S-S open form also started with disruptions of the α-helix. The free energy surfaces of the unfolding process revealed that the denaturation event of the S-S open form progresses through different sets of conformational ensembles. Subsequently, we compared the stabilities of individual β-sheet strands of both the S-S bonded and the S-S nonbonded forms of azurin. Further, we examined the contacts between individual residues for the central structures from the free energy surfaces of the S-S nonbonded form. The microscopic origin of the lowering in the denaturation temperature is further supplemented by thermodynamic analysis.

Role of Metal Cofactor in Enhanced Thermal Stability of Azurin

Metal cofactors are critical centers for different biochemical processes of metalloproteins, and often, this metal coordination renders additional structural stability. In this study, we explore the additional stability conferred by the copper ion on azurin by analyzing both the apo and holo forms using temperature replica exchange molecular dynamics (REMD) data. We find a 14 K decrease in denaturation temperature for apo (406 K) azurin relative to that of holo (420 K), indicating a copper ion-induced additional thermal stability for holo azurin. The unfolding of apo azurin begins with the melting of α-helix and β-sheet V, similar to that of holo form. β-Sheets IV, VII, and VIII are comparatively more stable than other β-strands and melt at higher temperatures. Similar to holo azurin, the strong hydrophobic interactions among the apolar residues in the protein core is the key factor that renders high stability to apo protein as well. We construct free energy surfaces at different temperatures to capture the major conformations along the unfolding basins of the protein. Using contact maps from different basins we show the changes in the interaction between different residues along the unfolding pathway. Furthermore, we compare the Cα root-mean-square fluctuations (Cα-RMSF) and B-factor of all residues of apo and holo forms to understand the flexibility of different regions. The concerted displacement of α-helix and β-sheets V and VI from the protein core is another distinction we observe for apo compared to the holo form, where β-sheet VI was relatively stable.

Folding of copper proteins: role of the metal?

Copper is a redox-active transition metal ion required for the function of many essential human proteins. For biosynthesis of proteins coordinating copper, the metal may bind before, during or after folding of the polypeptide. If the metal binds to unfolded or partially folded structures of the protein, such coordination may modulate the folding reaction. The molecular understanding of how copper is incorporated into proteins requires descriptions of chemical, thermodynamic, kinetic and structural parameters involved in the formation of protein-metal complexes. Because free copper ions are toxic, living systems have elaborate copper-transport systems that include particular proteins that facilitate efficient and specific delivery of copper ions to target proteins. Therefore, these pathways become an integral part of copper protein folding in vivo. This review summarizes biophysical-molecular in vitro work assessing the role of copper in folding and stability of copper-binding proteins as well as protein-protein copper exchange reactions between human copper transport proteins. We also describe some recent findings about the participation of copper ions and copper proteins in protein misfolding and aggregation reactions in vitro.

Quantification of excluded volume effects on the folding landscape of Pseudomonas aeruginosa apoazurin in vitro

Proteins fold and function inside cells that are crowded with macromolecules. Here, we address the role of the resulting excluded volume effects by in vitro spectroscopic studies of Pseudomonas aeruginosa apoazurin stability (thermal and chemical perturbations) and folding kinetics (chemical perturbation) as a function of increasing levels of crowding agents dextran (sizes 20, 40, and 70 kDa) and Ficoll 70. We find that excluded volume theory derived by Minton quantitatively captures the experimental effects when crowding agents are modeled as arrays of rods. This finding demonstrates that synthetic crowding agents are useful for studies of excluded volume effects. Moreover, thermal and chemical perturbations result in free energy effects by the presence of crowding agents that are identical, which shows that the unfolded state is energetically the same regardless of method of unfolding. This also underscores the two-state approximation for apoazurin's unfolding reaction and suggests that thermal and chemical unfolding experiments can be used in an interchangeable way. Finally, we observe increased folding speed and invariant unfolding speed for apoazurin in the presence of macromolecular crowding agents, a result that points to unfolded-state perturbations. Although the absolute magnitude of excluded volume effects on apoazurin is only on the order of 1-3 kJ/mol, differences of this scale may be biologically significant.

Role of metal in folding and stability of copper proteins in vitro

Metal coordination is required for function of many proteins. For biosynthesis of proteins coordinating a metal, the question arises if the metal binds before, during or after folding of the polypeptide. Moreover, when the metal is bound to the protein, how does its coordination affect biophysical properties such as stability and dynamics? Understanding how metals are utilized by proteins in cells on a molecular level requires accurate descriptions of the thermodynamic and kinetic parameters involved in protein-metal complexes. Copper is one of the essential transition metals found in the active sites of many key proteins. To avoid toxicity of free copper ions, living systems have developed elaborate copper-transport systems that involve dedicated proteins that facilitate efficient and specific delivery of copper to target proteins. This review describes in vitro and in silico biophysical work assessing the role of copper in folding and stability of copper-binding proteins. Examples of proteins discussed are: a blue-copper protein (Pseudomonas aeruginosa azurin), members of copper-transport systems (bacterial CopZ, human Atox1 and ATP7B domains) and multi-copper ferroxidases (yeast Fet3p and human ceruloplasmin). The consequences of interactions between copper proteins and platinum-complexes are also discussed. This article is part of a Special Issue entitled: Cell Biology of Metals.

A heme fusion tag for protein affinity purification and quantification

We report a novel affinity-based purification method for proteins expressed in Escherichia coli that uses the coordination of a heme tag to an L-histidine-immobilized sepharose (HIS) resin. This approach provides an affinity purification tag visible to the eye, facilitating tracking of the protein. We show that azurin and maltose binding protein are readily purified from cell lysate using the heme tag and HIS resin. Mild conditions are used; heme-tagged proteins are bound to the HIS resin in phosphate buffer, pH 7.0, and eluted by adding 200-500 mM imidazole or binding buffer at pH 5 or 8. The HIS resin exhibits a low level of nonspecific binding of untagged cellular proteins for the systems studied here. An additional advantage of the heme tag-HIS method for purification is that the heme tag can be used for protein quantification by using the pyridine hemochrome absorbance method for heme concentration determination.

Investigation of an anomalously accelerating substitution in the folding of a prototypical two-state protein

The folding rates of two-state single-domain proteins are generally resistant to small-scale changes in amino acid sequence. For example, having surveyed here over 700 single-residue substitutions in 24 well-characterized two-state proteins, we find that the majority (55%) of these substitutions affect folding rates by less than a factor of 2, and that only 9% affect folding rates by more than a factor of 8. Among those substitutions that significantly affect folding rates, we find that accelerating substitutions are an order of magnitude less common than those that decelerate the process. One of the most extreme outliers in this data set, an arginine-to-phenylalanine substitution at position 48 (R48F) of chymotrypsin inhibitor 2 (CI2), accelerates the protein's folding rate by a factor of 36 relative to that of the wild-type protein and is the most accelerating substitution reported to date in a two-state protein. In order to better understand the origins of this anomalous behavior, we have characterized the kinetics of multiple additional substitutions at this position. We find that substitutions at position 48 in CI2 fall into two distinct classes. The first, comprising residues that ablate the charge of the wild-type arginine but retain the hydrophobicity of its alkane chain, accelerate folding by at least 10-fold. The second class, comprising all other residues, produces folding rates within a factor of two of the wild-type rate. A significant positive correlation between hydrophobicity and folding rate across all of the residues we have characterized at this position suggests that the hydrophobic methylene units of the wild-type arginine play a significant role in stabilizing the folding transition state. Likewise, studies of the pH dependence of the histidine substitution indicate a strong correlation between folding rate and charge state. Thus, mutations that ablate the arginine's positive charge while retaining the hydrophobic contacts of its methylene units tend to dramatically accelerate folding. Previous studies have suggested that arginine 48 plays an important functional role in CI2, which may explain why it is highly conserved despite the anomalously large deceleration it produces in the folding of this protein.

Cavity-creating mutations in Pseudomonas aeruginosa azurin: effects on protein dynamics and stability

Changes in flexibility and structural stability of Pseudomonas aeruginosa azurin in response to cavity-creating mutations were probed by the phosphorescence emission of Trp-48, which was deeply buried in the compact hydrophobic core of the macromolecule, and by measurements of guanidinum hydrochloride unfolding, respectively. Replacement of the bulky side chains Phe-110, Phe-29, and Tyr-108 with the smaller Ala introduced cavities at different distances from the hydrophobic core. The phosphorescence lifetime (tau(0)) of Trp-48, buried inside the protein core, and the acrylamide quenching rate constant (k(q)) were used to monitor local and global flexibility changes induced by the introduction of the cavity. The results of this work demonstrate the following: 1), the effect on core flexibility of the insertion of cavities is not correlated readily to the distance of the cavity from the core; 2), the protein global flexibility results are related to the cavity distance from the packed core of the macromolecule; and 3), the increase in protein flexibility does not correspond necessarily to a comparable destabilizing effect of some mutations.

Reduction potential tuning of the blue copper center in Pseudomonas aeruginosa azurin by the axial methionine as probed by unnatural amino acids

The conserved axial ligand methionine 121 from Pseudomonas aeruginosa azurin (Az) has been replaced by isostructural unnatural amino acid analogues, oxomethionine (OxM), difluoromethionine (DFM), trifluoromethionine (TFM), selenomethionine (SeM), and norleucine (Nle) using expressed protein ligation. The replacements resulted in < 6 nm shifts in the S(Cys)-Cu charge transfer (CT) band in the electronic absorption spectra and < 8 gauss changes in the copper hyperfine coupling constants (AII) in the X-band electron paramagnetic resonance spectra, suggesting that isostructural replacement of Met resulted in minimal structural perturbation of the copper center. The slight blue shifts of the CT band follow the trend of stronger electronegativity of the ligands. This trend is supported by 19F NMR studies of the fluorinated methionine analogues. However, the order of AII differs, suggesting additional factors influencing AII. In contrast to the small changes in the UV-vis and EPR spectra, a large variation of > 227 mV in reduction potential was observed for the series of variants reported here. Additionally, a linear correlation was established between the reduction potentials and hydrophobicity of the variants. Extension of this analysis to other type 1 copper-containing proteins reveals a linear correlation between change in hydrophobicity and change in reduction potential, independent of the protein scaffold, experimental conditions, measurement techniques, and steric modifications. This analysis has also revealed for the first time high and low potential states for type 1 centers, and the difference may be attributable to destabilization of the protein fold by disruption of hydrophobic or hydrogen bonding interactions that stabilize the type 1 center.

Correlation between Protein Stability Cores and Protein Folding Kinetics: A Case Study on Pseudomonas aeruginosa Apo-Azurin

This paper reports a combined computational and experimental study of the correlation between protein stability cores and folding kinetics. An empirical potential function was developed, and it was used for analyzing interaction energies among secondary structure elements. Studies on a beta sandwich protein, Pseudomonas aeruginosa azurin, showed that the computationally identified substructure with the strongest interactions in the native state is identical to the "interlocked pair" of beta strands, an invariant motif found in most sandwich-like proteins. Moreover, previous and new in vitro folding results revealed that the identified substructure harbors most residues that form native-like interactions in the folding transition state. These observations demonstrate that the potential function is effective in revealing the relative strength of interactions among various protein parts; they also strengthen the suggestion that the most stable regions in native proteins favor stable interactions early during folding.

Phi-value analysis of apo-azurin folding: comparison between experiment and theory

Pseudomonas aeruginosa azurin is a 128-residue beta-sandwich metalloprotein; in vitro kinetic experiments have shown that it folds in a two-state reaction. Here, we used a variational free energy functional to calculate the characteristics of the transition state ensemble (TSE) for folding of the apo-form of P. aeruginosa azurin and investigate how it responds to thermal and mutational changes. The variational method directly yields predicted chevron plots for wild-type and mutant apo-forms of azurin. In parallel, we performed in vitro kinetic-folding experiments on the same set of azurin variants using chemical perturbation. Like the wild-type protein, all apo-variants fold in apparent two-state reactions both in calculations and in stopped-flow mixing experiments. Comparisons of phi (phi) values determined from the experimental and theoretical chevron parameters reveal an excellent agreement for most positions, indicating a polarized, highly structured TSE for folding of P. aeruginosa apo-azurin. We also demonstrate that careful analysis of side-chain interactions is necessary for appropriate theoretical description of core mutants.

Solvation of the folding-transition state in Pseudomonas aeruginosa azurin is modulated by metal: Solvation of azurin's folding nucleus

The role of water in protein folding, specifically its presence or not in the transition-state structure, is an unsolved question. There are two common classes of folding-transition states: diffuse transition states, in which almost all side chains have similar, rather low phi (phi) values, and polarized transition states, which instead display distinct substructures with very high phi-values. Apo-and zinc-forms of Pseudomonas aeruginosa azurin both fold in two-state equilibrium and kinetic reactions; while the apo-form exhibits a polarized transition state, the zinc form entails a diffuse, moving transition state. To examine the presence of water in these two types of folding-transition states, we probed the equilibrium and kinetic consequences of replacing core valines with isosteric threonines at six positions in azurin. In contrast to regular hydrophobic-to-alanine phi-value analysis, valine-to-threonine mutations do not disrupt the core packing but stabilize the unfolded state and can be used to assess the degree of solvation in the folding-transition state upon combination with regular phi-values. We find that the transition state for folding of apo-azurin appears completely dry, while that for zinc-azurin involves partially formed interactions that engage water molecules. This distinct difference between the apo-and holo-folding nuclei can be rationalized in terms of the shape of the free-energy barrier.

Unique complex between bacterial azurin and tumor-suppressor protein p53

The tumor-suppressor protein p53 is a major player in regulation of cell growth, genomic stability, and cell death. Recent work suggests that Pseudomonas aeruginosa azurin, as the only bacterial protein known to date, can enter cancer cells and interact with p53 promoting cell death. For the first time, here we demonstrate and characterize this proposed complex using purified proteins in vitro. We find that azurin binds to p53 with nanomolar affinity in a four-to-one stoichiometry (pH 7.5, 25 degrees C). Upon azurin binding, secondary structure is induced and tryptophan fluorescence is quenched, implying that interactions occur in the N-terminal p53 domain which is also the binding site for many oncogenes. Further biophysical studies may assist the design of novel cancer treatments that are based on azurin.

Snapshots of a Dynamic Folding Nucleus in Zinc-Substituted Pseudomonas aeruginosa Azurin

Zinc-substituted Pseudomonas aeruginosa azurin folds in two-state equilibrium and kinetic reactions. In the unfolded state, the zinc ion remains bound to the unfolded polypeptide via two native-state ligands (His117 and Cys112). The significantly curved Chevron plot for zinc-substituted azurin was earlier ascribed to movement of the folding-transition state. At low concentrations of denaturant, the transition state occurs early in the folding reaction (low Tanford beta-value), whereas at high-denaturant concentration, it moves closer to the native structure (high Tanford beta-value). Here, we use this movement to track the formation and growth of zinc-substituted azurin's folding nucleus with atomic resolution using protein engineering. The average phi (phi) value for 17 positions (covering all secondary-structure elements) goes from 0.25 in 0 M GuHCl (beta approximately 0.46) to 0.76 in 4 M GuHCl (beta approximately 0.86); a phi-value of 1 or 0 indicates native-like or unfolded-like interactions, respectively. Analysis of individual phi-values reveals a delocalized nucleus where structure condenses around a leading density centered on Leu50 in the core. The diffuse moving transition state for zinc-substituted azurin is in sharp contrast to the fixed polarized folding nucleus observed for apo-azurin. The dramatic difference in apparent kinetic behavior for the two forms of azurin can be rationalized as a minor alteration on a common free-energy profile that exhibits a broad activation barrier.

Role of structural determinants in folding of the sandwich-like protein Pseudomonas aeruginosa azurin

An invariant substructure that forms two interlocked pairs of neighboring beta-strands occurs in essentially all known sandwich-like proteins. Eight conserved positions in these strands were recently shown to act as structural determinants. To test whether the residues at these invariant positions are conserved for mechanistic (i.e., part of folding nucleus) or energetic (i.e., governing native-state stability) reasons, we characterized the folding behavior of eight point-mutated variants of the sandwich-like protein Pseudomonas aeruginosa apo-azurin. We find a simple relationship among the conserved positions: half of the residues form native-like interactions in the folding transition state, whereas the others do not participate in the folding nucleus but govern high native-state stability. Thus, evolutionary preservation of these specific positions gives both mechanistic and energetic advantages to members of the sandwich-like protein family.

Protein folding: defining a "standard" set of experimental conditions and a preliminary kinetic data set of two-state proteins

Recent years have seen the publication of both empirical and theoretical relationships predicting the rates with which proteins fold. Our ability to test and refine these relationships has been limited, however, by a variety of difficulties associated with the comparison of folding and unfolding rates, thermodynamics, and structure across diverse sets of proteins. These difficulties include the wide, potentially confounding range of experimental conditions and methods employed to date and the difficulty of obtaining correct and complete sequence and structural details for the characterized constructs. The lack of a single approach to data analysis and error estimation, or even of a common set of units and reporting standards, further hinders comparative studies of folding. In an effort to overcome these problems, we define here a "consensus" set of experimental conditions (25 degrees C at pH 7.0, 50 mM buffer), data analysis methods, and data reporting standards that we hope will provide a benchmark for experimental studies. We take the first step in this initiative by describing the folding kinetics of 30 apparently two-state proteins or protein domains under the consensus conditions. The goal of our efforts is to set uniform standards for the experimental community and to initiate an accumulating, self-consistent data set that will aid ongoing efforts to understand the folding process.

Role of Cofactors in Folding of the Blue-Copper Protein Azurin

Many proteins in living cells coordinate cofactors, such as metal ions, to attain their activity. Since the cofactors in such cases often can interact with their corresponding unfolded polypeptides in vitro, it is important to unravel how cofactors modulate protein folding. In this review, I will discuss the role of cofactors in folding of the blue-copper protein Pseudomonas aeruginosa azurin. In the case of both copper (Cu(II) and Cu(I)) and zinc (Zn(II)), the metal can bind to unfolded azurin. The residues involved in copper (Cu(II) and Cu(I)) coordination in the unfolded state have been identified as Cys112, His117, and Met121. The affinities of Cu(II), Cu(I), and Zn(II) are all higher for the folded than for the unfolded azurin polypeptide, resulting in metal stabilization of the native state as compared to the stability of apo-azurin. Cu(II), Zn(II), and several apo forms of azurin all fold in two-state kinetic reactions with roughly identical polypeptide-folding speeds. This suggests that the native-state beta-barrel topology, not cofactor interactions or thermodynamic stability, determines azurin's folding barrier. Nonetheless, copper binds much more rapidly (i.e., 4 orders of magnitude) to unfolded azurin than to folded azurin. Therefore, the fastest route to functional azurin is through copper binding before polypeptide folding; this sequence of events may be the relevant biological pathway.

Effects of cavity-forming mutations on the internal dynamics of azurin

The effects of two single-point cavity-forming mutations, F110S and I7S, on the internal dynamics of azurin from Pseudomonas aeruginosa were probed by the phosphorescence emission of Trp-48, deeply buried in the compact hydrophobic core of the macromolecule. Changes in flexibility of the protein matrix around the chromophore were monitored by the intrinsic phosphorescence lifetime (tau(0)) whereas more general effects on structural fluctuations were deduced from the phosphorescence acrylamide quenching rate constant (k(q)), which measures the diffusion of the solute through the protein fold. The results show a spectacular, 4-5 orders of magnitude, increase of k(q) emphasizing that large amplitude structural fluctuations permitting acrylamide migration to the protein core have been drastically enhanced in each azurin mutant. The large, 12-15 kcal/mol, decrease in the activation enthalpy associated to k(q) suggests that the rate enhancement is caused, rather than through a generalized increase of protein flexibility, by the elimination of an inner barrier to the diffusion process. According to tau(0) the chromophore environment is more fluid with I7S but strikingly more rigid with F110S, demonstrating that when internal cavities are formed local effects on the mobility at the mutation site are unpredictable. Both tau(0) and k(q) reveal a structure tightening role of bound Cd(2+) that correlates with the increase in stability from apo- to holo-azurin. While these alterations in internal dynamics of azurin do not seem to play a role on electron transfer through the central region, the enhanced migration of acrylamide emphasizes that cavities may be critical for the rapid diffusion of substrates to buried, solvent inaccessible sites of enzymes.

Methionine-121 coordination determines metal specificity in unfolded Pseudomonas aeruginosa azurin

Pseudomonas aeruginosa azurin binds copper so tightly that it remains bound even upon polypeptide unfolding. Copper can be substituted with zinc without change in protein structure, and also in this complex the metal remains bound upon protein unfolding. Previous work has shown that native-state copper ligands Cys112 and His117 are two of at least three metal ligands in the unfolded state. In this study we use isothermal titration calorimetry and spectroscopic methods to test if the native-state ligand Met121 remains a metal ligand upon unfolding. From studies on a point-mutated version of azurin (Met121Ala) and a set of model peptides spanning the copper-binding C-terminal part (including Cys112, His117 and Met121), we conclude that Met121 is a metal ligand in unfolded copper-azurin but not in the case of unfolded zinc-azurin. Combination of unfolding and metal-titration data allow for determination of copper (Cu(II) and Cu(I)) and zinc affinities for folded and unfolded azurin polypeptides, respectively.

Approaching the speed limit for Greek Key β-barrel formation: transition-state movement tunes folding rate of zinc-substituted azurin

Azurin is a blue-copper protein with a beta-barrel structure of Greek Key topology. In vitro, copper can be substituted with zinc without change in protein structure. We here analyze the kinetic folding behavior of zinc-substituted Pseudomonas aeruginosa azurin. Our findings can be summarized in three key conclusions: first, zinc remains strongly bound to the polypeptide upon unfolding, suggesting that the cofactor may bind to the protein before polypeptide folding in vivo. Second, the semi-logarithmic plot of folding and unfolding rates for zinc-substituted azurin as a function of denaturant concentration exhibits curvature due to a changing transition-state structure. Third, the extrapolated folding speed in water for zinc-substituted azurin is similar to that of other proteins with the same topology, implying that there is a speed limit that can be modulated by stability-driven transition-state movement for formation of beta-barrel structures with Greek Key topology.

If space is provided, bulky modification on the rim of azurin's beta-barrel results in folded protein

Pseudomonas aeruginosa azurin is a blue-copper protein with a beta-barrel fold. Here we report that, at conditions where thermal unfolding of apo-azurin is reversible, the reaction occurs in a single step with a transition midpoint (T(m)) of 69 degrees C (pH 7). The active-site mutation His117Gly creates a cavity in the beta-barrel near the surface but does not perturb the overall fold (T(m) of 64 degrees C, pH 7). Oxidation of the active-site cysteine (Cysteine-112) in wild-type azurin, which occurs readily at higher temperatures, results in a modified protein that cannot adopt a native-like structure. In sharp contrast, Cysteine-112 oxidation in His117Gly azurin yields a modified apo-azurin that appears folded and displays cooperative, reversible unfolding (T(m) approximately 55 degrees C, pH 7). We conclude that azurin's beta-barrel is a rigid structural element that constrains the structure of its surface; a bulky modification can only be accommodated if complementary space is provided.