Correlation of carboxylic acid pKa to protein binding and antibacterial activity of a novel class of bacterial translation inhibitors

Previously we reported the discovery and initial optimization of a novel anthranilic acid derived class of antibacterial agents which suffered from extensive protein binding. This report describes efforts directed toward understanding the relationship of the acidity of the carboxylic acid with the extent of protein binding. The pK(a) of the acid was modified via the synthesis of a number of anthranilic acid analogs which vary the aromatic ring substituent at the 4-position. The pK(a) and HSA binding constants have been determined for each of the analogs. Our results indicate a correlation between pK(a) and HSA K(d). The physical properties and antibacterial activities will be discussed as well as how these results help address the protein binding issue with this series of compounds.     

An important role of G638 in the cis-cleavage reaction of the Neurospora VS ribozyme revealed by a novel nucleotide analog incorporation method

We describe a chemical coupling procedure that allows joining of two RNAs, one of which contains a site-specific base analog substitution, in the absence of divalent ions. This method allows incorporation of nucleotide analogs at specific positions even into large, cis-cleaving ribozymes. Using this method we have studied the effects of substitution of G638 in the cleavage site loop of the VS ribozyme with a variety of purine analogs having different functional groups and pK(a) values. Cleavage rate versus pH profiles combined with kinetic solvent isotope experiments indicate an important role for G638 in proton transfer during the rate-limiting step of the cis-cleavage reaction.     

Study of the inhibition of neutral phosphatase from Pseudomonadaceae by derivatives of its substrates

The inhibition of neutral phosphatase isolated from the bacteria of the Pseudomonadaceae family by various fragments of the enzyme-hydrolyzed R-O-PO3H2 substrates, inorganic orthophosphate (KH2PO4) and its analogs as well as by adenine, adenosine, alcohols, sugars and amino acids, was studied. It was demonstrated that among other compounds tested only the orthophosphoric acid anions (H2PO4-) exhibit the properties of strong associative inhibitors (K1Vi = 4.35.10(-6)M of the enzyme. The pH dependence of the Michaelis constant [pKm0 = f(pH)] and the inhibition constant for phosphatase by potassium orthophosphate [pK1Vi(KH2PO4) = f(pH)] was studied. The presence in the enzyme active center of a carboxylic (pK = 4.3 +/- 0.1) (presumably, glutamine) and an imidazole (pK = 7.15 +/- 0.1) amino acid residues was postulated. The data obtained were compared to those for neutral, alkaline and acid phosphatases.     

DAT/SERT selectivity of flexible GBR 12909 analogs modeled using 3D-QSAR methods

The dopamine reuptake inhibitor GBR 12909 (1-{2-[bis(4-fluorophenyl)methoxy]ethyl}-4-(3-phenylpropyl)piperazine, 1) and its analogs have been developed as tools to test the hypothesis that selective dopamine transporter (DAT) inhibitors will be useful therapeutics for cocaine addiction. This 3D-QSAR study focuses on the effect of substitutions in the phenylpropyl region of 1. CoMFA and CoMSIA techniques were used to determine a predictive and stable model for the DAT/serotonin transporter (SERT) selectivity (represented by pK(i) (DAT/SERT)) of a set of flexible analogs of 1, most of which have eight rotatable bonds. In the absence of a rigid analog to use as a 3D-QSAR template, six conformational families of analogs were constructed from six pairs of piperazine and piperidine template conformers identified by hierarchical clustering as representative molecular conformations. Three models stable to y-value scrambling were identified after a comprehensive CoMFA and CoMSIA survey with Region Focusing. Test set correlation validation led to an acceptable model, with q(2)=0.508, standard error of prediction=0.601, two components, r(2)=0.685, standard error of estimate=0.481, F value=39, percent steric contribution=65, and percent electrostatic contribution=35. A CoMFA contour map identified areas of the molecule that affect pK(i) (DAT/SERT). This work outlines a protocol for deriving a stable and predictive model of the biological activity of a set of very flexible molecules.     

The pH dependence of red cell membrane transport of titratable anions studied by NMR spectroscopy

The effects of varying extracellular pH on the rates of uptake of titratable anions by human erythrocytes under conditions of constant intracellular pH have been determined for a series of highly related anions, the phosphate "analogs." These compounds are simply substituted phosphorus oxyacids, differing in the number and acidity of titratable protons: phosphate (HPO4(2-), pKa 6.8); phosphite (HPO3(2-), pKa 6.4); hypophosphite (H2PO2-); methylphosphonate ((CH3)PO3(2-), pKa 7.4); dimethylphosphinate ((CH3)2PO2-); fluorophosphate [PO3F2-, pKa 4.7); and thiophosphate (HSPO3(2-), pKa 5.5). Suspensions of intact, Cl(-)-loaded erythrocytes (intracellular pH, 7.2) were incubated at 37 degrees C in isotonic buffers (pH 4-8) containing 60 mM phosphate analog for specified time intervals, whereupon influx was halted by the addition of 1 mM 4-acetamido-4'-isothiocyanatostilbene-2,2'-disulfonic acid (SITS), an inhibitor of anion exchange. The intracellular anion concentrations were determined from 31P or 19F nuclear magnetic resonance spectra from the erythrocyte suspensions. The influx rates for the titratable phosphate analogs exhibited bimodal pH dependence, reaching maximal levels at pH values that increased with increasing anion pK. This pH-dependent behavior is consistent with a transport channel that contains a titratable regulatory site which interacts with the translocated anion. Based upon the Henderson-Hasselbalch equation, the probability that a titratable anion will have an electric charge of equal magnitude to that of the titratable carrier is highest at a pH value exactly midway between the pK of the regulatory site and that of the anion. The pH maxima observed for the phosphate analogs indicate a pK for this site of 5.5 at 37 degrees C. Intracellular pH changes associated with influx indicated that transport of the "fast" anion phosphite is largely in monoionized form. Intracellular pH changes associated with transport of slow anions were predominantly determined by partial ionic equilibrium effects and did not indicate the ionization state of the transported anion.     

pH dependent conformational changes in the T- and R-states of insulin in solution: circular dichroic studies in the pH range of 6 to 10.

Zinc insulin hexamer has been shown to undergo a phenol-induced T6 to R6 conformational transition in solution. Our circular dichroic (CD) studies demonstrate that insulin undergoes pH-dependent conformational changes over the pH range of 6-10 in the T-state and in the R- state. In order to determine which specific amino acid residues may be responsible for these pH-dependent changes, a series of insulin analogs were utilized. In the T-state, the pH dependent CD changes monitored in the far UV region have a pK of 8.2 and appear to be related to the titration of the A1-Gly amino group. Using the near UV CD a second pH-dependent conformational change was detected with a pK of 7.5 in the T-state. 1H N.M.R. studies suggest that B5-His may be responsible for this conformational transition. In the presence of m-cresol (R-state), the pK value was found to be 6.9. During this titration, the increased ellipticity for the R-state is diminishing as pH decreases from pH 8 to 6, and no difference in ellipticity was observed at 255 nm between T- and R-states at pH 6. Therefore, this may be due to the transition from the R back to the T-state.     

Binding of substrate analogs to hen egg-white lysozyme with an ester linkage between Glu 35 and Trp 108.

The interactions of the substrate analogs beta-methyl-GlcNAc, (GlcNAc)2, and (GlcNAc)3 with hen egg-white lysozyme [EC 3.2.1.17] in which an ester linkage had been formed between Glu 35 and Trp 108 (108 ester lysozyme), were studied by the circular dichroic and fluorescence techniques, and were compared with those for intact lysozyme. The binding constants of beta-methyl-GlcNAc and (GlcNAc)2 to 108 ester lysozyme were essentially the same as those for intact lysozyme in the pH range of 1 to 5. Above pH 5, the binding constants of these saccharides to 108 ester lysozyme did not change with pH, while the binding constants to intact lysozyme decreased. This indicates that Glu 35 (pK 6.0 in intact lysozyme) participates in the binding of these saccharides. The extent and direction of the pK shifts of Asp 52 (pK 3.5), Asp 48 (pK 4.4), and Asp 66 (pK 1.3) observed when beta-methyl-GlcNAc is bound to 108 ester lysozyme were the same as those for intact lysozyme. The participation of Asp 101 and Asp 66 in the binding of (GlcNAc)2 to 108 ester lysozyme was also the same as that for intact lysozyme. These findings indicate that the conformations of subsites B and C are not changed by the formation of the ester linkage. On the other hand, the binding constants of (GlcNAc)3 to 108 ester lysozyme were higher than those for intact lysozyme at all pH values studied. This result is interpreted in terms of an increase in the affinity for a GlcNAc residue of subsite D, which is situated near the esterified Glu 35.     

Proton transfer dynamics of GART: the pH-dependent catalytic mechanism examined by electrostatic calculations

The enzyme glycinamide ribonucleotide transformylase (GART) catalyzes the transfer of a formyl group from formyl tetrahydrofolate (fTHF) to glycinamide ribonucleotide (GAR), a process that is pH-dependent with pK(a) of approximately 8. Experimental studies of pH-rate profiles of wild-type and site-directed mutants of GART have led to the proposal that His108, Asp144, and GAR are involved in catalysis, with His108 being an acid catalyst, while forming a salt bridge with Asp144, and GAR being a nucleophile to attack the formyl group of fTHF. This model implied a protonated histidine with pK(a) of 9.7 and a neutral GAR with pK(a) of 6.8. These proposed unusual pK(a)s have led us to investigate the electrostatic environment of the active site of GART. We have used Poisson-Boltzmann-based electrostatic methods to calculate the pK(a)s of all ionizable groups, using the crystallographic structure of a ternary complex of GART involving the pseudosubstrate 5-deaza-5,6,7,8-THF (5dTHF) and substrate GAR. Theoretical mutation and deletion analogs have been constructed to elucidate pairwise electrostatic interactions between key ionizable sites within the catalytic site. Also, a construct of a more realistic catalytic site including a reconstructed pseudocofactor with an attached formyl group, in an environment with optimal local van der Waals interactions (locally minimized) that imitates closely the catalytic reactants, has been used for pK(a) calculations. Strong electrostatic coupling among catalytic residues His108, Asp144, and substrate GAR was observed, which is extremely sensitive to the initial protonation and imidazole ring flip state of His108 and small structural changes. We show that a proton can be exchanged between GAR and His108, depending on their relative geometry and their distance to Asp144, and when the proton is attached on His108, catalysis could be possible. Using the formylated locally minimized construct of GART, a high pK(a) for His108 was calculated, indicating a protonated histidine, and a low pK(a) for GAR(NH(2)) was calculated, indicating that GAR is in neutral form. Our results are in qualitative agreement with the current mechanistic picture of the catalytic process of GART deduced from the experimental data, but they do not reproduce the absolute magnitude of the pK(a)s extracted from fits of k(cat)-pH profiles, possibly because the static time-averaged crystallographic structure does not describe adequately the dynamic nature of the catalytic site during binding and catalysis. In addition, a strong effect on the pK(a) of GAR(NH(2)) is produced by the theoretical mutations of His108Ala and Asp144Ala, which is not in agreement with the observed insensitivity of the pK(a) of GAR(NH(2)) modeled from the experimental data using similar mutations. Finally, we show that important three-way electrostatic interactions between highly conserved His137, with His108 and Asp144, are responsible for stabilizing the electrostatic microenvironment of the catalytic site. In conclusion, our data suggest that further detailed computational and experimental work is necessary.     

Role of an active site adenine in hairpin ribozyme catalysis

The hairpin ribozyme is a small catalytic RNA that accelerates reversible cleavage of a phosphodiester bond. Structural and mechanistic studies suggest that divalent metals stabilize the functional structure but do not participate directly in catalysis. Instead, two active site nucleobases, G8 and A38, appear to participate in catalytic chemistry. The features of A38 that are important for active site structure and chemistry were investigated by comparing cleavage and ligation reactions of ribozyme variants with A38 modifications. An abasic substitution of A38 reduced cleavage and ligation activity by 14,000-fold and 370,000-fold, respectively, highlighting the critical role of this nucleobase in ribozyme function. Cleavage and ligation activity of unmodified ribozymes increased with increasing pH, evidence that deprotonation of some functional group with an apparent pK(a) value near 6 is important for activity. The pH-dependent transition in activity shifted by several pH units in the basic direction when A38 was substituted with an abasic residue, or with nucleobase analogs with very high or low pK(a) values that are expected to retain the same protonation state throughout the experimental pH range. Certain exogenous nucleobases that share the amidine group of adenine restored activity to abasic ribozyme variants that lack A38. The pH dependence of chemical rescue reactions also changed according to the intrinsic basicity of the rescuing nucleobase, providing further evidence that the protonation state of the N1 position of purine analogs is important for rescue activity. These results are consistent with models of the hairpin ribozyme catalytic mechanism in which interactions with A38 provide electrostatic stabilization to the transition state.     

Structure-activity relationship of ibogaine analogs interacting with nicotinic acetylcholine receptors in different conformational states

The interaction of ibogaine analogs with nicotinic acetylcholine receptors (AChRs) in different conformational states was studied by functional and structural approaches. The results established that ibogaine analogs: (a) inhibit (±)-epibatidine-induced Ca2? influx in human embryonic muscle AChRs with the following potency sequence (IC(50) in μM): (±)-18-methylaminocoronaridine (5.9±0.3)～(±)-18-methoxycoronaridine (18-MC) (6.8±0.8)>(-)-ibogaine (17±3)～(+)-catharanthine (20±1)>(±)-albifloranine (46±13), (b) bind to the [3H]TCP binding site with higher affinity when the Torpedo AChR is in the desensitized state compared to that in the resting state. Similar results were obtained using [3H]18-MC. These and docking results suggest a steric interaction between TCP and ibogaine analogs for the same site, (c) enhance [3H]cytisine binding to resting but not to desensitized AChRs, with desensitizing potencies (apparent EC??) that correlate very well with the pK(i) values in the desensitized state, and (d) there are good bilinear correlations between the ligand molecular volumes and their affinities in the desensitized and resting states, with an optimal volume of ～345 ?3 for the ibogaine site. These results indicate that the size of the binding sites for ibogaine analogs, located between the serine and nonpolar rings and shared with TCP, is an important structural feature for binding and for inducing desensitization.     

Multiple sugar binding sites in alpha-glucosidase

Twenty-five analogs of D-glucose were examined as reversible inhibitors of yeast alpha-glucosidase (EC 3.2.1.20). The K(i) values range from 0.38 mM for 6-deoxy-D-glucose (quinovose) to 1.0 M for D-lyxose at pH=6.3 (0.1 M NaCl, 25 degrees ). All the monosaccharides and the three disaccharides (maltose, isomaltose and alpha,alpha-trehalose) were found to be linear competitive inhibitors with respect to alpha-p-nitrophenyl glucoside (pNPG) hydrolysis. Multiple inhibition analysis reveals that there are at least three monosaccharide binding sites on the enzyme. One of these can be occupied by glucose [K(i)=1.8(+/-0.1) mM], one by D-galactose [K(i)=164(+/-11) mM] and one by D-mannose [K(i)=120(+/-9) mM]. The pH dependence for glucose binding closely follows that of V/K [pK(a1)=5.55(+/-0.15), pK(a2)=6.79(+/-0.15)], but the binding of mannose does not. Although the glucose subsite can be occupied simultaneously with the mannose or galactose subsites in the enzyme-product complex, no transglucosylation can be detected between pNPG and either mannose or galactose. This suggests that neither of these nonglucose subsites can be occupied in a productive manner in the covalent glucosyl-enzyme intermediate.     

pK values of histidine residues in ribonuclease Sa: effect of salt and net charge

The primary goal of this study was to gain a better understanding of the effect of environment and ionic strength on the pK values of histidine residues in proteins. The salt-dependence of pK values for two histidine residues in ribonuclease Sa (RNase Sa) (pI=3.5) and a variant in which five acidic amino acids have been changed to lysine (5K) (pI=10.2) was measured and compared to pK values of model histidine-containing peptides. The pK of His53 is elevated by two pH units (pK=8.61) in RNase Sa and by nearly one pH unit (pK=7.39) in 5K at low salt relative to the pK of histidine in the model peptides (pK=6.6). The pK for His53 remains elevated in 1.5M NaCl (pK=7.89). The elevated pK for His53 is a result of screenable electrostatic interactions, particularly with Glu74, and a non-screenable hydrogen bond interaction with water. The pK of His85 in RNase Sa and 5K is slightly below the model pK at low salt and merges with this value at 1.5M NaCl. The pK of His85 reflects mainly effects of long-range Coulombic interactions that are screenable by salt. The tautomeric states of the neutral histidine residues are changed by charge reversal. The histidine pK values in RNase Sa are always higher than the pK values in the 5K variant. These results emphasize that the net charge of the protein influences the pK values of the histidine residues. Structure-based pK calculations capture the salt-dependence relatively well but are unable to predict absolute histidine pK values.     

Deacylation kinetics of gamma-chymotrypsin in solution and in the crystal

The rate of catalyzed hydrolysis of the acyl enzyme analogs indolacryloyl-gamma-chymotrypsin and furylacryloyl-gamma-chymotrypsin in the crystal has been measured by single crystal microspectrophotometry and compared with the rate of catalyzed hydrolysis of acyl-gamma-chymotrypsin in solution and acyl-alpha-chymotrypsin both in solution and in the crystal. The maximal deacylation rate is the same for both species and independent of the physical state. However, the pH dependence of the deacylation rate of crystalline acyl-gamma-chymotrypsin shows a 0.9 unit shift in the pK of the catalytic system which is unique and probably consequent to specific lattice interactions.     

Structure-activity relationship of [Nphe1]-NC-(1-13)-NH2, a pure and selective nociceptin/orphanin FQ receptor antagonist.

A series of analogs of the ORL1 receptor antagonist [Nphe1]-NC(1-13)-NH2 was prepared and tested for agonistic and antagonistic activities in the mouse vas deferens, a preparation that shows high sensitivity to nociceptin and related peptides. The purpose of this study was to determine the role of the aromatic residue at the N-terminal for antagonism and eventually identify compounds with improved potency. Results indicated that all 23 compounds are inactive as agonists, and the antagonistic potency of the initial template [Nphe1]-NC(1-13)-NH2 is high (pKB 6.43) compared with those of all other compounds except [(S)(betaMe)Nphe1]NC(1-13)-NH2 (pK(B) 6.48). The other 22 compounds can be divided into two groups: 10 show antagonistic potencies (pK(B)) ranging from 5.30 to 5.86, whereas the other 12 compounds are inactive. This study clearly shows that the aromatic ring of Nphe is very critical for the interaction with the ORL1 receptor and can not be enlarged or sterically modified without significant loss of antagonistic potency.     

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.     

Lead optimization of 2-(piperidin-3-yl)-1H-benzimidazoles: identification of 2-morpholin- and 2-thiomorpholin-2-yl-1H-benzimidazoles as selective and CNS penetrating H₁-antihistamines for insomnia

The structure-activity relationships of 2-(piperidin-3-yl)-1H-benzimidazoles, 2-morpholine and 2-thiomorpholin-2-yl-1H-benzimidazoles are described. In the lead optimization process, the pK(a) and/or logP of benzimidazole analogs were reduced either by attachment of polar substituents to the piperidine nitrogen or incorporation of heteroatoms into the piperidine heterocycle. Compounds 9a and 9b in the morpholine series and 10g in the thiomorpholine series demonstrated improved selectivity and CNS profiles compared to lead compound 2 and these are potential candidates for evaluation as sedative hypnotics.     

Thermodynamic analysis of ionizable groups involved in the catalytic mechanism of human matrix metalloproteinase 7 (MMP-7)

Human matrix metalloproteinase 7 (MMP-7) exhibits a broad bell-shaped pH-dependence with the acidic and alkaline pK(e) (pK(e1) and pK(e2)) values of about 4 and 10. In this study, we estimated the ionizable groups involved in its catalytic mechanism by thermodynamic analysis. pK(a) of side chains of L-Asp, L-Glu, L-His, L-Cys, L-Tyr, L-Lys, and L-Arg at 25-45°C were determined by the pH titration of amino-acid solutions, from which their enthalpy changes, ?H°, of deprotonation were calculated. pK(e1) and pK(e2) of MMP-7 at 15-45°C were determined in the hydrolysis of (7-methoxycoumarin-4-yl)acetyl-L-Pro-L-Leu-Gly-L-Leu-[N(3)-(2,4-dinitrophenyl)-L-2,3-diaminopropionyl]-L-Ala-L-Arg-NH(2), from which ?H(o) for pK(e1) and pK(e2) was calculated. The ?H(o) for pK(e1) (-20.6±6.1kJmol(-1)) was similar to that for L-Glu (-23.6±5.8kJmol(-1)), and the ?H(o) for pK(e2) (89.9±4.0kJmol(-1)) was similar to those for L-Arg (87.6±5.5kJmol(-1)) and L-Lys (70.4±4.4kJmol(-1)). The mutation of the active-site residue Glu198 into Ala completely abolished the activity, suggesting that Glu198 is the ionizable group for pK(e1). On the other hand, no arginine or lysine residues are found in the active site of MMP-7. We proposed a possibility that a protein-bound water is the ionizable group for pK(e2).     

Kinetic properties and inhibition of the dimeric dUTPase-dUDPase from Leishmania major

Kinetic properties of the dimeric enzyme dUTPase from Leishmania major were studied using a continuous spectrophotometric method. dUTP was the natural substrate and dUMP and PPi the products of the hydrolysis. The trypanosomatid enzyme exhibited a low K(m) value for dUTP (2.11 microM), a k(cat) of 49 s(-1), strict Michaelis-Menten kinetics and is a potent catalyst of dUDP hydrolysis, whereas in other dUTPases described, this compound acts as a competitive inhibitor. Discrimination is achieved for the base and sugar moiety showing specificity constants for different dNTPs similar to those of bacterial, viral, and human enzymes. In the alkaline range, the K(m) for dUTP increases with the dissociation of ionizable groups showing pK(a) values of 8.8, identified as the uracil moiety of dUTP and 10, whereas in the acidic range, K(m) is regulated by an enzyme residue exhibiting a pK(a) of 7.1. Activity is strongly inhibited by the nucleoside triphosphate analog alpha-beta-imido-dUTP, indicating that the enzyme can bind triphosphate analogs. The existence of specific inhibition and the apparent structural and kinetic differences (reflected in different binding strength of dNTPs) with other eukaryotic dUTPases suggest that the present enzyme might be exploited as a target for new drugs against leishmaniasis.     

Self-association of glutamic acid-rich fusion peptide analogs of influenza hemagglutinin in the membrane-mimic environments: effects of positional difference of glutamic acids on side chain ionization constant and intra- and inter-peptide interactions deduced from NMR and gel electrophoresis measurements

Two glutamic acid-rich fusion peptide analogs of influenza hemagglutinin were synthesized to study the organization of the charged peptides in the membranous media. Fluorescence and gel electrophoresis experiments suggested a loose association between the monomers in the vesicles. A model was built which showed that a positional difference of 3, 7 and 4, 8 results in the exposure of Glu3 and Glu7 side chains to the apolar lipidic core. Supportive results include: first, pK(a) values of two pH units higher than reference value in aqueous medium for Glu3 and Glu7 CgammaH, whereas the deviation of pK(a) from the reference value for Glu4 and Glu8 CgammaH is substantially smaller; second, Hill coefficients of titration shift of these protons indicate anti-cooperativity for Glu3 and Glu7 side chain protons but less so for Glu4 and Glu8, implying a strong electrostatic interaction between Glu3 and Glu7 possibly resulting from their localization in an apolar environment; third, positive and larger titration shift for NH of Glu3 is observed compared to that of Glu4, suggesting stronger hydrogen bond between the NH and the carboxylic group of Glu3 than that of Glu4, consistent with higher degree of exposure to hydrophobic medium for the side chain of Glu3.     

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.