Rapid Calculation of Protein pK(a) Values Using Rosetta

We developed a Rosetta-based Monte Carlo method to calculate the pK(a) values of protein residues that commonly exhibit variable protonation states (Asp, Glu, Lys, His, and Tyr). We tested the technique by calculating pK(a) values for 264 residues from 34 proteins. The standard Rosetta score function, which is independent of any environmental conditions, failed to capture pK(a) shifts. After incorporating a Coulomb electrostatic potential and optimizing the solvation reference energies for pK(a) calculations, we employed a method that allowed side-chain flexibility and achieved a root mean-square deviation (RMSD) of 0.83 from experimental values (0.68 after discounting 11 predictions with an error over 2 pH units). Additional degrees of side-chain conformational freedom for the proximal residues facilitated the capture of charge-charge interactions in a few cases, resulting in an overall RMSD of 0.85 pH units. The addition of backbone flexibility increased the overall RMSD to 0.93 pH units but improved relative pK(a) predictions for proximal catalytic residues. The method also captures large pK(a) shifts of lysine and some glutamate point mutations in staphylococcal nuclease. Thus, a simple and fast method based on the Rosetta score function and limited conformational sampling produces pK(a) values that will be useful when rapid estimation is essential, such as in docking, design, and folding.     

Theoretical prediction of relative and absolute pKa values of aminopyridines

This work presents a study aimed at the theoretical prediction of pK(a) values of aminopyridines, as a factor responsible for the activity of these compounds as blockers of the voltage-dependent K(+) channels. To cover a large range of pK(a) values, a total of seven substituted pyridines is considered as a calibration set: pyridine, 2-aminopyridine, 3-aminopyridine, 4-aminopyridine, 2-chloropyridine, 3-chloropyridine, and 4-methylpirydine. Using ab initio G1, G2 and G3 extrapolation methods, and the CPCM variant of the Polarizable Continuum Model for solvation, we calculate gas phase and solvation free energies. pK(a) values are obtained from these data using a thermodynamic cycle for describing protonation in aqueous and gas phases. The results show that the relatively inexpensive G1 level of theory is the most accurate at predicting pK(a) values in aminopyridines. The highest standard deviation with respect to the experimental data is 0.69 pK(a) units for absolute values calculations. The difference increases slightly to 0.74 pK(a) units when the pK(a) is computed relative to the pyridine molecule. Considering only compounds at least as basic as pyridine (the values of interest for bioactive aminopyridines) the error falls to 0.10 and 0.12 pK(a) units for the absolute and relative computations, respectively. The technique can be used to predict the effect of electronegative substituents in the pK(a) of 4-AP, the most active aminopyridine considered in this work. Thus, 2-chloro and 3-chloro-4-aminopyridine are taken into account. The results show a decrease of the pK(a), suggesting that these compounds are less active than 4-AP at blocking the K(+) channel.     

Development, validation, and application of adapted PEOE charges to estimate pKa values of functional groups in protein-ligand complexes

For routine pK(a) calculations of protein-ligand complexes in drug design, the PEOE method to compute partial charges was modified. The new method is applicable to a large scope of proteins and ligands. The adapted charges were parameterized using experimental free energies of solvation of amino acids and small organic ligands. For a data set of 80 small organic molecules, a correlation coefficient of r(2) = 0.78 between calculated and experimental solvation free energies was obtained. Continuum electrostatics pK(a) calculations based on the Poisson-Boltzmann equation were carried out on a validation set of nine proteins for which 132 experimental pK(a) values are known. In total, an overall RMSD of 0.88 log units between calculated and experimentally determined data is achieved. In particular, the predictions of significantly shifted pK(a) values are satisfactory, and reasonable estimates of protonation states in the active sites of lysozyme and xylanase could be obtained. Application of the charge-assignment and pK(a)-calculation procedure to protein-ligand complexes provides clear structural interpretations of experimentally observed changes of protonation states of functional groups upon complex formation. This information is essential for the interpretation of thermodynamic data of protein-ligand complex formation and provides the basis for the reliable factorization of the free energy of binding in enthalpic and entropic contributions. The modified charge-assignment procedure forms the basis for future automated pK(a) calculations of protein-ligand complexes.     

Conformational preferences and pK(a) value of selenocysteine residue

The conformational preferences of the L-selenocysteine (Sec) dipeptides with selenol and selenolate groups (Ac-Sec-NHMe and Ac-Sec(-) -NHMe, respectively) and the apparent (i.e., macroscopic) pK(a) value of the Sec residue have been studied using the dispersion-corrected density functionals M06-2X and B2PLYP-D with the implicit solvation method in the gas phase and in water. In the gas phase, the backbone-to-backbone and/or side chain-to-backbone hydrogen bonds are found to contribute in stabilizing the most preferred conformations for the Sec and Sec(-) residues, as seen for the Cys and Cys(-) residues. However, the polyproline II-like conformations prevail over the conformations with the backbone-to-backbone hydrogen bonds in water because of the weakened hydrogen bonds by the favorable direct interactions between the backbone C?O and H?N groups and water molecules. The Sec and Sec(-) residues are found to adopt more various conformations than the Cys and Cys(-) residues in water, although the most preferred conformations of the neutral and/or anionic forms of the two residues are similar each other in the gas phase and in water. Using the statistically weighted free energies of the Sec and Sec(-) dipeptides in the gas phase and their solvation free energies, the pK(a) value of the Sec residue is estimated to be 5.47 at 25°C, which is in good agreement with the experimental value of 5.43 ± 0.02. It is found that the lower pK(a) value of the selenol side chain for the Sec residue by ～3 units than the thiol side chain for the Cys residue is ascribed to the higher gas-phase acidity of the Sec residue.     

Electrostatic effects in highly charged proteins: salt sensitivity of pKa values of histidines in staphylococcal nuclease

The pK(a) values of most histidines in small peptides and in myoglobin increase on average by 0.30 unit between 0.02 and 1.5 M NaCl [Kao et al. (2000) Biophys. J. 79, 1637]. The DeltapK(a) values reflect primarily the ionic strength dependence of the solvation energy; screening of Coulombic interactions contributes only in a minor way. This implies that Coulombic interactions are weak, or that attractive and repulsive contributions to the pK(a) values are balanced. To distinguish experimentally between these two possibilities, and to further characterize the magnitude and salt sensitivity of surface electrostatic interactions in proteins, the salt dependence of pK(a) values of histidines in staphylococcal nuclease was measured by (1)H NMR spectroscopy. Three of the four histidines titrated with significantly depressed pK(a) values, and the salt sensitivity of all histidine pK(a) values was substantial. In three cases, the pK(a) values increased by a full unit between 0.01 and 1.5 M KCl. Anion-specific effects were found; the pK(a) values measured under equivalent ionic strengths in SCN(-) and SO(4)(2-) were higher than in Cl(-); the order of the sensitivity of pK(a) values to anions was SCN(-) > Cl(-) > SO(4)(2-). Structure-based pK(a) calculations with continuum methods were performed to interpret the measured effects structurally and to test their ability to capture the experimental behavior. Calculations in which the protein interior was treated empirically with a dielectric constant of 20 reproduced the pK(a) values and their dependence on the concentration of Cl(-). According to the calculations, the pK(a) values are depressed because of unfavorable self-energies and repulsive Coulombic interactions. Their striking salt sensitivity reflects screening of weak, repulsive, Coulombic interactions among charges separated by more than 10 A. Long-range Coulombic interactions on the surfaces of proteins are weak, but they can add up to produce substantial electrostatic effects when positive and negative charges are not balanced.     

Peripheral factors in the mediation of cholecystokinin-induced satiety as assessed by comparative potencies of cholecystokinin antagonists.

Cholecystokinin COOH-terminal octapeptide (CCK-8) produces a satiating effect in the rat and other animals upon peripheral administration. Although it has been demonstrated that the receptors which mediate this action are located in the periphery and are of the CCK-A subtype, their anatomical location has not been firmly established. A dense population of CCK receptors in the pyloric sphincter has been suggested as a candidate. We here quantify the potency of several CCK antagonists to inhibit the contractile effect of CCK-8 on the rat pyloric sphincter in vitro. The potent and selective antagonist MK-329 has a Schild pK of 8.85; the less potent but selective antagonist lorglumide (CR-1409) a pK of 6.37; the related antagonist phenoxyacetylproglumide (phi oAc proglumide) a pK of 5.1; and the weak parent compound proglumide a pK of about 3.3. These data can be compared with the potencies of these compounds to inhibit the actions of CCK-8 to produce satiety in the rat; this comparison supports the contention that CCK receptors of the rat pyloric sphincter could in part mediate the satiety effect produced by exogenous CCK-8.     

The barrier for proton transport in aquaporins as a challenge for electrostatic models: the role of protein relaxation in mutational calculations

The origin of the barrier for proton transport through the aquaporin channel is a problem of general interest. It is becoming increasingly clear that this barrier is not attributable to the orientation of the water molecules across the channel but rather to the electrostatic penalty for moving the proton charge to the center of the channel. However, the reason for the high electrostatic barrier is still rather controversial. It has been argued by some workers that the barrier is due to the so-called NPA motif and/or to the helix macrodipole or to other specific elements. However, our works indicated that the main reason for the high barrier is the loss of the generalized solvation upon moving the proton charge from the bulk to the center of the channel and that this does not reflect a specific repulsive electrostatic interaction but the absence of sufficient electrostatic stabilization. At this stage it seems that the elucidation and clarification of the origin of the electrostatic barrier can serve as an instructive test case for electrostatic models. Thus, we reexamine the free-energy surface for proton transport in aquaporins using the microscopic free-energy perturbation/umbrella sampling (FEP/US) and the empirical valence bond/umbrella sampling (EVB/US) methods as well as the semimacroscopic protein dipole Langevin dipole model in its linear response approximation version (the PDLD/S-LRA). These extensive studies help to clarify the nature of the barrier and to establish the "reduced solvation effect" as the primary source of this barrier. That is, it is found that the barrier is associated with the loss of the generalized solvation energy (which includes of course all electrostatic effects) upon moving the proton charge from the bulk solvent to the center of the channel. It is also demonstrated that the residues in the NPA region and the helix dipole cannot be considered as the main reasons for the electrostatic barrier. Furthermore, our microscopic and semimacroscopic studies clarify the problems with incomplete alternative calculations, illustrating that the effects of various electrostatic elements are drastically overestimated by macroscopic calculations that use a low dielectric constant and do not consider the protein reorganization. Similarly, it is pointed out that microscopic potential of mean force calculations that do not evaluate the electrostatic barrier relative to the bulk water cannot be used to establish the origin of the electrostatic barrier. The relationship between the present study and calculations of pK(a)s in protein interiors is clarified, pointing out that approaches that are applied to study the aquaporin barrier should be validated by pK(a)s calculations. Such calculations also help to clarify the crucial role of solvation energies in establishing the barrier in aquaporins.     

Examining the promiscuous phosphatase activity of Pseudomonas aeruginosa arylsulfatase: a comparison to analogous phosphatases

Pseudomonas aeruginosa arylsulfatase (PAS) is a bacterial sulfatase capable of hydrolyzing a range of sulfate esters. Recently, it has been demonstrated to also show very high proficiency for phosphate ester hydrolysis. Such proficient catalytic promiscuity is significant, as promiscuity has been suggested to play an important role in enzyme evolution. Additionally, a comparative study of the hydrolyses of the p-nitrophenyl phosphate and sulfate monoesters in aqueous solution has demonstrated that despite superficial similarities, the two reactions proceed through markedly different transition states with very different solvation effects, indicating that the requirements for the efficient catalysis of the two reactions by an enzyme will also be very different (and yet they are both catalyzed by the same active site). This work explores the promiscuous phosphomonoesterase activity of PAS. Specifically, we have investigated the identity of the most likely base for the initial activation of the unusual formylglycine hydrate nucleophile (which is common to many sulfatases), and demonstrate that a concerted substrate-as-base mechanism is fully consistent with the experimentally observed data. This is very similar to other related systems, and suggests that, as far as the phosphomonoesterase activity of PAS is concerned, the sulfatase behaves like a "classical" phosphatase, despite the fact that such a mechanism is unlikely to be available to the native substrate (based on pK(a) considerations and studies of model systems). Understanding such catalytic versatility can be used to design novel artificial enzymes that are far more proficient than the current generation of designer enzymes.     

Physico-chemical properties of amphoteric, isoelectric, macroreticulate buffers

We report here the properties of a new family of resins possessing an amphoteric character and able to strongly buffer at their pI values. They have been adopted as carriers for growth of cells in tissue culture and for hydroponics (Righetti et al. 1991; J. Biotechnol. 17, 169-176) but it is to be expected that such resins could have interesting chromatographic applications. It has been found that such beads [made by incorporating a pK 6.2 weak acrylamido base and a pK 4.6 weak acrylamido acid in a 2:1 ratio (thus with a pI of 6.2) into a neutral polyacrylamide backbone], independently from their initial conditioning (acid- or base-washed), spontaneously seek their equilibrium position (pI value) upon washing off excess titrant. Thus, upon potentiometric titration, they are seen to buffer in both directions of the pH scale (contrary to the behaviour of a pure carboxyl or a pure amino surface, which will exhibit only unidirectional buffering power). From the behaviour of these amphoteric beads when polymerized in the absence or in the presence of salts (0.2 M NaCl), it is hypothesized that, for exerting buffering power, both the buffering ion and its counterion must be incorporated non-randomly in the chain, but as a couple or in close proximity. Upon random incorporation of the two ions, buffering power is lost.     

Calculation of pKas in RNA: on the structural origins and functional roles of protonated nucleotides

pK(a) calculations based on the Poisson-Boltzmann equation have been widely used to study proteins and, more recently, DNA. However, much less attention has been paid to the calculation of pK(a) shifts in RNA. There is accumulating evidence that protonated nucleotides can stabilize RNA structure and participate in enzyme catalysis within ribozymes. Here, we calculate the pK(a) shifts of nucleotides in RNA structures using numerical solutions to the Poisson-Boltzmann equation. We find that significant shifts are predicted for several nucleotides in two catalytic RNAs, the hairpin ribozyme and the hepatitis delta virus ribozyme, and that the shifts are likely to be related to their functions. We explore how different structural environments shift the pK(a)s of nucleotides from their solution values. RNA structures appear to use two basic strategies to shift pK(a)s: (a) the formation of compact structural motifs with structurally-conserved, electrostatic interactions; and (b) the arrangement of the phosphodiester backbone to focus negative electrostatic potential in specific regions.     

Recent advances in implicit solvent-based methods for biomolecular simulations

Implicit solvent-based methods play an increasingly important role in molecular modeling of biomolecular structure and dynamics. Recent methodological developments have mainly focused on the extension of the generalized Born (GB) formalism for variable dielectric environments and accurate treatment of nonpolar solvation. Extensive efforts in parameterization of GB models and implicit solvent force fields have enabled ab initio simulation of protein folding to native or near-native structures. Another exciting area that has benefited from the advances in implicit solvent models is the development of constant pH molecular dynamics methods, which have recently been applied to the calculations of protein pK(a) values and the studies of pH-dependent peptide and protein folding.     

The catalytic role of glutamate 151 in the leucine aminopeptidase from Aeromonas proteolytica

Glutamate 151 has been proposed to act as the general acid/base during the peptide hydrolysis reaction catalyzed by the co-catalytic metallohydrolase from Aeromonas proteolytica (AAP). However, to date, no direct evidence has been reported for the role of Glu-151 during catalytic turnover by AAP. In order to elucidate the catalytic role of Glu-151, altered AAP enzymes have been prepared in which Glu-151 has been substituted with a glutamine, an alanine, and an aspartate. The Michaelis constant (K(m)) does not change upon substitution to aspartate or glutamine, but the rate of the reaction changes drastically in the following order: glutamate (100% activity), aspartate (0.05%), glutamine (0.004%), and alanine (0%). Examination of the pH dependence of the kinetic constants k(cat) and K(m) revealed a change in the pK(a) of a group that ionizes at pH 4.8 in recombinant leucine aminopeptidase (rAAP) to 4.2 for E151D-AAP. The remaining pK(a) values at 5.2, 7.5, and 9.9 do not change. Proton inventory studies indicate that one proton is transferred in the rate-limiting step of the reaction at pH 10.50 for both rAAP and E151D-AAP, but at pH 6.50 two protons and general solvation effects are responsible for the observed effects in the reaction catalyzed by rAAP and E151D-AAP, respectively. Based on these data, Glu-151 is intrinsically involved in the peptide hydrolysis reaction catalyzed by AAP and can be assigned the role of a general acid and base.     

The pKa Cooperative: a collaborative effort to advance structure-based calculations of pKa values and electrostatic effects in proteins

The pK(a) Cooperative (http://www.pkacoop.org) was organized to advance development of accurate and useful computational methods for structure-based calculation of pK(a) values and electrostatic energies in proteins. The Cooperative brings together laboratories with expertise and interest in theoretical, computational, and experimental studies of protein electrostatics. To improve structure-based energy calculations, it is necessary to better understand the physical character and molecular determinants of electrostatic effects. Thus, the Cooperative intends to foment experimental research into fundamental aspects of proteins that depend on electrostatic interactions. It will maintain a depository for experimental data useful for critical assessment of methods for structure-based electrostatics calculations. To help guide the development of computational methods, the Cooperative will organize blind prediction exercises. As a first step, computational laboratories were invited to reproduce an unpublished set of experimental pK(a) values of acidic and basic residues introduced in the interior of staphylococcal nuclease by site-directed mutagenesis. The pK(a) values of these groups are unique and challenging to simulate owing to the large magnitude of their shifts relative to normal pK(a) values in water. Many computational methods were tested in this first Blind Prediction Challenge and critical assessment exercise. A workshop was organized in the Telluride Science Research Center to objectively assess the performance of many computational methods tested on this one extensive data set. This volume of Proteins: Structure, Function, and Bioinformatics introduces the pK(a) Cooperative, presents reports submitted by participants in the Blind Prediction Challenge, and highlights some of the problems in structure-based calculations identified during this exercise.     

Calculation of pK(a) in proteins with the microenvironment modulated-screened coulomb potential

The MM-SCP has been applied to predict pK(a) values of titratable residues in wild type and mutants of staphylococcal nuclease (SNase). The calculations were based on crystal structures made available by the Garcia-Moreno Laboratory. In the mutants, mostly deeply buried hydrophobic residues were replaced with ionizable residues, and thus their pK(a) values could be measured and calculated using various methods. The data set used here consisted of a set of WT SNase for which His pK(a) at several ionic strengths had been measured, a set of mutants for which measured pK(a) were available and a set of 11 mutants for which the measured pK(a) were not known at the time of calculation. For this latter set, blind predictions were submitted to the protein pK(a) cooperative, 2009 workshop at Telluride, where the results of the blind predictions were discussed (the RMSD of the submitted set was 1.10 pH units). The calculations on the structures with known pK(a) indicated that in addition to weaknesses of the method, structural issues were observed that led to larger errors (>1) in pK(a) predictions. For example, different crystallographic conditions or steric clashes can lead to differences in the local environment around the titratable residue, which can produce large differences in the calculated pK(a) . To gain further insight into the reliability of the MM-SCP, pK(a) of an extended set of 54 proteins belonging to several structural classes were carried out. Here some initial results from this study are reported to help place the SNase results in the appropriate context.     

Charge-charge interactions are key determinants of the pK values of ionizable groups in ribonuclease Sa (pI=3.5) and a basic variant (pI=10.2)

The pK values of the titratable groups in ribonuclease Sa (RNase Sa) (pI=3.5), and a charge-reversed variant with five carboxyl to lysine substitutions, 5K RNase Sa (pI=10.2), have been determined by NMR at 20 degrees C in 0.1M NaCl. In RNase Sa, 18 pK values and in 5K, 11 pK values were measured. The carboxyl group of Asp33, which is buried and forms three intramolecular hydrogen bonds in RNase Sa, has the lowest pK (2.4), whereas Asp79, which is also buried but does not form hydrogen bonds, has the most elevated pK (7.4). These results highlight the importance of desolvation and charge-dipole interactions in perturbing pK values of buried groups. Alkaline titration revealed that the terminal amine of RNase Sa and all eight tyrosine residues have significantly increased pK values relative to model compounds.A primary objective in this study was to investigate the influence of charge-charge interactions on the pK values by comparing results from RNase Sa with those from the 5K variant. The solution structures of the two proteins are very similar as revealed by NMR and other spectroscopic data, with only small changes at the N terminus and in the alpha-helix. Consequently, the ionizable groups will have similar environments in the two variants and desolvation and charge-dipole interactions will have comparable effects on the pK values of both. Their pK differences, therefore, are expected to be chiefly due to the different charge-charge interactions. As anticipated from its higher net charge, all measured pK values in 5K RNase are lowered relative to wild-type RNase Sa, with the largest decrease being 2.2 pH units for Glu14. The pK differences (pK(Sa)-pK(5K)) calculated using a simple model based on Coulomb's Law and a dielectric constant of 45 agree well with the experimental values. This demonstrates that the pK differences between wild-type and 5K RNase Sa are mainly due to changes in the electrostatic interactions between the ionizable groups. pK values calculated using Coulomb's Law also showed a good correlation (R=0.83) with experimental values. The more complex model based on a finite-difference solution to the Poisson-Boltzmann equation, which considers desolvation and charge-dipole interactions in addition to charge-charge interactions, was also used to calculate pK values. Surprisingly, these values are more poorly correlated (R=0.65) with the values from experiment. Taken together, the results are evidence that charge-charge interactions are the chief perturbant of the pK values of ionizable groups on the protein surface, which is where the majority of the ionizable groups are positioned in proteins.     

Progress in the prediction of pKa values in proteins

The pK(a) -cooperative aims to provide a forum for experimental and theoretical researchers interested in protein pK(a) values and protein electrostatics in general. The first round of the pK(a) -cooperative, which challenged computational labs to carry out blind predictions against pK(a) s experimentally determined in the laboratory of Bertrand Garcia-Moreno, was completed and results discussed at the Telluride meeting (July 6-10, 2009). This article serves as an introduction to the reports submitted by the blind prediction participants that will be published in a special issue of PROTEINS: Structure, Function and Bioinformatics. Here, we briefly outline existing approaches for pK(a) calculations, emphasizing methods that were used by the participants in calculating the blind pK(a) values in the first round of the cooperative. We then point out some of the difficulties encountered by the participating groups in making their blind predictions, and finally try to provide some insights for future developments aimed at improving the accuracy of pK(a) calculations.     

Acidic-basic properties of three alanine-based peptides containing acidic and basic side chains: comparison between theory and experiment

The purpose of this work was to evaluate the effect of the nature of the ionizable end groups, and the solvent, on their acid-base properties in alanine-based peptides. Hence, the acid-base properties of three alanine-based peptides: Ac-KK-(A)(7)-KK-NH(2) (KAK), Ac-OO-(A)(7)-DD-NH(2) (OAD), Ac-KK-(A)(7)-EE-NH(2) (KAE), where A, D, E, K, and O denote alanine, aspartic acid, glutamic acid, lysine, and ornithine, respectively, were determined in water and in methanol by potentiometry. With the availability of these data, the ability of two theoretical methods to simulate pH-metric titration of those peptides was assessed: (i) the electrostatically driven Monte Carlo method with the ECEPP/3 force field and the Poisson-Boltzmann approach to compute solvation energy (EDMC/PB/pH), and (ii) the molecular dynamics method with the AMBER force field and the Generalized Born model (MD/GB/pH). For OAD and KAE, pK(a1) and pK(a2) correspond to the acidic side chains. For all three compounds in both solvents, the pK(a1) value is remarkably lower than the pK(a) of a compound modeling the respective isolated side chain, which can be explained by the influence of the electrostatic field from positively charged ornithine or lysine side chains. The experimental titration curves are reproduced well by the MD/GB/pH approach, the agreement being better if restraints derived from NMR measurements are incorporated in the conformational search. Poorer agreement is achieved by the EDMC/PB/pH method.     

A low pKa cysteine at the active site of mouse methionine sulfoxide reductase A

Methionine sulfoxide reductase A is an essential enzyme in the antioxidant system which scavenges reactive oxygen species through cyclic oxidation and reduction of methionine and methionine sulfoxide. Recently it has also been shown to catalyze the reverse reaction, oxidizing methionine residues to methionine sulfoxide. A cysteine at the active site of the enzyme is essential for both reductase and oxidase activities. This cysteine has been reported to have a pK(a) of 9.5 in the absence of substrate, decreasing to 5.7 upon binding of substrate. Using three independent methods, we show that the pK(a) of the active site cysteine of mouse methionine sulfoxide reductase is 7.2 even in the absence of substrate. The primary mechanism by which the pK(a) is lowered is hydrogen bonding of the active site Cys-72 to protonated Glu-115. The low pK(a) renders the active site cysteine susceptible to oxidation to sulfenic acid by micromolar concentrations of hydrogen peroxide. This characteristic supports a role for methionine sulfoxide reductase in redox signaling.     

Contribution of aromatic-aromatic interactions to the anomalous pK(a) of tyrosine-9 and the C-terminal dynamics of glutathione S-transferase A1-1

Most cytosolic glutathione S-transferases (GSTs) exploit a hydrogen bond between an active site Tyr and the bound glutathione (GSH) cofactor to lower the pK(a) of the GSH and generate the nucleophilic thiolate anion, GS(-). In human (hGSTA1-1) and rat (rGSTA1-1) homologues, the active site Tyr-9 has a low pK(a) of 8.1-8.3, for which the functional significance is unknown. Crystal structures of GSTA1-1 suggest that weakly polar interactions between the electropositive ring edge of Phe-10 and the pi-cloud of Tyr-9, in the apoenzyme, could stabilize the tyrosinate anion and also modulate the pK(a) of GSH. Upon binding a product GSH conjugate, Phe-10 moves away from Tyr-9, allowing the highly dynamic C-terminus to "close" over the active site. To explore the role of Phe-10 in modulating the Tyr-9 pK(a) and in ligand binding, rGSTA1-1 mutants F10Y, F10L, and F10A were characterized. The pK(a)s of Tyr-9 in the apoenzymes were 8.2 +/- 0.2, 8.7 +/- 0.2, and 9.3 +/- 0.1, respectively, for F10Y, F10L, and F10A, compared to 8.3 +/- 0.2 for the "wild type". The experimentally determined pK(a)s qualitatively paralleled the energies required to remove a proton predicted by ab initio calculations using model compounds constrained to the coordinates of rGSTA1-1. The pK(a) of GSH in the binary complex was significantly less affected by these substitutions. In contrast, F220I and F220Y C-terminal mutations caused the pK(a) of Tyr-9 to decrease modestly. For the binary complex with S-hexyl-GSH, which induces the "closed" conformation, Tyr-9 retains a low pK(a) and the Phe-10 substitutions have significant effects. Presumably, Phe-10 plays a critical structural role in stabilizing the closed conformation. The mutations F10L and F10A also slowed the rate of GSH conjugate binding by 10-20-fold, as measured by stopped-flow fluorescence. The effects of Phe-10 substitution were large for both steps of the biphasic binding reaction, suggesting the importance of aromatic interactions throughout the reaction coordinate. A unified view of the C-terminal dynamics of GSTA1-1 is discussed, which emphasizes the coupling between Tyr-9 ionization, active site solvation, and C-terminal dynamics.     

UV Raman evidence of a tyrosine in apo-human serum transferrin with a low pK(a) that is elevated upon binding of sulphate

The binding of sulphate to human serum apo-transferrin has been examined by ultraviolet absorption and ultraviolet resonance Raman difference spectroscopies between pH 6.0 and 9.0. The ultraviolet absorption data reveals a negative feature at 245 nm that increases in magnitude with pH, with an apparent pK(a) of 7.57, which the Raman difference data reveals to be due to tyrosine. The pK(a) of this tyrosine is unusually low and is measured at 7.84 by the Raman difference method and is elevated to greater than 9.0 upon addition of sulphate. Previous studies on the N-lobe imply that Tyr 188 is the tyrosine with a low pK(a) and also that Arg 124 is the primary binding site for the sulphate. The functional relevance may be that with sulphate bound, both carbonate binding and the deprotonation of Tyr will be disfavoured, and as a result so is iron binding.