Disulfide formation and stability of a cysteine-rich repeat protein from Helicobacter pylori.

Helicobacter pylori cysteine-rich proteins (Hcps) are disulfide-containing repeat proteins. The repeating unit is a 36-residue, disulfide-bridged, helix-loop-helix motif. We use the protein HcpB, which has four repeats and four disulfide bridges arrayed in tandem, as a model to determine the thermodynamic stability of a disulfide-rich repeat protein and to study the formation and the contribution to stability of the disulfide bonds. When the disulfide bonds are intact, the chemical unfolding of HcpB at pH 5 is cooperative and can be described by a two-state reaction. Thermal unfolding is reversible between pH 2 and 5 and irreversible at higher pH 5. Differential scanning calorimetry shows noncooperative structural changes preceding the main thermal unfolding transition. Unfolding of the oxidized protein is not an all-or-none two-state process, and the disulfide bonds prevent complete unfolding of the polypeptide chain. The reduced protein is significantly less stable and does not unfold in a cooperative way. During oxidative refolding of the fully reduced protein, all the possible disulfide intermediates with a correct disulfide bond are formed. Formation of "wrong" (non-native) disulfide bonds could not be demonstrated, indicating that the reduced protein already has some partial repeating structure. There is a major folding intermediate with disulfides in the second, third, and fourth repeat and reduced cysteines in the first repeat. Disulfide formation in the first repeat limits the overall rate of oxidative refolding and contributes about half of the thermodynamic stability to native HcpB, estimated as 27 kJ mol(-1) at 25 degrees C and pH 7. The high contribution to stability of the first repeat may be explained by the repeat acting as a cap to protect the hydrophobic interior of the molecule.

Enthalpic and entropic effects of salt and polyol osmolytes on site-specific protein-DNA association: the integrase Tn916-DNA complex.

The effect of low molecular-weight compounds on the equilibrium constant K(A) can be used to explore the energetics and molecular mechanism of protein-DNA interactions. Here we use the complex composed of the integrase Tn916 DNA-binding domain and its target DNA duplex to investigate the effects of salt and the nonionic osmolytes glycerol and sorbitol on sequence-specific protein-DNA association. Increasing Na(+) concentration from 0.12 to 0.32 M weakens the binding affinity by a factor of 20. The decrease of affinity is dominated by a large loss of binding enthalpy but only a small loss of binding entropy. This contrasts the concept that the salt-induced weakening of protein-DNA binding is mainly entropic. The large enthalpy loss is discussed in the light of recent views about the nature of the general salt effect. Addition of up to 2.5 M sorbitol and up to 3.3 M glycerol causes a slight increase of the binding affinity. However, both osmolytes lead to a large enthalpy gain and a similarly large entropy loss. This intriguing enthalpy-entropy compensation can be explained in part by an enthalpic chelate effect: The osmolyte tightens the structure of the protein-DNA complex whereby the formation of enthalpically favorable noncovalent interactions is promoted at the entropic cost of a more rigid complex. The results were obtained by isothermal titration calorimetry. They are supported by kinetic experiments showing that the rate of formation of the complex is reduced by salt, but the rate of complex dissociation is not. Glycerol and sorbitol reduce both rates in line with an only small effect on complex stability. This work clarifies the thermodynamic and kinetic response of a novel protein-DNA complex to increased salt and the presence of two common, nonionic osmolytes.

Folding of a designed simple ankyrin repeat protein.

Ankyrin repeats (AR) are 33-residue motifs containing a beta-turn, followed by two alpha-helices connected by a loop. AR occur in tandem arrangements and stack side-by-side to form elongated domains involved in very different cellular tasks. Recently, consensus libraries of AR repeats were constructed. Protein E1_5 represents a member of the shortest library, and consists of only a single consensus repeat flanked by designed N- and C-terminal capping repeats. Here we present a biophysical characterization of this AR domain. The protein is compactly folded, as judged from the heat capacity of the native state and from the specific unfolding enthalpy and entropy. From spectroscopic data, thermal and urea-induced unfolding can be modeled by a two-state transition. However, scanning calorimetry experiments reveal a deviation from the two-state behavior at elevated temperatures. Folding and unfolding at 5 degrees C both follow monoexponential kinetics with k(folding) = 28 sec(-1) and k(unfolding) = 0.9 sec(-1). Kinetic and equilibrium unfolding parameters at 5 degrees C agree very well. We conclude that E1_5 folds in a simple two-state manner at low temperatures while equilibrium intermediates become populated at higher temperatures. A chevron-plot analysis indicates that the protein traverses a very compact transition state along the folding/unfolding pathway. This work demonstrates that a designed minimal ankyrin repeat protein has the thermodynamic and kinetic properties of a compactly folded protein, and explains the favorable properties of the consensus framework.

Inverse electrostatic effect: electrostatic repulsion in the unfolded state stabilizes a leucine zipper.

The pH-dependent stability of a protein is strongly affected by electrostatic interactions between ionizable residues in the folded as well as unfolded state. Here we characterize the individual contributions of charged Glu and His residues to stability and determine the NMR structure of the designed, heterodimeric leucine zipper AB consisting of an acidic A chain and a basic B chain. Thermodynamic parameters are compared with those of the homologous leucine zipper AB(SS) in which the A and B chains are disulfide-linked. NMR structures of AB based on (1)H NMR data collected at 600 MHz converge, and formation of the same six interchain salt bridges found previously in disulfide-linked AB(SS) [Marti, D. N., and Bosshard, H. R. (2003) J. Mol. Biol. 330, 621-637] is indicated. While the structures of AB and AB(SS) are very similar, their pH-dependent relative stabilities are strikingly different. The stability of AB peaks at pH approximately 4.5 and is higher at pH 8 than at pH 2. In contrast, AB(SS) is most stable at acidic pH where no interhelical salt bridges are formed. The different energetic contributions of charged Glu and His residues to stability of the two coiled coil structures were evaluated from pK(a) shifts induced by folding. The six charged Glu residues involved in salt bridges stabilize leucine zipper AB by 4.5 kJ/mol yet destabilize disulfide-linked AB(SS) by -1.1 kJ/mol. Two non-ion-paired Glu charges destabilize AB by only -1.8 kJ/mol but AB(SS) by -5.6 kJ/mol. The higher relative stability of AB at neutral pH is not caused by more favorable electrostatic interactions in the folded leucine zipper. It is due mainly to unfavorable electrostatic interactions in the unfolded A and B chains and may therefore be called an inverse electrostatic effect. This study illustrates the importance of residual interactions in the unfolded state and how the energetics of the unfolded state affect the stability of the folded protein.

Protein stabilization by salt bridges: concepts, experimental approaches and clarification of some misunderstandings.

Salt bridges in proteins are bonds between oppositely charged residues that are sufficiently close to each other to experience electrostatic attraction. They contribute to protein structure and to the specificity of interaction of proteins with other biomolecules, but in doing so they need not necessarily increase a protein's free energy of unfolding. The net electrostatic free energy of a salt bridge can be partitioned into three components: charge-charge interactions, interactions of charges with permanent dipoles, and desolvation of charges. Energetically favorable Coulombic charge-charge interaction is opposed by often unfavorable desolvation of interacting charges. As a consequence, salt bridges may destabilize the structure of the folded protein. There are two ways to estimate the free energy contribution of salt bridges by experiment: the pK(a) approach and the mutation approach. In the pK(a) approach, the contribution of charges to the free energy of unfolding of a protein is obtained from the change of pK(a) of ionizable groups caused by altered electrostatic interactions upon folding of the protein. The pK(a) approach provides the relative free energy gained or lost when ionizable groups are being charged. In the mutation approach, the coupling free energy between interacting charges is obtained from a double mutant cycle. The coupling free energy is an indirect and approximate measure of the free energy of charge-charge interaction. Neither the pK(a) approach nor the mutation approach can provide the net free energy of a salt bridge. Currently, this is obtained only by computational methods which, however, are often prone to large uncertainties due to simplifying assumptions and insufficient structural information on which calculations are based. This state of affairs makes the precise thermodynamic quantification of salt bridge energies very difficult. This review is focused on concepts and on the assessment of experimental methods and does not cover the vast literature.

Monomeric and dimeric bZIP transcription factor GCN4 bind at the same rate to their target DNA site.

Basic leucine zipper (bZIP) transcription factors are dimeric proteins that recognize dyadic and mostly palindromic DNA sites. Dimerization of bZIP transcription factor GCN4 is linked to the folding of its C-terminal leucine zipper domain. However, monomeric GCN4, lacking a folded leucine zipper, also recognizes the DNA site with dimerization taking place on the DNA. Here we report the kinetics of DNA recognition by unfolded monomeric and folded dimeric derivatives of GCN4 using a 19 bp double-stranded DNA containing a palindromic CRE site. The rate of DNA binding of both monomeric and dimeric GCN4 has a bimolecular rate constant of 3-5 x 10(8) M(-1) s(-1), which is near the diffusion limit. Because the rate of dimerization of GCN4 is slower (1.7 x 10(7) M(-1) s(-1)) than the rate of DNA association, the formation of the dimeric GCN4-DNA complex through consecutive binding of two monomers (monomer pathway) is faster when starting from free monomers. Thus, the results presented here support facilitated and rapid target recognition by the monomeric transcription factor. However, DNA binding of preformed folded dimeric GCN4 is as rapid as complex formation through the monomer pathway. Therefore, the monomer and dimer pathways are kinetically equivalent if monomeric and dimeric GCN4 are at equilibrium. Hence, the dimer pathway may also have a role under in vivo conditions. The lower affinity of GCN4 in which two DNA contacting residues have been mutated is due exclusively to the faster dissociation of the mutant protein-DNA complex and not to slower complex formation.

Fast folding of the HIV-1 and SIV gp41 six-helix bundles.

Human (HIV-1) and simian (SIV) immunodeficiency virus fusion with the host cell is promoted by the receptor-triggered refolding of the gp41 envelope protein into a stable trimer-of-hairpins structure that brings viral and cellular membranes into close proximity. The core of this hairpin structure is a six-helix bundle in which an inner homotrimeric coiled coil is buttressed by three antiparallel outer HR2 helices. We have used stopped-flow circular dichroism spectroscopy to characterize the unfolding and refolding kinetics of the six-helix bundle using the HIV-1 and SIV N34(L6)C28 polypeptides. In each case, the time-course of ellipticity changes in refolding experiments is well described by a simple two-state model involving the native trimer and the unfolded monomers. The unfolding free energy of the HIV-1 and SIV trimers and their urea dependence calculated from kinetic data are in very good agreement with data measured directly by isothermal unfolding experiments. Thus, formation of the gp41 six-helix bundle structure involves no detectable population of stable, partly folded intermediates. Folding of HIV-1 N34(L6)C28 is five orders of magnitudes faster than folding of its SIV counterpart in aqueous buffer: k(on),(HIV-1)=1.3 x 10(15)M(-2)s(-1) versus k(on),(SIV)=1.1 x 10(10)M(-2)s(-1). The unfolding rates are similar: k(off),(HIV-1)=1.1 x 10(-5)s(-1) versus k(off),(SIV=)5.7 x 10(-4)s(-1). Kinetic m-values indicate that the transition state for folding of the HIV-1 protein is significantly more compact than the transition state of the SIV protein. Replacement of a single SIV threonine by isoleucine corresponding to position 573 in the HIV-1 sequence significantly stabilizes the protein and renders the folding rate close to that of the HIV-1 protein yet without making the transition state of the mutant as compact as that of the HIV-1 protein. Therefore, the overall reduction of surface exposure in the high-energy transition state seems not to account for different folding rates. While the available biological evidence suggests that refolding of the gp41 protein is slow, our study implies that structural elements outside the trimer-of-hairpins limit the rate of HIV-1 fusion kinetics.

Electrostatic interactions in leucine zippers: thermodynamic analysis of the contributions of Glu and His residues and the effect of mutating salt bridges.

Electrostatic interactions play a complex role in stabilizing proteins. Here, we present a rigorous thermodynamic analysis of the contribution of individual Glu and His residues to the relative pH-dependent stability of the designed disulfide-linked leucine zipper AB(SS). The contribution of an ionized side-chain to the pH-dependent stability is related to the shift of the pK(a) induced by folding of the coiled coil structure. pK(a)(F) values of ten Glu and two His side-chains in folded AB(SS) and the corresponding pK(a)(U) values in unfolded peptides with partial sequences of AB(SS) were determined by 1H NMR spectroscopy: of four Glu residues not involved in ion pairing, two are destabilizing (-5.6 kJ mol(-1)) and two are interacting with the positive alpha-helix dipoles and are thus stabilizing (+3.8 kJ mol(-1)) in charged form. The two His residues positioned in the C-terminal moiety of AB(SS) interact with the negative alpha-helix dipoles resulting in net stabilization of the coiled coil conformation carrying charged His (-2.6 kJ mol(-1)). Of the six Glu residues involved in inter-helical salt bridges, three are destabilizing and three are stabilizing in charged form, the net contribution of salt-bridged Glu side-chains being destabilizing (-1.1 kJ mol(-1)). The sum of the individual contributions of protonated Glu and His to the higher stability of AB(SS) at acidic pH (-5.4 kJ mol(-1)) agrees with the difference in stability determined by thermal unfolding at pH 8 and pH 2 (-5.3 kJ mol(-1)). To confirm salt bridge formation, the positive charge of the basic partner residue of one stabilizing and one destabilizing Glu was removed by isosteric mutations (Lys-->norleucine, Arg-->norvaline). Both mutations destabilize the coiled coil conformation at neutral pH and increase the pK(a) of the formerly ion-paired Glu side-chain, verifying the formation of a salt bridge even in the case where a charged side-chain is destabilizing. Because removing charges by a double mutation cycle mainly discloses the immediate charge-charge effect, mutational analysis tends to overestimate the overall energetic contribution of salt bridges to protein stability.

Energetics of sequence-specific protein-DNA association: binding of integrase Tn916 to its target DNA.

The DNA binding domain of the transposon Tn916 integrase (INT-DBD) binds to its DNA target site by positioning the face of a three-stranded antiparallel beta-sheet within the major groove. Binding of INT-DBD to a 13 base pair duplex DNA target site was studied by isothermal titration calorimetry, differential scanning calorimetry, thermal melting followed by circular dichroism spectroscopy, and fluorescence spectroscopy. The observed heat capacity change accompanying the association reaction (DeltaC(p)) is temperature-dependent, decreasing from -1.4 kJ K(-1) mol(-1) at 4 degrees C to -2.9 kJ K(-1) mol(-1) at 30 degrees C. The reason is that the partial molar heat capacities of the free protein, the free DNA duplex, and the protein-DNA complex are not changing in parallel when the temperature increases and that thermal motions of the protein and the DNA are restricted in the complex. After correction for this effect, DeltaC(p) is -1.8 kJ K(-1) mol(-1) and temperature-independent. However, this value is still higher than DeltaC(p) of -1.2 kJ K(-1) mol(-1) estimated by semiempirical methods from dehydration of surface area buried at the complex interface. We propose that the discrepancy between the measured and the structure-based prediction of binding energetics is caused by incomplete dehydration of polar groups in the complex. In support, we identify cavities at the interface that are large enough to accommodate approximately 10 water molecules. Our results highlight the difficulties of structure-based prediction of DeltaC(p) (and other thermodynamic parameters) and emphasize how important it is to consider changes of thermal motions and soft vibrational modi in protein-DNA association reactions. This requires not only a detailed investigation of the energetics of the complex but also of the folding thermodynamics of the protein and the DNA alone, which are described in the accompanying paper [Milev et al. (2003) Biochemistry 42, 3492-3502].

Energetics of sequence-specific protein-DNA association: conformational stability of the DNA binding domain of integrase Tn916 and its cognate DNA duplex.

Sequence-specific DNA recognition by bacterial integrase Tn916 involves structural rearrangements of both the protein and the DNA duplex. Energetic contributions from changes of conformation, thermal motions and soft vibrational modi of the protein, the DNA, and the complex significantly influence the energetic profile of protein-DNA association. Understanding the energetics of such a complicated system requires not only a detailed calorimetric investigation of the association reaction but also of the components in isolation. Here we report on the conformational stability of the integrase Tn916 DNA binding domain and its cognate 13 base pair target DNA duplex. Using a combination of temperature and denaturant induced unfolding experiments, we find that the 74-residue DNA binding domain is compact and unfolds cooperatively with only small deviation from two-state behavior. Scanning calorimetry reveals an increase of the heat capacity of the native protein attributable to increased thermal fluctuations. From the combined calorimetric and spectroscopic experiments, the parameters of protein unfolding are T(m) = 43.8 +/- 0.3 degrees C, DeltaH(m) = 255 +/- 18 kJ mol(-1), DeltaS(m) = 0.80 +/- 0.06 kJ mol(-1), and DeltaC(p) = 5.0 +/- 0.8 kJ K(-1) mol(-1). The DNA target duplex displays a thermodynamic signature typical of short oligonucleotide duplexes: significant heat absorption due to end fraying and twisting precedes cooperative unfolding and dissociation. The parameters for DNA unfolding and dissociation are DeltaH(m) = 335 +/- 4 kJ mol(-1) and DeltaC(p) = 2.7 +/- 0.9 kJ K(-(1) mol(-1). The results reported here have been instrumental in interpreting the thermodynamic features of the association reaction of the integrase with its 13 base pair target DNA duplex reported in the accompanying paper [Milev et al. (2003) Biochemistry 42, 3481-3491].

Salt bridges destabilize a leucine zipper designed for maximized ion pairing between helices.

Interhelical salt bridges are common in leucine zippers and are thought to stabilize the coiled coil conformation. Here we present a detailed thermodynamic investigation of the designed, disulfide-linked leucine zipper AB(SS) whose high-resolution NMR structure shows six interhelical ion pairs between heptad positions g of one helix and e' of the other helix but no ion pairing within single helices. The average pK(a) value of the Glu side chain carboxyl groups of AB(SS) is slightly higher than the pK(a) of a freely accessible Glu in an unfolded peptide [Marti, D. N., Jelesarov, I., and Bosshard, H. R. (2000) Biochemistry 39, 12804-12818]. This indicates that the salt bridges are destabilizing, a prediction we now have confirmed by determining the pH +/- stability profile of AB(SS). Circular dichroism-monitored unfolding by urea and by heating and differential scanning calorimetry show that the coiled coil conformation is approximately 5 kJ/mol more stable when salt bridges are broken by protonation of the carboxyl side chains. Using guanidinium chloride as the denaturant, the increase in the free energy of unfolding on protonation of the carboxyl side chains is larger, approximately 17 kJ/mol. The discrepancy between urea and guanidinium chloride unfolding can be ascribed to the ionic nature of guanidinium chloride, which screens charge-charge interactions. This work demonstrates the difficulty of predicting the energetic contribution of salt bridges from structural data alone even in a case where the ion pairs are seen in high-resolution NMR structures. The reason is that the contribution to stability results from a fine balance between energetically favorable Coulombic attractions and unfavorable desolvation of charges and conformational constraints of the residues involved in ion pairing. The apparent discrepancy between the results presented here and mutational studies indicating stabilization by salt bridges is discussed and resolved. An explanation is proposed for why interhelical salt bridges are frequently found in natural coiled coils despite evidence that they do not directly contribute to stability.

Calculation of protein ionization equilibria with conformational sampling: pK(a) of a model leucine zipper, GCN4 and barnase.

The use of conformational ensembles provided by nuclear magnetic resonance (NMR) experiments or generated by molecular dynamics (MD) simulations has been regarded as a useful approach to account for protein motions in the context of pK(a) calculations, yet the idea has been tested occasionally. This is the first report of systematic comparison of pK(a) estimates computed from long multiple MD simulations and NMR ensembles. As model systems, a synthetic leucine zipper, the naturally occurring coiled coil GCN4, and barnase were used. A variety of conformational averaging and titration curve-averaging techniques, or combination thereof, was adopted and/or modified to investigate the effect of extensive global conformational sampling on the accuracy of pK(a) calculations. Clustering of coordinates is proposed as an approach to reduce the vast diversity of MD ensembles to a few structures representative of the average electrostatic properties of the system in solution. Remarkable improvement of the accuracy of pK(a) predictions was achieved by the use of multiple MD simulations. By using multiple trajectories the absolute error in pK(a) predictions for the model leucine zipper was reduced to as low as approximately 0.25 pK(a) units. The validity, advantages, and limitations of explicit conformational sampling by MD, compared with the use of an average structure and a high internal protein dielectric value as means to improve the accuracy of pK(a) calculations, are discussed.