RESUMEN
Urea and guanidine-hydrochloride (GdnHCl) are frequently used for protein denaturation in order to determine the Gibbs free energy of folding and kinetic folding/unfolding parameters. Constant pH value is applied in the folding/unfolding experiments at different denaturant concentrations and steady protonation state of titratable groups is assumed in the folded and unfolded protein, respectively. The apparent side-chain pKa values of Asp, Glu, His and Lys in the absence and presence of 6 M urea and GdnHCl, respectively, have been determined by 1H-NMR. pKa values of all four residues are up-shifted by 0.3-0.5 pH units in presence of 6 M urea by comparison with pKa values of the residues dissolved in water. In the presence of 6 M GdnHCl, pKa values are down-shifted by 0.2-0.3 pH units in the case of acidic and up-shifted by 0.3-0.5 pH units in the case of basic residues. Shifted pKa values in the presence of denaturant may have a pronounced effect on the outcome of the protein stability obtained from denaturant unfolding experiments.
Asunto(s)
Proteínas/química , Ácido Aspártico/química , Fenómenos Químicos , Química Física , Ácido Glutámico/química , Guanidina/química , Histidina/química , Concentración de Iones de Hidrógeno , Lisina/química , Espectroscopía de Resonancia Magnética/métodos , Pliegue de Proteína , Sensibilidad y Especificidad , Volumetría/métodos , Urea/químicaRESUMEN
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.
Asunto(s)
Leucina Zippers , Pliegue de Proteína , Termodinámica , Secuencia de Aminoácidos , Arginina/química , Dicroismo Circular , Cristalografía por Rayos X , Dimerización , Disulfuros/química , Ácido Glutámico/química , Calor , Concentración de Iones de Hidrógeno , Datos de Secuencia Molecular , Resonancia Magnética Nuclear Biomolecular/métodos , Fragmentos de Péptidos/química , Sales (Química)/química , Electricidad EstáticaRESUMEN
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.
Asunto(s)
Proteínas/química , Proteínas/metabolismo , Electricidad Estática , Animales , Conformación Proteica , Desnaturalización Proteica , Pliegue de Proteína , TermodinámicaRESUMEN
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.
Asunto(s)
Proteína gp41 de Envoltorio del VIH/química , VIH-1/química , Glicoproteínas de Membrana/química , Pliegue de Proteína , Estructura Secundaria de Proteína , Proteínas de los Retroviridae/química , Animales , Dicroismo Circular , Proteína gp41 de Envoltorio del VIH/metabolismo , Humanos , Glicoproteínas de Membrana/metabolismo , Modelos Moleculares , Desnaturalización Proteica , Proteínas de los Retroviridae/metabolismo , Virus de la Inmunodeficiencia de los Simios/químicaRESUMEN
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.
Asunto(s)
Ácido Glutámico/química , Histidina/química , Leucina Zippers/fisiología , Electricidad Estática , Termodinámica , Secuencia de Aminoácidos , Arginina/genética , Proteínas de Unión al ADN/química , Concentración de Iones de Hidrógeno , Leucina Zippers/genética , Lisina/genética , Espectroscopía de Resonancia Magnética , Modelos Moleculares , Datos de Secuencia Molecular , Mutagénesis Sitio-Dirigida , Norleucina/genética , Péptidos/química , Pliegue de Proteína , Proteínas Quinasas/química , Proteínas de Saccharomyces cerevisiae/química , Sales (Química)/químicaRESUMEN
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.
Asunto(s)
Leucina Zippers , Estructura Secundaria de Proteína , Secuencia de Aminoácidos , Ácido Glutámico/química , Concentración de Iones de Hidrógeno , Iones/química , Leucina Zippers/fisiología , Datos de Secuencia Molecular , Péptidos/química , Desnaturalización Proteica , Pliegue de Proteína , Sales (Química)/química , TemperaturaRESUMEN
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.