Nanometric design of extraordinary hydrophobic-induced pKa shifts for aspartic acid: relevance to protein mechanisms

Strategies for tailoring pH performances of glycoside hydrolases

Glycoside hydrolases (GHs) exhibit high activity and stability under harsh conditions, such as high temperatures and extreme pHs, given their wide use in industrial biotechnology. However, strategies for improving the acidophilic and alkalophilic adaptations of GHs are poorly summarized due to the complexity of the mechanisms of these adaptations. This review not only highlights the adaptation mechanisms of acidophilic and alkalophilic GHs under extreme pH conditions, but also summarizes the recent advances in engineering the pH performances of GHs with a focus on four strategies of protein engineering, enzyme immobilization, chemical modification, and medium engineering (additives). The examples described here summarize the methods used in modulating the pH performances of GHs and indicate that methods integrated in different protein engineering techniques or methods are efficient to generate industrial biocatalysts with the desired pH performance and other adapted enzyme properties.

Self-Assembly of an Amino Acid Derivative into an Antimicrobial Hydrogel Biomaterial

N-(4-Nitrobenzoyl)-Phe self-assembled into a transparent supramolecular hydrogel, which displayed high fibroblast and keratinocyte cell viability. The compound showed a mild antimicrobial activity against E. coli both as a hydrogel and in solution. Single-crystal XRD data revealed packing details, including protonation of the C-terminus due to an apparent pK_a shift, as confirmed by pH titrations. MicroRaman analysis revealed almost identical features between the gel and crystal states, although more disorder in the former. The hydrogel is thermoreversible and disassembles within a range of temperatures that can be fine-tuned by experimental conditions, such as gelator concentration. At the minimum gelling concentration of 0.63 wt %, the hydrogel disassembles in a physiological temperature range of 39-42 °C, thus opening the way to its potential use as a biomaterial.

Effect of sequence on the ionization behavior of a series of amphiphilic polypeptides

The behavior of five polypeptides made of hydrophilic and pH-responsive aspartic acid (Asp) and hydrophobic phenylalanine (Phe), which had been prepared by stitching together short well-defined sequences of Asp and Phe, was studied as a function of pH. The effect of pH on these polypeptides referred to as (Asp3Phe1)n, (Asp2Phe1)n, (Asp1Phe1)n, (Asp1Phe2)n, and (Asp1Phe3)n varied dramatically depending on their constituting sequence. The more hydrophobic polypeptides (Asp1Phe2)n and (Asp1Phe3)n behaved as if the Asp's were isolated from each other and showed an apparent pKa (pKa(app)) that remained constant with level of ionization (α = [Asp(-)]/[Asp]total) and equaled 5.4 and 6.4, respectively. The more hydrophilic polypeptides (Asp3Phe1)n and (Asp2Phe1)n behaved like weak polyacids showing a linear increase in pKa(app) with increasing α. The pKa(app) of (Asp1Phe1)n showed a trend as a function of α intermediate between the Asp-rich and Phe-rich polypeptides, behaving as if the Asp's were isolated at low α values (<0.35) but acting as a weak polyacid for large α values (>0.35). The effect that α, and thus the charge density of the polypeptides, had on the collapse and aggregation of the polypeptides was characterized by conducting static light scattering and fluorescence measurements. Static light scattering measurements demonstrated that all polypeptides precipitated and aggregated in solution at a critical charge density of 0.2. Fluorescence measurements with pyrene indicated that this behavior was due to the formation of Phe aggregates in water. Together, these experiments provide a complete description of how pH affects the behavior of a series of unique amphiphilic polypeptides designed with a well-defined sequence.

Effect of glycine substitution on Fmoc-diphenylalanine self-assembly and gelation properties

We have investigated the self-assembly behavior of fluorenyl-9-methoxycarbonyl (Fmoc)-FG, Fmoc-GG, and Fmoc-GF and compared it to that of Fmoc-FF using potentiometry, fluorescence and infrared spectroscopy, transmission electron microscopy, wide-angle X-ray scattering, and oscillatory rheometry. Titration experiments revealed a substantially shifted apparent pK(a) transition for Fmoc-FG, Fmoc-GG, and Fmoc-GF. The apparent pK(a) values observed correlated with the hydrophobicity (log P) of the Fmoc-dipeptide molecules. Fmoc-GG and Fmoc-GF were found to self-assemble only in their protonated form (below their apparent pK(a)), while Fmoc-FG formed self-assembled structures above and below its apparent pK(a). Fmoc-GG and Fmoc-FG were found to form hydrogels below their apparent pK(a) transitions in agreement with the entangled fibers morphologies revealed by TEM. Unlike Fmoc-FF and Fmoc-GG, Fmoc-FG showed unusual gelation behavior as gels were found to form upon heating. Fmoc-GF formed precipitates instead of a hydrogel below its apparent pK(a) in agreement with the formation of micrometer scale sheetlike structures observed by TEM. The fact that all four Fmoc-dipeptides were found to self-assemble suggests that the main driving force behind the self-assembly process is a combination of the hydrophobic and π-π interactions of the fluorenyl moieties with a secondary role for hydrogen bonding of the peptidic components. The nature of the peptidic tail was found to have a pronounced effect on the type of self-assembled structure formed. This work indicates that the substitution of phenylalanine by glycine significantly impacts on the mode of assembly and illustrates the versatility of aromatic peptide amphiphiles in the formation of structurally diverse nanostructures.

The optimization of polymalic acid peptide copolymers for endosomolytic drug delivery

Membranolytic macromolecules are promising vehicles for cytoplasmic drug delivery, but their efficiency and safety remains primary concerns. To address those concerns, membranolytic properties of various poly(β-L-malic acid) (PMLA) copolymers were extensively investigated as a function of concentration and pH. PMLA, a naturally occurring biodegradable polymer, acquires membranolytic activities after substitution of pendent carboxylates with hydrophobic amino acid derivatives. Ruled by hydrophobization and charge neutralization, membranolysis of PMLA copolymers increased as a function of polymer molecular weight and demonstrated a maximum with 50% substitution of carboxylates. Charge neutralization was achieved either conditionally by pH-dependent protonation or permanently by masking carboxylates. Membranolysis of PMLA copolymers containing tripeptides of leucine, tryptophan and phenylalanine were pH-dependent in contrast to pH-independent copolymers of Leucine ethyl ester and Leu-Leu-Leu-NH(2) with permanent charge neutralization. PMLA and tripeptides seemed a unique combination for pH-dependent membranolysis. In contrast to nontoxic pH-dependent PMLA copolymers, pH-independent copolymers were found toxic at high concentration, which is ascribed to their nonspecific disruption of plasma membrane at physiological pH. pH-Dependent copolymers were membranolytically active only at acidic pH typical of maturating endosomes, and are thus devoid of cytotoxicity. The PMLA tripeptide copolymers are useful for safe and efficient cytoplasmic delivery routed through endosome.

Elastic-contractile model proteins: Physical chemistry, protein function and drug design and delivery

This review presents the structure and physico-chemical properties of ECMPs, elastic-contractile model proteins using sparse design modifications of elastic (GVGVP)(n); it describes the capacity of ECMP to perform the energy conversions that sustain living organisms; it arrives at the hydration thermodynamics of ECMP in terms of the change in Gibbs free energy of hydrophobic association, ΔG(HA), and the apolar-polar repulsive free energy of hydration, ΔG(ap); it applies ΔG(HA), ΔG(ap), and the nature of elasticity to describe the function of basic diverse proteins, namely - the F₁-motor of ATP synthase, Complex III of mitochondria, the KscA potassium-channel, and the molecular chaperonin, GroEL/ES; it applies ΔG(HA) and ΔG(ap) to describe the function of ABC exporter proteins that confer multi-drug resistance (MDR) on micro-organisms and human carcinomas and suggests drug modifications with which to overcome MDR. Using ECMP, means are demonstrated, for quantifying drug hydrophobicity with which to combat MDR and for preparing ECMP drug delivery nanoparticles, ECMPddnp, decorated with synthetic antigen-binding fragments, Fab1 and Fab2, with which to target specific up-regulated receptors, characteristic of human carcinoma cells, for binding and localized drug release.

Determination of ligand-binding sites on proteins using long-range hydrophobic potential

Here we developed a new program, HydrophoBicity On a Protein (HBOP), to find the ligand-binding site of a protein using the long-range hydrophobic-potential function estimated from the experimental data of Israelachvili and Pashley. We calculated the hydrophobic-potential energies at each grid point of a lattice around a protein using the potential function. The hydrophobic potential was evaluated using the carbon atoms of the hydrophobic residues, with the exception of those of the amide groups. We tested HBOP on 26 types of protein (72 protein-ligand complexes), the three-dimensional structures of which were determined experimentally. Although only one hydrophobic function was used, HBOP could successfully identify the binding sites in all of the proteins tested. Moreover, in 24 of the proteins, the binding sites were located in the most hydrophobic region. Surprisingly, the binding sites on sugar binding proteins were the most hydrophobic sites. It implies that the hydrophobic interaction plays an important role in the formation of protein-ligand complexes.

Molecular biophysics of elastin structure, function and pathology

Owing to the presence of the recurring sequence XPGX' (where X and X' are hydrophobic residues), the molecular structure of the sequences between cross-links in elastin is viewed primarily as a series of beta-turns which become helically ordered by hydrophobic folding into beta-spirals, which in turn assemble hydrophobically into twisted filaments. Both hydrophobic folding and assembly occur when the temperature is raised above Tt, the onset of an inverse temperature transition. Using poly[fv(VPGVG),fx(VPGXG)] (where fv and fx are mole fractions with fv + fx = 1 and X is now any of the naturally occurring amino acid residues), plots of fx versus Tt result in a new hydrophobicity scale based directly on the hydrophobic folding and assembly processes of interest. With the reference values chosen at fx = 1, the most hydrophobic residues of elastin, Tyr (Y) and Phe (F), have low values of Tt, -55 and -30 degrees C, respectively, and the most hydrophilic residues, Glu (E-), Asp (D-) and Lys (K+), have high values of 250, 170 and 120 degrees C, respectively. Raising the average value of Tt for a chain or chain segment from below to above physiological temperature drives hydrophobic unfolding and disassembly; lowering Tt does the reverse. This delta Tt mechanism has been used reversibly to interconvert many energy forms and is used here to explain initiating events of elastogenesis, pulmonary emphysema, solar elastosis and the paucity of elastic fibres in scar tissue. In general, oxidation and/or photolysis convert(s) hydrophobic residues into polar residues with the consequences of irreversibly raising Tt to above 37 degrees C, hydrophobic unfolding and disassembly (fibre swelling), and greater susceptibility to proteolysis.

Topology of the porin MspA in the outer membrane of Mycobacterium smegmatis

MspA is the major porin of Mycobacterium smegmatis mediating the exchange of hydrophilic solutes across the outer membrane (OM). It is the prototype of a new family of octameric porins with a single central channel of 9.6 nm in length and consists of two hydrophobic beta-barrels of 3.7 nm in length and a more hydrophilic, globular rim domain. The length of the hydrophobic domain of MspA does not match the thicknesses of mycobacterial OMs of 5-12 nm as derived from electron micrographs. Further, the membrane topology of MspA is unknown as it is for any other mycobacterial OM protein. We used MspA as a molecular ruler to define the boundaries of the OM of M. smegmatis by surface labeling of single cysteine mutants. Seventeen mutants covered the surface of the rim domain and were biotinylated with a membrane-impermeable reagent. The label efficiencies in vitro were remarkably similar to the predicted accessibilities of the cysteines. By contrast, six of these mutants were protected from biotinylation in M. smegmatis cells. Tryptophan 21 defines a horizontal plane that dissects the surface-exposed versus the membrane-protected residues of MspA. The 8 phenylalanines at position 99 form a ring at the periplasmic end of the hydrophobic beta-barrel domain. These results indicated that (i) the membrane boundaries of MspA are defined by aromatic girdles as in porins of Gram-negative bacteria and (ii) loops and a 3.4-nm long part of the hydrophilic rim domain are embedded into the OM of M. smegmatis. This is the first report suggesting that elements other than hydrophobic alpha-helices or beta-sheets are integrated into a lipid membrane.

Stabilization of internal charges in a protein: water penetration or conformational change?

The ionizable amino acid side chains of proteins are usually located at the surface. However, in some proteins an ionizable group is embedded in an apolar internal region. Such buried ionizable groups destabilize the protein and may trigger conformational changes in response to pH variations. Because of the prohibitive energetic cost of transferring a charged group from water to an apolar medium, other stabilizing factors must be invoked, such as ionization-induced water penetration or structural changes. To examine the role of water penetration, we have measured the 17O and 2H magnetic relaxation dispersions (MRD) for the V66E and V66K mutants of staphylococcal nuclease, where glutamic acid and lysine residues are buried in predominantly apolar environments. At neutral pH, where these residues are uncharged, we find no evidence of buried water molecules near the mutation site. This contrasts with a previous cryogenic crystal structure of the V66E mutant, but is consistent with the room-temperature crystal structure reported here. MRD measurements at different pH values show that ionization of Glu-66 or Lys-66 is not accompanied by penetration of long-lived water molecules. On the other hand, the MRD data are consistent with a local conformational change in response to ionization of the internal residues.

Protein engineering of bacterial alpha-amylases

alpha-Amylases constitute a very diverse family of glycosyl hydrolases that cleave alpha1-->4 linkages in amylose and related polymers. Recent structural and mutagenic studies of archeael, mammalian and bacterial alpha-amylases have resulted in a wealth of information on the catalytic mechanism and on the structural features of this enzyme class. Because of their high thermo-stability, the Bacillus alpha-amylases have found widespread use in industrial processes, and much attention has been devoted to optimising these enzymes for the very harsh conditions encountered there. Stability has been a major area of focus in this respect, and several remarkably stable bacterial alpha-amylases have been produced by bioengineering techniques. Protein engineering studies of pH-activity profiles and of substrate specificities have also been initiated, although without much success. In the coming years it is likely, however, that the focus of alpha-amylase engineering will shift from engineering stability to these new areas.

The ionization of a buried glutamic acid is thermodynamically linked to the stability of Leishmania mexicana triose phosphate isomerase

The amino acid sequence of Leishmania mexicana triose phosphate isomerase is unique in having at position 65 a glutamic acid instead of a glutamine. The stability properties of LmTIM and the E65Q mutant were investigated by pH and guanidinium chloride-induced unfolding. The crystal structure of E65Q was determined. Three important observations were made: (a) there are no structural rearrangements as the result of the substitution; (b) the mutant is more stable than the wild-type; and (c) the stability of the wild-type enzyme shows strong pH dependence, which can be attributed to the ionization of Glu65. Burying of the Glu65 side chain in the uncharged environment of the dimer interface results in a shift in pKa of more than 3 units. The pH-dependent decrease in overall stability is due to weakening of the monomer-monomer interactions (in the dimer). The E65Q substitution causes an increase in stability as the result of the formation of an additional hydrogen bond in each subunit (DeltaDeltaG degrees of 2 kcal.mol-1 per monomer) and the elimination of a charged group in the dimer interface (DeltaDeltaG degrees of at least 9 kcal.mol-1 per dimer). The computated shift in pKa and the stability of the dimer calculated from the charge distribution in the protein structure agree closely with the experimental results. The guanidinium chloride dependence of the unfolding constant was smaller than expected from studies involving monomeric model proteins. No intermediates could be identified in the unfolding equilibrium by combining fluorescence and CD measurements. Study of a stable monomeric triose phosphate isomerase variant confirmed that the phenomenon persists in the monomer.

Streptavidin-binding and -dimerizing ligands discovered by phage display, topochemistry, and structure-based design

Structural and mechanistic determinants of affinity of streptavidin-binding peptide ligands discovered by phage display are reviewed along with the use of streptavidin as a paradigm for structure-based design. A novel way of producing protein-dimerizing ligands in the streptavidin model system is discussed, in which crystal packing topochemically mediates or even catalyzes dimerization of adjacent bound ligands whose reactive ligating groups are presented toward one another in productive orientations in the crystal lattice. Finally, through crystallography on a set of streptavidin complexes with small molecule and peptide ligands at multiple pHs in two space groups, the mechanism by which ligands enhance intersubunit stabilization of the streptavidin tetramer is probed.

Electrostatics in the active site of an alpha-amylase

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

Electrophilic assistance by Asp-99 of 3-oxo-Delta 5-steroid isomerase

3-Oxo-Delta 5-steroid isomerase (Delta 5-3-ketosteroid isomerase, KSI; EC 5.3.3.1) catalyzes the conversion of a variety of beta, gamma-unsaturated 3-oxosteroids to their corresponding alpha, beta-unsaturated isomers at rates that approach the diffusion limit for specific substrates. The reaction proceeds through a dienolate intermediate, with two amino acid residues (Asp-38 and Tyr-14) known to be involved in catalysis. When the complete three-dimensional structure of KSI was determined recently by NMR methods, an additional polar residue (Asp-99) was found in the active site and this group was shown to be important for catalytic activity. In this work, we examine the properties of several mutant KSIs to determine the nature of catalysis by Asp-99 of KSI. The electrophoretic mobilities of wild-type (WT) KSI and several mutants (D99A, D99N, D38N, and D38N/D99A) on native gels were determined at pH values ranging from 6.0 to 8.5. The results demonstrate that the pKa of Asp-99 is >8.5 in wild-type KSI. The pH-rate profiles for the D99A, D99N, and D38H/D99A mutants of KSI were also determined. For all three mutants, kcat and kcat/KM do not decrease at high pH, in contrast to those for WT and D38H, which lose activity above pH 9 and 8, respectively. Mutation of Asp-99 to Asn decreases kcat for the substrate 5-androstene-3,17-dione by 27-fold and kcat/Km by 23-fold, substantially less than the loss of activity (3000-fold in kcat and 2200-fold in kcat/Km) observed when Asp-99 is mutated to Ala, consistent with a hydrogen bonding role for Asp-99. Taken together, these results provide evidence that Asp-99 participates in catalysis in its protonated form, with a pKa of >9 in WT and approximately 8.5 in the D38H mutant. Asp-99 likely donates a hydrogen bond to O-3 of the steroid, helping to stabilize the transition state(s) of the KSI-catalyzed reaction.

Engineering protein-based machines to emulate key steps of metabolism (biological energy conversion)

Metabolism is the conversion of available energy sources to those energy forms required for sustaining and propagating living organisms; this is simply biological energy conversion. Proteins are the machines of metabolism; they are the engines of motility and the other machines that interconvert energy forms not involving motion. Accordingly, metabolic engineering becomes the use of natural protein-based machines for the good of society. In addition, metabolic engineering can utilize the principles, whereby proteins function, to design new protein-based machines to fulfill roles for society that proteins have never been called upon throughout evolution to fulfill. This article presents arguments for a universal mechanism whereby proteins perform their diverse energy conversions; it begins with background information, and then asserts a set of five axioms for protein folding, assembly, and function and for protein engineering. The key process is the hydrophobic folding and assembly transition exhibited by properly balanced amphiphilic protein sequences. The fundamental molecular process is the competition for hydration between hydrophobic and polar, e.g., charged, residues. This competition determines Tt, the onset temperature for the hydrophobic folding and assembly transition, Nhh, the numbers of waters of hydrophobic hydration, and the pKa of ionizable functions. Reported acid-base titrations and pH dependence of microwave dielectric relaxation data simultaneously demonstrate the interdependence of Tt, Nhh and the pKa using a series of microbially prepared protein-based poly(30mers) with one glutamic acid residue per 30mer and with an increasing number of more hydrophobic phenylalanine residues replacing valine residues. Also, reduction of nicotinamides and flavins is shown to lower Tt, i.e., to increase hydrophobicity. Furthermore, the argument is presented, and related to an extended Henderson-Hasselbalch equation, wherein reduction of nicotinamides represents an increase in hydrophobicity and resulting hydrophobic-induced pKa shifts become the basis for understanding a primary energy conversion (proton transport) process of mitochondria. Copyright 1998 John Wiley & Sons, Inc.

In crystals of complexes of streptavidin with peptide ligands containing the HPQ sequence the pKa of the peptide histidine is less than 3.0

The pH dependences of the affinities for streptavidin of linear and cyclic peptide ligands containing the HPQ sequence discovered by phage display were determined by plasmon resonance measurements. At pH values ranging from 3.0 to 9.0, the Kd values for Ac-AEFSHPQNTIEGRK-NH2, cyclo-Ac-AE[CHPQGPPC]IEGRK-NH2, and cyclo-Ac-AE[CHPQFC]IEGRK-NH2, were determined by competition, and those for cyclo-[5-S-valeramide-HPQGPPC]K-NH2 were determined directly by equilibrium affinity measurements. The Kd values of the ligands increase by an average factor of 3.0 +/- 0.8 per decrease in pH unit between pH approximately 4.5 and pH approximately 6.3. Below pH approximately 4.5 there is a smaller increase in Kd values, and above pH approximately 6.3 the Kd values become relatively pH-independent. We determined the crystal structures of complexes of streptavidin with cyclo-[5-S-valeramide-HPQGPPC]K-NH2 at pH 1.5, 2.5, 3.0, and 3.5, with cyclo-Ac-[CHPQFC]-NH2 at pH 2.0, 3.0, 3.6, 4.2, 4.8, and 11.8, with cyclo-Ac-[CHPQGPPC]-NH2 at pH 2.5, 2.9, and 3.7, and with FSHPQNT at pH 4.0 and compared the structures with one another and with those previously determined at other pH values. At pH values from 3.0 to 11.8, the electron density for the peptide His side chain is strong, flat, and well defined. A hydrogen bond between the Ndelta1 atom of the His and the peptide Gln amide group indicates the His of the bound peptide in the crystals is uncharged at pH >/= 3.0. By determining selected structures in two different space groups, I222 with two crystallographically inequivalent ligand sites and I4122 with one site, we show that below pH approximately 3.0, the pKa of the bound peptide His in the crystals is influenced by crystal packing interactions. The presence of the Ndelta1His-NGln hydrogen bond along with pH dependences of the peptide affinities suggest that deprotonation of the peptide His is required for high affinity binding of HPQ-containing peptides to streptavidin both in the crystals and in solution.

Scanning mutagenesis of the putative transmembrane segments of Kir2.1, an inward rectifier potassium channel

Structural models of inward rectifier K+ channels incorporate four identical or homologous subunits, each of which has two hydrophobic segments (M1 and M2) which are predicted to span the membrane as alpha helices. Since hydrophobic interactions between proteins and membrane lipids are thought to be generally of a nonspecific nature, we attempted to identify lipid-contacting residues in Kir2.1 as those which tolerate mutation to tryptophan, which has a large hydrophobic side chain. Tolerated mutations were defined as those which produced measurable inwardly rectifying currents in Xenopus oocytes. To distinguish between water-accessible positions and positions adjacent to membrane lipids or within the protein interior we also mutated residues in M1 and M2 individually to aspartate, since an amino acid with a charged side chain should not be tolerated at lipid-facing or interior positions, due to the energy cost of burying a charge in a hydrophobic environment. Surprisingly, 17 out of 20 and 17 out of 22 non-tryptophan residues in M1 and M2, respectively, tolerated being mutated to tryptophan. Moreover, aspartate was tolerated at 15 out of 22 and 15 out of 21 non-aspartate M1 and M2 positions respectively. Periodicity in the pattern of tolerated vs. nontolerated mutations consistent with alpha helices or beta strands did not emerge convincingly from these data. We consider the possibility that parts of M1 and M2 may be in contact with water.

Self-consistent field approach to protein structure and stability. I: pH dependence of electrostatic contribution

Starting from the simple case of an external field acting on noninteracting particles, a formulation of the self-consistent field theory for treating proteins and unfolded protein chains with multiple interacting titratable groups is given. Electrostatic interactions between the titratable groups are approximated by a Debye-Huckel expression. Amino acid residues are treated as polarizable bodies with a single dielectric constant. Dielectric properties of protein molecules are described in terms of local dielectric constants determined by the space distribution of residue volume density around each ionized residue. Calculations are based on average charges of titratable groups, distance of separation between them, on their pKa's, residue volumes and on the local dielectric constant. A set of different residue volumes is used to analyze the influence of the permanent dipole of polar parts of the residue on calculated titration curves, electrostatic contribution to the free energy of protein stability, and pK shifts. Calculations with zero volumes--which means that charged portions of protein molecules are viewed as part of the high dielectric medium--give good agreement with experimental data. The theory was tested against most accurate approaches currently available for the calculation of the pKa's of ionizable groups based upon finite difference solutions of the Poisson-Boltzmann equation (FDPB). For 70 theoretically calculated pKa's in a total of six proteins the accuracy of the approach presented here is assessed by comparison of computed pKa's with that measured. The overall root-mean-square error is 0.79, compared to the value 0.89 obtained by FDPB approach given in the paper of Antosiewicz et al. (J. Mol. Biol. 238:415-436, 1994). The test of Debye-Huckel approximation for the electrostatic pair interactions shows that it is in excellent agreement with experimental data as well as the calculations of the FDPB and Tanford-Kirkwood methods on the pK shifts of His64 in the active site of subtilisin over the whole range of ionic strengths. (Gilson and Honig, Proteins 3:32-52, 1988; Russell et al., J.Mol.Biol. 193:803-813, 1987). The theory was also analytically and numerically tested on a simple models where the exact statistical mechanical treatment is still simple (Yang et al., Proteins 15:252-265, 1993; Bashford and karplus, J. Phys. Chem. 95:9556-9561, 1991).

Experimental measurement of the effective dielectric in the hydrophobic core of a protein

The dielectric inside a protein is a key physical determinant of the magnitude of electrostatic interactions in proteins. We have measured this dielectric phenomenologically, in terms of the dielectric that needs to be used with the Born equation in order to reproduce the observed pKa shifts induced by burial of an ionizable group in the hydrophobic core of a protein. Mutants of staphylococcal nuclease with a buried lysine residue at position 66 were engineered for this purpose. The pKa values of buried lysines were measured by difference potentiometry. The extent of coupling between the pKa and the global stability of the protein was evaluated by measuring pKa values in hyperstable forms of nuclease engineered to be 3.3 or 6.5 kcal mol-1 more stable than the wild type. The crystallographic structure of one mutant was determined to describe the environment of the buried lysine. The dielectrics that were measured range from 10 to 12. Published pKa values of buried ionizable residues in other proteins were analyzed in a similar fashion and the dielectrics obtained from these values are consistent with the ones measured in nuclease. These results argue strongly against the prevalent use of dielectrics of 4 or lower to describe the dielectric effect inside a protein in structure-based calculations of electrostatic energies with continuum dielectric models.

A large photolysis-induced pKa increase of the chromophore counterion in bacteriorhodopsin: implications for ion transport mechanisms of retinal proteins

The proton-pumping mechanism of bacteriorhodopsin is dependent on a photolysis-induced transfer of a proton from the retinylidene Schiff base chromophore to the aspartate-85 counterion. Up until now, this transfer was ascribed to a > 7-unit decrease in the pKa of the protonated Schiff base caused by photoisomerization of the retinal. However, a comparably large increase in the pKa of the Asp-85 acceptor also plays a role, as we show here with infrared measurements. Furthermore, the shifted vibrational frequency of the Asp-85 COOH group indicates a transient drop in the effective dielectric constant around Asp-85 to approximately 2 in the M photointermediate. This dielectric decrease would cause a > 40 kJ-mol-1 increase in free energy of the anionic form of Asp-85, fully explaining the observed pK alpha increase. An analogous photolysis-induced destabilization of the Schiff base counterion could initiate anion transport in the related protein, halorhodopsin, in which aspartate-85 is replaced by Cl- and the Schiff base proton is consequently never transferred.

Theory of electrostatic interactions in macromolecules

In the past year, substantial progress has been made in the modeling of electrostatic interactions in biomolecules. This review highlights advances in the following areas: first, the efficient computation of long-range electrostatic interactions in detailed molecular simulations; second, the application of the Poisson-Boltzmann electrostatic model in conformational analysis; third, the application of the Poisson-Boltzmann model in quantum chemistry calculations; fourth, the development of atomic parameters; and finally, the modeling of ionization equilibria in proteins.