Carboxyl pKa Values and Acid Denaturation of BBL

The protein BBL undergoes structural transitions and acid denaturation between pH 1.2 and 8.0. Using NMR spectroscopy, we measured the pK_a values of all the carboxylic residues in this pH range. We employed ¹³C direct-detection two-dimensional IPAP (in-phase antiphase) CACO NMR spectroscopy to monitor the ionization state of different carboxylic groups and demonstrated its advantages over other NMR techniques in measuring pK_a values of carboxylic residues. The two residues Glu161 and Asp162 had significantly lowered pK_a values, showing that these residues are involved in a network of stabilizing electrostatic interactions, as is His166. The other carboxylates had unperturbed values. The pH dependence of the free energy of denaturation was described quantitatively by the ionizations of those three residues of perturbed pK_a, and, using thermodynamic cycles, we could calculate their pK_as in the native and denatured states as well as the equilibrium constants for denaturation of the different protonation states. We also measured ¹³C^α chemical shifts of individual residues as a function of pH. These shifts sense structural transitions rather than ionizations, and they titrated with pH consistent with the change in equilibrium constant for denaturation. Kinetic measurements of the folding of BBL E161Q indicated that, at pH 7, the stabilizing interactions with Glu161 are formed mainly in the transition state. We also found that local interactions still exist in the acid-denatured state of BBL, which attenuate somewhat the flexibility of the acid-denatured state.     

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.     

The pKa Values of Acidic and Basic Residues Buried at the Same Internal Location in a Protein Are Governed by Different Factors

The pK_a values of internal ionizable groups are usually very different from the normal pK_a values of ionizable groups in water. To examine the molecular determinants of pK_a values of internal groups, we compared the properties of Lys, Asp, and Glu at internal position 38 in staphylococcal nuclease. Lys38 titrates with a normal or elevated pK_a, whereas Asp38 and Glu38 titrate with elevated pK_a values of 7.0 and 7.2, respectively. In the structure of the L38K variant, the buried amino group of the Lys38 side chain makes an ion pair with Glu122, whereas in the structure of the L38E variant, the buried carboxyl group of Glu38 interacts with two backbone amides and has several nearby carboxyl oxygen atoms. Previously, we showed that the pK_a of Lys38 is normal owing to structural reorganization and water penetration concomitant with ionization of the Lys side chain. In contrast, the pK_a values of Asp38 and Glu38 are perturbed significantly owing to an imbalance between favorable polar interactions and unfavorable contributions from dehydration and from Coulomb interactions with surface carboxylic groups. Their ionization is also coupled to subtle structural reorganization. These results illustrate the complex interplay between local polarity, Coulomb interactions, and structural reorganization as determinants of pK_a values of internal groups in proteins. This study suggests that improvements to computational methods for pK_a calculations will require explicit treatment of the conformational reorganization that can occur when internal groups ionize.     

Electrostatic effects in a network of polar and ionizable groups in staphylococcal nuclease

His121 and His124 are embedded in a network of polar and ionizable groups on the surface of staphylococcal nuclease. To examine how membership in a network affects the electrostatic properties of ionizable groups, the tautomeric state and the pK_a values of these histidines were measured with NMR spectroscopy in the wild-type nuclease and in 13 variants designed to disrupt the network. In the background protein, His121 and His124 titrate with pK_a values of 5.2 and 5.6, respectively. In the variants, where the network was disrupted, the pK_a values range from 4.03 to 6.46 for His121, and 5.04 to 5.99 for His124. The largest decrease in a pK_a was observed when the favorable Coulomb interaction between His121 and Glu75 was eliminated; the largest increase was observed when Tyr91 or Tyr93 was substituted with Ala or Phe. In all variants, the dominant tautomeric state at neutral pH was the N^ε2 state. At one level the network behaves as a rigid unit that does not readily reorganize when disrupted: crystal structures of the E75A or E75Q variants show that even when the pivotal Glu75 is removed, the overall configuration of the network was unaffected. On the other hand, a few key hydrogen bonds appear to govern the conformation of the network, and when these bonds are disrupted the network reorganizes. Coulomb interactions within the network report an effective dielectric constant of 20, whereas a dielectric constant of 80 is more consistent with the magnitude of medium to long-range Coulomb interactions in this protein. The data demonstrate that when structures are treated as static, rigid bodies, structure-based pK_a calculations with continuum electrostatics method are not useful to treat ionizable groups in cases where pK_a values are governed by short-range polar and Coulomb interactions.     

Uncovering Specific Electrostatic Interactions in the Denatured States of Proteins

The stability and folding of proteins are modulated by energetically significant interactions in the denatured state that is in equilibrium with the native state. These interactions remain largely invisible to current experimental techniques, however, due to the sparse population and conformational heterogeneity of the denatured-state ensemble under folding conditions. Molecular dynamics simulations using physics-based force fields can in principle offer atomistic details of the denatured state. However, practical applications are plagued with the lack of rigorous means to validate microscopic information and deficiencies in force fields and solvent models. This study presents a method based on coupled titration and molecular dynamics sampling of the denatured state starting from the extended sequence under native conditions. The resulting denatured-state pK_as allow for the prediction of experimental observables such as pH- and mutation-induced stability changes. I show the capability and use of the method by investigating the electrostatic interactions in the denatured states of wild-type and K12M mutant of NTL9 protein. This study shows that the major errors in electrostatics can be identified by validating the titration properties of the fragment peptides derived from the sequence of the intact protein. Consistent with experimental evidence, our simulations show a significantly depressed pK_a for Asp⁸ in the denatured state of wild-type, which is due to a nonnative interaction between Asp⁸ and Lys¹². Interestingly, the simulation also shows a nonnative interaction between Asp⁸ and Glu⁴⁸ in the denatured state of the mutant. I believe the presented method is general and can be applied to extract and validate microscopic electrostatics of the entire folding energy landscape.     

Catalytic role of the calcium ion in GH97 inverting glycoside hydrolase

The role of calcium ion in the active site of the inverting glycoside hydrolase family 97 enzyme, BtGH97a, was investigated through structural and kinetic studies. The calcium ion was likely directly involved in the catalytic reaction. The pH dependence of k_cat/K_m values in the presence or absence of calcium ion indicated that the calcium ion lowered the pK_a of the base catalyst. The significant decreases in k_cat/K_m for hydrolysis of substrates with basic leaving groups in the absence of calcium ion confirmed that the calcium ion facilitated the leaving group departure.     

Intrinsically disordered inhibitor of glutamine synthetase is a functional protein with random‐coil‐like pKa values

The sequential action of glutamine synthetase (GS) and glutamate synthase (GOGAT) in cyanobacteria allows the incorporation of ammonium into carbon skeletons. In the cyanobacterium Synechocystis sp. PCC 6803, the activity of GS is modulated by the interaction with proteins, which include a 65‐residue‐long intrinsically disordered protein (IDP), the inactivating factor IF7. This interaction is regulated by the presence of charged residues in both IF7 and GS. To understand how charged amino acids can affect the binding of an IDP with its target and to provide clues on electrostatic interactions in disordered states of proteins, we measured the pK_a values of all IF7 acidic groups (Glu32, Glu36, Glu38, Asp40, Asp58, and Ser65, the backbone C‐terminus) at 100 mM NaCl concentration, by using NMR spectroscopy. We also obtained solution structures of IF7 through molecular dynamics simulation, validated them on the basis of previous experiments, and used them to obtain theoretical estimates of the pK_a values. Titration values for the two Asp and three Glu residues of IF7 were similar to those reported for random‐coil models, suggesting the lack of electrostatic interactions around these residues. Furthermore, our results suggest the presence of helical structure at the N‐terminus of the protein and of conformational changes at acidic pH values. The overall experimental and in silico findings suggest that local interactions and conformational equilibria do not play a role in determining the electrostatic features of the acidic residues of IF7.     

Electrostatic analysis of the hepatitis C virus NS3 helicase reveals both active and allosteric site locations

Multi-conformation continuum electrostatics (MCCE) was used to analyze various structures of the NS3 RNA helicase from the hepatitis C virus in order to determine the ionization state of amino acid side chains and their pK_as. In MCCE analyses of HCV helicase structures that lacked ligands, several active site residues were identified to have perturbed pK_as in both the nucleic acid binding site and in the distant ATP-binding site, which regulates helicase movement. In all HCV helicase structures, Glu493 was unusually basic and His369 was abnormally acidic. Both these residues are part of the HCV helicase nucleic acid binding site, and their roles were analyzed by examining the pH profiles of site-directed mutants. Data support the accuracy of MCCE predicted pK_a values, and reveal that Glu493 is critical for low pH enzyme activation. Several key residues, which were previously shown to be involved in helicase-catalyzed ATP hydrolysis, were also identified to have perturbed pK_as including Lys210 in the Walker-A motif and the DExD/H-box motif residues Asp290 and His293. When DNA was present in the structure, the calculated pK_as shifted for both Lys210 and Asp290, demonstrating how DNA binding might lead to electrostatic changes that stimulate ATP hydrolysis.     

On the pH-dependence of the light-induced hydrogen ion gradient in spinach chloroplasts

The dependence of the light-induced H⁺ gradient in chloroplasts (ΔpH) on external pH was examined using the distribution of aniline, an amine of low pK_a. ΔpH was essentially independent of pH over the range of 7–8. It was previously reported that ΔpH, determined from the distribution of relatively polar amines of high pK_a, decreased as the pH was lowered below 8. It is suggested that, in the case of amines of high pK_a, ΔpH values determined at low external pH values are too low because the permeability of chloroplasts to the amine cation relative to that of the unprotonated form may be significant.     

Minimum-Energy Path for a U6 RNA Conformational Change Involving Protonation, Base-Pair Rearrangement and Base Flipping

The U6 RNA internal stem-loop (U6 ISL) is a highly conserved domain of the spliceosome that is important for pre-mRNA splicing. The U6 ISL contains an internal loop that is in equilibrium between two conformations controlled by the protonation state of an adenine (pK_a = 6.5). Lower pH favors formation of a protonated C-A⁺ wobble pair and base flipping of the adjacent uracil. Higher pH favors stacking of the uracil and allows an essential metal ion to bind at this position. Here, we define the minimal-energy path for this conformational transition. To do this, we solved the U6 ISL structure at higher pH (8.0) in order to eliminate interference from the low-pH conformer. This structure reveals disruption of the protonated C-A⁺ pair and formation of a new C-U pair, which explains the preference for a stacked uracil at higher pH. Next, we used nudged elastic band molecular dynamics simulations to calculate the minimum-energy path between the two conformations. Our results indicate that the C-U pair is dynamic, which allows formation of the more stable C-A⁺ pair upon adenine protonation. After formation of the C-A⁺ pair, the unpaired uracil follows a minor-groove base-flipping pathway. Molecular dynamics simulations suggest that the extrahelical uracil is stabilized by contacts with the adjacent helix.     

Mutation of non-conserved amino acids surrounding catalytic site to shift pH optimum of Bacillus circulans xylanase

Improvement of enzyme function by engineering pH dependence of enzymatic activity is of importance for industrial application of Bacillus circulans xylanases. Target mutation sites were selected by structural alignment between B. circulans xylanase and other xylanases having different pH optima. We selected non-conserved mutant sites within 8 Å from the catalytic residues, to see whether these residues have some role in modulating pK_as of the catalytic residues. We hypothesized that the non-conserved residues which may not have any role in enzyme catalysis might perturb pK_as of the catalytic residues. Change in pK_a of a titratable group due to change in electrostatic potential of a mutation was calculated and the change in pH optimum was predicted from the change in pK_a of the catalytic residues. Our strategy is proved to be useful in selection of promising mutants to shift the pH optimum of the xylanases towards desired side.     

The ionic environment determines ribozyme cleavage rate by modulation of nucleobase pK a

Several small ribozymes employ general acid–base catalysis as a mechanism to enhance site-specific RNA cleavage, even though the functional groups on the ribonucleoside building blocks of RNA have pK_a values far removed from physiological pH. The rate of the cleavage reaction is strongly affected by the identity of the metal cation present in the reaction solution; however, the mechanism(s) by which different cations contribute to rate enhancement has not been determined. Using the Neurospora VS ribozyme, we provide evidence that different cations confer particular shifts in the apparent pK_a values of the catalytic nucleobases, which in turn determines the fraction of RNA in the protonation state competent for general acid–base catalysis at a given pH, which determines the observed rate of the cleavage reaction. Despite large differences in observed rates of cleavage in different cations, mathematical models of general acid–base catalysis indicate that k₁, the intrinsic rate of the bond-breaking step, is essentially constant irrespective of the identity of the cation(s) in the reaction solution. Thus, in contrast to models that invoke unique roles for metal ions in ribozyme chemical mechanisms, we find that most, and possibly all, of the ion-specific rate enhancement in the VS ribozyme can be explained solely by the effect of the ions on nucleobase pK_a. The inference that k₁ is essentially constant suggests a resolution of the problem of kinetic ambiguity in favor of a model in which the lower pK_a is that of the general acid and the higher pK_a is that of the general base.     

Kinetic evidence for surface residues influencing the active site of Coprinus cinereus peroxidase: analysis of the pH dependence of G154E,P90H and P90H–G154E substrate entrance mutants

Three mutants of Coprinus cinereus peroxidase (CIP) were made to mimic the substrate entrance histidine 82–glutamic acid 146 pair of the substrate channel in lignin peroxidase (LIP). Compound I formation of LIP has a low pH optimum around pH 3, while optimal formation of CIP compound I is obtained at pH 6–11. The mutants were glycine 154→glutamic acid (G154E), proline 90→histidine (P90H) and the double mutant P90H–G154E. All three showed kinetics of compound I formation similar to that of wt CIP between pH 3 and 9. However, the stability of compound I was strongly affected by these mutations. In wt CIP compound I is stable for approximately 30 min, while compound I of the mutants were stable for 5 s or less. The P90H and P90H–G154E mutants showed pK_a values for the alkaline transition at least one pH unit lower than for wt CIP and the G154E mutant. We suggest that the changed electrostatic field results in destabilisation of the oxidised heme in compound I and II and that the P90H residue increases the electrostatic potential in the distal cavity thereby decreasing the pK_a for the alkaline transition.     

Electrostatic control of enzyme reactions: The mechanism of inhibition of glucose oxidase by putrescine

The interaction of putrescine dihydrochloride with glucose oxidase is reported. At pH 7.65 glucose oxidase is strongly anionic (Z = ?80). The pK_a of an essential acidic group on the reduced form of the enzyme is extremely sensitive to ionic strength, as predicted by simple electrostatic theory [J. G. Voet, J. Coe, J. Epstein, V. Matossian, and T. Shipley (1981), Biochemistry, 20, 7182–7185]. Putrescine dihydrochloride was found to inhibit glucose oxidase at pH 7.65 at a constant ionic strength of 0.05. The kinetics do not obey simple competitive inhibition, however. The data can best be explained by a model in which change in the electrostatic potential of the enzyme on putrescine binding changes the observed pK_a of the essential acidic group. The pH dependence of putrescine inhibition supports this interpretation. At I = 0.05, 5 mM putrescine was found to change the pK_a of the essential acidic group from 7.6 to 7.1. The shift in the pK_a as a function of putrescine concentration at pH 7.7 and I = 0.05 also supports the model presented. The K_a for putrescine to the active form of the enzyme was calculated to be 4.2 mm.     

Combined computational and biochemical study reveals the importance of electrostatic interactions between the “pH sensor” and the cation binding site of the sodium/proton antiporter NhaA of Escherichia coli

Sodium proton antiporters are essential enzymes that catalyze the exchange of sodium ions for protons across biological membranes. The crystal structure of NhaA has provided a basis to explore the mechanism of ion exchange and its unique regulation by pH. Here, the mechanism of the pH activation of the antiporter is investigated through functional and computational studies of several variants with mutations in the ion‐binding site (D163, D164). The most significant difference found computationally between the wild type antiporter and the active site variants, D163E and D164N, are low pK_a values of Glu78 making them insensitive to pH. Although in the variant D163N the pK_a of Glu78 is comparable to the physiological one, this variant cannot demonstrate the long‐range electrostatic effect of Glu78 on the pH‐dependent structural reorganization of trans‐membrane helix X and, hence, is proposed to be inactive. In marked contrast, variant D164E remains sensitive to pH and can be activated by alkaline pH shift. Remarkably, as expected computationally and discovered here biochemically, D164E is viable and active in Na⁺/H⁺ exchange albeit with increased apparent K_M. Our results unravel the unique electrostatic network of NhaA that connect the coupled clusters of the “pH sensor” with the binding site, which is crucial for pH activation of NhaA. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Structures and free-energy landscapes of the wild type and mutants of the Abeta(21-30) peptide are determined by an interplay between intrapeptide electrostatic and hydrophobic interactions

The initial events in protein aggregation involve fluctuations that populate monomer conformations, which lead to oligomerization and fibril assembly. The highly populated structures, driven by a balance between hydrophobic and electrostatic interactions in the protease-resistant wild-type Aβ_21-30 peptide and mutants E22Q (Dutch), D23N (Iowa), and K28N, are analyzed using molecular dynamics simulations. Intrapeptide electrostatic interactions were connected to calculated pK_a values that compare well with the experimental estimates. The pK_a values of the titratable residues show that E22 and D23 side chains form salt bridges only infrequently with the K28 side chain. Contacts between E22-K28 are more probable in “dried” salt bridges, whereas D23-K28 contacts are more probable in solvated salt bridges. The strength of the intrapeptide hydrophobic interactions increases as D23N < WT < E22Q < K28A. Free-energy profiles and disconnectivity representation of the energy landscapes show that the monomer structures partition into four distinct basins. The hydrophobic interactions cluster the Aβ_21-30 peptide into two basins, differentiated by the relative position of the DVG(23-25) and GSN(25-27) fragments about the G25 residue. The E22Q mutation increases the population with intact VGSN turn compared to the wild-type (WT) peptide. The increase in the population of the structures in the aggregation-prone Basin I in E22Q, which occurs solely due to the difference in charge states between the Dutch mutant and the WT, gives a structural explanation of the somewhat larger aggregation rate in the mutant. The D23N mutation dramatically reduces the intrapeptide interactions. The K28A mutation increases the intrapeptide hydrophobic interactions that promote population of structures in Basin I and Basin II whose structures are characterized by hydrophobic interaction between V24 and K28 side chains but with well-separated ends of the backbone atoms in the VGSN turn. The intrapeptide electrostatic interactions in the WT and E22Q peptides roughen the free-energy surface compared to the K28A peptide. The D23N mutation has a flat free-energy surface, corresponding to an increased population of random coil-like structures with weak hydrophobic and electrostatic interactions. We propose that mutations or sequences that enhance the probability of occupying Basin I would promote aggregation of Aβ peptides.     

Arginine: Its pKa value revisited

Using complementary approaches of potentiometry and NMR spectroscopy, we have determined that the equilibrium acid dissociation constant (pK_a value) of the arginine guanidinium group is 13.8 ± 0.1. This is substantially higher than that of ∼12 often used in structure-based electrostatics calculations and cited in biochemistry textbooks. The revised intrinsic pK_a value helps explains why arginine side chains in proteins are always predominantly charged, even at pH values as great as 10. The high pK_a value also reinforces the observation that arginine side chains are invariably protonated under physiological conditions of near neutral pH. This occurs even when the guanidinium moiety is buried in a hydrophobic micro-environment, such as that inside a protein or a lipid membrane, thought to be incompatible with the presence of a charged group.     

Steady-state kinetics of electron transfer through cytochrome chain of uncoupled submitochondrial particles. II. Influence of pH on kinetics of electron transfer

pH Dependences of steady-state kinetic parameters of cytochrome chains of submitochondrial particles have been studied. It has been shown that the lifetimes of activated states (τ) of the pairs of cytochromes b → c₁ and a → a₃ have different pH dependences; those for the c₁ → c and c → a cytochrome pairs being similar. The rate constants for the non-activated state of the respiratory chains decreased for the b → c₁ pair and increased for the a → a₃ pair when the pH value was increased.The values of pK calculated from these dependences for the pairs b → c₁ and a → a₃ were 7.2 and 8.9, respectively. It has been supposed that the ratio of activated to non-activated electron carriers may be controlled by the local pH value in the mitochondrial membrane, the latter being dependent upon the rate of electron transfer. The kinetic model based on this assumption allows one to explain the experimental dependences on pH of the rate constants for cytochromes b → c, and a → a₃.The values of the diffusion rate constants for H⁺ and OH^? ions in the mitochondrial membrane estimated from these kinetic data obtained in this study weree 10⁴–10⁵ s^?1 and 10²–10³ s^?1, respectively.