Structure-based pKa prediction provides a thermodynamic basis for the role of histidines in pH-induced conformational transitions in dengue virus

pH-induced conformational changes in dengue virus (DENV) are critical to its ability to infect host cells. The envelope protein heterodimers that make up the viral envelope shift from a dimer to a trimer conformation at low-pH during membrane fusion. Previous studies have suggested that the ionization of histidine residues at low-pH is central to this pH-induced conformational change. We sought out to use molecular modeling with structure-based pKa prediction to provide a quantitative basis for the role of histidines in pH-induced conformational changes and identify which histidine residues were primarily responsible for this transition. We combined existing crystallographic and cryo-electron microscopy data to construct templates of the dimer and trimer conformations for the mature and immature virus. We then generated homology models for the four DENV serotypes and carried out structure-based pKa prediction using Rosetta. Our results showed that the pKa values of a subset of conserved histidines in DENV successfully capture the thermodynamics necessary to drive pH-induced conformational changes during fusion. Here, we identified the structural determinants underlying these pKa values and compare our findings with previous experimental results.     

Electrostatics in the active site of an alpha-amylase.

The importance of electrostatics in catalysis has been emphasized in the literature for a large number of enzymes. We examined this hypothesis for the Bacillus licheniformis alpha-amylase by constructing site-directed mutants that were predicted to change the pKa values of the catalytic residues and thus change the pH-activity profile of the enzyme. To change the pKa of the catalytic residues in the active site, we constructed mutations that altered the hydrogen bonding network, mutations that changed the solvent accessibility, and mutations that altered the net charge of the molecule. The results show that changing the hydrogen bonding network near an active site residue or changing the solvent accessibility of an active site residue will very likely result in an enzyme with drastically reduced activity. The differences in the pH-activity profiles for these mutants were modest. pH-activity profiles of mutants which change the net charge on the molecule were significantly different from the wild-type pH-activity profile. The differences were, however, difficult to correlate with the electrostatic field changes calculated. In several cases we observed that pH-activity profiles shifted in the opposite direction compared to the shift predicted from electrostatic calculations. This strongly suggests that electrostatic effects cannot be solely responsible for the pH-activity profile of the B. licheniformis alpha-amylase.     

DNA minor groove recognition by bis-benzimidazole analogues of Hoechst 33258: insights into structure-DNA affinity relationships assessed by fluorescence titration measurements

Fluorescence titration measurements have been used to examine the binding interaction of a number of analogues of the bis -benzimidazole DNA minor groove binding agent Hoechst 33258 with the decamer duplex d(GCAAATTTGC)2. The method of continuous variation in ligand concentration (Job plot analysis) reveals a 1:1 binding stoichiometry for all four analogues; binding constants are independent of drug concentration (in the range [ligand] = 0.1-5 microM). The four analogues studied were chosen in order to gain some insight into the relative importance of a number of key structural features for minor groove recognition, namely (i) steric bulk of the N -methylpiperazine ring, (ii) ligand hydrophobicity, (iii) isohelicity with the DNA minor groove and (iv) net ligand charge. This was achieved, first, by replacing the bulky, non-planar N -methylpiperazine ring with a less bulky planar charged imidazole ring permitting binding to a narrower groove, secondly, by linking the N -methylpiperazine ring to the phenyl end of the molecule to give the molecule a more linear, less isohelical conformation and, finally, by introducing a charged imidazole ring in place of the phenolic OH making it dicationic, enabling the contribution of the additional electrostatic interaction and extended conformation to be assessed. Delta G values were measured at 20 degrees C in the range -47.6 to -37.5 kJ mol-1 and at a number of pH values between 5.0 and 7.2. We find a very poor correlation between Delta G values determined by fluorescence titration and effects of ligand binding on DNA melting temperatures, concluding that isothermal titration methods provide the most reliable method of determining binding affinities. Our results indicate that the bulky N -methylpiperazine ring imparts a large favourable binding interaction, despite its apparent requirement for a wider minor groove, which others have suggested arises in a large part from the hydrophobic effect. The binding constant appears to be insensitive to the isohelical arrangement of the constituent rings which in these analogues gives the same register of hydrogen bonding interactions with the floor of the groove.     

Molecular mechanisms of pH-driven conformational transitions of proteins: insights from continuum electrostatics calculations of acid unfolding

The acid unfolding of staphylococcal nuclease (SNase) is very cooperative (Whitten and García-Moreno, Biochemistry 2000;39:14292-14304). As many as seven hydrogen ions (H+) are bound preferentially by the acid-unfolded state relative to the native (N) state in the pH range 3.2-3.9. To investigate the mechanism of acid unfolding, structure-based pKa calculations were performed with a variety of continuum electrostatic methods. The calculations reproduced successfully the H+ binding properties of the N state between pH 5 and 9, but they systematically overestimated the number of H+ bound upon acid unfolding. The calculated pKa values of all carboxylic residues in the N state were more depressed than they should be. The discrepancy between the observed and the calculated H+ uptake upon acid unfolding was not improved by using high protein dielectric constants, structures relaxed with molecular dynamics, or other empirical modifications implemented previously by others to maximize agreement between measured and calculated pKa values. This suggests an important role for conformational fluctuations of the backbone as important determinants of pKa values of carboxylic groups. Because no global or subglobal conformational changes have been observed previously for SNase under acidic conditions above the acid-unfolding region, these fluctuations must be local. The acid unfolding of SNase does not seem to involve the disruption of the N state by accruement of intramolecular repulsive interactions, nor the protonation of key ion paired carboxylic residues. It is more consistent with modest contributions from many H+ binding groups, with an important role for local conformational fluctuations in the coupling between H+ binding and the global structural transition.     

Side chain electrostatic interactions and pH‐dependent expansion of the intrinsically disordered,highly acidic carboxyl‐terminus of γ‐tubulin

Intramolecular electrostatic attraction and repulsion strongly influence the conformational sampling of intrinsically disordered proteins and domains (IDPs). In order to better understand this complex relationship, we have used nuclear magnetic resonance to measure side chain pK_a values and pH‐dependent translational diffusion coefficients for the unstructured and highly acidic carboxyl‐terminus of γ‐tubulin (γ‐CT), providing insight into how the net charge of an IDP relates to overall expansion or collapse of the conformational ensemble. Many of the pK_a values in the γ‐CT are shifted upward by 0.3–0.4 units and exhibit negatively cooperative ionization pH profiles, likely due to the large net negative charge that accumulates on the molecule as the pH is raised. pK_a shifts of this magnitude correspond to electrostatic interaction energies between the affected residues and the rest of the charged molecule that are each on the order of 1 kcal mol^?1. Diffusion of the γ‐CT slowed with increasing net charge, indicative of an expanding hydrodynamic radius (r_H). The degree of expansion agreed quantitatively with what has been seen from comparisons of IDPs with different charge content, yielding the general trend that every 0.1 increase in relative charge (|Q|/res) produces a roughly 5% increase in r_H. While γ‐CT pH titration data followed this trend nearly perfectly, there were substantially larger deviations for the database of different IDP sequences. This suggests that other aspects of an IDP's primary amino acid sequence beyond net charge influence the sensitivity of r_H to electrostatic interactions.     

Generation of the O630 photointermediate of bacteriorhodopsin is controlled by the state of protonation of several protein residues

The last stages of the photocycle of the photosynthetic pigment all-trans bacteriorhodopsin (bR570), as well as its proton pump mechanism, are markedly pH dependent. We have measured the relative amount of the accumulated O630 intermediate (Phir), as well as its rise and decay rate constants (kr and kd, respectively), over a wide pH range. The experiments were carried out in deionized membrane suspensions to which varying concentrations of metal cations and of large organic cations were added. The observed pH dependencies, s-shaped curves in the case of Phir and bell-shaped curves for kr and kd, are interpreted in terms of the titration of three protein residues denoted as R1, R2, and R3. The R1 titration is responsible for the increase in Phir, kr, and kd upon lowering the pH from pH approximately 9.5 to 7. At low pH Phir exhibits a secondary rise which is attributed to the titration of a low pKa group, R2. After reaching a maximum at pH approximately 7, kr and kd undergo a decrease upon decreasing the pH, which is attributed to the titration of R3. All three titrations exhibit pKa values which decrease upon increasing the salt concentration. As in the case of the Purple (bR570) if Blue (bR605) equilibrium, divalent cations are substantially more effective than monovalent cations in shifting the pKa values. Moreover, bulky organic cations are as effective as small metal cations. It is concluded that analogously to the Purple if Blue equilibrium, the salt binding sites which control the pKa values of R1, R2, and R3 are located on, or close to, the membrane surface. Possible identifications of the three protein residues are considered. Experiments with the E204Q mutant show that the mutation has markedly affected the R2 (Phir) titration, suggesting that R2 should be identified with Glu-204 or with a group whose pKa is affected by Glu-204. The relation between the R1, R2 and R3 titrations and the proton pump mechanism is discussed. It is evident that the pH dependence of Phir is unrelated to the measured pKa of the group (XH) which releases the proton to the extracellular medium during the photocycle. However, since the same residue may exhibit different pKa values at different stages of the photocycle, it cannot be excluded that R2 or R3 may be identified with XH.     

Electrostatic contributions to residue-specific protonation equilibria and proton binding capacitance for a small protein

Charge-charge interactions in proteins are important in a host of biological processes. Here we use 13C NMR chemical shift data for individual aspartate and glutamate side chain carboxylate groups to accurately detect site-specific protonation equilibria in a variant of the B1 domain of protein G (PGB1-QDD). Carbon chemical shifts are dominated by changes in the electron distribution within the side chain and therefore excellent reporters of the charge state of individual groups, and the data are of high precision. We demonstrate that it is possible to detect local charge interactions within this small protein domain that stretch and skew the chemical shift titration curves away from "ideal" behavior and introduce a framework for the analysis of such convoluted data to study local charge-charge interactions and electrostatic coupling. It is found that, due to changes in electrostatic potential, the proton binding affinity, Ka, of each carboxyl group changes throughout the titration process and results in a linearly pH dependent pKa value. This result could be readily explained by calculations of direct charge-charge interactions based on Coulomb's law. In addition, the slope of pKa versus pH was dependent on screening by salt, and this dependence allowed the selective study of charge-charge interactions. For PGB1-QDD, it was established that mainly differences in self-energy, and not direct charge-charge interactions, are responsible for shifted pKa values within the protein environment.     

The titrations of Asp-85 and of the cation binding residues in bacteriorhodopsin are not coupled

An outstanding problem relating to the structure and function of bacteriorhodopsin (bR), which is the only protein in the purple membrane of the photosynthetic microorganism Halobacterium salinarium, is the relation between the titration of Asp-85 and the binding/unbinding of metal cations. An extensively accepted working hypothesis has been that the two titrations are coupled, namely, protonation of Asp-85 (located in the vicinity of the retinal chromophore) and cation unbinding occur concurrently. We have carried out a series of experiments in which the purple blue equilibrium and the binding of Mn2+ ions (monitored by electron spin resonance) were followed as a function of pH for several (1-4) R = [Mn2+]/[bR] molar ratios. Data were obtained for native bR, bR mutants, artificial bR and chemically modified bR. We find that in the native pigment the two titrations are separated by more than a pKa unit [delta pKa = pKa(P/B)-pKa(Mn2+) = (4.2-2.8) = 1.4]. In the non-native systems, delta pKa values as high as 5 units, as well as negative delta pKas, are observed. We conclude that the pH titration of cation binding residues in bR is not directly related to the titration of Asp-85. This conclusion is relevant to the nature of the high affinity cation sites in bR and to their role in the photosynthetic function of the pigment.     

Titration_DB: Storage and analysis of NMR‐monitored protein pH titration curves

NMR‐monitored pH titration experiments are routinely used to measure site‐specific protein pKa values. Accurate experimental pKa values are essential in dissecting enzyme catalysis, in studying the pH‐dependence of protein stability and ligand binding, in benchmarking pKa prediction algorithms, and ultimately in understanding electrostatic effects in proteins. However, due to the complex ways in which pH‐dependent electrostatic and structural changes manifest themselves in NMR spectra, reported apparent pKa values are often dependent on the way that NMR pH‐titration curves are analyzed. It is therefore important to retain the raw NMR spectroscopic data to allow for documentation and possible re‐interpretation. We have constructed a database of primary NMR pH‐titration data, which is accessible via a web interface. Here, we report statistics of the database contents and analyze the data with a global perspective to provide guidelines on best practice for fitting NMR titration curves. Titration_DB is available at http://enzyme.ucd.ie/Titration_DB . Proteins 2010. © 2009 Wiley‐Liss, Inc.     

Conformation, stability, and active-site cysteine titrations of Escherichia coli D26A thioredoxin probed by Raman spectroscopy.

The active-site cysteines (Cys 32 and Cys 35) of Escherichia coli thioredoxin are oxidized to a disulfide bridge when the protein mediates substrate reduction. In reduced thioredoxin, Cys 32 and Cys 35 are characterized by abnormally low pKa values. A conserved side chain, Asp 26, which is sterically accessible to the active site, is also essential to oxidoreductase activity. pKa values governing cysteine thiol-thiolate equilibria in the mutant thioredoxin, D26A, have been determined by direct Raman spectrophotometric measurement of sulfhydryl ionizations. The results indicate that, in D26A thioredoxin, both sulfhydryls titrate with apparent pKa values of 7.5+/-0.2, close to values measured previously for wild-type thioredoxin. Sulfhydryl Raman markers of D26A and wild-type thioredoxin also exhibit similar band shapes, consistent with minimal differences in respective cysteine side-chain conformations and sulfhydryl interactions. The results imply that neither the Cys 32 nor Cys 35 SH donor is hydrogen bonded directly to Asp 26 in the wild-type protein. Additionally, the thioredoxin main-chain conformation is largely conserved with D26A mutation. Conversely, the mutation perturbs Raman bands diagnostic of tryptophan (Trp 28 and Trp 31) orientations and leads to differences in their pH dependencies, implying local conformational differences near the active site. We conclude that, although the carboxyl side chain of Asp 26 neither interacts directly with active-site cysteines nor is responsible for their abnormally low pKa values, the aspartate side chain may play a role in determining the conformation of the enzyme active site.     

A pH-dependent conformational change in the coat protein subunits from potato virus X.

Both the circular dichroism and fluorescence spectra of the dissociated coat protein subunits from potato virus X changed substantially over the pH range 8 to 4, irreversible changes resulted below pH 4, with tyrosyl and tryptophanyl residues affected most. The titration curves show a pKa of about 5.6 and do not require cooperative interactions between the coat protein subunits, thus they are in marked contrast to titrations of tobacco mosaic virus A-protein. The spectra of the intact virus were little changed between pH 8 and 4 and suggested that the coat protein was locked into a conformation similar to that of the subunits in solution at pH 7. It is proposed that the pH induced conformational change is responsible for determining the acidic branch of the pH profile for reconstitution of potato virus X from its dissociated coat protein subunits and RNA.     

Investigation of the haem-nicotinate interaction in leghaemoglobin. Role of hydrogen bonding.

A strategic assessment of the contributions of two active-site hydrogen bonds in the binding of nicotinate to recombinant ferric soybean leghaemoglobin a (rLb) was carried out by mutagenic replacement of the hydrogen-bonding residues (H61A and Y30A variants) and by complementary chemical substitution of the carboxylate functionality on the nicotinate ligand. Dissociation constants, Kd (pH 5.5, mu = 0.10 M, 25.0 +/- 0.1 degrees C), for binding of nicotinate to ferric rLb, H61A and Y30A were 1.4 +/- 0.3 microM, 19 +/- 1 microM and 11 +/- 1 microM, respectively; dissociation constants for binding of nicotinamide were, respectively, 38 +/- 1 mM, 50 +/- 2 mM and 12 +/- 1 mM, and for binding of pyridine were 260 +/- 50 microM, 4.5 +/- 0.5 microM and 66 +/- 8 microM, respectively. Binding of cyanide and azide to the H61A and Y30A variants was unaffected by the mutations. The pH-dependence of nicotinate binding for rLb and Y30A was consistent with a single titration process (pKa values 6.9 +/- 0.1 and 6.7 +/- 0.2, respectively); binding of nicotinate to H61A was independent of pH. Reduction potentials for the rLb and rLb-nicotinate derivatives were 29 +/- 2 mV (pH 5.40, 25.0 degrees C, mu = 0.10 M) and - 65 +/- 2 mV (pH 5.42, 25.0 degrees C, mu = 0.10 M), respectively. The experiments provide a quantitative assessment of the role of individual hydrogen bonds in the binding process, together with a definitive determination of the pKa of His61 and unambiguous evidence that titration of His61 controls binding in the neutral to alkaline region.     

Nuclear magnetic resonance titration curves of histidine ring protons. Human metmyoglobin and the effects of azide on human, horse, and sperm whale metmyoglobins.

Four titrating histidine ring C2 and C4 proton resonances are observed in 220 MHz proton NMR spectra of human metmyoglobin as a function of pH. Values of ionization constants determined from the NMR titration data using an equation describing a simple proton association-dissociation equilibrium are curves (1) 6.6, (2) 7.0, (3) 5.8, and (4) 7.4. Four histidine residues have also been found to be solvent-accessible in human metmyoglobin by carboxymethylation studies (Harris, C.M., and Hill, R.L. (1969) J. Biol. Chem. 244, 2195-2203). Two of the titration curves (3 and 4) deviate significantly from the chemical shift values normally observed for histidine C2 proton resonances. Curve 3, with a low pKa, is shifted downfield at high values of pH and also exhibits a second minor inflection with a pKa value of 8.8. On the other hand, the high pKa curve, 4, is shifted upfield at all values of pH. The characteristics of the NMR titration curves with the lowest and highest pKa values (3 and4) are very similar to curves observed previously with sperm whale and horse metmyoglobins (Cohen, J.S., Hagenmaier, H., Pollard, H., and Schechter, A.N. (1972) J. Mol. Biol. 71, 513-519). These results indicate that the histidine residues from which these curves are derived have unusual and characteristic environments in this series of homologous proteins. The NMR spectra of all three metmyoglobins are changed extensively as a result of azide ion binding, indicating conformational changes affecting the environments of several imidazole side chains. The presence of azide ion causes a selective downfield chemical shift for the low pKa curve and a selective upfield chemical shift for the high pKa curve in all three proteins. Azide also abolishes the second inflection seen in the low pKa curve at high pH. In addition to these effects, the presence of azide ion permits the observation of two additional titrating proton resonances for all three metmyoglobins. Increasing the azide to protein ratio at several fixed values of pH yields results which show that a slow exchange process is occurring with each of the metmyoglobins. In the azide titration studies the maximum changes in the NMR spectra occurred at approximately equimolar concentrations. The NMR results for these proteins in the absence and presence of azide ion are related to x-ray crystallographic studies of sperm whale metmyoglobin and the known alkylation properties of the histidine residues. Tentative assignments of the titrating resonances observed are suggested.     

Remeasuring HEWL pK(a) values by NMR spectroscopy: methods, analysis, accuracy, and implications for theoretical pK(a) calculations

Site-specific pK(a) values measured by NMR spectroscopy provide essential information on protein electrostatics, the pH-dependence of protein structure, dynamics and function, and constitute an important benchmark for protein pK(a) calculation algorithms. Titration curves can be measured by tracking the NMR chemical shifts of several reporter nuclei versus sample pH. However, careful analysis of these curves is needed to extract residue-specific pK(a) values since pH-dependent chemical shift changes can arise from many sources, including through-bond inductive effects, through-space electric field effects, and conformational changes. We have re-measured titration curves for all carboxylates and His 15 in Hen Egg White Lysozyme (HEWL) by recording the pH-dependent chemical shifts of all backbone amide nitrogens and protons, Asp/Glu side chain protons and carboxyl carbons, and imidazole protonated carbons and protons in this protein. We extracted pK(a) values from the resulting titration curves using standard fitting methods, and compared these values to each other, and with those measured previously by 1H NMR (Bartik et al., Biophys J 1994;66:1180–1184). This analysis gives insights into the true accuracy associated with experimentally measured pK(a) values. We find that apparent pK(a) values frequently differ by 0.5–1.0 units depending upon the nuclei monitored, and that larger differences occasionally can be observed. The variation in measured pK(a) values, which reflects the difficulty in fitting and assigning pH-dependent chemical shifts to specific ionization equilibria, has significant implications for the experimental procedures used for measuring protein pK(a) values, for the benchmarking of protein pK(a) calculation algorithms, and for the understanding of protein electrostatics in general.     

Free energy of proton binding in proteins

In this article we use literature data on the titration of denatured ribonuclease to test the accuracy of proton-binding distributions obtained using our recent approach employing moments. We find that using only the local slope of the titration curve at a small number of points (five, for example) we can reproduce the detailed proton-binding distribution at all pH values. Our method gives the complete proton-binding polynomial for a given protein and each coefficient in this polynomial in turn yields the free energy for binding a given number of protons in all ways to the protein. Using these net free energies, we can then compute the average proton-binding free energy per proton as a function of the fraction of protons bound. We find that this function is remarkably similar for different proteins, even for proteins that exhibit quite different titration behavior. For the special case of binding to independent sites, we obtain simple relations for the first and last terms in the free energy per-proton function. For this special case we also can calculate the distribution functions giving the probability that a molecule has a given number of positive or negative charges and the joint distribution that a molecule simultaneously has a given number of positive and negative charge.     

pH-induced structural changes in bacteriorhodopsin studied by Fourier transform infrared spectroscopy.

Previous C13-NMR studies showed that two of the four internal aspartic acid residues (Asp-96 and Asp-115) of bacteriorhodopsin (bR) are protonated up to pH = 10, but no accurate pKa of these residues has been determined. In this work, infrared spectroscopy with the attenuated total reflection technique was used to characterize pH-dependent structural changes of ground-state, dark-adapted wild-type bacteriorhodopsin and its mutant (D96N) with aspartic acid-96 replaced by asparagine. Data indicated deprotonation of Asp-96 at high pH (pKa = 11.4 +/- 0.1), but no Asp-115 titration was observed. The analysis of the whole spectral region characteristic to complex conformational changes in the protein showed a more complicated titration with an additional pKa value (pKa1 = 9.3 +/- 0.3 and pKa2 = 11.5 +/- 0.2). Comparison of results obtained for bR and the D96N mutant of bR shows that the pKa approximately 11.5 characterizes not a direct titration of Asp-96 but a protein conformational change that makes Asp-96 accessible to the external medium.     

pH dependence of the interaction of hirudin with thrombin.

The kinetics of the inhibition of human alpha-thrombin by recombinant hirudin have been studied over the pH range from 6 to 10. The association rate constant for hirudin did not vary significantly over this pH range. The dissociation constant of hirudin depended on the ionization state of groups with pKa values of about 7.1, 8.4, and 9.2. Optimal binding of hirudin to thrombin occurred when the groups with pKa values of 8.4 and 9.0 were protonated and the other group with a pKa of 7.1 was deprotonated. The pH kinetics of genetically engineered forms of hirudin were examined in an attempt to assign these pKa values to particular groups. By using this approach, it was possible to show that protonation His51 and ionization of acidic residues in the C-terminal region of hirudin were not responsible for the observed pKa values. In contrast, the pKa value of 8.4 was not observed when a form of hirudin with an acetylated alpha-amino group was examined, and, thus, this pKa value was assigned to the alpha-amino group of hirudin. The requirement for this group to be protonated for optimal binding to thrombin is discussed in terms of the crystal structure of the thrombin-hirudin complex. Examination of this structure allowed the other pKa values of 7.1 and 9.2 to be tentatively attributed to His57 and the alpha-amino group of Ile16 of thrombin.     

Influence of anions and pH on the conformational change of horse liver alcohol dehydrogenase induced by binding of oxidized nicotinamide adenine dinucleotide: binding of chloride to the catalytic metal ion

The conformational change of horse liver alcohol dehydrogenase induced by binding of NAD+ was studied by electronic absorption spectroscopy using cobalt as a spectroscopic probe in the active site. The complex of the enzyme with NAD+ exists in an acidic and an alkaline form. The transition between the two forms proceeds through several intermediates and is controlled by an apparent pKa of 6.9. Only at pH values below this pKa can a complex between enzyme, NAD+, and Cl- be formed. The spectral changes indicate that chloride displaces the cobalt-bound water molecule in a tetracoordinate structure. We conclude that a negative charge at the active site is necessary to stabilize the closed conformation of the enzyme in the presence of NAD+. Spectral correlations are given which strongly support the postulation of a metal-bound alkoxide in the closed structure of the enzyme as an essential feature of the catalytic mechanism of horse liver alcohol dehydrogenase.     

Calcium-binding properties of cardiac and skeletal troponin C as determined by circular dichroism and ultraviolet difference spectroscopy.

Calcium titration of the conformational change in cardiac and skeletal troponin C (TN-C) was followed by circular dichroism (CD) at pH values in the range from 5.2 to 7.4. Computer analysis was used to resolve the contributions from the different classes of Ca2+ -binding sites. At pH 6.94 in skeletal TN-C, apparent affinity constants for calcium of 1.8 x 10(7) and 4.5 x 10(5) M-1 were determined for the two classes of binding sites. The more sophisticated computer analysis of the data has revealed a substantial CD contribution from the low-affinity sites (approximately 30% of the high affinity contribution at pH 6.94) and suggests that skeletal TN-C with Ca2+ bound at the low-affinity sites is in a different conformation from that when just the high-affinity sites are occupied, in agreement with a recent nuclear magnetic resonance (NMR) study on this system (Seaman, K. B., Hartshorne, D. J. & Bothener-By, A. A. (1977) Biochemistry 16,4039-4046). With the cardiac protein at pH 7.07, an apparent affinity constant for calcium of 2.0 x 10(7) M-1 was calculated while no low-affinity site at this pH was detected by CD. On the other hand, at lower pH values, such as 6.05, a CD contribution from the cardiac low-affinity Ca2+ -binding site is detected with an apparent binding constant of 3.7 +/- 0.7 x 10(4) M-1. At the lower pH values, protonation of a class of carboxyl groups in each protein which possesses a high pKa (6.2-6.3) elicits the conformational change at the high-affinity sites with a corresponding decrease in the overall magnitude of the Ca2+ -evoked changes. The expression of a conformational change upon Ca2+ binding at the level of the low-affinity sites is enchanced by protonation of a class of carboxyls with a pKa of 6.3 in cardiac TN-C and 6.7-6.8 with the skeletal homologue. In both cases, this contribution is reduced upon protonation of carboxyls with pKa less than or equal to 5.5. It was also observed that the low-affinity sites of skeletal TN-C have a much larger role to play in the total conformational change than the low-affinity sites of cardiac TN-C, a finding probably related to the inability of site 1 in the cardiac protein to bind calcium. In the cardiac protein, the Ca2+ -induced tyrosine difference-spectrum maximum is reduced from deltaepsilonM,287nm =330M-1.cm-1 to 20M-1.cm-1 by protonation of a class of groups with a pKa of 6.4, presumably the same carboxyl groups as those invoved in the CD conformational contribution from the high-affinity binding sites. No such effect was observed for the skeletal protein where deltaepsilonM,287nm was constant at 110M-1 .cm-1 over the pH range studied. The dramatic alterations in the tyrosine environment of cardiac TN-C with pH are attributed to either or both of the tyrosines located in the two high-affinity Ca2+ -binding sites (sites 3 and 4)...     

Structure, mechanism, and conformational dynamics of O-acetylserine sulfhydrylase from Salmonella typhimurium: comparison of A and B isozymes

O-Acetylserine sulfhydrylase is a pyridoxal 5'-phosphate-dependent enzyme that catalyzes the final step in the cysteine biosynthetic pathway in enteric bacteria and plants, the replacement of the beta-acetoxy group of O-acetyl-l-serine by a thiol to give l-cysteine. Two isozymes are found in Salmonella typhimurium, with the A-isozyme expressed under aerobic and the B-isozyme expressed under anaerobic conditions. The structure of O-acetylserine sulfhydrylase B has been solved to 2.3 A and exhibits overall a fold very similar to that of the A-isozyme. The main difference between the two isozymes is the more hydrophilic active site of the B-isozyme with two ionizable residues, C280 and D281, replacing the neutral residues S300 and P299, respectively, in the A-isozyme. D281 is above the re face of the cofactor and is within hydrogen-bonding distance to Y286, while C280 is located about 3.4 A from the pyridine nitrogen (N1) of the internal Schiff base. The B-isozyme has a turnover number (V/Et) 12.5-fold higher than the A-isozyme and an approximately 10-fold lower Km for O-acetyl-l-serine. Studies of the first half-reaction by rapid-scanning stopped-flow indicate a first-order conversion of the internal Schiff base to the alpha-aminoacrylate intermediate at any concentration of O-acetyl-l-serine. The Kd values for formation of the external Schiff base with cysteine and serine, obtained by spectral titration, are pH dependent and exhibit a pKa of 7.0-7.5 (for a group that must be unprotonated for optimum binding) with values, above pH 8.0, of about 3.0 and 30.0 mM, respectively. In both cases the neutral enolimine is favored at high pH. Failure to observe the pKa for the alpha-amines of cysteine and serine in the pKESB vs pH profile suggests a compensatory effect resulting from titration of a group on the enzyme with a pKa in the vicinity of the alpha-amine's pKa. The pH dependence of the first-order rate constant for decay of the alpha-aminoacrylate intermediate to give pyruvate and ammonia gives a pKa of about 9 for the active site lysine (K41), a pH unit higher than that of the A-isozyme. The difference in pH dependence of the pKESB for cysteine and serine, the higher pKa for K41, and the preference for the neutral species at high pH compared to the A-isozyme can be explained by titration of C280 to give the thiolate. Subtle conformational differences between O-acetylserine sulfhydrylase A and O-acetylserine sulfhydrylase B are detected by comparing the absorption and emission spectra of the internal aldimine in the absence and presence of the product acetate and of the external aldimine with l-serine. The two isozymes show a different equilibrium distribution of the enolimine and ketoenamine tautomers, likely as a result of a more polar active site for O-acetylserine sulfhydrylase B. The distribution of cofactor tautomers is dramatically affected by the ligation state of the enzyme. In the presence of acetate, which occupies the alpha-carboxylate subsite, the equilibrium between tautomers is shifted toward the ketoenamine tautomer, as a result of a conformational change affecting the structure of the active site. This finding, in agreement with structural data, suggests for the O-acetylserine sulfhydrylase B-isozyme a higher degree of conformational flexibility linked to catalysis.