Structural origins of high apparent dielectric constants experienced by ionizable groups in the hydrophobic core of a protein

The side chains of Lys66, Asp66, and Glu66 in staphylococcal nuclease are fully buried and surrounded mainly by hydrophobic matter, except for internal water molecules associated with carboxylic oxygen atoms. These ionizable side chains titrate with pK_a values of 5.7, 8.8, and 8.9, respectively. To reproduce these pK_a values with continuum electrostatics calculations, we treated the protein with high dielectric constants. We have examined the structural origins of these high apparent dielectric constants by using NMR spectroscopy to characterize the structural response to the ionization of these internal side chains. Substitution of Val66 with Lys66 and Asp66 led to increased conformational fluctuations of the microenvironments surrounding these groups, even under pH conditions where Lys66 and Asp66 are neutral. When Lys66, Asp66, and Glu66 are charged, the proteins remain almost fully folded, but resonances for a few backbone amides adjacent to the internal ionizable residues are broadened. This suggests that the ionization of the internal groups promotes a local increase in dynamics on the intermediate timescale, consistent with either partial unfolding or increased backbone fluctuations of helix 1 near residue 66, or, less likely, with increased fluctuations of the charged side chains at position 66. These experiments confirm that the high apparent dielectric constants reported by internal Lys66, Asp66, and Glu66 reflect localized changes in conformational fluctuations without incurring detectable global structural reorganization. To improve structure-based pK_a calculations in proteins, we will need to learn how to treat this coupling between ionization of internal groups and local changes in conformational fluctuations explicitly.     

Are acidic and basic groups in buried proteins predicted to be ionized?

Ionizable residues play essential roles in proteins, modulating protein stability, fold and function. Asp, Glu, Arg, and Lys make up about a quarter of the residues in an average protein. Multi-conformation continuum electrostatic (MCCE) calculations were used to predict the ionization states of all acidic and basic residues in 490 proteins. Of all 36,192 ionizable residues, 93.5% were predicted to be ionized. Thirty-five percent have lost 4.08 kcal/mol solvation energy (DeltaDeltaG(rxn)) sufficient to shift a pK(a) by three pH units in the absence of other interactions and 17% have DeltaDeltaG(rxn) sufficient to shift pK(a) by five pH units. Overall 85% of these buried residues (DeltaDeltaG(rxn)>5DeltapK units) are ionized, including 92% of the Arg, 86% of the Asp, 77% of the Glu, and 75% of the Lys. Ion-pair interactions stabilize the ionization of both acids and bases. The backbone dipoles stabilize anions more than cations. The interactions with polar side-chains are also different for acids and bases. Asn and Gln stabilize all charges, Ser and Thr stabilize only acids while Tyr rarely stabilize Lys. Thus, hydroxyls are better hydrogen bond donors than acceptors. Buried ionized residues are more likely to be conserved than those on the surface. There are 3.95 residues buried per 100 residues in an average protein.     

Carbanions from decarboxylation of orotate analogs: stability and mechanistic implications

The pKa's of the 6-CH groups of 1,3-dimethyluracil, N-methyl-2-pyridone, and N-methyl-4-pyridone were determined through their reactions with bases derived from carbon acids with known pKa and the reactions of their corresponding carbanions with the carbon acids. No correlation between the stability of the carbanions and the rate of decarboxylation of corresponding carboxylic acids was found.     

Structural and mechanistic insights into unusual thiol disulfide oxidoreductase

Cytoplasmic desulfothioredoxin (Dtrx) from the anaerobe Desulfovibrio vulgaris Hildenborough has been identified as a new member of the thiol disulfide oxidoreductase family. The active site of Dtrx contains a particular consensus sequence, CPHC, never seen in the cytoplasmic thioredoxins and generally found in periplasmic oxidases. Unlike canonical thioredoxins (Trx), Dtrx does not present any disulfide reductase activity, but it presents instead an unusual disulfide isomerase activity. We have used NMR spectroscopy to gain insights into the structure and the catalytic mechanism of this unusual Dtrx. The redox potential of Dtrx (-181 mV) is significantly less reducing than that of canonical Trx. A pH dependence study allowed the determination of the pK(a) of all protonable residues, including the cysteine and histidine residues. Thus, the pK(a) values for the thiol group of Cys(31) and Cys(34) are 4.8 and 11.3, respectively. The His(33) pK(a) value, experimentally determined for the first time, differs notably as a function of the redox states, 7.2 for the reduced state and 4.6 for the oxidized state. These data suggest an important role for His(33) in the molecular mechanism of Dtrx catalysis that is confirmed by the properties of mutant DtrxH33G protein. The NMR structure of Dtrx shows a different charge repartition compared with canonical Trx. The results presented are likely indicative of the involvement of this protein in the catalysis of substrates specific of the anaerobe cytoplasm of DvH. The study of Dtrx is an important step toward revealing the molecular details of the thiol-disulfide oxidoreductase catalytic mechanism.     

Protonation state of Asp30 exerts crucial influence over surface loop rearrangements responsible for NO release in nitrophorin 4

pK(a) values of ionizable residues were calculated for the crystal structures describing the pH and NO binding dependant conformations of nitrophorin 4, a pH sensitive NO carrier heme protein. Comparison of resultant H-bonding patterns allowed the identification of the amino acids that take part in signaling pH change. We carried out MD simulations to show that the protonation state of Asp30, buried in the closed conformation, is crucial for maintaining the tight packed conformation of the closed form of the complex - presenting a model for the functional decrease of NO binding affinity of nitrophorins at physiological pH.     

Roles of N-terminal pyroglutamate in maintaining structural integrity and pKa values of catalytic histidine residues in bullfrog ribonuclease 3

Many proteins and bioactive peptides contain an N-terminal pyroglutamate residue (Pyr1). This residue reduces the susceptibility of the protein to aminopeptidases and often has important functional roles. The antitumor ribonuclease RC-RNase 3 (RNase 3) from oocytes of Rana catesbeiana (bullfrog) is one such protein. We have produced recombinant RNase 3 containing the N-terminal Pyr1 (pRNase 3) and found it to be indistinguishable from the native RNase 3 by mass spectrometry and a variety of other biochemical and immunological criteria. We demonstrated by NMR analysis that the Pyr1 of pRNase 3 forms hydrogen bonds with Lys9 and Ile96 and stabilizes the N-terminal alpha-helix in a rigid conformation. In contrast, the N-terminal alpha-helix becomes flexible and the pKa values of the catalytic residues His10 and His97 altered when Pyr1 formation is blocked by an extra methionine at the N terminus in the recombinant mqRNase 3. Thus, our results provide a mechanistic explanation on the essential role of Pyr1 in maintaining the structural integrity, especially at the N-terminal alpha-helix, and in providing the proper environment for the ionization of His10 and His97 residues for catalysis and cytotoxicity against HeLa cells.     

Spectroscopic characterization of a fruit-specific metallothionein: M. acuminata MT3

Metallothioneins (MTs) are ubiquitous low molecular mass, cysteine-rich proteins with the ability to bind d¹⁰ metal ions in the form of metal-thiolate clusters. In contrast to the vertebrate forms, knowledge about the properties of members of the plant metallothionein family is still scarce. The amino acid sequences of plant MTs are distinctively different to the sequences of other MT species. The protein under investigation, Musa acuminata (banana) MT3, belongs to the plant MT fruit-specific p3 subfamily. With a total of 10 cysteine residues, MT3 features a cysteine content and percentage that is more comparable to fungal and prokaryotic MTs than to the well characterized mammalian iso-forms. The gene sequence encoding MT3 was cloned into a suitable vector and the protein was recombinantly overexpressed in Escherichia coli. MT3 is able to coordinate a maximum of four divalent d¹⁰ metal ions under the formation of metal-thiolate clusters. The hitherto unknown spectroscopic behavior of MT3 in combination with the metal ions Zn²⁺, Cd²⁺, Pb²⁺, and Hg²⁺ will be presented and gives rise to the existence of a weaker metal ion coordination site. The pH stability of the investigated zinc and cadmium clusters is comparable to the values found for other plant metallothioneins though significantly lower than for the mammalian iso-forms. Possible metal-thiolate cluster structures will additionally be discussed.     

Spectroscopic characterization of Cicer arietinum metallothionein 1

The plant metallothioneins differ distinctively from other metallothionein families with respect to the cysteine distribution patterns, the presence of aromatic amino acids in most and histidine in some forms, as well as long cysteine-free amino acid stretches between cysteine-rich regions. Although known for more than 25 years, research activity on plant metallothioneins has been low increasing only in the past few years. In the following, we will present the first characterization of Cicer arietinum (chickpea) MT1. In this root-specific protein two cysteine-rich regions with six cysteine residues each are separated by a 42 amino acids long linker region. A synthetic gene encoding MT1 was designed, cloned into a suitable vector, and the protein was over-expressed in Escherichia coli. We could show, that MT1 has the ability to coordinate up to five Zn²⁺ or Cd²⁺ ions and even higher amounts of Hg²⁺. According to titration experiments pH-dependent zinc- and cadmium-thiolate cluster stability in MT1 is considerably lower than in vertebrate metallothioneins. The approximate contribution of secondary structural elements to the overall structure was assessed with circular dichroism and infrared spectroscopy. Hypothetical metal-thiolate cluster structures will be presented.     

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

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

Mutational analysis demonstrates that specific electrostatic interactions can play a key role in the denatured state ensemble of proteins

The nature of the denatured state ensemble has been controversial for decades owing, in large part, to the difficulty in characterizing the structure and energetics of denatured state interactions. There is increasing evidence for relatively non-specific hydrophobic clustering in the denatured states of some proteins but other types of interactions are much less well characterized. Here, we report the characterization of highly specific electrostatic interactions in the denatured state of a small alpha-beta protein, the N-terminal domain of the ribosomal protein L9 (NTL9). Mutation of Lys12 to Met has been shown to increase the stability of NTL9 significantly through the disruption of denatured state interactions. Here, we describe the analysis of the pH-dependent stability of 13 mutants designed to probe the nature of the Lys12 denatured state interaction. Lys12 is located in a lysine-rich region of the protein but analysis of a set of Lys to Met mutants shows that it plays a unique role in the denatured state. Analysis of mutants of all of the acidic residues in NTL9 shows that Lys12 forms a specific non-native electrostatic interaction with Asp8 in the denatured state ensemble. Thus the distribution of charge-charge interactions in the denatured state ensemble of NTL9 appears to be biased by few key interactions and is very different from that expected in a random coil. We propose that these interactions are not encoded by local sequence effects but rather reflect interactions among residues more distant in sequence. These results demonstrate that electrostatic as well as hydrophobic interactions can play an important role in the denatured state ensemble.     

Novel bisubstrate analog inhibitors of serotonin N-acetyltransferase: the importance of being neutral

Linker modified novel bisubstrate analog inhibitors 4-7 for serotonin N-acetyltransferase (arylalkylamine N-acetyltransferase, AANAT) have been designed and synthesized. Examination of these inhibitors with AANAT in vitro suggested that: (i) linker hydrogen bonding makes only modest contributions to the affinity of bisubstrate analog inhibitors studied; (ii) greater than or equal to four methylene groups between the indole and the coenzyme A (CoASH) moieties are required for a bisubstrate analog inhibitor to achieve strong AANAT inhibition; (iii) the AANAT active site appears not to accommodate positively charged linkers as well as neutral ones; and (iv) substrate amine pKa depression may constitute one strategy for AANAT substrate recognition and catalysis. The results reported here have enhanced our understanding of AANAT substrate recognition/catalysis, and are important for novel inhibitor design. Since AANAT belongs to the GCN5-related N-acetyltransferase (GNAT) superfamily, our experimental strategies should find applications for other acetyltransferases.     

Probing the sweet determinants of brazzein: wild-type brazzein and a tasteless variant, brazzein-ins(R18a-I18b), exhibit different pH-dependent NMR chemical shifts

Brazzein is a small, intensely sweet protein. As a probe of the functional properties of its solvent-exposed loop, two residues (Arg-Ile) were inserted between Leu18 and Ala19 of brazzein. Psychophysical testing demonstrated that this mutant is totally tasteless. NMR chemical shift mapping of differences between this mutant and brazzein indicated that residues affected by the insertion are localized to the mutated loop, the region of the single alpha-helix, and around the Cys16-Cys37 disulfide bond. Residues unaffected by this mutation included those near the C-terminus and in the loop connecting the alpha-helix and the second beta-strand. In particular, several residues of brazzein previously shown to be essential for its sweetness (His31, Arg33, Glu41, Arg43, Asp50, and Tyr54) exhibited negligible chemical shift changes. Moreover, the pH dependence of the chemical shifts of His31, Glu41, Asp50, and Tyr54 were unaltered by the insertion. The insertion led to large chemical shift and pKa perturbation of Glu36, a residue shown previously to be important for brazzein's sweetness. These results serve to refine the known sweetness determinants of brazzein and lend further support to the idea that the protein interacts with a sweet-taste receptor through a multi-site interaction mechanism, as has been postulated for brazzein and other sweet proteins (monellin and thaumatin).     

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.     

Electrostatic interactions in leucine zippers: thermodynamic analysis of the contributions of Glu and His residues and the effect of mutating salt bridges

Electrostatic interactions play a complex role in stabilizing proteins. Here, we present a rigorous thermodynamic analysis of the contribution of individual Glu and His residues to the relative pH-dependent stability of the designed disulfide-linked leucine zipper AB(SS). The contribution of an ionized side-chain to the pH-dependent stability is related to the shift of the pK(a) induced by folding of the coiled coil structure. pK(a)(F) values of ten Glu and two His side-chains in folded AB(SS) and the corresponding pK(a)(U) values in unfolded peptides with partial sequences of AB(SS) were determined by 1H NMR spectroscopy: of four Glu residues not involved in ion pairing, two are destabilizing (-5.6 kJ mol(-1)) and two are interacting with the positive alpha-helix dipoles and are thus stabilizing (+3.8 kJ mol(-1)) in charged form. The two His residues positioned in the C-terminal moiety of AB(SS) interact with the negative alpha-helix dipoles resulting in net stabilization of the coiled coil conformation carrying charged His (-2.6 kJ mol(-1)). Of the six Glu residues involved in inter-helical salt bridges, three are destabilizing and three are stabilizing in charged form, the net contribution of salt-bridged Glu side-chains being destabilizing (-1.1 kJ mol(-1)). The sum of the individual contributions of protonated Glu and His to the higher stability of AB(SS) at acidic pH (-5.4 kJ mol(-1)) agrees with the difference in stability determined by thermal unfolding at pH 8 and pH 2 (-5.3 kJ mol(-1)). To confirm salt bridge formation, the positive charge of the basic partner residue of one stabilizing and one destabilizing Glu was removed by isosteric mutations (Lys-->norleucine, Arg-->norvaline). Both mutations destabilize the coiled coil conformation at neutral pH and increase the pK(a) of the formerly ion-paired Glu side-chain, verifying the formation of a salt bridge even in the case where a charged side-chain is destabilizing. Because removing charges by a double mutation cycle mainly discloses the immediate charge-charge effect, mutational analysis tends to overestimate the overall energetic contribution of salt bridges to protein stability.     

Charge equilibria and pK(a) of malvidin-3-glucoside by electrophoresis

Paper electrophoresis has been used over the pH range 1.2 to 10.4 to measure apparent pK(a) values for malvidin-3-O-glucoside of pK(a(1)) 1.76+/-0.07, pK(a(2)) 5.36+/-0.04, and pK(a(3)) 8.39+/-0.07. Using solvent partitioning between buffered aqueous solutions and n-octanol, several micro-pK(a) constants for malvidin-3-O-glucoside were also identified, highlighting the complex nature of malvidin-3-glucoside equilibria. As a nonspectrophotometric procedure, the charge-dependent electrophoretic mobility method provided independent information on the net charge and color of anthocyanin species at wine pH (ca. 3.6). At this pH, the color of malvidin-3-glucoside in red wines is consistent only with the uncharged quinonoidal base as a major colored component of the equilibria.     

RFP tags for labeling secretory pathway proteins

Red fluorescent proteins (RFPs) are useful tools for live cell and multi-color imaging in biological studies. However, when labeling proteins in secretory pathway, many RFPs are prone to form artificial puncta, which may severely impede their further uses. Here we report a fast and easy method to evaluate RFPs fusion properties by attaching RFPs to an environment sensitive membrane protein Orai1. In addition, we revealed that intracellular artificial puncta are actually colocalized with lysosome, thus besides monomeric properties, pKa value of RFPs is also a key factor for forming intracellular artificial puncta. In summary, our current study provides a useful guide for choosing appropriate RFP for labeling secretory membrane proteins. Among RFPs tested, mOrange2 is highly recommended based on excellent monomeric property, appropriate pKa and high brightness.     

Sensing pH via p-cyanophenylalanine fluorescence: Application to determine peptide pKa and membrane penetration kinetics

We expand the spectroscopic utility of a well-known infrared and fluorescence probe, p-cyanophenylalanine, by showing that it can also serve as a pH sensor. This new application is based on the notion that the fluorescence quantum yield of this unnatural amino acid, when placed at or near the N-terminal end of a polypeptide, depends on the protonation status of the N-terminal amino group of the peptide. Using this pH sensor, we are able to determine the N-terminal pK_a values of nine tripeptides and also the membrane penetration kinetics of a cell-penetrating peptide. Taken together, these examples demonstrate the applicability of using this unnatural amino acid fluorophore to study pH-dependent biological processes or events that accompany a pH change.