Atomic size packing defects in proteins

Conserved buried water molecules enable the β-trefoil architecture

Available high-resolution crystal structures for the family of β-trefoil proteins in the structural databank were queried for buried waters. Such waters were classified as either: (a) unique to a particular domain, family, or superfamily or (b) conserved among all β-trefoil folds. Three buried waters conserved among all β-trefoil folds were identified. These waters are related by the threefold rotational pseudosymmetry characteristic of this protein architecture (representing three instances of an identical structural environment within each repeating trefoil-fold motif). The structural properties of this buried water are remarkable and include: residing in a cavity space no larger than a single water molecule, exhibiting a positional uncertainty (i.e., normalized B-factor) substantially lower than the average Cα atom, providing essentially ideal H-bonding geometry with three solvent-inaccessible main chain groups, simultaneously serving as a bridging H-bond for three different β-strands at a point of secondary structure divergence, and orienting conserved hydrophobic side chains to form a nascent core-packing group. Other published work supports an interpretation that these interactions are key to the formation of an efficient folding nucleus and folded thermostability. The fundamental threefold symmetric structural element of the β-trefoil fold is therefore, surprisingly, a buried water molecule.

Relationship between the magnitude of IgE production in mice and conformational stability of the house dust mite allergen, Der p 2

BACKGROUND

Protein antigens are degraded by endosomal protease in antigen presentation cell. T cells recognize peptides derived from antigen proteins bound to class II major histocompatibility complex molecules. We previously reported that an increase in the conformational stability of an antigen depressed its immunogenicity. However, there is little information on antigens with differences in molecular properties such as net charges and molecular weight.

METHODS

Denaturation experiments against guanidine hydrochloride. The serum IgE levels to protein antigens at 35days after the first immunization analyzed using ELISA.

RESULTS

The Der p 2 mutations in which Ile13 is mutated to Ala (I13A) and Ala122 is mutated to Ile (A122I) were shown to have lower and higher conformational stability than the wild-type, respectively, by denaturation experiments. The amount of IgE production by the less stable I13A mutant was higher and that of the stable A122I mutant was lower than that of the wild-type.

CONCLUSION

Our results suggest that the increased conformational stability of Der p 2 depressed the IgE production in mice.

GENERAL SIGNIFICANCE

These findings should provide a milestone for the engineering of allergen vaccines.

Statistical survey of the buried waters in the Protein Data Bank

The structures of buried water molecules were studied in an ensemble of high-quality and non-redundant protein crystal structures. Buried water molecules were clustered and classified in lake-like clusters, which are completely isolated from the bulk solvent, and bay-like clusters, which are in contact with the bulk solvent through a surface water molecule. Buried water molecules are extremely common: lake-like clusters are found in 89 % of the protein crystal structures and bay-like clusters in 93 %. Clusters with only one water molecule are much more common than larger clusters. Both cluster types incline to be surrounded by loop residues, and to a minor extent by residues in extended secondary structure. Helical residues on the contrary do not tend to surround clusters of buried water molecules. One buried water molecule is found every 30-50 amino acid residues, depending on the secondary structures that are more abundant in the protein. Both main- and side-chain atoms are in contact with buried waters; they form four hydrogen bonds with the first water and 1-1.5 additional hydrogen bond for each additional water in the cluster. Consequently, buried water molecules appear to be firmly packed and rigid like the protein atoms. In this regard, it is remarkable to observe that prolines often surround water molecules buried in the protein interior. Interestingly, clusters of buried water molecules tend to be just beneath the protein surface. Moreover, water molecules tend to form a one-dimensional wire rather than more compact arrangements. This agrees with recent evidence of the mechanisms of solvent exchange between internal cavities and bulk solvent.

How structure defines affinity in protein-protein interactions

Protein-protein interactions (PPI) in nature are conveyed by a multitude of binding modes involving various surfaces, secondary structure elements and intermolecular interactions. This diversity results in PPI binding affinities that span more than nine orders of magnitude. Several early studies attempted to correlate PPI binding affinities to various structure-derived features with limited success. The growing number of high-resolution structures, the appearance of more precise methods for measuring binding affinities and the development of new computational algorithms enable more thorough investigations in this direction. Here, we use a large dataset of PPI structures with the documented binding affinities to calculate a number of structure-based features that could potentially define binding energetics. We explore how well each calculated biophysical feature alone correlates with binding affinity and determine the features that could be used to distinguish between high-, medium- and low- affinity PPIs. Furthermore, we test how various combinations of features could be applied to predict binding affinity and observe a slow improvement in correlation as more features are incorporated into the equation. In addition, we observe a considerable improvement in predictions if we exclude from our analysis low-resolution and NMR structures, revealing the importance of capturing exact intermolecular interactions in our calculations. Our analysis should facilitate prediction of new interactions on the genome scale, better characterization of signaling networks and design of novel binding partners for various target proteins.

Molecular interactions within the halophilic, thermophilic, and mesophilic prokaryotic ribosomal complexes: clues to environmental adaptation

Using the available crystal structures of 50S ribosomal subunits from three prokaryotic species: Escherichia coli (mesophilic), Thermus thermophilus (thermophilic), and Haloarcula marismortui (halophilic), we have analyzed different structural features of ribosomal RNAs (rRNAs), proteins, and of their interfaces. We have correlated these structural features with the environmental adaptation strategies of the corresponding species. While dense intra-rRNA packing is observed in thermophilic, loose intra-rRNA packing is observed in halophilic (both compared to mesophilic). Interestingly, protein-rRNA interfaces of both the extremophiles are densely packed compared to that of the mesophilic. The intersubunit bridge regions are almost devoid of cavities, probably ensuring the proper formation of each bridge (by not allowing any loosely packed region nearby). During rRNA binding, the ribosomal proteins experience some structural transitions. Here, we have analyzed the intrinsically disordered and ordered regions of the ribosomal proteins, which are subjected to such transitions. The intrinsically disordered and disorder-to-order transition sites of the thermophilic and mesophilic ribosomal proteins are simultaneously (i) highly conserved and (ii) slowly evolving compared to rest of the protein structure. Although high conservation is observed at such sites of halophilic ribosomal proteins, but slow rate of evolution is absent. Such differences between thermophilic, mesophilic, and halophilic can be explained from their environmental adaptation strategy. Interestingly, a universal biophysical principle evident by a linear relationship between the free energy of interface formation, interface area, and structural changes of r-proteins during assembly is always maintained, irrespective of the environmental conditions.

A comparison of structural and evolutionary attributes of Escherichia coli and Thermus thermophilus small ribosomal subunits: signatures of thermal adaptation

Here we compare the structural and evolutionary attributes of Thermus thermophilus and Escherichia coli small ribosomal subunits (SSU). Our results indicate that with few exceptions, thermophilic 16S ribosomal RNA (16S rRNA) is densely packed compared to that of mesophilic at most of the analogous spatial regions. In addition, we have located species-specific cavity clusters (SSCCs) in both species. E. coli SSCCs are numerous and larger compared to T. thermophilus SSCCs, which again indicates densely packed thermophilic 16S rRNA. Thermophilic ribosomal proteins (r-proteins) have longer disordered regions than their mesophilic homologs and they experience larger disorder-to-order transitions during SSU-assembly. This is reflected in the predicted higher conformational changes of thermophilic r-proteins compared to their mesophilic homologs during SSU-assembly. This high conformational change of thermophilic r-proteins may help them to associate with the 16S ribosomal RNA with high complementary interfaces, larger interface areas, and denser molecular contacts, compared to those of mesophilic. Thus, thermophilic protein-rRNA interfaces are tightly associated with 16S rRNA than their mesophilic homologs. Densely packed 16S rRNA interior and tight protein-rRNA binding of T. thermophilus (compared to those of E. coli) are likely the signatures of its thermal adaptation. We have found a linear correlation between the free energy of protein-RNA interface formation, interface size, and square of conformational changes, which is followed in both prokaryotic and eukaryotic SSU. Disorder is associated with high protein-RNA interface polarity. We have found an evolutionary tendency to maintain high polarity (thereby disorder) at protein-rRNA interfaces, than that at rest of the protein structures. However, some proteins exhibit exceptions to this general trend.

Cavities in protein-DNA and protein-RNA interfaces

An analysis of cavities present in protein–DNA and protein–RNA complexes is presented. In terms of the number of cavities and their total volume, the interfaces formed in these complexes are akin to those in transient protein–protein heterocomplexes. With homodimeric proteins protein–DNA interfaces may contain cavities involving both the protein subunits and DNA, and these are more than twice as large as cavities involving a single protein subunit and DNA. A parameter, cavity index, measuring the degree of surface complementarity, indicates that the packing of atoms in protein–protein/DNA/RNA is very similar, but it is about two times less efficient in the permanent interfaces formed between subunits in homodimers. As within the tertiary structure and protein–protein interfaces, protein–DNA interfaces have a higher inclination to be lined by β-sheet residues; from the DNA side, base atoms, in particular those in minor grooves, have a higher tendency to be located in cavities. The larger cavities tend to be less spherical and solvated. A small fraction of water molecules are found to mediate hydrogen-bond interactions with both the components, suggesting their primary role is to fill in the void left due to the local non-complementary nature of the surface patches.

Dielectric relaxation of cytochrome c oxidase: Comparison of the microscopic and continuum models

We have studied a charge-insertion process that models the deprotonation of a histidine side chain in the active site of cytochrome c oxidase (CcO) using both the continuum electrostatic calculations and the microscopic simulations. The group of interest is a ligand to Cu(B) center of CcO, which has been previously suggested to play the role of the proton pumping element in the enzyme; the group is located near a large internal water cavity in the protein. Using the nonpolarizable Amber-99 force field in molecular dynamics (MD) simulations, we have calculated the nuclear part of the reaction-field energy of charging of the His group and combined it with the electronic part, which we estimated in terms of the electronic continuum (EC) model, to obtain the total reaction-field energy of charging. The total free energy obtained in this MDEC approach was then compared with that calculated using pure continuum electrostatic model with variable dielectric parameters. The dielectric constant for the "dry" protein and that of the internal water cavity of CcO were determined as those parameters that provide best agreement between the continuum and microscopic MDEC model. The nuclear (MD) polarization alone (without electronic part) of a dry protein was found to correspond to an unphysically low dielectric constant of only about 1.3, whereas the inclusion of electronic polarizability increases the protein dielectric constant to 2.6-2.8. A detailed analysis is presented as to how the protein structure should be selected for the continuum calculations, as well as which probe and atomic radii should be used for cavity definition. The dielectric constant of the internal water cavity was found to be 80 or even higher using "standard" parameters of water probe radius, 1.4 A, and protein atomic radii from the MD force field for cavity description; such high values are ascribed to the fact that the standard procedure produces unphysically small cavities. Using x-ray data for internal water in CcO, we have explored optimization of the parameters and the algorithm of cavity description. For Amber radii, the optimal probe size was found to be 1.25 A; the dielectric of water cavity in this case is in the range of 10-16. The most satisfactory cavity description, however, was achieved with ProtOr atomic radii, while keeping the probe radius to be standard 1.4 A. In this case, the value of cavity dielectric constant was found to be in the range of 3-6. The obtained results are discussed in the context of recent calculations and experimental measurements of dielectric properties of proteins.

Cavities and atomic packing in protein structures and interfaces

A comparative analysis of cavities enclosed in a tertiary structure of proteins and interfaces formed by the interaction of two protein subunits in obligate and non-obligate categories (represented by homodimeric molecules and heterocomplexes, respectively) is presented. The total volume of cavities increases with the size of the protein (or the interface), though the exact relationship may vary in different cases. Likewise, for individual cavities also there is quantitative dependence of the volume on the number of atoms (or residues) lining the cavity. The larger cavities tend to be less spherical, solvated, and the interfaces are enriched in these. On average 15 A(3) of cavity volume is found to accommodate single water, with another 40-45 A(3) needed for each additional solvent molecule. Polar atoms/residues have a higher propensity to line solvated cavities. Relative to the frequency of occurrence in the whole structure (or interface), residues in beta-strands are found more often lining the cavities, and those in turn and loop the least. Any depression in one chain not complemented by a protrusion in the other results in a cavity in the protein-protein interface. Through the use of the Voronoi volume, the packing of residues involved in protein-protein interaction has been compared to that in the protein interior. For a comparable number of atoms the interface has about twice the number of cavities relative to the tertiary structure.

Pressure- and temperature-induced unfolding studies: thermodynamics of core hydrophobicity and packing of ribonuclease A

In this work we demonstrate that heat and pressure induce only slightly different energetic changes in the unfolded state of RNase A. Using pressure and temperature as denaturants on a significant number of variants, and by determining the free energy of unfolding at different temperatures, we estimated the stability of variants unable to complete the unfolding transition owing to the experimental conditions required for pressure experiments. The overall set of results allowed us to map the contributions to stability of the hydrophobic core residues of RNase A, with the positions most critical for stability being V54, V57, I106 and V108. We also show that the stability differences can be attributed to both hydrophobic interactions and packing density with an equivalent energetic magnitude. The main hydrophobic core of RNase A is tightly packed, as shown by the small-to-large and isosteric substitutions. In addition, we found that large changes in the number of methylene groups have non-additive positive stability interaction energies that are consistent with exquisite tight core packing and rearrangements of van der Waals' interactions in the protein interior, even after drastic deleterious substitutions.

Standard atomic volumes in double-stranded DNA and packing in protein--DNA interfaces

Standard volumes for atoms in double-stranded B-DNA are derived using high resolution crystal structures from the Nucleic Acid Database (NDB) and compared with corresponding values derived from crystal structures of small organic compounds in the Cambridge Structural Database (CSD). Two different methods are used to compute these volumes: the classical Voronoi method, which does not depend on the size of atoms, and the related Radical Planes method which does. Results show that atomic groups buried in the interior of double-stranded DNA are, on average, more tightly packed than in related small molecules in the CSD. The packing efficiency of DNA atoms at the interfaces of 25 high resolution protein-DNA complexes is determined by computing the ratios between the volumes of interfacial DNA atoms and the corresponding standard volumes. These ratios are found to be close to unity, indicating that the DNA atoms at protein-DNA interfaces are as closely packed as in crystals of B-DNA. Analogous volume ratios, computed for buried protein atoms, are also near unity, confirming our earlier conclusions that the packing efficiency of these atoms is similar to that in the protein interior. In addition, we examine the number, volume and solvent occupation of cavities located at the protein-DNA interfaces and compared them with those in the protein interior. Cavities are found to be ubiquitous in the interfaces as well as inside the protein moieties. The frequency of solvent occupation of cavities is however higher in the interfaces, indicating that those are more hydrated than protein interiors. Lastly, we compare our results with those obtained using two different measures of shape complementarity of the analysed interfaces, and find that the correlation between our volume ratios and these measures, as well as between the measures themselves, is weak. Our results indicate that a tightly packed environment made up of DNA, protein and solvent atoms plays a significant role in protein-DNA recognition.

Stabilization of hen egg white lysozyme by a cavity-filling mutation

Stabilization of a protein using cavity-filling strategy has hardly been successful because of unfavorable van der Waals contacts. We succeeded in stabilizing lysozymes by cavity-filling mutations. The mutations were checked by a simple energy minimization in advance. It was shown clearly that the sum of free energy change caused by the hydrophobicity and the cavity size was correlated very well with protein stability. We also considered the aromatic-aromatic interaction. It is reconfirmed that the cavity-filling mutation in a hydrophobic core is a very useful method to stabilize a protein when the mutation candidate is selected carefully.

The protein state of matter

Three ways are generally used to visualize proteins: (1) a static model in which the atomic positions are defined, (2) a dynamic model taking into account fluctuations, and (3) a reactive model that reflects the internal and external electric fields of the molecule. The properties of chromophoric prosthetic groups can be probed by optical spectroscopy, and when high resolution techniques are used, the results reveal information about the local electric fields in proteins, as influenced and determined by atomic positions and dynamics.

Conformational stabilities of Escherichia coli RNase HI variants with a series of amino acid substitutions at a cavity within the hydrophobic core

Escherichia coli ribonuclease HI has a cavity within the hydrophobic core. Two core residues, Ala52 and Val74, resided at both ends of this cavity. We have constructed a series of single mutant proteins at Ala52, and double mutant proteins, in which Ala52 was replaced by Gly, Val, Ile, Leu, or Phe, and Val74 was replaced by Ala or Leu. All of these mutant proteins, except for A52W, A52R, and A52G/V74A, were overproduced and purified. Measurement of the thermal denaturations of the proteins at pH 3.2 by CD suggests that the cavity is large enough to accommodate three methyl or methylene groups without creating serious strains. A correlation was observed between the protein stability and the hydrophobicity of the substituted residue. As a result, a number of the mutant proteins were more stable than the wild-type protein. The stabilities of the mutant proteins with charged or extremely bulky residues at the cavity were lower than those expected from the hydrophobicities of the substituted residues, suggesting that considerable strains are created at the mutation sites in these mutant proteins. However, examination of the far- and near-UV CD spectra and the enzymatic activities suggest that all of the mutant proteins have structures similar to that of the wild-type protein. These results suggest that the cavity in the hydrophobic core of E. coli RNase HI is conformationally fairly stable. This may be the reason why the cavity-filling mutations effectively increase the thermal stability of this protein.

NMR identification of hydrophobic cavities with low water occupancies in protein structures using small gas molecules

Magnetization transfer through dipole-dipole interactions (nuclear Overhauser effects, NOEs) between water protons and the protons lining two small hydrophobic cavities in hen egg-white lysozyme demonstrates the presence of water molecules with occupancies of approximately 10-50%. Similarly, NOEs were observed between the cavity protons and the protons of hydrogen, methane, ethylene or cyclopropane applied at 1-200 bar pressure. These gases can thus be used as general NMR indicators of empty or partially hydrated hydrophobic cavities in proteins. All gases reside in the cavities for longer than 1 ns in marked contrast to common belief that gas diffusion in proteins is not much slower than in water. Binding to otherwise empty cavities may be a major aspect of the anesthetic effect of small organic gas molecules.

Properties of the protein matrix revealed by the free energy of cavity formation

BACKGROUND

The classical picture of the hydrophobic stabilization of proteins invokes a resemblance between the protein interior and nonpolar solvents, but the extent to which this is the case has often been questioned. The protein interior is believed to be at least as tightly packed as organic crystals, and was shown to have very low compressibility. There is also evidence that these properties are not uniform throughout the protein, and conflicting views exist on the nature of sidechain packing and on its influence on the properties of the protein.

RESULTS

In order to probe the physical properties of the protein, the free energy associated with the formation of empty cavities has been evaluated for two proteins: barnase and T4 lysozyme. To this end, the likelihood of encountering such cavities was computed from room temperature molecular dynamics trajectories of these proteins in water. The free energy was evaluated in each protein taken as a whole and in submolecular regions. The computed free energies yielded information on the manner in which empty space is distributed in the system, while the latter undergoes thermal motion, a property hitherto not analyzed in heterogeneous media such as proteins. Our results showed that the free energy of cavity formation is higher in proteins than in both water and hexane, providing direct evidence that the native protein medium differs in fundamental ways from the two liquids. Furthermore, although the packing density was found to be higher in nonpolar regions of the protein than in polar ones, the free energy cost of forming atomic size cavities is significantly lower in nonpolar regions, implying that these regions contain larger chunks of empty space, thereby increasing the likelihood of containing atomic size packing defects. These larger empty spaces occur preferentially where buried hydrophobic sidechains belonging to secondary structures meet one another. These particular locations also appear to be more compressible than other parts of the core or surface of the protein.

CONCLUSIONS

The cavity free energy calculations described here provide a much more detailed physical picture of the protein matrix than volume and packing calculations. According to this picture, the packing of hydrophobic sidechains is tight in the interior of the protein, but far from uniform. In particular, the packing is tighter in regions where the backbone forms less regular hydrogen-bonding interactions than at interfaces between secondary structure elements, where such interactions are fully developed. This may have important implications on the role of sidechain packing in protein folding and stability.

The structure of the C-terminal domain of methionine synthase: presenting S-adenosylmethionine for reductive methylation of B12

BACKGROUND

In both mammalian and microbial species, B12-dependent methionine synthase catalyzes methyl transfer from methyltetrahydrofolate (CH3-H4folate) to homocysteine. The B12 (cobalamin) cofactor plays an essential role in this reaction, accepting the methyl group from CH3-H4folate to form methylcob(III)alamin and in turn donating the methyl group to homocysteine to generate methionine and cob(I)alamin. Occasionally the highly reactive cob(I)alamin intermediate is oxidized to the catalytically inactive cob(II)alamin form. Reactivation to sustain enzyme activity is achieved by a reductive methylation, requiring S-adenosylmethionine (AdoMet) as the methyl donor and, in Esherichia coli, flavodoxin as an electron donor. The intact system is controlled and organized so that AdoMet, rather than methyltetrahydrofolate, is the methyl donor in the reactivation reaction. AdoMet is not wasted as a methyl donor in the catalytic cycle in which methionine is synthesized from homocysteine. The structures of the AdoMet binding site and the cobalamin-binding domains (previously determined) provide a starting point for understanding the methyl transfer reactions of methionine synthase.

RESULTS

We report the crystal structure of the 38 kDa C-terminal fragment of E.coli methionine synthase that comprises the AdoMet-binding site and is essential for reactivation. The structure, which includes residues 901-1227 of methionine synthase, is a C-shaped single domain whose central feature is a bent antiparallel betasheet. Database searches indicate that the observed polypeptide has no close relatives. AdoMet binds near the center of the inner surface of the domain and is held in place by both side chain and backbone interactions.

CONCLUSIONS

The conformation of bound AdoMet, and the interactions that determine its binding, differ from those found in other AdoMet-dependent enzymes. The sequence Arg-x-x-x-Gly-Tyr is critical for the binding of AdoMet to methionine synthase. The position of bound AdoMet suggests that large areas of the C-terminal and cobalamin-binding fragments must come in contact in order to transfer the methyl group of AdoMet to cobalamin. The catalytic and activation cycles may be turned off and on by alternating physical separation and approach of the reactants.

Evidence on close packing and cavities in proteins

The packing of a protein's constituent atoms and the attendant constraints placed upon them form the basis of many attempts to understand and predict protein structure, stability, folding and even function. Although the significance of packing is yet to be fully comprehended, recent experimental and theoretical investigations have increased our understanding through the description of mutational effects on structure and stability, determination of the limits of packing constraints for both protein folding and structure prediction, and delineation of packing guidelines on the basis of observed cavities in the native protein folds. These advances and allowing protein modellers, engineers and designers to tackle their problems from a more rational perspective.

Cavities and packing at protein interfaces

An analysis of internal packing defects or "cavities" (both empty and water-containing) within protein structures has been undertaken and includes 3 cavity classes: within domains, between domains, and between protein subunits. We confirm several basic features common to all cavity types but also find a number of new characteristics, including those that distinguish the classes. The total cavity volume remains only a small fraction of the total protein volume and yet increases with protein size. Water-filled "cavities" possess a more polar surface and are typically larger. Their constituent waters are necessary to satisfy the local hydrogen bonding potential. Cavity-surrounding atoms are observed to be, on average, less flexible than their environments. Intersubunit and interdomain cavities are on average larger than the intradomain cavities, occupy a larger fraction of their resident surfaces, and are more frequently water-filled. We observe increased cavity volume at domain-domain interfaces involved with shear type domain motions. The significance of interfacial cavities upon subunit and domain shape complementarity and the protein docking problem, as well as in their structural and functional role in oligomeric proteins, will be discussed. The results concerning cavity size, polarity, solvation, general abundance, and residue type constituency should provide useful guidelines for protein modeling and design.

Buried waters and internal cavities in monomeric proteins

We have analyzed the buried water molecules and internal cavities in a set of 75 high-resolution, nonhomologous, monomeric protein structures. The number of hydrogen bonds formed between each water molecule and the protein varies from 0 to 4, with 3 being most common. Nearly half of the water molecules are found in pairs or larger clusters. Approximately 90% are shown to be associated with large cavities within the protein, as determined by a novel program, PRO_ACT. The total volume of a protein's large cavities is proportional to its molecular weight and is not dependent on structural class. The largest cavities in proteins are generally elongated rather than globular. There are many more empty cavities than hydrated cavities. The likelihood of a cavity being occupied by a water molecule increases with cavity size and the number of available hydrogen bond partners, with each additional partner typically stabilizing the occupied state by 0.6 kcal/mol.

An analysis of packing in the protein folding problem

The number of globular proteins for which high resolution structures are available is rapidly increasing. In each case the particular sequence of the polypeptide appears to yield only a single, compact, biologically active structure. However, peptides with no obvious sequence similarity may form remarkably similar structures.

Prediction and analysis of structure, stability and unfolding of thermolysin-like proteases

Bacillus neutral proteases (NPs) form a group of well-characterized homologous enzymes, that exhibit large differences in thermostability. The three-dimensional (3D) structures of several of these enzymes have been modelled on the basis of the crystal structures of the NPs of B. thermoproteolyticus (thermolysin) and B. cereus. Several new techniques have been developed to improve the model-building procedures. Also a 'model-building by mutagenesis' strategy was used, in which mutants were designed just to shed light on parts of the structures that were particularly hard to model. The NP models have been used for the prediction of site-directed mutations aimed at improving the thermostability of the enzymes. Predictions were made using several novel computational techniques, such as position-specific rotamer searching, packing quality analysis and property-profile database searches. Many stabilizing mutations were predicted and produced: improvement of hydrogen bonding, exclusion of buried water molecules, capping helices, improvement of hydrophobic interactions and entropic stabilization have been applied successfully. At elevated temperatures NPs are irreversibly inactivated as a result of autolysis. It has been shown that this denaturation process is independent of the protease activity and concentration and that the inactivation follows first-order kinetics. From this it has been conjectured that local unfolding of (surface) loops, which renders the protein susceptible to autolysis, is the rate-limiting step. Despite the particular nature of the thermal denaturation process, normal rules for protein stability can be applied to NPs. However, rather than stabilizing the whole protein against global unfolding, only a small region has to be protected against local unfolding. In contrast to proteins in general, mutational effects in proteases are not additive and their magnitude is strongly dependent on the location of the mutation. Mutations that alter the stability of the NP by a large amount are located in a relatively weak region (or more precisely, they affect a local unfolding pathway with a relatively low free energy of activation). One weak region, that is supposedly important in the early steps of NP unfolding, has been determined in the NP of B. stearothermophilus. After eliminating this weakest link a drastic increase in thermostability was observed and the search for the second-weakest link, or the second-lowest energy local unfolding pathway is now in progress. Hopefully, this approach can be used to unravel the entire early phase of unfolding.

Shape distributions of protein topography

A method is presented for measuring protein surface shape. It is an improvement of an earlier method that intersects a sphere with the solvent-excluded volume of a protein molecule. The new method, called a shape distribution, produces a more sophisticated description of the region of the sphere inside the protein than is provided by simply measuring the region's area or solid angle. This method is applied to the prediction of molecular complexes in three systems: the hemoglobin nonallosteric interface, trypsin and trypsin inhibitor, and heme and apomyoglobin. It does not uniquely predict the correct structure, even though the individual structures are taken from the experimentally determined complex structure. However, it does provide a list of several hundred predicted complexes, one of which is correct, and from which the correct complex might be extracted by a subsequent chemical filter.

A cavity-containing mutant of T4 lysozyme is stabilized by buried benzene

The hydrophobic cores of proteins are generally well packed, with few cavities. Mutations in which a bulky buried residue such as leucine or phenylalanine is replaced with a small residue such as alanine can create cavities in the core of a protein (our unpublished results). The sizes and shapes of such cavities can vary substantially depending on factors such as local geometry, whether or not a cavity already exists at the site of substitution, and the degree to which the protein structure relaxes to occupy the space vacated by the substituted residue. We show by crystallographic and thermodynamic analysis that the cavity created by the replacement Leu 99----Ala in T4 lysozyme is large enough to bind benzene and that ligand binding increases the melting temperature of the protein by 6.0 degrees C at pH 3.0. Benzene does not, however, bind to the cavity created by the Phe 153----Ala replacement. The results show that cavities can be engineered in proteins and suggest that such cavities might be tailored to bind specific ligands. The binding of benzene at an internal site 7 A from the molecular surface also illustrates the dynamic nature of proteins, even in crystals.

Hydrophobic packing in T4 lysozyme probed by cavity-filling mutants

To probe the nature of the hydrophobic cores of proteins and to test potential ways of increasing protein thermostability, an attempt was made to improve the packing within T4 bacteriophage lysozyme by engineered amino acid replacements. Two mutations, Leu-133----Phe and Ala-129----Val, which were designed to fill the largest cavities that exist in the folded structure of the native protein, were constructed. The mutant proteins have normal activities and their thermal stabilities are marginally lower than that of wild-type lysozyme. Crystal structure analysis of the mutant proteins shows that the introduced amino acids are accommodated with very little perturbation of the three-dimensional structure. Incorporation of the more bulky hydrophobic residues within the core of the protein is expected to provide an increase in hydrophobic stabilization, but this is seen to be offset by the introduction of strain. Inspection of the mutant structures shows that in each case the introduced amino acid side chain is forced to adopt a non-optimal dihedral angle X1. Strain is also observed in the form of bond angle distortion and in unfavorable van der Waals contacts. The results illustrate how the observed core structures of proteins represent a compromise between the hydrophobic effect, which will tend to maximize the core packing density, and the strain energy that would be incurred in eliminating all packing defects. The results also suggest that mutations designed to increase protein stability by filling existing cavities may be effective in some cases but are unlikely to provide a general method for increasing protein stability.

Protein-ligand dynamics. A 96 picosecond simulation of a myoglobin-xenon complex

A 96 picosecond dynamics trajectory of myoglobin with five xenon-probe ligands in internal cavities is examined to study the effect of protein motions on ligand motion and internal cavity fluctuations. Average structural and energetic properties indicate that the simulation is well behaved. The average protein volume is similar to the volume of the X-ray model and the main-chain atom root-mean-square deviation between the X-ray model and the average dynamical structure is 1.25 A. The protein volume oscillates 3 to 4% around the volume of the X-ray structure. These fluctuations lead to changes in the internal free volume and in the size, shape and location of atom-sized cavity features. Transient cavities produced in the simulation have a crucial role in the movement of two of the ligands. One of the ligands escapes to the protein surface, whilst a second ligand travels through the protein interior. Complex gating processes involving several protein residues are responsible for producing the necessary pores through which the ligand passes between transient cavities or packing defects.

Criteria that discriminate between native proteins and incorrectly folded models

Various theoretical concepts, such as free energy potentials, electrostatic interaction potentials, atomic packing, solvent-exposed surface, and surface charge distribution, were tested for their ability to discriminate between native proteins and misfolded protein models. Misfolded models were constructed by introducing incorrect side chains onto polypeptide backbones: side chains of the alpha-helical hemerythrin were modeled on the beta-sheeted backbone of immunoglobulin VL domain, whereas those of the VL domain were similarly modeled on the hemerythrin backbone. CONGEN, a conformational space sampling program, was used to construct the side chains, in contrast to the previous work, where incorrect side chains were modeled in all trans conformations. Capability of the conformational search procedure to reproduce native conformations was gauged first by rebuilding (the correct) side chains in hemerythrin and the VL domain: constructs with r.m.s. differences from the x-ray side chains 2.2-2.4 A were produced, and many calculated conformations matched the native ones quite well. Incorrectly folded models were then constructed by the same conformational protocol applied to incorrect amino acid sequences. All CONGEN constructs, both correctly and incorrectly folded, were characterized by exceptionally small molecular surfaces and low potential energies. Surface charge density, atomic packing, and Coulomb formula-based electrostatic interactions of the misfolded structures and the correctly folded proteins were similar, and therefore of little interest for diagnosing incorrect folds. The following criteria clearly favored the native structures over the misfolded ones: 1) solvent-exposed side-chain nonpolar surface, 2) number of buried ionizable groups, and 3) empirical free energy functions that incorporate solvent effects.