Thermodynamic principles for the engineering of pH-driven conformational switches and acid insensitive proteins

The general thermodynamic principles behind pH driven conformational transitions of biological macromolecules are well understood. What is less obvious is how they can be used to engineer pH switches in proteins. The acid unfolding of staphylococcal nuclease (SNase) was used to illustrate different factors that can affect pH-driven conformational transitions. Acid unfolding is a structural transition driven by preferential H⁺ binding to the acid unfolded state (U) over the native (N) state of a protein. It is the result of carboxylic groups that titrate with more normal pK_a values in the U state than in the N state. Acid unfolding profiles of proteins reflect a balance between electrostatic and non-electrostatic contributions to stability. Several strategies were used in attempts to turn SNase into an acid insensitive protein: (1) enhancing global stability of the protein with mutagenesis or with osmolytes, (2) use of high salt concentrations to screen Coulomb interactions, (3) stabilizing the N state through specific anion effects, (4) removing Asp or Glu residues that titrate with depressed pK_a values in the N state, and (5) removing basic residues that might have strong repulsive interactions in the N state at low pH. The only effective way to engineer acid resistance in SNase is not through modulation of pK_a values of Asp/Glu but by enhancing the global stability of the protein. Modulation of pH-driven conformational transitions by selective manipulation of the electrostatic component of the switch is an extremely difficult undertaking.     

pH dependence of conformational fluctuations of the protein backbone

Proton binding equilibria (pK_a values) of ionizable groups in proteins are exquisitely sensitive to their microenvironments. Apparent pK_a values measured for individual ionizable residues with NMR spectroscopy are actually population‐weighted averages of the pK_a in different conformational microstates. NMR spectroscopy experiments with staphylococcal nuclease were used to test the hypothesis that pK_a values of surface Glu and Asp residues are affected by pH‐sensitive fluctuations of the backbone between folded and locally unfolded conformations. ¹⁵N spin relaxation studies showed that as the pH decreases from the neutral into the acidic range the amplitudes of backbone fluctuations in the ps‐ns timescale increase near carboxylic residues. Hydrogen exchange experiments suggested that backbone conformational fluctuations promoted by decreasing pH also reflect slower local or sub‐global unfolding near carboxylic groups. This study has implications for structure‐based pK_a calculations: (1) The timescale of the backbone's response to ionization events in proteins can range from ps to ms, and even longer; (2) pH‐sensitive fluctuations of the backbone can be localized to both the segment the ionizable residue is attached to or the one that occludes the ionizable group; (3) Structural perturbations are not necessarily propagated through Coulomb interactions; instead, local fluctuations appear to be coupled through the co‐operativity inherent to elements of secondary structure and to networks of hydrogen bonds. These results are consistent with the idea that local conformational fluctuations and stabilities are important determinants of apparent pK_a values of ionizable residues in proteins. Proteins 2014; 82:3132–3143. © 2014 Wiley Periodicals, Inc.     

Effect of a buried ion pair in the hydrophobic core of a protein: An insight from constant pH molecular dynamics study

Constant pH molecular dynamics (CpHMD) is a commonly used sampling method, which incorporates the coupling of conformational flexibility and protonation state of a protein during the simulation by using pH as an external parameter. The effects on the structure and stability of a hyperstable variant of staphylococcal nuclease (Δ+PHS) protein of an artificial charge pair buried in its hydrophobic core are investigated by applying both CpHMD and accelerated molecular dynamics coupled with constant pH (CpHaMD) methods. Generalized Born electrostatics is used to model the solvent water. Two sets of starting coordinates of V23E/L36K variant of Δ+PHS, namely, Maestro generated coordinates from Δ+PHS and crystal structure coordinates of the same are considered for detail investigations. On the basis of root mean square displacement (RMSD) and root mean square fluctuations (RMSF) calculations, it is observed that this variant is stable over a wide range of pH. The calculated pKa values for aspartate and glutamate residues based on both CpHMD and CpHaMD simulations are consistent with the reported experimental values (within ± 0.5 to ± 1.5 pH unit), which clearly indicates that the local chemical environment of the carboxylic acids in V23E/L36K variant are comparable to the parent form. The strong salt bridge interaction between the mutated pair, E23/K36 and additional hydrogen bonds formed in the V23E/L36K variant, may help to compensate for the unfavorable self‐energy experienced by the burial of these residues in the hydrophobic core. However, from RMSD, RMSF, and pKa analysis, no significant change in the global conformation of V23E/L36K variant with respect to the parent form, Δ+PHS is noticed. © 2014 Wiley Periodicals, Inc. Biopolymers 103: 148–157, 2015.     

Model for the differential stabilities of rhinovirus and poliovirus to mild acidic pH, based on electrostatics calculations.

Previous calculations of electrostatic interactions in the rhinovirus capsid have identified a subset of histidine residues, paired with lysine or arginine, that may be involved in pH-induced conformational changes related to viral uncoating. Further calculations with the finite difference method, accounting for the dielectric environment of the ionizable groups, suggest that charge burial in the crystal conformation will prevent protonation of these histidine residues in the pentamer-pentamer interface. Calculations with a modelled pentamer-pentamer interface in which three beta-strands are removed recover mildly acidic pKa values for the histidines. These results are discussed in the context of the structural interactions of these three beta-strands, which form a beta-sheet extension from the rest of the capsid, and with regard to the conformation of the homologous beta-sheet extension in poliovirus, which also possesses homologous histidine-lysine/arginine pairs. A model is developed in which the structural stability of the beta-sheet extension is related to the difference in acid stability of rhinovirus and poliovirus. It is suggested that, for poliovirus prior to cell receptor binding, the beta-sheet extension is stable at pH 3, the pentamer-pentamer interface histidines remain buried, and the virus is acid-stable. Cell receptor binding of poliovirus destabilizes the beta-sheet extension and the acid lability that is proposed to result could be involved in viral uncoating. For rhinovirus it is suggested that the observed conformational change in the absence of cell receptor binding involves a further acidic pH-activated process or conformational fluctuations that rearrange the beta-sheet extension and expose the pentamer-pentamer interface histidine residues to the acidic medium. Sequence analysis and electrostatics calculations reveal an aspartic acid in the beta-sheet extension that may have different pKa values in rhinovirus and poliovirus.     

The pH-induced release of iron from transferrin investigated with a continuum electrostatic model.

A reduction in pH induces the release of iron from transferrin in a process that involves a conformational change in the protein from a closed to an open form. Experimental evidence suggests that there must be changes in the protonation states of certain, as yet not clearly identified, residues in the protein accompanying this conformational change. Such changes in protonation states of residues and the consequent changes in electrostatic interactions are assumed to play a large part in the mechanism of release of iron from transferrin. Using the x-ray crystal structures of human ferri- and apo-lactoferrin, we calculated the pKa values of the titratable residues in both the closed (iron-loaded) and open (iron-free) conformations with a continuum electrostatic model. With the knowledge of a residue's pKa value, its most probable protonation state at any specified pH may be determined. The preliminary results presented here are in good agreement with the experimental observation that the binding of ferric iron and the synergistic anion bicarbonate/carbonate results in the release of approximately three H+ ions. It is suggested that the release of these three H+ ions may be accounted for, in most part, by the deprotonation of the bicarbonate and residues Tyr-92, Lys-243, Lys-282, and Lys-285 together with the protonation of residues Asp-217 and Lys-277.     

Characterizing a partially folded intermediate of the villin headpiece domain under non-denaturing conditions: contribution of His41 to the pH-dependent stability of the N-terminal subdomain

The contribution of interactions involving the imidazole ring of His41 to the pH-dependent stability of the villin headpiece (HP67) N-terminal subdomain has been investigated by nuclear magnetic resonance (NMR) spin relaxation. NMR-derived backbone N-H order parameters (S2) for wild-type (WT) HP67 and H41Y HP67 indicate that reduced conformational flexibility of the N-terminal subdomain in WT HP67 is due to intramolecular interactions with the His41 imidazole ring. These interactions, together with desolvation effects, contribute to significantly depress the pKa of the buried imidazole ring in the native state. 15N R1rho relaxation dispersion data indicate that WT HP67 populates a partially folded intermediate state that is 10.9 kJ mol(-1) higher in free energy than the native state under non-denaturing conditions at neutral pH. The partially folded intermediate is characterized as having an unfolded N-terminal subdomain while the C-terminal subdomain retains a native-like fold. Although the majority of the residues in the N-terminal subdomain sample a random-coil distribution of conformations, deviations of backbone amide 1H and 15N chemical shifts from canonical random-coil values for residues within 5A of the His41 imidazole ring indicate that a significant degree of residual structure is maintained in the partially folded ensemble. The pH-dependence of exchange broadening is consistent with a linear three-state exchange model whereby unfolding of the N-terminal subdomain is coupled to titration of His41 in the partially folded intermediate with a pKa,I=5.69+/-0.07. Although maintenance of residual interactions with the imidazole ring in the unfolded N-terminal subdomain appears to reduce pKa,I compared to model histidine compounds, protonation of His41 disrupts these interactions and reduces the difference in free energy between the native state and partially folded intermediate under acidic conditions. In addition, chemical shift changes for residues Lys70-Phe76 in the C-terminal subdomain suggest that the HP67 actin binding site is disrupted upon unfolding of the N-terminal subdomain, providing a potential mechanism for regulating the villin-dependent bundling of actin filaments.     

Acid-induced equilibrium folding intermediate of human platelet profilin

The acid-induced unfolding of human platelet profilin (HPP) can be minimally modeled as a three-state process. Equilibrium unfolding studies have been performed on human platelet profilin1 (HPP) and monitored by far-UV circular dichroism, tryptophan fluorescence, ANS binding, and NMR spectroscopy. Far-UV CD measurements obtained by acid titration demonstrate that HPP unfolds via a three-state mechanism (N --> I --> U), with a highly populated intermediate between pH 4 and 5. Approximately 80% of native helical secondary structural content remains at pH 4, as indicated by monitoring the CD signal at 222 nm. The stability (DeltaGH2O) of the native conformation at pH 7.0 (obtained by monitoring the change in tryptophan signal as a function of urea concentration) is 5.56 +/- 0.51 kcal mol-1; however, the DeltaGH2O for the intermediate species at pH 4 is 2.01 +/- 0.47 kcal mol-1. The calculated m-values for the pH 7.0 and pH 4.0 species were 1.64 +/- 0.15 and 1.34 +/- 0.17 kcal mol-1 M-1, respectively, which is an indication that the native and intermediate species are similarly compact. Additionally, translational diffusion measurements obtained by NMR spectroscopy and ANS binding studies are consistent with a globular and compact conformation at both pH 7.0 and 4.0. The pKa values for the two histidine (His) residues located on helix 4 of HPP were determined to be 5.6 and 5.7 pH units. These pKa values coincide with the midpoint of the far-UV CD acid titration curve and suggest that the protonation of one or both His residues may play a role in the formation of the unfolding intermediate. Stable intermediate species populate the 2D 1H-15N HSQC NMR spectra between pH 4 and 5. A number of backbone and side-chain resonances show significant perturbations relative to the native spectrum; however, considerable nativelike tertiary contacts remain. Interestingly, the residues on HPP that are significantly altered at low pH coincide with segments of the G-actin binding surface and poly-l-proline binding interface. The earlier reports that a decrease in pH below 6.0 induces structural alterations in profilin, favoring dissociation of the profilin-actin complex, corresponds with the structural alterations observed in the partially unfolded species. Our findings suggest that a novel mechanism for pH induced disruption of the profilin-G-actin complex involve a nativelike unfolding intermediate of profilin.     

pH-induced denaturation of proteins: a single salt bridge contributes 3-5 kcal/mol to the free energy of folding of T4 lysozyme

The energetics of a salt bridge formed between the side chains of aspartic acid 70 (Asp70) and histidine 31 (His31) of T4 lysozyme have been examined by nuclear magnetic resonance techniques. The pKa values of the residues in the native state are perturbed from their values in the unfolded protein such that His31 has a pKa value of 9.1 in the native state and 6.8 in the unfolded state at 10 degrees C in moderate salt. Similarly, the aspartate pKa is shifted to a value of about 0.5 in the native state from its value of 3.5-4.0 in the unfolded state. These shifts in pKa show that the salt bridge is stabilized 3-5 kcal/mol. This implies that the salt bridge stabilizes the native state by 3-5 kcal/mol as compared to the unfolded state. This is reflected in the thermodynamic stability of mutants of the protein in which Asp70, His31, or both are replaced by asparagine. These observations and consideration of the thermodynamic coupling of protonation state to folding of proteins suggest a mechanism of acid denaturation in which the unfolded state is progressively stabilized by protonation of its acid residues as pH is lowered below pH 4. The unfolded state is stabilized only if acidic groups in the folded state have lower pKa values than in the unfolded state. When the pH is sufficiently low, the acid groups of both the native and unfolded states are fully protonated, and the apparent unfolding equilibrium constant becomes pH independent. Similar arguments apply to base-induced unfolding.(ABSTRACT TRUNCATED AT 250 WORDS)     

pH-dependent conformational states of horse liver alcohol dehydrogenase.

The quenching of liver alcohol dehydrogenase protein fluorescence at alkaline pH indicates two conformational states of the enzyme with a pKa of 9.8+/-0.2, shifted to 10.6+/-0.2 in D2O. NAD+ and 2-p-toluidinonaphthalene-6-sulfonate, a fluorescent probe competitive with coenzyme, bind to the acid conformation of the enzyme. The pKa of the protein-fluorescence quenching curve is shifted toward 7.6 in the presence of NAD+, and the ternary complex formation with NAD+ and trifluoroethanol results in a pH-independent maximal quench. At pH (pD) 10.5, the rate constant for NAD+ binding was 2.6 times faster in D2O2 than in H2O due to the shift of the pKa. Based on these results, a scheme has been proposed in which the state of protonation of an enzyme functional group with a pKa of 9.8 controls the conformational state of the enzyme. NAD+ binds to the acid conformation and subsequently causes another conformational change resulting in the perturbation of the pKa to 7.6. Alcohol then binds to the unprotonated form of the functional group with a pKa of 7.6 in the binary enzyme-NAD+ complex and converts the enzyme to the alkaline conformation. Thus, at neutral pH liver alcohol dehydrogenase undergoes two conformational changes en route to the ternary complex in which hydride transfer occurs.     

Linkage of proton binding to the thermal unfolding of Sso7d from the hyperthermophilic archaebacterium Sulfolobus solfataricus.

In this study the pH dependence of the thermal stability of Sso7d from Sulfolobus solfataricus is analyzed. This small globular protein of 63 residues shows a very marked dependence of thermal stability on pH: the denaturation temperature passes from 65.2 degrees C at pH 2.5 to 97.9 degrees C at pH 4.5. Analysis of the data points out that the binding of at least two protons is coupled to the thermal unfolding. By linking the proton binding to the conformational unfolding equilibrium, a thermodynamic model, which is able to describe the dependence upon the solution pH of both the excess heat capacity function and the denaturation Gibbs energy change for Sso7d, is developed. The decreased stability in very acid conditions is due to the binding of two protons on identical and noninteracting sites of the unfolded state. Actually, such sites are two carboxyl groups possessing very low pKa values in the native structure, probably involved in salt-bridges on the protein surface.     

Two steps in the transition between the native and acid states of bovine α-lactalbumin detected by circular polarization of luminescence: Evidence for a premolten globule state?

A few studies indirectly support the existence of an intermediate in the transition of Ca(2+)-saturated bovine alpha-lactalbumin (alpha-LA) from the native (N) to the acidic (A) state, known as the molten globule state. However, direct experimental evidence for the appearance of this intermediate has not been obtained. The signal of circular polarization of luminescence (CPL) is sensitive to fine conformational transitions because of its susceptibility to changes in the environmental asymmetry of fluorescent chromophores in their excited electronic states. In the present study, CPL measurements were applied using the intrinsic tryptophan fluorescence of alpha-LA as well as the fluorescence of 8-anilino-1-naphthalenesulfonic acid (ANS) bound to alpha-LA. CPL of tryptophan and ANS was measured in the pH range of 2.5-6 in order to find direct experimental evidence for the proposed intermediate. CPL (characterized by the emission anisotropy factor, g(em)) depends on the asymmetry of the protein molecular structure in the environment of the tryptophan and the ANS chromophores in the excited electronic state. The pH dependence of both the gab, absorption anisotropy factor determined by CD, and the ANS steady state fluorescence, showed a single transition at pH 3-3.7 as already reported elsewhere. This transition was interpreted as being a result of a change of the alpha-LA tertiary structure, which resulted in a loss of asymmetry of the environment of both the tryptophan residues and the ANS hydrophobic binding sites. The pH dependence of the tryptophan and ANS g(em) showed an additional conformational transition at pH 4-5, which coincided with the pKa of Ca2+ dissociation (pKa 5), as predicted by Permyakov et al. (1981, Biochem Biophys Res Commun 100:191-197). The titration curve showed that there is a pH range between 3.7 and 4.1 in which alpha-LA exists in an intermediate state between the N- and A-state. We suggest that the intermediate is the premolten globule state characterized by a reduced Ca2+ binding to the alpha-LA, native-like tertiary structure, and reduced asymmetric fluctuation of the tertiary structure on the nanosecond time scale. This intermediate resembles the "critical activated state" theoretically deduced by Kuwajima et al. (1989, J Mol Biol 206:547-561). The present study demonstrates the power of CPL measurements for the investigation of folding/unfolding transitions in proteins.     

Effect of a specific hydrogen bond (N138ND2-Q106O) on conformational integrity, stability, and activity of staphylococcal nuclease

There are two hydrogen bonding interactions (N138ND2-Q106O and Y54OH-S141OG) between the C-terminal region and the main body of staphylococcal nuclease (SNase). To examine the role of these hydrogen bonds, SNase(141) and its three mutants, SNase(141)N138D, SNase(141)S141A, and SNase(141)N138D/S141A, were created. The N138D mutation has the N138ND2-Q106O interaction deleted and the S141A mutation has the Y54OH-S141OG and S141OG-N138O interactions deleted. The conformational features, stability, and activity of the proteins have been compared by using circular dichroism, intrinsic and ANS-binding fluorescence, GdnHCl-induced denaturation, and activity assay. The results clearly show that the N138D mutation significantly alters the secondary and tertiary structures of the protein, producing a partially unfolding state; in contrast, the S141A mutation has no such effect on structure. These results strongly suggest that the specific hydrogen bond, N138ND2-Q106O, plays an important role in maintaining the conformational integrity and stability of the nuclease.     

Apparent pKa shifts of titratable residues at high denaturant concentration and the impact on protein stability

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.     

pH dependent thermodynamic and amide exchange studies of the C-terminal domain of the ribosomal protein L9: implications for unfolded state structure

It is now recognized that unfolded states of globular proteins are not random coils but instead can contain significant amounts of residual structure. Here, we combine amide H/D exchange studies and thermodynamic measurements to probe pH dependent structure in the unfolded state of the small, mixed alpha-beta protein CTL9. The m value measured by urea denaturation is strongly dependent upon pD, increasing by 40% from pD 7.5 to 4.85. Likewise, the change in heat capacity upon unfolding, deltaCp(o), increases significantly from pD 7.5 to 5.5. These studies argue that the unfolded state contains interactions, presumably hydrophobic in nature, that lead to a more compact state at high pH. The expansion at lower pH correlates with the estimated unfolded state pKa values of the three histidines in CTL9 with additional contributions from acid side chains at the lower pH. Amide H/D exchange studies were conducted at pD 5.0, 6.0, and 7.0. At pD 5.0, the exchange rates could be measured for 44 residues, 29 of which exchanged by global unfolding. No evidence was found for any super protected sites, that is, sites that exchange at rates slower than those expected for global exchange. The estimated precision for the experiments limits detection to residues that are protected 2.3-fold above the intrinsic exchange rate. Thirty-seven residues could be followed at pD 6 and 27 residues at pD 7. Again no evidence for a significant super protected structure was observed. The properties of CTL9(11) are compared to other structured denatured states.     

Determinants of conductance of a bacterial voltage-gated sodium channel

Through molecular dynamics (MD) and free energy simulations in electric fields, we examine the factors influencing conductance of bacterial voltage-gated sodium channel Na_vMs. The channel utilizes four glutamic acid residues in the selectivity filter (SF). Previously, we have shown, through constant pH and free energy calculations of pKa values, that fully deprotonated, singly protonated, and doubly protonated states are all feasible at physiological pH, depending on how many ions are bound in the SF. With 173 MD simulations of 450 or 500 ns and additional free energy simulations, we determine that the conductance is highest for the deprotonated state and decreases with each additional proton bound. We also determine that the pKa value of the four glutamic residues for the transition between deprotonated and singly protonated states is close to the physiological pH and that there is a small voltage dependence. The pKa value and conductance trends are in agreement with experimental work on bacterial Na_v channels, which show a decrease in maximal conductance with lowering of pH, with pKa in the physiological range. We examine binding sites for Na⁺ in the SF, compare with previous work, and note a dependence on starting structures. We find that narrowing of the gate backbone to values lower than the crystal structure's backbone radius reduces the conductance, whereas increasing the gate radius further does not affect the conductance. Simulations with some amount of negatively charged lipids as opposed to purely neutral lipids increases the conductance, as do simulations at higher voltages.     

Asymmetric protonation of EmrE

The small multidrug resistance transporter EmrE is a homodimer that uses energy provided by the proton motive force to drive the efflux of drug substrates. The pKa values of its “active-site” residues—glutamate 14 (Glu14) from each subunit—must be poised around physiological pH values to efficiently couple proton import to drug export in vivo. To assess the protonation of EmrE, pH titrations were conducted with ¹H-¹⁵N TROSY-HSQC nuclear magnetic resonance (NMR) spectra. Analysis of these spectra indicates that the Glu14 residues have asymmetric pKa values of 7.0 ± 0.1 and 8.2 ± 0.3 at 45°C and 6.8 ± 0.1 and 8.5 ± 0.2 at 25°C. These pKa values are substantially increased compared with typical pKa values for solvent-exposed glutamates but are within the range of published Glu14 pKa values inferred from the pH dependence of substrate binding and transport assays. The active-site mutant, E14D-EmrE, has pKa values below the physiological pH range, consistent with its impaired transport activity. The NMR spectra demonstrate that the protonation states of the active-site Glu14 residues determine both the global structure and the rate of conformational exchange between inward- and outward-facing EmrE. Thus, the pKa values of the asymmetric active-site Glu14 residues are key for proper coupling of proton import to multidrug efflux. However, the results raise new questions regarding the coupling mechanism because they show that EmrE exists in a mixture of protonation states near neutral pH and can interconvert between inward- and outward-facing forms in multiple different protonation states.     

Predicting extreme pKa shifts in staphylococcal nuclease mutants with constant pH molecular dynamics

Accurate computational methods of determining protein and nucleic acid pK(a) values are vital to understanding pH-dependent processes in biological systems. In this article, we use the recently developed method constant pH molecular dynamics (CPHMD) to explore the calculation of highly perturbed pK(a) values in variants of staphylococcal nuclease (SNase). Simulations were performed using the replica exchange (REX) protocol for improved conformational sampling with eight temperature windows, and yielded converged proton populations in a total sampling time of 4 ns. Our REX-CPHMD simulations resulted in calculated pK(a) values with an average unsigned error (AUE) of 0.75 pK units for the acidic residues in Δ + PHS, a hyperstable variant of SNase. For highly pK(a)-perturbed SNase mutants with known crystal structures, our calculations yielded an AUE of 1.5 pK units and for those mutants based on modeled structures an AUE of 1.4 pK units was found. Although a systematic underestimate of pK shifts was observed in most of the cases for the highly perturbed pK mutants, correlations between conformational rearrangement and plasticity associated with the mutation and error in pK(a) prediction was not evident in the data. This study further extends the scope of electrostatic environments explored using the REX-CPHMD methodology and suggests that it is a reliable tool for rapidly characterizing ionizable amino acids within proteins even when modeled structures are employed.     

pH dependence of the hydrogen exchange in the SH3 domain of alpha-spectrin

Using nuclear magnetic resonance we have measured the hydrogen exchange (HX) in the Src homology region 3 (SH3) domain of alpha-spectrin as a function of pH*. At very acidic pH* values the exchange of most residues appears to occur via global unfolding, although several residues show abnormally large Gibbs energies of exchange, suggesting the presence of some residual structure in the unfolded state. At higher pH* HX occurs mainly via local or partial unfoldings. We have been able to characterize the coupling between the electrostatic interactions in this domain and the conformational fluctuations occurring under native conditions by analyzing the dependence upon pH* of the Gibbs energy of exchange. The SH3 domain seems to be composed of a central core, which requires large structural disruptions to become exposed to the solvent, surrounded by smaller subdomains, which fluctuate independently.     

13C]Methylated ribonuclease A. 13C NMR studies of the interaction of lysine 41 with active site ligands

A unique resonance in the 13C NMR spectrum of [13C]methylated ribonuclease A has been assigned to a N epsilon, N-dimethylated active site residue, lysine 41. The chemical shift of this resonance was studied over the pH range 3 to 11, and the titration curve showed two inflection points, at pH 5.7 and 9.0. The higher pKa, designated pKa1, was assigned to the ionization of the lysyl residue itself while the pKa of 5.7, designated pKa2, was assigned on the basis of its pKa to the ionization of a histidyl residue which is somehow coupled to lysine 41. Both pKa values are measurably perturbed by the binding of active site ligands including nucleotides, nucleosides, phosphate, and sulfate. In most cases, the alterations in pKa values induced by the ligands were larger for pKa2. The ligand-induced perturbations in pKa2 generally paralleled those reported for histidine 12, another active site residue (Griffin, J. H., Schechter, A. N., and Cohen, J. S. (1973) Ann. N. Y. Acad. Sci. 222, 693-708). The sensitivity of the N epsilon, N-dimethylated lysine 41 resonance to the histidyl ionization may result from a conformational change in the active site region of ribonuclease which is coupled to the histidyl ionization. This coupling between lysine 41 and another ribonuclease residue, which has not been documented previously, offers new insight into the interrelationship between residues in the active site of this well characterized enzyme.     

Conformation and hydrogen ion titration of proteins: a continuum electrostatic model with conformational flexibility.

A new method for including local conformational flexibility in calculations of the hydrogen ion titration of proteins using macroscopic electrostatic models is presented. Intrinsic pKa values and electrostatic interactions between titrating sites are calculated from an ensemble of conformers in which the positions of titrating side chains are systematically varied. The method is applied to the Asp, Glu, and Tyr residues of hen lysozyme. The effects of different minimization and/or sampling protocols for both single-conformer and multi-conformer calculations are studied. For single-conformer calculations it is found that the results are sensitive to the choice of all-hydrogen versus polar-hydrogen-only atomic models and to the minimization protocol chosen. The best overall agreement of single-conformer calculations with experiment is obtained with an all-hydrogen model and either a two-step minimization process or minimization using a high dielectric constant. Multi-conformational calculations give significantly improved agreement with experiment, slightly smaller shifts between model compound pKa values and calculated intrinsic pKa values, and reduced sensitivity of the intrinsic pKa calculations to the initial details of the structure compared to single-conformer calculations. The extent of these improvements depends on the type of minimization used during the generation of conformers, with more extensive minimization giving greater improvements. The ordering of the titrations of the active-site residues, Glu-35 and Asp-52, is particularly sensitive to the minimization and sampling protocols used. The balance of strong site-site interactions in the active site suggests a need for including site-site conformational correlations.