Structural and thermodynamic analysis of the binding of solvent at internal sites in T4 lysozyme.

To investigate the structural and thermodynamic basis of the binding of solvent at internal sites within proteins a number of mutations were constructed in T4 lysozyme. Some of these were designed to introduce new solvent-binding sites. Others were intended to displace solvent from preexisting sites. In one case Val-149 was replaced with alanine, serine, cysteine, threonine, isoleucine, and glycine. Crystallographic analysis shows that, with the exception of isoleucine, each of these substitutions results in the binding of solvent at a polar site that is sterically blocked in the wild-type enzyme. Mutations designed to perturb or displace a solvent molecule present in the native enzyme included the replacement of Thr-152 with alanine, serine, cysteine, valine, and isoleucine. Although the solvent molecule was moved in some cases by up to 1.7 A, in no case was it completely removed from the folded protein. The results suggest that hydrogen bonds from the protein to bound solvent are energy neutral. The binding of solvent to internal sites within proteins also appears to be energy neutral except insofar as the bound solvent may prevent a loss of energy due to potential hydrogen bonding groups that would otherwise be unsatisfied. The introduction of a solvent-binding site appears to require not only a cavity to accommodate the water molecule but also the presence of polar groups to help satisfy its hydrogen-bonding potential. It may be easier to design a site to accommodate two or more water molecules rather than one as the solvent molecules can then hydrogen-bond to each other. For similar reasons it is often difficult to design a point mutation that will displace a single solvent molecule from the core of a protein.

Size versus polarizability in protein-ligand interactions: binding of noble gases within engineered cavities in phage T4 lysozyme.

To investigate the relative importance of size and polarizability in ligand binding within proteins, we have determined the crystal structures of pseudo wild-type and cavity-containing mutant phage T4 lysozymes in the presence of argon, krypton, and xenon. These proteins provide a representative sample of predominantly apolar cavities of varying size and shape. Even though the volumes of these cavities range up to the equivalent of five xenon atoms, the noble gases bind preferentially at highly localized sites that appear to be defined by constrictions in the walls of the cavities, coupled with the relatively large radii of the noble gases. The cavities within pseudo wild-type and L121A lysozymes each bind only a single atom of noble gas, while the cavities within mutants L133A and F153A have two independent binding sites, and the L99A cavity has three interacting sites. The binding of noble gases within two double mutants was studied to characterize the additivity of binding at such sites. In general, when a cavity in a protein is created by a "large-to-small" substitution, the surrounding residues relax somewhat to reduce the volume of the cavity. The binding of xenon and, to a lesser degree, krypton and argon, tend to expand the volume of the cavity and to return it closer to what it would have been had no relaxation occurred. In nearly all cases, the extent of binding of the noble gases follows the trend xenon>krypton>argon. Pressure titrations of the L99A mutant have confirmed that the crystallographic occupancies accurately reflect fractional saturation of the binding sites. The trend in noble gas affinity can be understood in terms of the effects of size and polarizability on the intermolecular potential. The plasticity of the protein matrix permits repulsion due to increased ligand size to be more than compensated for by attraction due to increased ligand polarizability. These results have implications for the mechanism of general anesthesia, the migration of small ligands within proteins, the detection of water molecules within apolar cavities and the determination of crystallographic phases.

Accurate calculation of the density of proteins.

On the basis of theoretical calculations, Andersson & Hovmöller have recently suggested that the long-established value of 1.35 g cm(-3) for the mean density of proteins should be revised to 1.22 g cm(-3) [Andersson & Hovmöller (2000), Acta Cryst. D56, 789-790]. To substantiate their assertion, these authors used the Voronoi algorithm to calculate the mean atomic volume for 30 representative protein structures. The Voronoi procedure requires that atoms of interest be bounded on all sides by other atoms. Volume calculations for surface atoms that are not surrounded or are only sparsely surrounded by other atoms either are not possible or may be unreliable. In an attempt to circumvent this problem, Andersson & Hovmöller rejected atoms with calculated volumes that were indeterminate or were greater than 50 A(3). In the present study, it is shown that this criterion is not sufficiently restrictive to ensure accurate volume determinations. When only strictly buried atoms are included in the volume calculations using the Voronoi algorithm, the mean density is found to be 1.47 +/- 0.05 g cm(-3). In addition, an alternate procedure based on the Connolly algorithm that permits all protein atoms to be included in volume calculations gives 1.43 +/- 0.03 g cm(-3) for the mean density of the same set of proteins. The latter two calculated values are mutually consistent and are in better agreement with the experimental value.

Structural and functional effects of apolar mutations of the distal valine in myoglobin.

High-resolution structures of the aquomet, deoxy, and CO forms of Ala68, Ile68, Leu68, and Phe68 sperm whale myoglobins have been determined by X-ray crystallography. These 12 new structures, plus those of wild-type myoglobin, have been used to interpret the effects of mutations at position 68 and the effects of cobalt substitution on the kinetics of O2, CO, and NO binding. Molecular dynamics simulations based on crystal structures have provided information about the time-dependent behavior of photolyzed ligands for comparison with picosecond geminate recombination studies. The Val68-->Ala mutation has little effect on the structure and function of myoglobin. In Ala68 deoxymyoglobin, as in the wild-type protein, a water molecule hydrogen-bonded to the N epsilon atom of the distal histidine restricts ligand binding and appears to be more important in regulating the function of myoglobin than direct steric interactions between the ligand and the C gamma atoms of the native valine side-chain. This distal pocket water molecule is displaced by the larger side-chains at position 68 in the crystal structures of Leu68 and Ile68 deoxymyoglobins. The Leu68 side-chain can rotate about its C alpha-C beta and C beta-C gamma bonds to better accommodate bound ligands, resulting in net increases in overall association rate constants and affinities due to the absence of the distal pocket water molecule. However, the flexibility of Leu68 makes simulation of picosecond NO recombination difficult since multiple starting conformations are possible. In the case of Ile68, rotation of the substituted side-chain is restricted due to branching at the beta carbon, and as a result, the delta methyl group is located close to the iron atom in both the deoxy and liganded structures. The favorable effect of displacing the distal pocket water molecule is offset by direct steric hindrance between the bound ligand and the terminal carbon atom of the isoleucine side-chain, resulting in net decreases in affinity for all three ligands and inhibition of geminate recombination which is reproduced in the molecular dynamics simulations. In Phe68 myoglobin, the benzyl side-chain is pointed away from the ligand binding site, occupying a region in the back of the distal pocket. As in wild-type and Ala68 myoglobins, a well-defined water molecule is found hydrogen bonded to the distal histidine in Phe68 deoxymyoglobin. This water molecule, in combination with the large size of the benzyl side-chain, markedly reduces the speed and extent of ligand movement into the distal pocket. (ABSTRACT TRUNCATED AT 400 WORDS)

Structural determinants of the stretching frequency of CO bound to myoglobin.

In order to assess the relative importance of polar versus steric interactions, infrared spectra and overall CO binding properties were measured at room temperature for 41 different recombinant myoglobins containing mutations at His64(E7), Val68(E11), Phe43(CD1), Arg45(CD3), Phe46(CD4), and Leu29(B10). The results were compared to the crystal structures of wild-type, Phe29, Val46, Ala68, Phe68, Gln64, Leu64, and Gly64 sperm whale CO-myoglobin and that of Thr68 pig CO-myoglobin. As observed in several previous studies, replacement of the distal histidine (His64) with aliphatic amino acids results in the appearance of a single IR band in the 1960-1970-cm-1 region and in large increases in CO affinity (KCO). More complex behavior is observed for Gly, Ala, Gln, Met, and Trp substitutions at position 64, but in each case there is a net increase in the intensity of this high-frequency component. Replacement of Val68 with Ala, Leu, Ile, and Phe produces little effect on the IR spectrum, whereas these mutations cause 20-fold changes in KCO, presumably due to steric effects. Replacement of Val68 with Thr decreases KCO 4-5-fold, whereas the position of the major IR band increases from 1945 to 1961 cm-1. Replacement of Val68 with Asn also causes a large decrease in KCO, but in this case, the peak position of the major IR band decreases from 1945 to 1916 cm-1. Nine replacements were made in the CD corner at positions 43, 45, and 46. All of the resultant mutants show increased stretching frequencies that can be correlated with movement of the imidazole side chain of His64 away from the bound ligand. All five substitutions at position 29 cause changes in the IR spectra. The Leu29-->Phe mutation had the largest effect, producing a single band centered at 1932 cm-1. Together these data demonstrate that there is little direct correlation between affinity, vCO, and Fe-C-O geometry. The major factor governing vCO appears to be the electrostatic potential surrounding the bound ligand and not steric hindrance. The presence of positive charges from proton donors, such as N epsilon of His64 and N delta of Asn68, cause a decrease in the bond order and stretching frequency of bound CO. In contrast, the negative portion of the Thr68 dipole points directly toward the bound ligand and increases the C-O bond order and stretching frequency. Movement of His64 away from the bound ligand or replacement of this residue with aliphatic amino acids prevents hydrogen-bonding interactions, causing vCO to increase.(ABSTRACT TRUNCATED AT 250 WORDS)

His64(E7)-->Tyr apomyoglobin as a reagent for measuring rates of hemin dissociation.

To develop an assay for hemin dissociation, His64(E7) was replaced by Tyr in sperm whale myoglobin producing a holoprotein with a distinct green color due to an intense absorption band at 600 nm. Val68(E11) was replaced by Phe in the same protein to increase its stability. When excess Tyr64-Val68 apoglobin is mixed with either metmyoglobin or methemoglobin, the solution turns from brown to green, and the absorbance changes can be used to measure complete time courses for hemin dissociation from either holoprotein. This assay has been used to measure rates of hemin dissociation from native metmyoglobin, four myoglobin mutants (Ala64(E7), Ala68(E11), Phe68(E11), and Glu45(CD3)), native methemoglobin, valence hybrid hemoglobins, and two mutant hemoglobins ((alpha(Gly-E7)beta(native))2, and (alpha(native)beta(Gly-E7))2). Two kinetic phases were observed for hemin dissociation from native human hemoglobin at pH 7.0 and 37 degrees C. Valence and mutant hybrid hemoglobins were used to assign the faster phase (k = 7.8 +/- 2.0 h-1) to hemin dissociation from ferric beta subunits and the slower (k = 0.6 +/- 0.15 h-1) to dissociation from alpha subunits. The corresponding rate for wild-type metmyoglobin is 0.007 +/- 0.004 h-1.

High-resolution crystal structures of distal histidine mutants of sperm whale myoglobin.

The highly conserved distal histidine residue (His64) of sperm whale myoglobin modulates the affinity of ligands. In an effort to fully characterize the effects of mutating residue 64, we have determined the high-resolution crystal structures of the Gly64, Val64, Leu64, Thr64 and Gln64 mutants in several liganded forms. Metmyoglobins with hydrophobic substitutions at residue 64 (Val64 and Leu64) lack a water molecule at the sixth coordination position, while those with polar amino acid residues at this position (wild-type and Gln64) retain a covalently bound water molecule. In the Thr64 mutant, the bound water position is only partially occupied. In contrast, mutating the distal histidine residue to glycine does not cause loss of the coordinated water molecule, because the hydrogen bond from the imidazole side-chain is replaced by one from a well-ordered solvent water molecule. Differences in water structure around the distal pocket are apparent also in the structures of deoxymyoglobin mutants. The water molecule that is hydrogen-bonded to the N epsilon atom of histidine 64 in wild-type deoxymyoglobin is not found in any of the position 64 mutant structures that were determined. Comparison of the carbonmonoxy structures of wild-type, Gly64, Leu64 and Gln64 myoglobins in the P6 crystal form shows that the conformation of the Fe-C-O complex is nearly linear and is independent of the identity of the amino acid residue at position 64. However, the effect of CO binding on the conformation of residue 64 is striking. Superposition of deoxy and carbonmonoxy structures reveals significant displacements of the residue 64 side-chain in the wild-type and Gln64 myoglobins, but no displacement in the Leu64 mutant. These detailed structural studies provide key insights into the mechanisms of ligand binding and discrimination in myoglobin.

A novel site-directed mutant of myoglobin with an unusually high O2 affinity and low autooxidation rate.

Mutants of sperm whale myoglobin were constructed at position 29 (B10 in helix notation) to examine the effects of distal pocket size on the rates of ligand binding and autooxidation. Leu29 was replaced with Ala, Val, and Phe using the synthetic gene and Escherichia coli expression system of Springer and Sligar (Springer, B. A., and Sligar, S. G. (1987) Proc. Natl. Acad. Sci. U. S. A. 84, 8961-8965). Structures of the ferric forms of Val29 and Phe29, and the oxy form of Phe29 myoglobin were determined to 1.7 A by x-ray crystallography. The ferric mutant proteins are remarkably isomorphous with the wild type protein except in the immediate vicinity of residue 29. Thus, the protein structure in the distal pocket of myoglobin can accommodate either a large "hole" (i.e. Ala or Val) or a large side chain (i.e. Phe) at position 29 without perturbation of tertiary structure. Phe29 oxymyoglobin is also identical to the native oxy protein in terms of overall structure and interactions between the bound O2 and His64, Val68, Phe43, and Ile107. The distance between the nearest side chain atom of residue 29 and the second atom of the bound oxygen molecule is 3.2 A in the Phe29 protein and 4.9 A in native myoglobin. The equilibrium constants for O2 binding to Ala29, Val29, and Leu29 (native) myoglobin are the same, approximately 1.0 x 10(6) M-1 at 20 degrees C, whereas that for the Phe29 protein is markedly greater, 15 x 10(6) M-1. This increase in affinity is due primarily to a 10-fold decrease in the O2 dissociation rate constant for the Phe29 mutant and appears to be the result of stabilizing interactions between the negative portion of the bound O2 dipole and the partially positive edge of the phenyl ring. Increasing the size of residue 29 causes large decreases in the rate of autooxidation of myoglobin: k(ox) = 0.24, 0.23, 0.055, and 0.005 h-1 for Ala29, Val29, Leu29 (native), and Phe29 myoglobin, respectively, in air at 37 degrees C. Thus, the Leu29----Phe mutation produces a reduced protein that is remarkably stable and is expressed in E. coli as 100% MbO2. The selective pressure to conserve Leu29 at the B10 position probably represents a compromise between reducing the rate of autooxidation and maintaining a large enough O2 dissociation rate constant to allow rapid oxygen release during respiration.