Probing the Structural Basis for Enzyme-Substrate Recognition In Cu,Zn Superoxide Dismutase

A full understanding of enzyme-substrate interactions requires a detailed knowledge of their structural basis at atomic resolution. Crystallographic and biochemical data have been analyzed with coupled computational and computer graphic approaches to characterize the molecular basis for recognition of the superoxide anion substrate by Cu. Zn superoxide dismutase (SOD). Detailed analysis of the bovine SOD structure aligned with SOD sequences from 15 species provides new results concerning the significance and molecular basis for sequence conservation. Specific roles have been assigned for all 23 invariant residues and additional residues exhibiting functional equivalence. Sequence invariance is dominated by 15 residues that form the active site stcreochemistry. supporting a primary biological function of superoxide dismutation. Using data from crystallographic structures and site-directed mutants, we are testing the role of individual residues in the active site channel, including (in human SOD) Glu132, Glu133, Lys136, Thr137, and Arg 143. Electrostatic calculations incorporating molecular flexibility suggest that the region of positive electrostatic potential in and over the active site channel above the Cu ion sweeps through space during molecular motion to enhance the facilitated diffusion responsible for the enzyme's rapid catalytic rate.     

Electrostatic calculations and model-building suggest that DNA bound to CAP is sharply bent

Two observations suggest that DNA, upon binding to E. coli catabolite gene activator protein (CAP), is sharply bent by a total angle of at least 100-150 degrees: (1) The electrostatic potential field of CAP shows regions of positive potential that form a ramp on 3 sides of the protein. (2) The DNA binding site size as determined by DNA ethylation interference with binding, (Majors: "Control of the E. coli Lac Operon at the Molecular Level." Ph.D. Thesis, Harvard University, Cambridge, 1977) and by relative affinities of DNA fragments of various lengths (Liu-Johnson et al.: Cell 47:995-1005, 1986) requires severe bending of the DNA to maintain its favorable electrostatic contact with the protein.     

The effects of truncating long-range forces on protein dynamics

This paper considers the effects of truncating long-range forces on protein dynamics. Six methods of truncation that we investigate as a function of cutoff criterion of the long-range potentials are (1) a shifted potential; (2) a switching function; (3) simple atom-atom truncation based on distance; (4) simple atom-atom truncation based on a list which is updated periodically (every 25 steps); (5) simple group-group truncation based on distance; and (6) simple group-group truncation based on a list which is updated periodically (every 25 steps). Based on 70 calculations of carboxymyoglobin we show that the method and distance of long range cutoff have a dramatic effect on overall protein behavior. Evaluation of the different methods is based on comparison of a simulation's rms fluctuation about the average coordinates, the rms deviation from the average coordinates of a no cutoff simulation and from the X-ray structure of the protein. The simulations in which long-range forces are truncated by a shifted potential shows large rms deviations for cutoff criteria less than 14 A, and reasonable deviations and fluctuations at this cutoff distance or larger. Simulations using a switching function are investigated by varying the range over which electrostatic interactions are switched off. Results using a short switching function that switches off the potential over a short range of distances are poor for all cutoff distances. A switching function over a 5-9 A range gives reasonable results for a distance-dependent dielectric, but not using a constant dielectric. Both the atom-atom and group-group truncation methods based on distance shows large rms deviation and fluctuation for short cutoff distances, while for cutoff distances of 11 A or greater, reasonable results are achieved. Although comparison of these to distance-based truncation methods show surprisingly larger rms deviations for the group-group truncation, contrary to simulation studies of aqueous ionic solutions. The results of atom-atom or group-group list-based simulations generally appear to be less stable than the distance-based simulations, and require more frequent velocity scaling or stronger coupling to a heat bath.     

Do salt bridges stabilize proteins? A continuum electrostatic analysis

The electrostatic contribution to the free energy of folding was calculated for 21 salt bridges in 9 protein X-ray crystal structures using a continuum electrostatic approach with the DELPHI computer-program package. The majority (17) were found to be electrostatically destabilizing; the average free energy change, which is analogous to mutation of salt bridging side chains to hydrophobic isosteres, was calculated to be 3.5 kcal/mol. This is fundamentally different from stability measurements using pKa shifts, which effectively measure the strength of a salt bridge relative to 1 or more charged hydrogen bonds. The calculated effect was due to a large, unfavorable desolvation contribution that was not fully compensated by favorable interactions within the salt bridge and between salt-bridge partners and other polar and charged groups in the folded protein. Some of the salt bridges were studied in further detail to determine the effect of the choice of values for atomic radii, internal protein dielectric constant, and ionic strength used in the calculations. Increased ionic strength resulted in little or no change in calculated stability for 3 of 4 salt bridges over a range of 0.1-0.9 M. The results suggest that mutation of salt bridges, particularly those that are buried, to "hydrophobic bridges" (that pack at least as well as wild type) can result in proteins with increased stability. Due to the large penalty for burying uncompensated ionizable groups, salt bridges could help to limit the number of low free energy conformations of a molecule or complex and thus play a role in determining specificity (i.e., the uniqueness of a protein fold or protein-ligand binding geometry).     

Optimization of the electrostatic interactions in proteins of different functional and folding type

The 3-dimensional optimization of the electrostatic interactions between the charged amino acid residues was studied by Monte Carlo simulations on an extended representative set of 141 protein structures with known atomic coordinates. The proteins were classified by different functional and structural criteria, and the optimization of the electrostatic interactions was analyzed. The optimization parameters were obtained by comparison of the contribution of charge-charge interactions to the free energy of the native protein structures and for a large number of randomly distributed charge constellations obtained by the Monte Carlo technique. On the basis of the results obtained, one can conclude that the charge-charge interactions are better optimized in the enzymes than in the proteins without enzymatic functions. Proteins that belong to the mixed αβ folding type are electrostatically better optimized than pure α-helical or β-strand structures. Proteins that are stabilized by disulfide bonds show a lower degree of electrostatic optimization. The electrostatic interactions in a native protein are effectively optimized by rejection of the conformers that lead to repulsive charge-charge interactions. Particularly, the rejection of the repulsive contacts seems to be a major goal in the protein folding process. The dependence of the optimization parameters on the choice of the potential function was tested. The majority of the potential functions gave practically identical results.     

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 pK_as 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 pK_as in the protein, in contrast to multiple simulations for the free energy perturbation method. pK_int 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 pK_int shifts is discussed. © 1994 Wiley-Liss, Inc.     

Dominant role of local dipoles in stabilizing uncompensated charges on a sulfate sequestered in a periplasmic active transport protein.

Electrostatic interactions are among the key factors determining the structure and function of proteins. Here we report experimental results that illuminate the functional importance of local dipoles to these interactions. The refined 1.7-A X-ray structure of the liganded form of the sulfate-binding protein, a primary sulfate active transport receptor of Salmonella typhimurium, shows that the sulfate dianion is completely buried and bound by hydrogen bonds (mostly main-chain peptide NH groups) and van der Waals forces. The sulfate is also closely linked, via one of these peptide units, to a His residue. It is also adjacent to the N-termini of three alpha-helices, of which the two shortest have their C-termini "capped" by Arg residues. Site-directed mutagenesis of the recombinant Escherichia coli sulfate receptor had no effect on sulfate-binding activity when an Asn residue was substituted for the positively charged His and the two Arg (changed singly and together) residues. These results, combined with other observations, further solidify the idea that stabilization of uncompensated charges in a protein is a highly localized process that involves a collection of local dipoles, including those of peptide units confined to the first turns of helices. The contribution of helix macrodipoles appears insignificant.     

Comparison of the crystal structure of bacteriophage T4 lysozyme at low, medium, and high ionic strengths

Crystals of bacteriophage T4 lysozyme used for structural studies are routinely grown from concentrated phosphate solutions. It has been found that crystals in the same space group can also be grown from solutions containing 0.05 M imidazole chloride, 0.4 M sodium choride, and 30% polyethylene glycol 3500. These crystals, in addition, can also be equilibrated with a similar mother liquor in which the sodium chloride concentration is reduced to 0.025 M. The availability of these three crystal variants has permitted the structure of T4 lysozyme to be compared at low, medium, and high ionic strength. At the same time the X-ray structure of phage T4 lysozyme crystallized from phosphate solutions has been further refined against a new and improved X-ray diffraction data set. The structures of T4 lysozyme in the crystals grown with polyethylene glycol as a precipitant, regardless of the sodium chloride concentration, were very similar to the structure in crystals grown from concentrated phosphate solutions. The main differences are related to the formation of mixed disulfides between cysteine residues 54 and 97 and 2-mercaptoethanol, rather than to the differences in the salt concentration in the crystal mother liquor. Formation of the mixed disulfide at residue 54 resulted in the displacement of Arg-52 and the disruption of the salt bridge between this residue and Glu-62. Other than this change, no obvious alterations in existing salt bridges in T4 lysozyme were observed. Neither did the reduction in the ionic strength of the mother liquor result in the formation of new salt bridge interactions. These results are consistent with the ideas that a crystal structure determined at high salt concentrations is a good representation of the structure at lower ionic strengths, and that models of electrostatic interactions in proteins that are based on crystal structures determined at high salt concentrations are likely to be relevant at physiological ionic strengths.     

Electrostatic fields at the active site of ribulose-1,5-bisphosphate carboxylase.

A macroscopic approach has been employed to calculate the electrostatic potential field of nonactivated ribulose-1,5-bisphosphate carboxylase and of some complexes of the enzyme with activator and substrate. The overall electrostatic field of the L2-type enzyme from the photosynthetic bacterium Rhodospirillum rubrum shows that the core of the dimer, consisting of the two C-terminal domains, has a predominantly positive potential. These domains provide the binding sites for the negatively charged phosphate groups of the substrate. The two N-terminal domains have mainly negative potential. At the active site situated between the C-terminal domain of one subunit and the N-terminal domain of the second subunit, a large potential gradient at the substrate binding site is found. This might be important for polarization of chemical bonds of the substrate and the movement of protons during catalysis. The immediate surroundings of the activator lysine, K191, provide a positive potential area which might cause the pK value for this residue to be lowered. This observation suggests that the electrostatic field at the active site is responsible for the specific carbamylation of the epsilon-amino group of this lysine side chain during activation. Activation causes a shift in the electrostatic potential at the position of K166 to more positive values, which is reflected in the unusually low pK of K166 in the activated enzyme species. The overall shape of the electrostatic potential field in the L2 building block of the L8S8-type Rubisco from spinach is, despite only 30% amino acid homology for the L-chains, strikingly similar to that of the L2-type Rubisco from Rhodospirillum rubrum. A significant difference between the two species is that the potential is in general more positive in the higher plant Rubisco. In particular, the second phosphate binding site has a considerably more positive potential, which might be responsible for the higher affinity for the substrate of L8S8-type enzymes. The higher potential at this site might be due to two remote histidine residues, which are conserved in the plant enzymes.