Finite Difference Poisson-Boltzmann Method Research Articles

The Role of Electrostatic Interactions in the Regulation of the Membrane Association of G Protein βγ Heterodimers

In this paper we report calculations of electrostatic interactions between the transducin (G(t)) betagamma heterodimer (G(t)betagamma) and phospholipid membranes. Although membrane association of G(t)betagamma is due primarily to the hydrophobic penetration into the membrane interior of a farnesyl chain attached to the gamma subunit, structural studies have revealed that there is a prominent patch of basic residues on the surface of the beta subunit surrounding the site of farnesylation that is exposed upon dissociation from the G(t)alpha subunit. Moreover, phosducin, which produces dissociation of G(t)betagamma from membranes, interacts directly with G(t)betagamma and introduces a cluster of acidic residues into this region. The calculations, which are based on the finite difference Poisson-Boltzmann method, account for a number of experimental observations and suggest that charged residues play a role in mediating protein-membrane interactions. Specifically, the calculations predict the following. 1) Favorable electrostatic interactions enhance the membrane partitioning due to the farnesyl group by an order of magnitude although G(t)betagamma has a large net negative charge (-12). 2) This electrostatic attraction positions G(t)betagamma so that residues implicated in mediating the interaction of G(t)betagamma with its membrane-bound effectors are close to the membrane surface. 3) The binding of phosducin to G(t)betagamma diminishes the membrane partitioning of G(t)betagamma by an order of magnitude. 4) Lowering the ionic strength of the solution converts the electrostatic attraction into a repulsion. Sequence analysis and homology model building suggest that our conclusions may be generalized to other Gbetagamma and phosducin isoforms as well.

Salt Effects on Ionization Equilibria of Histidines in Myoglobin

The salt dependence of histidine p K a values in sperm whale and horse myoglobin and in histidine-containing peptides was measured by 1H-NMR spectroscopy. Structure-based p K a calculations were performed with continuum methods to test their ability to capture the effects of solution conditions on p K a values. The measured p K a of most histidines, whether in the protein or in model compounds, increased by 0.3 pH units or more between 0.02 M and 1.5 M NaCl. In myoglobin two histidines (His 48 and His 36) exhibited a shallower dependence than the average, and one (His 113) showed a steeper dependence. The 1H-NMR data suggested that the salt dependence of histidine p K a values in the protein was determined primarily by the preferential stabilization of the charged form of histidine with increasing salt concentrations rather than by screening of electrostatic interactions. The magnitude and salt dependence of interactions between ionizable groups were exaggerated in p K a calculations with the finite-difference Poisson-Boltzmann method applied to a static structure, even when the protein interior was treated with arbitrarily high dielectric constants. Improvements in continuum methods for calculating salt effects on p K a values will require explicit consideration of the salt dependence of model compound p K a values used for reference in the calculations.

Empirical calculation of the relative free energies of peptide binding to the molecular chaperone DnaK.

We describe a methodology to calculate the relative free energies of protein-peptide complex formation. The interaction energy was decomposed into nonpolar, electrostatic and entropic contributions. A free energy-surface area relationship served to calculate the nonpolar free energy term. The electrostatic free energy was calculated with the finite difference Poisson-Boltzmann method and the entropic contribution was estimated from the loss in the conformational entropy of the peptide side chains. We applied this methodology to a series of DnaK*peptide complexes. On the basis of the single known crystal structure of the peptide-binding domain of DnaK with a bound heptapeptide, we modeled ten other DnaK*heptapeptide complexes with experimentally measured K(d) values from 0.06 microM to 11 microM, using molecular dynamics to refine the structures of the complexes. Molecular dynamic trajectories, after equilibration, were used for calculating the energies with greater accuracy. The calculated relative binding free energies were compared with the experimentally determined free energies. Linear scaling of the calculated terms was applied to fit them to the experimental values. The calculated binding free energies were between -7.1 kcal/mol and - 9.4 kcal/mol with a correlation coefficient of 0.86. The calculated nonpolar contributions are mainly due to the central hydrophobic binding pocket of DnaK for three amino acid residues. Negative electrostatic fields generated by the protein increase the binding affinity for basic residues flanking the hydrophobic core of the peptide ligand. Analysis of the individual energy contributions indicated that the nonpolar contributions are predominant compared to the other energy terms even for peptides with low affinity and that inclusion of the change in conformational entropy of the peptide side chains does not improve the discriminative power of the calculation. The method seems to be useful for predicting relative binding energies of peptide ligands of DnaK and might be applicable to other protein-peptide systems, particularly if only the structure of one protein-ligand complex is available.

Conformational analysis of single-base bulges in A-form DNA and RNA using a hierarchical approach and energetic evaluation with a continuum solvent model

The analysis and prediction of non-canonical structural motifs in RNA is of great importance for an understanding of the function and design of RNA structures. A hierarchical method has been employed to generate a large variety of sterically possible conformations for a single-base adenine bulge structure in A-form DNA and RNA. A systematic conformational search was performed on the isolated bulge motif and neighboring nucleotides under the constraint to fit into a continuous helical structure. These substructures were recombined with double-stranded DNA or RNA. Energy minimization resulted in more than 300 distinct bulge conformations. Energetic evaluation using a solvation model based on the finite-difference Poisson-Boltzmann method identified three basic classes of low-energy structures. The three classes correspond to conformations with the bulge base stacked between flanking nucleotides (I), location of the bulge base in the minor groove (II) and conformations with a continuous stacking of the flanking helices and a looped out bulge base (III). For the looped out class, two subtypes (IIIa and IIIb) with different backbone geometries at the bulge site could be distinguished. The conformation of lowest calculated energy was a class I structure with backbone torsion angles close to those in standard A-form RNA. Conformations very close to the extra-helical looped out bulge structure determined by X-ray crystallography were also among the low-energy structures. In addition, topologies observed in other experimentally determined bulge structures have been found among low-energy conformers. The implicit solvent model was further tested by comparing an uridine and adenine bulge flanked by guanine:cytosine base-pairs, respectively. In agreement with the experimental observation, a looped out form was found as the energetically most favorable form for the uridine bulge and a stacked conformation in case of the adenine bulge. The inclusion of solvation effects especially electrostatic reaction field contributions turned out to be critically important in order to select realistic low-energy bulge structures from a large number of sterically possible conformations. The results indicate that the approach might be useful to model the three-dimensional structure of non-canonical motifs embedded in double-stranded RNA, in particular, to restrict the number of possible conformations to a manageable number of conformers with energies below a certain threshold.

Molecular modelling investigation of wild-type and the R528H mutated segment IIS4 of human L-type voltage-gated calcium channels.

A molecular modelling study was performed in order to investigate the pathologically modified properties of L-type voltage-gated calcium channels caused by the arginine-to-histidine mutation at position 528 (R528H) in segment IIS4. For an appropriate consideration of the ionization state the finite difference Poisson-Boltzmann method was applied to compute the apparent pKa values of all titratable residues using standard conditions and an explicit lipophilic environment, respectively. Restrained molecular dynamics simulations were carried out for the alpha-helical transmembrane segments of the wild-type and the R528H mutant to explore their conformational behaviour. While both structures showed almost the same side-chain flexibility around the conserved residues, only the mutant partially formed a hydrogen bond from H528 to R531 during dynamics simulations. This local interaction not only causes a lower mobility of the directly involved residues but also leads to a global distortion of all positively charged amino acids of the mutant. Mostly affected is the side-chain of R534 that is shifted about 61 degrees closer to R531 and about 5 A in the direction of the cytoplasm. Subsequent examination of the molecular characteristics of this putative voltage sensor of the channel revealed considerable variations with regard to hydrogen bonding and electronic properties. Most obvious are the dramatic loss of a strong positive molecular electrostatic potential and the reduced hydrogen donor activities around position 528 of the mutant. How these results may be interpreted in relation to an enhanced inactivation rate is discussed, considering earlier findings at homologous voltage-gated potassium and sodium channels.

Calculation of HyHel10-lysozyme binding free energy changes: Effect of ten point mutations

The change in free energy of binding of hen egg white lysozyme (HEL) to the antibody HyHel-10 arising from ten point mutations in HEL (D101K, D101G, K96M, K97D, K97G, K97G, R21E, R21K, W62Y, and W63Y) was calculated using a combination of the finite difference Poisson-Boltzmann method for the electrostatic contribution, a solvent accessible surface area term for the non-polar contribution, and rotamer counting for the sidechain entropy contribution. Comparison of experimental and calculated results indicate that because of pKa shifts in some of the mutated residues, primarily those involving Aspartate or Glutamate, proton uptake or release occurs in binding. When this effect was incorporated into the binding free energy calculations, the agreement with experiment improved significantly, and resulted in a mean error of about 1.9 kcal/mole. Thus these calculations predict that there should be a significant pH dependence to the change in binding caused by these mutations. The other major contributions to binding energy changes comes from solvation and charge charge interactions, which tend to oppose each other. Smaller contributions come from nonpolar interactions and sidechain entropy changes. The structures of the HyHel-10-HEL complexes with mutant HEL were obtained by modeling, and the effect of the modeled structure on the calculations was also examined. "Knowledge based" modeling and automatic generation of models using molecular mechanics produced comparable results.

Calculation of HyHel10‐lysozyme binding free energy changes: Effect of ten point mutations

The change in free energy of binding of hen egg white lysozyme (HEL) to the antibody HyHel-10 arising from ten point mutations in HEL (D101K, D101G, K96M, K97D, K97G, K97G, R21E, R21K, W62Y, and W63Y) was calculated using a combination of the finite difference Poisson-Boltzmann method for the electrostatic contribution, a solvent accessible surface area term for the non-polar contribution, and rotamer counting for the sidechain entropy contribution. Comparison of experimental and calculated results indicate that because of pKa shifts in some of the mutated residues, primarily those involving Aspartate or Glutamate, proton uptake or release occurs in binding. When this effect was incorporated into the binding free energy calculations, the agreement with experiment improved significantly, and resulted in a mean error of about 1.9 kcal/mole. Thus these calculations predict that there should be a significant pH dependence to the change in binding caused by these mutations. The other major contributions to binding energy changes comes from solvation and charge charge interactions, which tend to oppose each other. Smaller contributions come from nonpolar interactions and sidechain entropy changes. The structures of the HyHel-10-HEL complexes with mutant HEL were obtained by modeling, and the effect of the modeled structure on the calculations was also examined. “Knowledge based” modeling and automatic generation of models using molecular mechanics produced comparable results. Proteins 33:39–48, 1998. © 1998 Wiley-Liss, Inc.

Electrostatic Contributions to Heat Capacity Changes of DNA-Ligand Binding

Significant heat capacity changes ( ΔC p) often accompany protein unfolding, protein binding, and specific DNA-ligand binding reactions. Such changes are widely used to analyze contributions arising from hydrophobic and polar hydration. Current models relate the magnitude of ΔC p to the solvent accessible surface area (ASA) of the molecule. However, for many binding systems—particularly those involving non-peptide ligands—these models predict a ΔC p that is significantly different from the experimentally measured value. Electrostatic interactions provide a potential source of heat capacity changes and do not scale with ASA. Using finite-difference Poisson-Boltzmann methods (FDPB), we have determined the contribution of electrostatics to the ΔC p associated with binding for DNA binding reactions involving the ligands DAPI, netropsin, lexitropsin, and the λ repressor binding domain.

Investigating the dielectric effects of channel pore water on the electrostatic barriers of the permeation ion by the finite difference Poisson-Boltzmann method.

In this paper, the finite difference Poisson-Boltzmann (FDPB) method with four dielectric constants is developed to study the effect of dielectric saturation on the electrostatic barriers of the permeation ion. In this method, the inner shape of the channel pore is explicitly represented, and the fact that the dielectric constant inside the channel pore is different from that of bulk water is taken into account. A model channel system which is a righthanded twist bundle with four alpha-helical segments is provided for this study. From the FDPB calculations, it is found that the difference of the ionic electrostatic solvation energy for wider domains depends strongly on the pore radius in the vicinity of the ion when the pore dielectric constant is changed from 78 to 5. However, the electrostatic solvation energy of the permeation ion can not be significantly affected by the dielectric constant in regions with small pore radii. Our results indicate that the local electrostatic interactions inside the ion channel are of major importance for ion electrostatic solvation energies, and the effect of dielectric saturation on the electrostatic barriers is coupled to the interior channel dimensions.

Titration calculations of foot-and-mouth disease virus capsids and their stabilities as a function of pH

Foot-and-mouth disease virus (FMDV), a non-enveloped picornavirus, is sensitive to acidic conditions. At pH values below 7 the icosahedral virus capsid, formed from 60 copies of a protomer containing four polypeptides (VP1 to 4), dissociates into 12 pentamers, releasing the viral RNA. Evidence suggests that this acid lability may assist FMDV cell entry via an endosomal pathway. Calculations of titration curves and pH-stability profiles are presented for three different strains of FMDV, O1BFS, A1061 and A22 Iraq, and compared with experimental data for complete virions and empty capsids (which lack RNA). The finite difference Poisson-Boltzmann method was used for the calculation of electrostatic free energies with the solvent treated as a dielectric continuum. The inter-pentamer interface in the virus is formed by two protomers related by 2-fold icosahedral symmetry. As a simple model for inter-pentamer interactions, a dimer and two separate protomers were compared. The association free energy was computed by integrating the difference between the titration curves of the two species. The calculations reproduced the observed decrease in capsid stability at acidic pH but not the difference in pH sensitivities of the two type A viruses. It is shown that only residues within 15 Å of the interface play a significant role in determining acid lability. For the experimentally studied pH range (5 to 7.6), histidine residues were found to dominate the pH-dependence of the stability. Two histidine residues in VP3, H142 and H145, are shown to have the greatest effect by virtue of their interactions with many polar residues across the inter-pentamer interface; the interaction of H142 with an α-helix in the opposite pentamer contributes only a small proportion of the destabilization energy.

Molecular docking of superantigens with class II major histocompatibility complex proteins.

The molecular recognition of two superantigens with class II major histocompatibility complex molecules was simulated by using protein-protein docking. Superantigens studied were staphylococcal enterotoxin B (SEB) and toxic shock syndrome toxin-1 (TSST-1) in their crystallographic assemblies with HLA-DR1. Rigid-body docking was performed sampling configurational space of the interfacial surfaces by employing a strategy of partitioning the contact regions on HLA-DR1 into separate molecular recognition units. Scoring of docked conformations was based on an electrostatic continuum model evaluated with the finite-difference Poisson-Boltzmann method. Estimates of nonpolar contributions were derived from the buried molecular surface areas. We found for both superantigens that docking the HLA-DR1 surface complementary with the SEB and TSST-1 contact regions containing a homologous hydrophobic surface loop provided sufficient recognition for the reconstitution of native-like conformers exhibiting the highest-scoring free energies. For the SEB complex, the calculations were successful in reproducing the total association free energy. A comparison of the free-energy determinants of the conserved hydrophobic contact residue indicates functional similarity between the two proteins for this interface. Though both superantigens share a common global association mode, differences in binding topology distinguish the conformational specificities underlying recognition.

Simulation of electrostatic and hydrodynamic properties of Serratia endonuclease.

We analyze the electrostatic and hydrodynamic properties of a nuclease from the pathogenic gram-negative bacterium Serratia marcescens using finite-difference Poisson-Boltzmann methods for electrostatic calculations and a bead-model approach for diffusion coefficient calculations. Electrostatic properties are analyzed for the enzyme in monomeric and dimeric forms and also in the context of DNA binding by the nuclease. Our preliminary results show that binding of a double-stranded DNA dodecamer by nuclease causes an overall shift in the charge of the protein by approximately three units of elementary charge per monomer, resulting in a positively charged protein at physiologic pH. In these calculations, the free enzyme was found to have a negative (-1 e) charge per monomer at pH 7. The most dramatic shift in pKa involves His 89 whose pKa increases by three pH units upon DNA binding. This shift leads to a protonated residue at pH 7, in contrast to the unprotonated form in the free enzyme. DNA binding also leads to a decrease in the energetic distances between the most stable protonation states of the enzyme. Dimerization has no significant effect on the electrostatic properties of each of the monomers for both free enzyme and that bound to DNA. Results of hydrodynamic calculations are consistent with the dimeric form of the enzyme in solution. The computed translational diffusion coefficient for the dimer model of the enzyme is in very good agreement with measurements from light scattering experiments. Preliminary electrooptical calculations indicate that the dimer should possess a large dipole moment (approximately 600 Debye units) as well as substantial optical anisotropy (limiting reduced linear electric dichroism of about 0.3). Therefore, this system may serve as a good model for investigation of electric and hydrodynamic properties by relaxation electrooptical experiments.

Electrostatic interactions in hirudin-thrombin binding

Hirudin is a good anticoagulant owing to potent inhibition of the serine protease thrombin. An aspartate- and glutamate-rich portion of hirudin plays an important part in its tight binding to thrombin through a ladder of salt bridges, and these residues have previously been mutated to asparagine or glutamine. Detailed calculations of the electrostatic contribution to changes in binding from these mutations have been performed using the finite-difference Poisson-Boltzmann method which include charge-charge interactions, solvation interactions, the residual electrostatic interaction of mutant residues, p K a shifts, and ionic strength. Single mutant effects on binding energy were close to experimental values, except for the D55N mutant whose effect is overestimated, perhaps because of displacement of a bound chloride ion from the site where it binds. Multiple mutation values were generally overestimated. The effect of p K a shifts upon the binding is significant for one hirudin residue E58, but this appears to be due to a poor salt bridge with thrombin caused by crystal contacts. Electrostatic interaction between the acidic residues is unfavorable. However, analysis of experimental multiple mutation/single mutation data shows apparently negative interactions between these residues, from which it is concluded that structural changes can occur in the complex to relieve an unfavorable interaction when more than one acidic residue is mutated. In all cases, there is a loss in stability of the complex from mutations due to loss of favorable charge-charge interactions with thrombin, but this is largely compensated for by reduced unfavorable desolvation interactions, and by residual polar interactions in the Asn/Gln mutants.

A compound derived from prodocarpic acid inhibits influenza a virus replication in vitro: K.A. Staschke, S.D. Hatch, J.C. Tang, W.J. Hornback, S. Mauldin, J.E. Munroe, J.M. Colacino∗, and M.A. Muesing∗. Lilly Research Laboratories, Indianapolis, IN USA

The pH-induced conformational change of influenza virus hemagglutinin (HA) has been investigated by calculating the change of electrostatic energy of the fragment of HA2 upon pH change. The average charge and electrostatic free energy are calculated as a function of pH for the fusion peptide (residues 1–20 of HA2) and the polypeptide of residues 54–77 of HA2 by using the finite difference Poisson-Boltzmann method. It is found that as pH decreases from 8 to 5, the electrostatic free energy of the fusogenic state is lowered by ∼2 kcal/\ and the fusogenic state is less ionized compared to that of the native state for both polypeptides. For the fusion peptide at the fusogenic state, most of ionizable residues are neutral at acidic pH except Glu-11. For the polypeptide of residues 54–77 at the fusogenic state, most of residues except Glu-74 and His-64 are fully charged between pH 5 and pH 8.

Characterization of a new class of respiratory syncytial virus agents inhibiting the fusion process: B. O'Hara 1, G. Ellestad 1, J. Morin 1, A. Gazumyan 1, J. LaRocque 1, B. Mitsner 1, U. Raifeld 1, A. Nikitenko 1, W.-D. Ding 1, P. Wyde 2, B. Gilbert 2, and Y. Gluzman 1 1 Wyeth-Ayerst Research, Pearl River NY 10965, 2 Baylor

The pH-induced conformational change of influenza virus hemagglutinin (HA) has been investigated by calculating the change of electrostatic energy of the fragment of HA2 upon pH change. The average charge and electrostatic free energy are calculated as a function of pH for the fusion peptide (residues 1–20 of HA2) and the polypeptide of residues 54–77 of HA2 by using the finite difference Poisson-Boltzmann method. It is found that as pH decreases from 8 to 5, the electrostatic free energy of the fusogenic state is lowered by ∼2 kcal/\ and the fusogenic state is less ionized compared to that of the native state for both polypeptides. For the fusion peptide at the fusogenic state, most of ionizable residues are neutral at acidic pH except Glu-11. For the polypeptide of residues 54–77 at the fusogenic state, most of residues except Glu-74 and His-64 are fully charged between pH 5 and pH 8.

A free energy analysis of nucleic acid base stacking in aqueous solution

This paper reports a theoretical study of the free energy contributions to nucleic acid base stacking in aqueous solution. Electrostatic interactions are treated by using the finite difference Poisson-Boltzmann method and nonpolar effects are treated with explicit calculation of van der Waals interactions and/or free energy-surface area relationships. Although for some pairs of bases there is a favorable Coulombic interaction in the stacked conformation, generally the net effect of electrostatic interactions is to oppose stacking. This result is caused by the loss of favorable base-solvent electrostatic interactions, that accompany the partial removal of polar atoms from water in the stacked conformation. Nonpolar interactions, involving the hydrophobic effect and enhancement of van der Waals interactions caused by close-packing, drive stacking. The calculations qualitatively reproduce the experimental dependence of stacking free energy on purine-pyrimidine composition.

Computation of the dipole moments of proteins

A simple and computationally feasible procedure for the calculation of net charges and dipole moments of proteins at arbitrary pH and salt conditions is described. The method is intended to provide data that may be compared to the results of transient electric dichroism experiments on protein solutions. The procedure consists of three major steps: (i) calculation of self energies and interaction energies for ionizable groups in the protein by using the finite-difference Poisson-Boltzmann method, (ii) determination of the position of the center of diffusion (to which the calculated dipole moment refers) and the extinction coefficient tensor for the protein, and (iii) generation of the equilibrium distribution of protonation states of the protein by a Monte Carlo procedure, from which mean and root-mean-square dipole moments and optical anisotropies are calculated. The procedure is applied to 12 proteins. It is shown that it gives hydrodynamic and electrical parameters for proteins in good agreement with experimental data.

Free Energy Determinants of Secondary Structure Formation: I. α-Helices

The Zimm-Bragg parameters sand σ are calculated for the helix-coil transition of poly- l-alanine. The theoretical approach involves evaluating gas phase conformational energies for both coil and helical states using the CHARMM potential function and accounting for solvation effects with various continuum solvation models. Conformational free energies are then incorporated into a formalism developed by Go et al.for the calculation of s and σ. Calculated values for both s and σ as well as the enthalpy change associated with helix formation are in good agreement with experimental data when the Finite Difference Poisson-Boltzmann (FDPB) method is used to treat solvent effects. The driving force for the helix-coil transition is analyzed in terms of individual free energy components. Hydrogen bond formation is found to contribute little to helix stability because the internal hydrogen bonding energy is largely canceled by the large free energy cost associated with removing polar groups from water. The entropic cost associated with fixing backbone dihedral angles in the helical conformation is found to be ∼7 e.u./residue (about 2 kcal/mol at room temperature). The major driving force favoring helix formation can be associated with interactions including enhanced van der Waals interactions in the close-packed helix conformation and the hydrophobic effect. These contribute about 2 kcal/mol favoring the helical state. The differences in helical propensities between alanine and glycine are attributed primarily to hydrophobic and packing interactions involving the C βwith a smaller contribution arising from increased conformational freedom for glycine in the coil state. The description of helix formation presented here is consistent with previous conclusions regarding tertiary structure formation which suggest that hydrophobic and close-packed interactions provide stability while hydrogen bond formation constitutes a structural constraint imposed by the high free energy cost associated with burying unsatisfied hydrogen bonding groups. α-Helix formation may thus be viewed as a form of hydrophobic collapse constrained by the requirement that polar groups be either exposed to solvent or form hydrogen bonds. More generally it appears from this study that for a folding model to be a realistic, it must properly account for the chemical nature of the polypeptide chain, particularly the solvation energetics of amide groups.

A structural model for the prostate disease marker, human prostate‐specific antigen

Prostate-specific antigen (PSA) provides an excellent serum marker for prostate cancer, the most frequent form of cancer in American males. PSA is a 237-residue protease based on sequence homology to kallikrein-like enzymes. To predict the 3-dimensional structure of PSA, homology modeling studies were performed based on sequence and structural alignments with tonin, pancreatic kallikrein, chymotrypsin, and trypsin. The structurally conserved regions of the 4 reference X-ray proteins provided the core structure of PSA, whereas the loop structures were modeled on the loops of tonin and kallikrein. The unique "kallikrein loop" insert, between Ser 95b and Pro 95k of kallikrein, was constructed using molecular mechanics, dynamics, and electrostatics calculations. In the resulting PSA structure, the catalytic triad, involving residues His 57, Asp 102, and Ser 195, and hydrophobic and electrostatic interactions typical of serine proteases were extremely well conserved. Similarly, the 5-disulfide bonds of kallikrein were also conserved in PSA. These results, together with the fact that no major steric clashes arose during the modeling process, provide strong evidence for the validity of the PSA model. Calculation of the electrostatic potential contours of kallikrein and PSA was carried out using the finite difference Poisson-Boltzmann method. The calculations revealed matching areas of negative potential near the catalytic triad, but differences in the positive potential surrounding the active site. The PSA glycosylation site, Asn 61, is fully accessible to the solvent and is enclosed in a positive region of the isopotential map. The bottom of the substrate specificity pocket, residue S1, is a serine (Ser 189) as in chymotrypsin, rather than aspartate (Asp 189) as in tonin, kallikrein, and trypsin. This fact, plus other features of the S1 binding-pocket region, suggest that PSA would prefer substrates with hydrophobic residues at the P1 position. The location of a potential zinc ion binding site involving the side chain of histidines 91, 101, and 233 is also suggested. This PSA model should facilitate the understanding and prediction of structural and functional properties of this important cancer marker.

Intrinsic pKas of ionizable residues in proteins: an explicit solvent calculation for lysozyme.

Molecular dynamics simulations of triclinic hen egg white lysozyme in aqueous solution were performed to calculate the intrinsic pKas of 14 ionizable residues. An all-atom model was used for both solvent and solute, and a single 180 ps simulation in conjunction with a Gaussian fluctuation analysis method was used. An advantage of the Gaussian fluctuation method is that it only requires a single simulation of the system in a reference state to calculate all the pKas in the protein, in contrast to multiple simulations for the free energy perturbation method. pKint shifts with respect to reference titratable residues were evaluated and compared to results obtained using a finite difference Poisson-Boltzmann (FDPB) method with a continuum solvent model; overall agreement with the direction of the shifts was generally observed, though the magnitude of the shifts was typically larger with the explicit solvent model. The contribution of the first solvation shell to the total charging free energies of the titratable groups was explicitly evaluated and found to be significant. Dielectric shielding between pairs of titratable groups was examined and found to be smaller than expected. The effect of the approximations used to treat the long-range interactions on the pKint shifts is discussed.