A buried lysine that titrates with a normal pKa: role of conformational flexibility at the protein-water interface as a determinant of pKa values

Previously we reported that Lys, Asp, and Glu residues at positions 66 and 92 in staphylococcal nuclease (SNase) titrate with pK(a) values shifted by up to 5 pK(a) units in the direction that promotes the neutral state. In contrast, the internal Lys-38 in SNase titrates with a normal pK(a). The crystal structure of the L38K variant shows that the side chain of Lys-38 is buried. The ionizable moiety is approximately 7 A from solvent and ion paired with Glu-122. This suggests that the pK(a) value of Lys-38 is normal because the energetic penalty for dehydration is offset by a favorable Coulomb interaction. However, the pK(a) of Lys-38 was also normal when Glu-122 was replaced with Gln or with Ala. Continuum electrostatics calculations were unable to reproduce the pK(a) of Lys-38 unless the protein was treated with an artificially high dielectric constant, consistent with structural reorganization being responsible for the normal pK(a) value of Lys-38. This reorganization must be local because circular dichroism and NMR spectroscopy indicate that the L38K protein is native-like under all conditions studied. In molecular dynamics simulations, the ion pair between Lys-38 and Glu-122 is unstable. The simulations show that a minor rearrangement of a loop is sufficient to allow penetration of water to the amino moiety of Lys-38. This illustrates both the important roles of local flexibility and water penetration as determinants of pK(a) values of ionizable groups buried near the protein-water interface, and the challenges faced by structure-based pK(a) calculations in reproducing these effects.     

Role of a Helix B Lysine Residue in the Photoactive Site in Channelrhodopsins

In most studied microbial rhodopsins two conserved carboxylic acid residues (the homologs of Asp-85 and Asp-212 in bacteriorhodopsin) and an arginine residue (the homolog of Arg-82) form a complex counterion to the protonated retinylidene Schiff base, and neutralization of the negatively charged carboxylates causes red shifts of the absorption maximum. In contrast, the corresponding neutralizing mutations in some relatively low-efficiency channelrhodopsins (ChRs) result in blue shifts. These ChRs do not contain a lysine residue in the second helix, conserved in higher efficiency ChRs (Lys-132 in the crystallized ChR chimera). By action spectroscopy of photoinduced channel currents in HEK293 cells and absorption spectroscopy of detergent-purified pigments, we found that in tested ChRs the Lys-132 homolog controls the direction of spectral shifts in the mutants of the photoactive site carboxylic acid residues. Analysis of double mutants shows that red spectral shifts occur when this Lys is present, whether naturally or by mutagenesis, and blue shifts occur when it is replaced with a neutral residue. A neutralizing mutation of the Lys-132 homolog alone caused a red spectral shift in high-efficiency ChRs, whereas its introduction into low-efficiency ChR1 from Chlamydomonas augustae (CaChR1) caused a blue shift. Taking into account that the effective charge of the carboxylic acid residues is a key factor in microbial rhodopsin spectral tuning, these findings suggest that the Lys-132 homolog modulates their pK_a values. On the other hand, mutation of the Arg-82 homolog that fulfills this role in bacteriorhodopsin caused minimal spectral changes in the tested ChRs. Titration revealed that the pK_a of the Asp-85 homolog in CaChR1 lies in the alkaline region unlike in most studied microbial rhodopsins, but is substantially decreased by introduction of a Lys-132 homolog or neutralizing mutation of the Asp-212 homolog. In the three ChRs tested the Lys-132 homolog also alters channel current kinetics.     

The role of the active site residues in human galactokinase: implications for the mechanisms of GHMP kinases

Galactokinase catalyses the phosphorylation of galactose at the expense of ATP. Like other members of the GHMP family of kinases it is postulated to function through an active site base mechanism in which Asp-186 abstracts a proton from galactose. This asparate residue was altered to alanine and to asparagine by site-directed mutagenesis of the corresponding gene. This resulted in variant enzyme with no detectable galactokinase activity. Alteration of Arg-37, which lies adjacent to Asp-186 and is postulated to assist the catalytic base, to lysine resulted in an active enzyme. However, alteration of this residue to glutamate abolished activity. All the variant enzymes, except the arginine to lysine substitution, were structurally unstable (as judged by native gel electrophoresis in the presence of urea) compared to the wild type. This suggests that the lack of activity results from this structural instability, in addition to any direct effects on the catalytic mechanism. Computational estimations of the pK_a values of the arginine and aspartate residues, suggest that Arg-37 remains protonated throughout the catalytic cycle whereas Asp-186 has an abnormally high pK_a value (7.18). Quantum mechanics/molecular mechanics (QM/MM) calculations suggest that Asp-186 moves closer to the galactose molecule during catalysis. The experimental and theoretical studies presented here argue for a mechanism in which the C₁-OH bond in the sugar is weakened by the presence of Asp-186 thus facilitating nucleophilic attack by the oxygen atom on the γ-phosphorus of ATP.     

Role of a Helix B Lysine Residue in the Photoactive Site in Channelrhodopsins

In most studied microbial rhodopsins two conserved carboxylic acid residues (the homologs of Asp-85 and Asp-212 in bacteriorhodopsin) and an arginine residue (the homolog of Arg-82) form a complex counterion to the protonated retinylidene Schiff base, and neutralization of the negatively charged carboxylates causes red shifts of the absorption maximum. In contrast, the corresponding neutralizing mutations in some relatively low-efficiency channelrhodopsins (ChRs) result in blue shifts. These ChRs do not contain a lysine residue in the second helix, conserved in higher efficiency ChRs (Lys-132 in the crystallized ChR chimera). By action spectroscopy of photoinduced channel currents in HEK293 cells and absorption spectroscopy of detergent-purified pigments, we found that in tested ChRs the Lys-132 homolog controls the direction of spectral shifts in the mutants of the photoactive site carboxylic acid residues. Analysis of double mutants shows that red spectral shifts occur when this Lys is present, whether naturally or by mutagenesis, and blue shifts occur when it is replaced with a neutral residue. A neutralizing mutation of the Lys-132 homolog alone caused a red spectral shift in high-efficiency ChRs, whereas its introduction into low-efficiency ChR1 from Chlamydomonas augustae (CaChR1) caused a blue shift. Taking into account that the effective charge of the carboxylic acid residues is a key factor in microbial rhodopsin spectral tuning, these findings suggest that the Lys-132 homolog modulates their pK_a values. On the other hand, mutation of the Arg-82 homolog that fulfills this role in bacteriorhodopsin caused minimal spectral changes in the tested ChRs. Titration revealed that the pK_a of the Asp-85 homolog in CaChR1 lies in the alkaline region unlike in most studied microbial rhodopsins, but is substantially decreased by introduction of a Lys-132 homolog or neutralizing mutation of the Asp-212 homolog. In the three ChRs tested the Lys-132 homolog also alters channel current kinetics.     

Highly perturbed pKa values in the unfolded state of hen egg white lysozyme

The majority of pK_a values in protein unfolded states are close to the amino acid model pK_a values, thus reflecting the weak intramolecular interactions present in the unfolded ensemble of most proteins. We have carried out thermal denaturation measurements on the WT and eight mutants of HEWL from pH 1.5 to pH 11.0 to examine the unfolded state pK_a values and the pH dependence of protein stability for this enzyme. The availability of accurate pK_a values for the folded state of HEWL and separate measurements of mutant-induced effects on the folded state pK_a values, allows us to estimate the pK_a values of seven acidic residues in the unfolded state of HEWL. Asp-48 and Asp-66 display pK_a values of 2.9 and 3.1 in our analysis, thus representing the most depressed unfolded state pK_a values observed to date. We observe a strong correlation between the folded state pK_a values and the unfolded state pK_a values of HEWL, thus suggesting that the unfolded state of HEWL possesses a large degree of native state characteristics.     

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.     

Cloning and expression of <Emphasis Type="Italic">β</Emphasis>-glucosidase gene from <Emphasis Type="Italic">Bacillus licheniformis</Emphasis> into <Emphasis Type="Italic">E. coli</Emphasis> BL 21 (DE3)

A 1.4 Kb fragment of Bacillus licheniformis ATCC 14580 encoding β-glucosidase was cloned and expressed in Escherichia coli. β-Glucosidase expressed by E. coli harboring cloned gene was located entirely in the intracellular fraction. Recombinant β-glucosidase protein was purified to homogeneity level and the molecular weight was found to be 53 kDa using sodium dodecyl sulfate polyacrylamide gel electrophoresis analysis. It gave maximum activity at 50°C and pH 6. K _m and V _max were 0.206 mM and 1.26 U/mg, respectively, with p-nitrophenyl-β-D-glucopyranoside, while activation energy E_a, enthalpy of activation ?H and entropy of activation ΔS were found to be 66.31 kJ/mol, 64.04 kJ/mol and 48.28 J/mol/K, respectively. The pK_a1 and pK_a2 of the ionizable groups of active site residues involved in V_max were found to be 5.5 and 7.0, respectively. When the recombinant β-glucosidase protein was used as a member of consortium with endoglucanase and exoglucanase for the saccharification of wheat straw, 123% increase in saccharification was observed.     

Factors That Differentiate the H-bond Strengths of Water Near the Schiff Bases in Bacteriorhodopsin and Anabaena Sensory Rhodopsin

Bacteriorhodopsin (BR) functions as a light-driven proton pump, whereas Anabaena sensory rhodopsin (ASR) is believed to function as a photosensor despite the high similarity in their protein sequences. In Fourier transform infrared (FTIR) spectroscopic studies, the lowest O-D stretch for D₂O was observed at ∼2200 cm⁻¹ in BR but was significantly higher in ASR (>2500 cm⁻¹), which was previously attributed to a water molecule near the Schiff base (W402) that is H-bonded to Asp-85 in BR and Asp-75 in ASR. We investigated the factors that differentiate the lowest O-D stretches of W402 in BR and ASR. Quantum mechanical/molecular mechanical calculations reproduced the H-bond geometries of the crystal structures, and the calculated O-D stretching frequencies were corroborated by the FTIR band assignments. The potential energy profiles indicate that the smaller O-D stretching frequency in BR originates from the significantly higher pK_a(Asp-85) in BR relative to the pK_a(Asp-75) in ASR, which were calculated to be 1.5 and −5.1, respectively. The difference is mostly due to the influences of Ala-53, Arg-82, Glu-194–Glu-204, and Asp-212 on pK_a(Asp-85) in BR and the corresponding residues Ser-47, Arg-72, Ser-188-Asp-198, and Pro-206 on pK_a(Asp-75) in ASR. Because these residues participate in proton transfer pathways in BR but not in ASR, the presence of a strongly H-bonded water molecule near the Schiff base ultimately results from the proton-pumping activity in BR.     

Prediction of Secondary Ionization of the Phosphate Group in Phosphotyrosine Peptides

A computational approach, based on a continuum molecular electrostatics model, for the calculation of the pK_a values of secondary ionization of the phosphate group in phenyl phosphate derivatives is described. The method uses the ESP atomic charges of the mono-anionic and di-anionic forms of the ionizable phosphate group, computed with the use of the density functional method, and applies a new concept of the model group, being the reference state for the pK_a calculations. Both conformational flexibility and tautomeric degrees of freedom are taken into account in the calculations. The method was parameterized using experimentally available pK_a values of four derivatives of phenyl phosphates, and phosphotyrosine. Subsequently this parameterization was used to predict pK_a of the phosphate group in a short peptide Gly-Gly-Tyr(P)-Ala, and in a longer peptide consisting of 12 residues, the latter in water, and in a complex with a protein—phospholipase. The agreement between the computed and the experimental pK_a values is better than ±0.3 pH units for the optimized solute dielectric constant of 11–13. This approach is promising and its extension to other phospho-amino acids is in progress.     

P100, a transcriptional coactivator,is a human homologue of staphylococcal nuclease.

Staphylococcus aureus nuclease (SNase) homologues, previously thought to be restricted to bacteria and archaea, are demonstrated by sequence analysis to be present also in eukaryotes. The human cellular coactivator p100 is shown to contain four repeats, each of which is a SNase homologue. Surprisingly, these repeats are unlikely to possess SNase-like activities as each lacks equivalent SNase catalytic residues, yet they may mediate p100''s single-stranded DNA-binding function. Products of Corydalis sempervirens and Saccharomyces cerevisiae open reading frames are predicted to adopt the same fold and possess similar functions as SNase. Five additional hypothetical proteins of bacterial origin are also predicted to be active SNase-like nucleases, including one that appears to be C-terminally truncated in a manner analogous to an engineered active SNase variant. Conservation of Asp-19 and Asp-83 among these homologues suggests a re-evaluation of the roles of these residues in Ca(2+)-binding and/or catalysis.     

Experimental Evidence for the Existence of a Stable Half-Barrel Subdomain in the (β/α)8-Barrel Fold

The (β/α)₈-barrel is one of the most common folds functioning as enzymes. The emergence of two (β/α)₈-barrel enzymes involved in histidine biosynthesis, each of which has a twofold symmetric structure, has been proposed to be a consequence of tandem duplication and fusion of a (β/α)₄-half-barrel. However, little evidence has been found for the existence of an ancestral half-barrel in the evolution of other (β/α)₈-barrel proteins. In order to detect remnants of an ancestral half-barrel in the (β/α)₈-barrel structure of Escherichia coli N-(5′-phosphoribosyl)anthranilate isomerase, we engineered three potential half-barrel units, (β/α)_1-4, (β/α)_3-6, and (β/α)_5-8. Among these three arrangements, only (β/α)_3-6 is stable; it exists in equilibrium between monomeric and dimeric forms. Thus, the central segment of N-(5′-phosphoribosyl)anthranilate isomerase from E. coli can serve as a half-barrel precursor. A tandem duplication of (β/α)_3-6 yielded predominantly monomeric structures that were quite stable. This result exemplified that the structural characteristics of noncovalently assembled half-barrels could be improved by duplication and fusion. Moreover, our results may provide information regarding the local structural units that encompass interactions important for the early folding events of this ubiquitous protein conformation.     

Crystallographic study of hydration of an internal cavity in engineered proteins with buried polar or ionizable groups

Although internal water molecules are essential for the structure and function of many proteins, the structural and physical factors that govern internal hydration are poorly understood. We have examined the molecular determinants of internal hydration systematically, by solving the crystal structures of variants of staphylococcal nuclease with Gln-66, Asn-66, and Tyr-66 at cryo (100 K) and room (298 K) temperatures, and comparing them with existing cryo and room temperature structures of variants with Glu-66, Asp-66, Lys-66, Glu-92 or Lys-92 obtained under conditions of pH where the internal ionizable groups are in the neutral state. At cryogenic temperatures the polar moieties of all these internal side chains are hydrated except in the cases of Lys-66 and Lys-92. At room temperature the internal water molecules were observed only in variants with Glu-66 and Tyr-66; water molecules in the other variants are probably present but they are disordered and therefore undetectable crystallographically. Each internal water molecule establishes between 3 and 5 hydrogen bonds with the protein or with other internal water molecules. The strength of interactions between internal polar side chains and water molecules seems to decrease from carboxylic acids to amides to amines. Low temperature, low cavity volume, and the presence of oxygen atoms in the cavity increase the positional stability of internal water molecules. This set of structures and the physical insight they contribute into internal hydration will be useful for the development and benchmarking of computational methods for artificial hydration of pockets, cavities, and active sites in proteins.     

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

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

Conserved residues in the aromatic acid/H symporter family are important for benzoate uptake by NCgl2325 in Corynebacterium glutamicum

Corynebacterium glutamicum is a model organism for genetic and physiological studies in Gram-positive bacteria. NCgl2325 in C. glutamicum, a transporter belonging to the aromatic acid/H⁺ symporter family, has previously been reported to be involved in benzoate assimilation. Here, we showed that this transporter, fused with GFP, was associated with the cell membrane in Escherichia coli and C. glutamicum. Uptake assays with [¹⁴C]-labeled benzoate demonstrated that NCgl2325 transported benzoate into the cells at a V_max of 0.19 ± 0.01 nmol/min/mg of dry weight, and the K_m value was determined to be 1.11 ± 0.24 ??M. Among the competing substrates tested, hydroxyl-substituted benzoates resulted in significant inhibition (>50%) of benzoate uptake. Site-directed mutagenesis of conserved residues in the hydrophilic cytoplasmic loops (Gly-80, Asp-84 and Asp-312) and the hydrophobic transmembrane regions (Asp-35, Arg-119, Glu-139 and Arg-386) resulted in loss of benzoate transport activity. This is the first study to investigate the molecular basis of benzoate transport.     

Self‐guided Langevin dynamics study of regulatory interactions in NtrC

Multiple self‐guided Langevin dynamics (SGLD) simulations were performed to examine structural and dynamical properties of the receiver domain of nitrogen regulatory protein C (NtrC^r). SGLD and MD simulations of the phosphorylated active form structure suggest a mostly stable but broad structural ensemble of this protein. The finite difference Poisson–Boltzmann calculations of the pK_a values of the active site residues suggest an increase in the pK_a of His‐84 on phosphorylation of Asp‐54. In SGLD simulations of the phosphorylated active form with charged His‐84, the average position of the regulatory helix α4 is found closer to the starting structure than in simulations with the neutral His‐84. To model the transition pathway, the phosphate group was removed from the simulations. After 7 ns of simulations, the regulatory helix α4 was found approximately halfway between positions in the NMR structures of the active and inactive forms. Removal of the phosphate group stimulated loss of helix α4, suggesting that the pathway of conformational transition may involve partial unfolding mechanism. The study illustrates the potential utility of the SGLD method in studies of the coupling between ligand binding and conformational transitions. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Induced β-Barrel Formation of the Alzheimer's Aβ25-35 Oligomers on Carbon Nanotube Surfaces: Implication for Amyloid Fibril Inhibition

Recent experimental studies show that carbon nanotubes impact the aggregation process of proteins associated with neurodegenerative diseases. However, the details of molecular interactions between proteins and carbon nanotubes are still not well understood. In this study, we investigate the initial adsorption features and dynamics of the Alzheimer's amyloid-β peptide spanning residues 25-35 (Aβ25-35) on a single-walled carbon nanotube (SWNT) surface using fully atomic molecular dynamics simulations (MD) in explicit solvent. The initial configurations of the Aβ25-35 peptides consist of two preformed bilayer β-sheets, each with four or five β-strands in parallel or mixed antiparallel-parallel orientations. Our simulations show, for what we believe is the first time, that two disjointed Aβ25-35 β-sheets with mixed antiparallel-parallel strands can assemble into β-barrels wrapping the SWNT. In contrast, both simulations of Aβ25-35 without SWNT, and simulations of SWNT−Aβ25-35 with purely parallel β-strands, lead to disordered aggregates. We find that Aβ25-35 β-barrel formation involves at least two steps: i), curving of the Aβ25-35 β-sheets as a result of strong hydrophobic interactions with carbon nanotube concomitantly with dehydration of the SWNT-peptide interface; and ii), intersheet backbone hydrogen bond formation with fluctuating intrasheet hydrogen bonds. Detailed analysis of the conversion shows that β-barrel formation on SWNT surface results from the interplay of dehydration and peptide-SWNT/peptide-peptide interactions. Implications of our results on amyloid fibril inhibition are discussed.     

Hydrogen Ion Equilibria of Cajanus cajan Lectin

Hydrogen ion titration of an affinity-purified mannose/glucose-specific lectin from Cajanus cajan pulse was carried out at 30°C and ionic strength of 0.15 by a discontinuous method. The titration was reversible in the pH range 2–12.0. The numbers of different ionizable groups per 39,000 g of the lectin were 43 carboxyl groups (pK_int = 3.93), 10 imidazole groups, 21 ε-amino groups, 12.8 phenoxyl groups (pK_int = 10.0), and 5 guanidyl groups. Only seven tyrosine residues of the lectin were dissociated under native conditions. The remaining six tyrosines became available for titration upon denaturation of the lectin in 9 M urea.     

Graphical analysis of pH-dependent properties of proteins predicted using PROPKA

Background  

Charge states of ionizable residues in proteins determine their pH-dependent properties through their pK_a values. Thus, various theoretical methods to determine ionization constants of residues in biological systems have been developed. One of the more widely used approaches for predicting pK_a values in proteins is the PROPKA program, which provides convenient structural rationalization of the predicted pK_a values without any additional calculations.